The use of various interplanetary scintillation indices within geomagnetic forecasts

Interplanetary scintillation (IPS), the twinkling of small angular diameter radio sources, is caused by the interaction of the signal with small-scale plasma irregularities in the solar wind. The technique may be used to sense remotely the near-Earth heliosphere and observations of a sufficiently large number of sources may be used to track large-scale disturbances as they propagate from close to the Sun to the Earth. Therefore, such observations have potential for use within geomagnetic forecasts. We use daily data from the Mullard Radio Astronomy Observatory, made available through the World Data Centre, to test the success of geomagnetic forecasts based on IPS observations. The approach discussed here was based on the reduction of the information in a map to a single number or series of numbers. The advantages of an index of this nature are that it may be produced routinely and that it could ideally forecast both the occurrence and intensity of geomagnetic activity. We start from an index that has already been described in the literature, INDEX35. On the basis of visual examination of the data in a full skymap format modifications were made to the way in which the index was calculated. It was hoped that these would lead to an improvement in its forecasting ability. Here we assess the forecasting potential of the index using the value of the correlation coefficient between daily Ap and the IPS index, with IPS leading by 1 day. We also compare the forecast based on the IPS index with forecasts of Ap currently released by the Space Environment Services Center (SESC). Although we find that the maximum improvement achieved is small, and does not represent a significant advance in forecasting ability, the IPS forecasts at this phase of the solar cycle are of a similar quality to those made by SESC.     

Proposed major format change to geomagnetic activity reports and forecasts produced by the SESC,Boulder, Colorado,U.S.A.

The Space Environment Services Center (SESC) provides real-time information on the state of the solar-terrestrial environment. In particular, the geomagnetic field is closely monitored. Summaries of geomagnetic activity are issued daily, as are forecasts for each of the next three days. Customers are alerted when activity rises above selected thresholds of theK- andA-indices. This present geomagnetic reporting and forecasting format is no longer the most appropriate for the expanding number of agencies and industries which the SESC serves. These customers need more detailed information about the geographic distribution of geomagnetic activity stated in plainer language. The proposed new format would quantitatively describe geomagnetic activity in terms of 6 qualitative categories: quiet, unsettled, active, minor storm, major storm, and severe storm. Summaries and forecasts would be made for 2 specific latitude zones: mid-latitudes (20–50 deg geomagnetic) and high-latitudes (50–80 deg geomagnetic). The format adapts well to longer-range (27-day) forecasts.     

Non-linear predictions ofAp by activity class and numerical value

Two neural network algorithms are applied to the short-term,1 to 3 days, prediction of theAp geomagnetic index. A multi-layer, back-propagation (MBP) network is used to implement a self-prediction filter forAp and this provides a forecast of the numerical value of the index. A probabilistic neural network (PNN) is used to estimate the probability distribution of theAp index, in six activity classes, and to provide a forecast of the single most likely activity class for each day. BothAp and an index of solar activity, based on the daily reports issued by the Space Environment Services Centre (Boulder), are input to the probabilistic net. It is found that the numerical forecasts of the MBP filter are most accurate at low, non-storm, levels of activity. This non-linear method provides quantitatively better estimates of activity than are produced by an existing linear prediction filter, particularly with increasing forward forecasting lag. At high levels of the solar activity index the PNN is found to anticipate storm classAp with around 60% accuracy in 1992 and 1993. Some details of the algorithms and implementation issues are described. It is concluded that interplanetary field and solar wind data will be significant components of any of the possible future developments which are discussed.     

大别山库区降水预报性能评估及应用对策

对降水预报进行性能评估及应用对策研究可以更好地发挥降水预报在水库调度中的决策支持作用.基于大别山库区近10 a汛期(2007—2016年5月1日—9月30日)24~168 h共7个预见期降水预报和地面降水观测资料,采用正确率、TS评分、概率统计、ROC曲线以及CTS等方法评估大别山库区降水预报性能,并以响洪甸水库为重点研究区域分析降水预报在水库调度中的应用对策.结果表明:1)大别山库区各量级的降水预报都有正预报技巧;24~72 h预见期降水预报的TS评分较高且空报率、漏报率也较低,具有较高的预报性能;但96 h及以上预见期降水预报性能明显下降,中雨以上量级空报率、漏报率较大,特别是对大暴雨及其以上量级的降水预报性能显著下降.2)大别山库区预报降水量级与实况降水量级基本符合,预报降水量级大于等于实况降水量级的概率超过75%;虽然降水预报量级上呈现出过度预报的现象,但降水过程预报对水库调度仍有较好的应用价值,应用时要考虑到降水预报量级可能存在偏差.3)转折性天气预报96 h及以上预见期CTS评分较低,但72 h以内预见期的性能明显改进,尤其是24 h预见期CTS评分也提高到了38.2%;水库调度可从长预见期的降水预报获取降水过程及其可能发生转折的信息,根据短预见期的降水预报进行调度方案调整.     

