Abstract: A method of generating scintillation prediction data comprises: •a) for a plurality of satellites (12A), and reference stations (10A, 10B, 10C, 10D) measuring phase scintillation data, satellite by satellite for each reference station during multiple epochs (t?2, t?1, t, t+k); •b) forecasting an expected phase scintillation value for each satellite and reference station for a period of at least 24 hours based on a cyclical prediction model; •c) for a given user location (20) and a given satellite, spatially interpolating the expected phase scintillation values of the plurality of reference stations to determine a predicted phase scintillation index; •d) repeating step c) for further satellites visible from the user location.

Abstract: A method of generating scintillation prediction data comprises: •a) for a plurality of satellites (12A), and reference stations (10A, 10B, 10C, 10D) measuring phase scintillation data, satellite by satellite for each reference station during multiple epochs (t?2, t?1, t, t+k); •b) forecasting an expected phase scintillation value for each satellite and reference station for a period of at least 24 hours based on a cyclical prediction model; •c) for a given user location (20) and a given satellite, spatially interpolating the expected phase scintillation values of the plurality of reference stations to determine a predicted phase scintillation index; •d) repeating step c) for further satellites visible from the user location.