Relation of foF2 to parameters of the stratosphere: Specification of the dependence on magnetic activity

The dependence of the correlation coefficient r(h, fo) between the stratospheric parameter h(100) and the critical frequency foF2 on geomagnetic activity is studied. In one of the previous publications of the authors, a general conclusion was drawn that if one chooses for the analysis only the days with geomagnetic index Ap < 12, this coefficient grows. In this paper we present a more complicated and statistically substantiated analysis of the effect of the r(h, fo) dependence on the limitations on Ap imposed on the analyzed days. Wide range of limitations on Ap from Ap = 4 to Ap = 40, five stations, and years of high (1980) and low (1987) solar activity are considered.     

关于利用行星际闪烁（IPS）观测进行地球磁暴急始预报的初步探讨

本文在分析行星际激波传播的空间特性的基础上,提出了一种物理模式加统计的预报方法。对1975-1982年间58个耀斑-IPS激波事件是否引起地球磁暴急始进行预报的结果如下:有无磁暴急始发生的报准率是77.6％,虚报率为22.4％;磁暴急始发生时间的预报准确性,大多在±10小时内。这表明利用IPS观测进行磁暴急始预报有很大的潜力和现实可能性。     

The Impact of Weather Forecasts of Various Lead Times on Snowmaking Decisions Made for the 2010 Vancouver Olympic Winter Games

The Vancouver 2010 Winter Olympics were held from 12 to 28 February 2010, and the Paralympic events followed 2 weeks later. During the Games, the weather posed a grave threat to the viability of one venue and created significant complications for the event schedule at others. Forecasts of weather with lead times ranging from minutes to days helped organizers minimize disruptions to sporting events and helped ensure all medal events were successfully completed. Of comparable importance, however, were the scenarios and forecasts of probable weather for the winter in advance of the Games. Forecasts of mild conditions at the time of the Games helped the Games' organizers mitigate what would have been very serious potential consequences for at least one venue. Snowmaking was one strategy employed well in advance of the Games to prepare for the expected conditions. This short study will focus on how operational decisions were made by the Games' organizers on the basis of both climatological and snowmaking forecasts during the pre-Games winter. An attempt will be made to quantify, economically, the value of some of the snowmaking forecasts made for the Games' operators. The results obtained indicate that although the economic value of the snowmaking forecast was difficult to determine, the Games' organizers valued the forecast information greatly. This suggests that further development of probabilistic forecasts for applications like pre-Games snowmaking would be worthwhile.     

Correlation between the near-Earth solar wind parameters and the source surface magnetic field

The solar magnetic field B _s at the Earth’s projection onto the solar-wind source surface has been calculated for each day over a long time interval (1976–2004). These data have been compared with the daily mean solar wind (SW) velocities and various components of the interplanetary magnetic field (IMF) near the Earth. The statistical analysis has revealed a rather close relationship between the solar-wind parameters near the Sun and near the Earth in the periods without significant sporadic solar and interplanetary disturbances. Empirical numerical models have been proposed for calculating the solar-wind velocity, IMF intensity, and IMF longitudinal and B _z components from the solar magnetic data. In all these models, the B _s value plays the main role. It is shown that, under quiet or weakly disturbed conditions, the variations in the geomagnetic activity index Ap can be forecasted for 3–5 days ahead on the basis of solar magnetic observations. Such a forecast proves to be more reliable than the forecasts based on the traditional methods.     

Probability gains of an epidemic-type aftershock sequence model in retrospective forecasting of <Emphasis Type="Italic">M</Emphasis> ≥ 5 earthquakes in Italy

A stochastic triggering (epidemic) model incorporating short-term clustering was fitted to the instrumental earthquake catalog of Italy for event with local magnitudes 2.6 and greater to optimize its ability to retrospectively forecast 33 target events of magnitude 5.0 and greater that occurred in the period 1990–2006. To obtain an unbiased evaluation of the information value of the model, forecasts of each event use parameter values obtained from data up to the end of the year preceding the target event. The results of the test are given in terms of the probability gain of the epidemic-type aftershock sequence (ETAS) model relative to a time-invariant Poisson model for each of the 33 target events. These probability gains range from 0.93 to 32000, with ten of the target events yielding a probability gain of at least 10. As the forecasting capability of the ETAS model is based on seismic activity recorded prior to the target earthquakes, the highest probability gains are associated with the occurrence of secondary mainshocks during seismic sequences. However, in nine of these cases, the largest mainshock of the sequence was marked by a probability gain larger than 50, having been preceded by previous smaller magnitude earthquakes. The overall evaluation of the performance of the epidemic model has been carried out by means of four popular statistical criteria: the relative operating characteristic diagram, the R score, the probability gain, and the log-likelihood ratio. These tests confirm the superior performance of the method with respect to a spatially varying, time-invariant Poisson model. Nevertheless, this method is characterized by a high false alarm rate, which would make its application in real circumstances problematic.     

Verification of SPE probability forecasts at the Space Environment Prediction Center (SEPC)

In space weather forecasting, forecast verification is necessary so that the forecast quality can be assessed. This paper provides an example of how to choose and devise verification methods and techniques according to different space weather forecast products. Solar proton events (SPEs) are hazardous space weather events, and forecasting them is one of the major tasks of the Space Environment Prediction Center (SEPC) at the National Space Science Center of the Chinese Academy of Sciences. Through analyzing SPE occurrence characteristics, SPE forecast properties, and verification requirements at SEPC, verification methods for SPE probability forecasts are identified, and verification results obtained. Overall, SPE probability forecasts at SEPC exhibit good accuracy, reliability, and discrimination. Compared with climatology and persistence forecasts, the SPE forecasts are more accurate. However, the forecasts for SPE onset days are substantially underestimated and need to be considerably improved.     

Forecasting geomagnetic activity

A study of the daily, monthly and annual values of the geomagnetic activity indexAp for 1932 onwards and of the annualaa indices for 1868 onwards indicated that their variations had large random components. Long-term predictions were not possible from time-series extrapolations.     

Geoeffectiveness of XRA events associated with RSP II and/or RSP IV estimated using the artificial neural network

A forecasting scheme of geomagnetic activity is presented, based on the analysis of the geoeffectiveness of X-ray flares, accompanied by Type II and/or Type IV radio bursts (RSP) observed on the solar disc in the years 1996–2004. The neural network was used to construct this scheme enabling us to determine the probability, with which flares will be followed by a geomagnetic response of a particular intensity. The successfulness of forecasts produced after the fact depended on the flare class and on the combination of radio-burst types. In the case of RSP IV, 58% of the geomagnetic responses of X-ray flares of at least B class were successful. If only RSP II was observed, the forecast was successful only for flares of the X class (67% of successful forecasts). In the second step, a strong geomagnetic response was correctly forecast after geoeffective flares in 58% of the cases. The results are in a good agreement with recent papers based on physical modelling. fridrich@geomag.sk ph@ig.cas.cz, jboch@ig.cas.cz     

Manifestation of planetary wave-type oscillations in variations in the critical frequencies of the ionospheric <Emphasis Type="Italic">F</Emphasis>2 layer in the Asian region

Results of studies of the wave structure of the critical frequencies of the ionospheric F2 layer with periods of planetary waves for two Asian stations—Irkutsk and Wuhan (China)—are presented. Estimates of the appearance frequency, amplitudes, and the lifetime of oscillations with periods typical of planetary waves (2–25 days) are obtained. It is shown that these characteristics depend on the season and place of observation. The appearance of joint periodicities in the critical frequencies at both stations, as well as in the planetary index of geomagnetic activity Ap, is noted.     

Solar origin of geomagnetic storms and prediction of storms with the use of neural networks

This review deals with how the changes of the large-scale solar magnetic fields are related to the occurrence of solar phenomena, which are associated with geomagnetic storms. The review also describes how artificial neural networks have been used to forecast geomagnetic storms either from daily solar input data or from hourly solar wind data. With solar data as input predictions 1–3 days or a month in advance are possible, while using solar wind data as input predictions about an hour in advance are possible. The predictions one hour ahead of the geomagnetic storm indexD _st from only solar wind input data have reached such high accuracy, that they are of practical use in combination with real-time solar wind observations at L1. However, the predictions days and a month ahead need to be much improved in order to be of real practical use.     

Relation between daytime and nighttime values of the critical frequency foF2

The relation between the daytime in the nighttime values of the critical frequencies (foF2) of the ionospheric F ₂ layer is considered. The correlation coefficient of foF2 measured at 1400 and 0200 LT of the same day is considered in various seasons of years with different solar activity (during the complete cycle of solar activity in 1979–1989). Special accent is made on the dependencies of the above mentioned correlation on a choice of magnetically quiet days with various limitations on maximal values of geomagnetic index Ap. It has been obtained that a statistically significant negative correlation between the foF2(1400) and foF2(0200) is more pronounced in the periods of high solar activity. The effect increases with increasing limitation of the considered days on value of Ap: the largest values of the correlation coefficient are observed if only very quiet days are considered (Ap < 6). There are preliminary indications that the considered relation between daytime and nighttime foF2 values depends on latitude.     

Seismic pattern recognition techniques to predict large eruptions at the Popocatépetl,Mexico, volcano

Using pattern recognition techniques, we formulate a simple prediction rule for a retrospective prediction of the three last largest eruptions of the Popocatépetl, Mexico, volcano that occurred on 23 April–30 June 1997 (Eruption 1; VEI ~ 2–3); 11 December 2000–23 January 2001 (Eruption 2; VEI ~ 3–4) and 7 June–4 September 2002 (Eruption 3; explosive dome extrusion and destruction phase). Times of Increased Probability (TIP) were estimated from the seismicity recorded by the local seismic network from 1 January 1995 to 31 December 2005. A TIP is issued when a cluster of seismic events occurs under our algorithm considerations in a temporal window several days (or weeks) prior to large volcanic activity providing sufficient time to organize an effective alert strategy. The best predictions of the three analyzed eruptions were obtained when averaging seismicity rate over a 5-day window with a threshold value of 12 events and declaring an alarm for 45 days. A TIP was issued about six weeks before Eruption 1. TIPs were detected about one and four weeks before Eruptions 2 and 3, respectively. According to our objectives, in all cases, the observed TIPs would have allowed the development of an effective civil protection strategy. Although, under our model considerations the three eruptive events were successfully predicted, one false alarm was also issued by our algorithm. An analysis of the epicentral and depth distribution of the local seismicity used by our prediction rule reveals that successful TIPs were issued from microearthquakes that took place below and towards SE of the crater. On the contrary, the seismicity that issued the observed false alarm was concentrated below the summit of the volcano. We conclude that recording of precursory seismicity below and SE of the crater together with detection of TIPs as described here, could become an important tool to predict future large eruptions at Popocatépetl. Although our model worked well for events that occurred in the past, it is necessary to verify the real capability of the model for future eruptive events.     

华北地区中强以上地震与前兆震群的相关分析

应用“八五”攻关推广的有关震群参数计算的软件 ,对华北地区 (33°～ 43°N,1 1 4°～ 1 2 4°E)自 1 970年以来所发生的 40次震群事件的各项参数 (b,U,H,K,F,ρ等值 )进行了计算 ,旨在依据前兆震群的判定指标 ,寻找前兆震群与中强以上地震 (MS≥ 434)的对应关系。结果表明 ,无论是单项指标或综合指标的报准率均在 65%左右 ,仅有个别达到 75%以上 ,并且有约30 %具有前兆意义的震群属于漏报和虚报。所以 ,如何判定前兆震群和充分利用前兆震群预报地震和监视未来地震活动有待更深入的研究和探讨。     

Probabilistic radar rainfall nowcasts using empirical and theoretical uncertainty models

Weather radar has a potential to provide accurate short‐term (0–3 h) forecasts of rainfall (i.e. radar nowcasts), which are of great importance in warnings and risk management for hydro‐meteorological events. However, radar nowcasts are affected by large uncertainties, which are not only linked to limitations in the forecast method but also because of errors in the radar rainfall measurement. The probabilistic quantitative precipitation nowcasting approach attempts to quantify these uncertainties by delivering the forecasts in a probabilistic form. This study implements two forms of probabilistic quantitative precipitation nowcasting for a hilly area in the south of Manchester, namely, the theoretically based scheme [ensemble rainfall forecasts (ERF)‐TN] and the empirically based scheme (ERF‐EM), and explores which one exhibits higher predictive skill. The ERF‐TN scheme generates ensemble forecasts of rainfall in which each ensemble member is determined by the stochastic realisation of a theoretical noise component. The so‐called ERF‐EM scheme proposed and applied for the first time in this study, aims to use an empirically based error model to measure and quantify the combined effect of all the error sources in the radar rainfall forecasts. The essence of the error model is formulated into an empirical relation between the radar rainfall forecasts and the corresponding ‘ground truth’ represented by the rainfall field from rain gauges measurements. The ensemble members generated by the two schemes have been compared with the rain gauge rainfall. The hit rate and the false alarm rate statistics have been computed and combined into relative operating characteristic curves. The comparison of the performance scores for the two schemes shows that the ERF‐EM achieves better performance than the ERF‐TN at 1‐h lead time. The predictive skills of both schemes are almost identical when the lead time increases to 2 h. In addition, the relation between uncertainty in the radar rainfall forecasts and lead time is also investigated by computing the dispersion of the generated ensemble members. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of Creep Events as Potential Earthquake Precursors: Application to the Creeping Section of the San Andreas Fault, California

—We report the analysis of over 16 years of fault creep and seismicity data from part of the creeping section of the San Andreas fault to examine and assess the temporal association between creep events and subsequent earthquakes. The goal is to make a long-term evaluation of creep events as a potential earthquake precursor. We constructed a catalog of creep events from available digital creepmeter data and compared it to a declustered seismicity catalog for the area between San Juan Bautista and San Benito, California, for 1980 to 1996. For magnitude thresholds of 3.8 and above and time windows of 5 to 10 days, we find relatively high success rates (40% to 55% 'hits') but also very high false alarm rates (generally above 90%). These success rates are statistically significant (0.0007 < P < 0.04). We also tested the actual creep event catalog against two different types of synthetic seismicity catalogs, and found that creep events are followed closely in time by earthquakes from the real catalog far more frequently than the average for the synthetic catalogs, generally by more than two standard deviations. We find no identifiable spatial pattern between the creep events and earthquakes that are hit or missed. We conclude that there is a significant temporal correlation between creep events and subsequent small to moderate earthquakes, however that additional information (such as from other potential precursory phenomena) is required to reduce the false alarm rate to an acceptable level.     

Earthquake Forecasting in Diverse Tectonic Zones of the Globe

We present a simple method for long- and short-term earthquake forecasting (estimating earthquake rate per unit area, time, and magnitude). For illustration we apply the method to the Pacific plate boundary region and the Mediterranean area surrounding Italy and Greece. Our ultimate goal is to develop forecasting and testing methods to validate or falsify common assumptions regarding earthquake potential. Our immediate purpose is to extend the forecasts we made starting in 1999 for the northwest and southwest Pacific to include somewhat smaller earthquakes and then adapt the methods to apply in other areas. The previous forecasts used the CMT earthquake catalog to forecast magnitude 5.8 and larger earthquakes. Like our previous forecasts, the new ones here are based on smoothed maps of past seismicity and assume spatial clustering. Our short-term forecasts also assume temporal clustering. An important adaptation in the new forecasts is to abandon the use of tensor focal mechanisms. This permits use of earthquake catalogs that reliably report many smaller quakes with no such mechanism estimates. The result is that we can forecast earthquakes at higher spatial resolution and down to a magnitude threshold of 4.7. The new forecasts can be tested far more quickly because smaller events are considerably more frequent. Also, our previous method used the focal mechanisms of past earthquakes to estimate the preferred directions of earthquake clustering, however the method made assumptions that generally hold in subduction zones only. The new approach escapes those assumptions. In the northwest Pacific the new method gives estimated earthquake rate density very similar to that of the previous forecast.