Modeling geomagnetically induced electric field and currents by combining a global MHD model with a local one-dimensional method
|
|
- Catherine Chase
- 6 years ago
- Views:
Transcription
1 SPACE WEATHER, VOL. 10,, doi: /2012sw000772, 2012 Modeling geomagnetically induced electric field and currents by combining a global MHD model with a local one-dimensional method J. J. Zhang, 1,2 C. Wang, 1 and B. B. Tang 1 Received 31 January 2012; revised 10 April 2012; accepted 15 April 2012; published 30 May [1] Geomagnetically induced currents (GIC) flowing in long conductor systems on the ground are a well-known space weather hazard. We develop a new approach to simulating GICs by combining a global MHD model with a local one-dimensional method. As an example, we apply this approach to model the GIC at the Pirttikoski 400 kv transformer of the Finnish high-voltage power system during the space weather event of September The modeled results can capture the main observational features, and the model performances is better than two GIC persistence models, which demonstrates this promising new approach in GIC forecasting. Citation: Zhang, J. J., C. Wang, and B. B. Tang (2012), Modeling geomagnetically induced electric field and currents by combining a global MHD model with a local one-dimensional method, Space Weather, 10,, doi: /2012sw Introduction [2] Induced by temporal variation of the external magnetic field, the geoelectric field occurs at the Earth s surface producing geomagnetically induced currents (GIC) in technological facilities. It may bring about severe electric blackout, communication outage, and corrosion of oil and gas pipelines [Boteler and Pirjola, 1998], leading to economic losses in the society. It is always of interest to model the geomagnetically induced electric field and currents both in basic understanding of the the solar wind-magnetosphereionosphere-earth interaction chain [Viljanen et al., 2004] and in practical application to mitigate or avoid economic losses. [3] Generally speaking, the calculation of GIC in a conductor system is performed in two independent steps. The first step is referred to as geophysical step, which is to determine the geoelectric field in the area of the conductor system under study; the second step is the engineering step, which obtains the induced currents due to the given electric field. Once the induced electric field is obtained, GIC in technological systems can be easily computed. The geophysical part is much more difficult to solve, since the geoelectric field not only is affected by spatially and temporally varying ionospheric currents, but also depends on 1 State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing, China. 2 Also at Graduate University of the Chinese Academy of Sciences, Beijing, China. Corresponding author: J. J. Zhang, State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing , China. (jjzhang@spaceweather.ac.cn) the Earth s conductivity. Theoretical modeling of the geomagnetically induced electric field have been carried out since 1980s. Häkkinen and Pirjola [1986] derived exact formulas for getting the electromagnetic field at the surface of a layered Earth produced by a three-dimensional auroral current system. But the calculations of the method are computationally demanding and thus impractical. The calculation of the surface fields is greatly simplified and accelerated when the complex image method (CIM) is applied [Boteler and Pirjola, 1998; Pirjola and Viljanen, 1998]. It considers an image of the primary ionospheric source at a complex depth to represent the induced telluric current when the Earth has a layered conductivity structure. Viljanen et al. [2004] pointed out that although CIM is a convenient tool for theoretical modeling purposes, it is not the most optimal choice for more operational applications. They showed that the CIM can be replaced by a more convenient local 1-D method. In the local 1-D method, the geoelectric field is determined from ground-based geomagnetic recording based on the plane wave relation between the horizontal electric and magnetic fields at the Earth s surface, a local 1-D Earth s conductivity model with no lateral variations is assumed in the method. Based on the local 1-D method, nowcast research of GIC has been carried out in some high latitude countries to monitor the impact of the geomagnetic perturbation on manmade technological systems. The GIC Now! server in Finland ( has been providing nowcast estimates of GIC in the Finnish natural gas pipeline system since June 2003 [Viljanen et al., 2006]. Boteler et al. [2007] developed a system named Real- Time GIC Simulator to simulate the GIC flowing in a power system in Canada. Copyright 2012 by the American Geophysical Union 1of11
2 [4] Some works have been done to estimate the level of the GIC fluctuations by using the information of the upstream solar wind as input. Weigel et al. [2003] and Wintoft [2005] have reproduced some features of proxies of GIC relatively successfully by means of empirical models. On the basis of elementary physical principles, global magnetohydrodynamics (MHD) models have played a significant role in space weather studies in recent decades. And the ability of the models are being extended to modeling GIC. Directly related to the induced electric field and further to GIC that affect power transmission systems and pipelines, the time derivative of the geomagnetic field has been modeled by using MHD models by Raeder et al. [2001] and Yu and Ridley [2008]. Instead of simulating proxies of GIC, Pulkkinen et al. [2007] took the first attempt to modeling GIC directly based on MHD model, although the GIC sites in their study are fictitious. In their work, they computed the induced electric field on ground from the output of the threedimensional global MHD code BATS-R-US [Powell et al., 1999] by using the complex image method (CIM). The GIC at more than thirty magnetometer stations of IMAGE [Lühr, 1994] and Greenland [Friis-Christensen et al., 1985] magnetometer arrays were simulated during the geomagnetically active period of October The magnetometer stations were considered as GIC sites. Due to lack of the real measured GIC data, they calculated measured GIC from the measured geomagnetic recording at each magnetometer station. At last, comparison between simulated GIC and measured GIC at the geomagnetic stations was made. They found that for the studied event the model captures some central features of the measured GIC, although the modeled peak values are smaller. The GIC modeling approach they proposed has been used to the Solar Shield project, which is launched to establish an experimental system for forecasting and mitigating space weather effects on high-voltage power transmission systems [Pulkkinen et al., 2010a]. [5] As Viljanen et al. [2004] pointed out that the application of the local 1-D method is more convenient than the CIM method which has been used in Pulkkinen et al. s GIC modeling approach [Pulkkinen et al., 2007] in determination of the geoelectric field, we introduce a new approach to model the GIC flowing in long conductor systems on ground by combining a global MHD model with the local 1-D method. This will provide a new modeling process to extend the framework of the MHD-based space weather modeling down to the surface of the Earth and enable us test the capabilities of the MHD model to simulating space weather events from the GIC viewpoint. Moreover, rather than model GIC fluctuations at some fictitious GIC sites as what Pulkkinen et al. [2007] did, we will simulate GIC flowing in a real power system in this study and compare the modeled GIC with the GIC recording both qualitatively and quantitatively. [6] To satisfy the requirement of this study, a specific space weather event and power system are selected. During September 1999, the Solar Terrestrial Energy Program- Research, Applications and Modeling Phase(S-RAMP) Group of the Scientific Committee for Solar Terrestrial Physics (SCOSTEP) coordinated a Space Weather Month campaign to follow the progress of space weather events from their initiation on the Sun to their impacts at the Earth, including their effects on space-based and groundbased technological systems and the assessment of the accuracy of specification and forecasting techniques. GIC measured at two 400 kv transformers of the Finnish highvoltage power system are publicly available for the period of the campaign. The transformers are located at Pirttikoski (PIR) and Rauma (RAU), respectively. It provides us with an opportunity to make comparison between modeled and real measured GIC. Taking the space weather event during September 1999 as an example, we will carry out the modeling process to simulate the GIC flowing along the neutral lead of the PIR 400 kv transformer to the Earth during the event. [7] In the remainder of this paper we will first describe the methodology of the modeling in section 2. In section 3, we present the simulation results of GIC at the PIR 400 kv transformer during the space weather event of September 1999, and make comparison of the simulated and the measured GIC qualitatively. In section 4, two kinds of metrics will be introduced to estimate the performance of the model quantitatively. Finally, we will give discussion and summary in section Methodology 2.1. Numerical Model [8] Many three-dimensional (3-D) MHD codes have been developed to simulate the global behavior of the solar wind-magnetosphere-ionosphere system in literature, for instance, the Lyon-Fedder-Mobarry (LFM) code [Lyon et al., 1998], the BATS-R-US code of the Michigan group [Powell et al., 1999], the European GUMICS code [Janhunen, 1996], and so on. The global MHD code utilized in this study is the so-called piecewise parabolic method with a Lagrangian remap (PPMLR)-MHD code which is developed by Hu et al. [2005, 2007]. The code is an extension of the Lagrangian version of the piecewise parabolic method (PPMLR) [Colella and Woodward, 1984] to MHD. It possesses a formal accuracy of the third order in space and the second order in time, and a low numerical dissipation. [9] The PPMLR-MHD code solves ideal MHD equations for the solar wind-magnetosphere-ionosphere coupling system. The coupling process between magnetosphere and ionosphere has been described in detail by Wang et al. [2011]; here we just give a brief description. An inner magnetosphere boundary is set at r =3R E, and an electrostatic ionosphere is set at r = R E. A magnetosphereionosphere electrostatic coupling model is imbedded between inner boundary and ionosphere to drive the inner boundary convection. From the inner boundary the fieldaligned currents (FAC) are mapped along the Earth s dipole magnetic field lines to the ionosphere where they are used as source term of a two-dimensional Poisson equation for electric potential. Once the potential is 2of11
3 obtained, they are mapped back to the inner boundary to calculate convection velocity. To calculate the conductance tensor of the ionosphere, two models are applied together. For the contribution from the solar EUV radiation, we use a model in which the conductance depends on the solar flux F 10.7 and solar zenith angle c [Moen and Brekke, 1993]. For the auroral region, the model developed by Ahn et al. [1998] is used, in which the conductance is empirically derived from the geomagnetic disturbance. The Hall and Pedersen conductance over the dark polar cap region and subauroral region are given to be constants 2.0 S and 1.0 S, respectively, regardless of the magnetic perturbation level Calculation of Geomagnetic Variation [10] To calculate the geomagnetic disturbance, we derive the equivalent current system (ECS) first, detailed description of the method to obtain the ECS in simulation can be found in Raeder et al. [2001] and Wang et al. [2011]. The total height-integrated horizontal ionospheric current can be generally split into two parts: the poloidal current J P and toroidal current J T [Kamide et al., 1981;Raeder et al., 2001]. Fukushima [1976] pointed out that the geomagnetic disturbance produced by poloidal ionospheric current and fieldaligned current (FAC) cancel under the assumption of radial field lines. And the contribution from the magnetospheric currents could be ignored compared with the ionospheric currents for the high latitude region we concerned, so the ionospheric toroidal current J T could be considered as the equivalent current. Thus we can simply obtain the geomagnetic perturbations on ground by using Biot-Savart s law with J T dbðrþ ¼ m Z 0 R r J T ðr Þ 4p s jr r j 3 d2 r ; ð1þ where r denotes the position of the current source, R is the location of the point of interest on ground. The integration domain is limited to about 1000 km both in latitudinal and longitudinal directions overhead, which is enough to produce most fraction of horizontal perturbation approximately in auroral region [Yu and Ridley, 2008]. The effect of the induced (telluric) current on the geomagnetic variations at the Earth s surface is difficult to quantify for its effect is depend on the ground conductivity and the frequency of the external magnetic field. For simplicity, we neglect the effect of it. This is reasonable since the ratio of telluric current contribution to total horizontal magnetic variations close to strong currents is typically 10 20% [Tanskanen et al., 2001], though the ratio can reach to 40% during very rapid changes such as substorm onsets with the duration of a few minutes, the magnetic field variations on ground are still mainly of ionospheric origin [Pulkkinen et al., 2003]. Neglecting the effect of telluric current means here that the amplitude of the geomagnetic field on ground is slightly underestimated. Further this can cause slightly underestimation of the ground electric field and GIC according to the relations between ground magnetic field and electric field and between the electric field and GIC, which will be introduced later Local 1-D Method [11] Based on the plane wave method [Cagniard, 1953] and the Earth s local 1-D conductivity assumption, i.e., the conductivity depends only on depth, the induced electric field is related to the magnetic field in spectral domain as follows [Viljanen et al., 2004]: Z ~E x ðwþ ¼ ~ 0 ðwþ By ~ ðwþ; E ~ Z y ðwþ ¼ ~ 0 ðwþ Bx ~ ðwþ; m 0 where m 0 is the vacuum permeability, ~ Z 0 (w) is the surface impedance dependent of the thicknesses and conductivities of the layered Earth. The subscripts x and y denote the north and the east in geographic coordinates, respectively. This is the so-called local 1-D method. [12] The time series of geomagnetic variations at the concerned point on ground derived from the MHD model first are transformed from geomagnetic coordinates to geographic coordinates. The computations of the local 1-D method are carried out in the spectral domain, transformations between the temporal and the spectral domain are conducted using Fourier transformation. With the information about the Earth s 1-D local conductivity model and simulated geomagnetic variations, the geoelectric field is thus obtained Calculation of GIC From Electric Field [13] The calculation of GIC in a power system is a straightforward task after the electric field is given. If we take the simplest assumption that the electric field is spatially uniform over the region of the whole system, GIC can obtained from m 0 GICðÞ¼aE t x ðþþbe t y ðþ; t where E x, E y are the modeled values close to or at the GIC site; a and b are constant coefficients, which are determined by the topology and electrical properties of the system of interest. They can be obtained theoretically when knowing the power grid information including the network topology, station coordinates, transformer resistances, transmission line resistances, and station earthing resistances [Lehtinen and Pirjola, 1985]. If the electric field and the corresponding GIC are known they also can be deduced empirically [Pulkkinen et al., 2006; Wik et al., 2008]. But it should be noted that the empirical a and b are only valid for the network configuration during the period when the coefficients are determined. If changes occur in the network, new empirical coefficients a and b need be determined. Equation (3) can also be applied for continuously earthed systems like buried pipelines. Figure 1 gives the whole modeling process schematically. 3. Example of the Usage of the Method 3.1. The September 1999 Event [14] During the Space Weather Month campaign of September 1999, an interplanetary magnetic cloud (IMC) ð2þ ð3þ 3of11
4 3.2. Finnish High-Voltage Power System [15] We obtained GIC data during the event of September 1999 at the PIR and RAU 400 kv transformers of the Finnish high-voltage power system (Figure 4). The locations of the two transformers PIR and RAU, geomagnetic station SGO (Sodankylä Geophysical Observatory) which is the nearest magnetometer station from PIR, Figure 1. Schematic diagram of the modeling process. extending from 22 to 23 September triggered a large storm. It is considered that the source of this event was a faint but complete halo coronal mass ejection (CME) event observed at 0606 UT on 20 September [Lam et al., 2002], it reached Earth 63 hours later with a required travel speed of 660 km/s. Figure 2 shows the solar wind and interplanetary magnetic field (IMF) conditions at the Earth s bow shock nose for 33 hours from 0100 UT on 22 September to 1000 UT on 23 September, including the three components of IMF and solar wind velocity, solar wind temperature and plasma density. Following a forward shock at about 1200 UT on 22 September, an IMC event extends from 2000 UT on 22 September to 0700 UT on 23 September with maximum southward IMF of 22 nt, and the IMF rotated from southward to northward during the passage of the IMC. The minimum SYM H index of the triggered storm is 166 nt, and the minimum AL index is 1465 nt during the period. The indices are plotted in Figure 3. The data presented in Figure 2 and Figure 3 are 1 min High Resolution OMNI (HRO) data from NASA OMNI database ( Figure 2. The solar wind and interplanetary magnetic field (IMF) conditions at the Earth s bow shock nose from 0100 UT on 22 September to 1000 UT on 23 September From top to bottom, plots show the three components of the IMF and solar wind velocity, the solar wind plasma density, and the temperature. 4of11
5 Table 1. Coordinates of GIC Sites PIR and RAU and Geomagnetic Observatories SGO and NGO Station Geographic Latitude Geographic Longitude Geomagnetic Latitude Geomagnetic Longitude PIR SGO RAU NGO Figure 3. AL/AU index and SYM H index from 0100 UT on 22 September to 1000 UT on 23 September and NGO (Nurmijärvi Geophysical Observatory) which is the nearest geomagnetic observatory from RAU are marked in Figure 4. The coordinates of the transformers and geomagnetic observatories are shown in Table 1. In this study we simulate the GIC at the PIR 400 kv transformer only, because the magnetic latitude of the RAU 400 kv transformer is beyond the lowest boundary of the ionosphere that can be included in the MHD-based model. This limitation of the model will be discussed later. [16] In order to calculate GIC, the parameters a and b in equation (3) should be determined for the PIR 400 kv transformer. Without the power grid data, we determine the two parameters empirically. The geomagnetic variations recorded at SGO and a 1-D conductivity model (Table 2) [Viljanen et al., 1999] are used to determine the electric field at PIR by using local 1-D method. The chosen ground conductivity model is a resistive conductivity model which is the typical model for the ground of central Finland. We crudely assume the model is appropriate for the whole power grid in this study. The distance between the PIR 400 kv transformer and SGO is about 100 km and it is sufficiently short to determine the electric field at the GIC site by using the magnetic recordings at SGO, in Finnish GIC studies sufficiently close has been found to be of the order of 200 km [Pulkkinen et al., 2006]. On the basis of equation (3), we obtain a and b by a least squares fit using measured GIC at PIR and the calculated geoelectric field. We choose the most disturbed 4 hours from 2000 UT to 2400 UT on 22 September for the fitting process, and get a = 9.6 Akm/V, b = 5.9 Akm/V, the cross-correlation coefficient between measured and fitted GIC is The fitted GIC and measured GIC during the entire period of 33 hours are shown in Figure 5. Good agreement between the fitted and measured GIC indicates the chosen ground conductivity model is reasonable. Same ground conductivity model and coefficients a, b are applied in the simulation part, the key point of the whole simulation process lies in whether or not the MHD model can reproduce realistic geomagnetic variations at the site of interest, from which the GIC is derived at last. Thus the approach to obtain a and b and the crude assumption of the ground conductivity model in this section can be considered acceptable. [17] The geomagnetic data at SGO is obtained from the IMAGE magnetometer website ( reqform/dataform.html). The time resolution of magnetic field data and GIC data are 10 s originally and then they are averaged to 1 min in order to keep accord with the 1 min resolution of the simulation result. Table 2. Ground Conductivity Model Figure 4. Finnish high-voltage power system. The Pirttikoski (PIR) 400 kv transformer, the Rauma (RAU) 400 kv transformer, and the geomagnetic observatories SGO and NGO are marked on the map. Depth (km) Resistivity (Wm) , , , ,000 > of11
6 Figure 5. Measured (grey) and fitted (black) geomagnetically induced current (GIC) at the PIR 400 kv transformer for the investigated period from 0100 UT on 22 September to 1000 UT on 23 September Simulation Results [18] We simulate the GIC at the PIR 400 kv transformer when the IMC swept across the Earth using the modeling process introduced in section 2. The 1 min solar wind and IMF data which are shown in Figure 2 drive the global PPMLR-MHD model at the inflow boundary without the correction of the time difference for the solar wind propagation from 30 R E to the bow shock nose, which is about 3 min here. While the IMF B x is fixed to keep divergencefree condition in the simulation. First, the geomagnetic variations at SGO are modeled to compute further the induced geoelectric field by using the local 1-D method. One should note that in practical GIC simulation, the geomagnetic variations at the GIC site under study can be calculated directly from the output of the MHD model to obtain the modeled geoelectric field. Here we calculate the ground magnetic disturbance at SGO rather than PIR in order to make the comparison of the measured and simulated magnetic field fluctuations. Moreover, the distance between SGO and PIR is short enough, the magnetic field variations at SGO can be used to determine the geoelectric field at the GIC site as mentioned in section 3.2. [19] In Figure 6 we present the x-component (geographic north) and y-component (geographic east) of the measured and simulated geomagnetic field perturbation at SGO. The simulation generally reproduces the intense magnetic fluctuation caused by the IMC. The measured and simulated time derivative of B x and B y are also presented in the figure. The timing of main activities and the basic fluctuation from the simulation are similar to the observation. The geoelectric field computed from the simulated and measured geomagnetic field at SGO are presented in Figure 7, the 1-D ground conductivity model of Table 1 is adopted in calculation. Similar fluctuation of the main active period between the two results can be seen. The modeled and real measured GIC are shown in Figure 8. We can see that the modeled GIC curve is very similar to the measured GIC curve though the amplitudes of the modeled peak values are lower than those of the measured peak values, and the timing of the major activities between simulation and observation is very close. 4. Estimation of the Model Performance [20] In order to estimate the performance of the model during the investigated event quantitatively, two metrics are used in this section Log-Spectral Distance [21] The spectral power determines the level of disturbance in various frequency ranges that can be generated by the model compared with the observed level of fluctuations. In order to quantify the models capability to reproduce the GIC-related features of the magnetic field power spectra, the metric which is called log-spectral distance was developed by Pulkkinen et al. [2010b], and this part of the study builds on the work of them. The logspectral distance M s is a single number measuring the discrepancy between the observed spectral distribution and that obtained by a model in a given spectral range. The derivation of the M s is described as follows. [22] In spectral domain, equation (3) can be written as G ~ ICðwÞ ¼ a E ~ x þ b E ~ y ¼ a Z0 ~ By ~ b Z0 ~ Bx ~ : m 0 m 0 The relation (4) takes the form of a triangle inequality for power spectra G ~ ICðwÞ a j jj Z m ~ 0 jj B ~ y jþj b jj Z 0 m ~ 0 jj B ~ x j: 0 To calculate the log-spectral distance, the logarithm of the ratio of the right-hand side of equation (5) between observed and simulated variables is defined as m s " jaj Bx ~ m s ðwþ ¼ log obs þjbj B ~ # y obs jaj B ~ x sim þjbj B ~ : ð6þ y sim ð4þ ð5þ 6of11
7 Figure 6. The x-component, y-component, and time derivative of the measured (grey) and the simulated (black) geomagnetic field fluctuation at SGO from 0100 UT on 22 September to 1000 UT on 23 September Here, rather than assume a = b as what Pulkkinen et al. [2010b] did originally, we put a = 9.6 Akm/V and b = 5.9 Akm/V derived in section 3.2 into equation (6). Finally, log-spectral distance M s is computed sffiffiffiffiffiffiffiffiffiffiffiffi 1 X M s ¼ m 2 s N ; ð7þ where N is the number of frequencies. As a dimensionless quantity, M s can be equal or larger than zero, M s = 0 denotes the perfect model performance. To carry out spectral analysis, first we divide the time series of measured and simulated magnetic fluctuations of the whole period we investigated into 2 hours long segments with w 50% overlap between adjacent segments. After being multiplied with the Hann window function, data within each segment are transformed to the spectral domain by Fourier transformation. Spectra from all segments are averaged to form the spectra of measured and simulated magnetic variation in equation (6), and the sum in equation (7) is conducted over the periods min. Finally we get M s = 6.18, the value indicates the performance of the model is not bad according to the work of Pulkkinen et al. [2011], which provided a benchmark for monitoring the performance of space weather models Utility Metric [23] From the user s perspective, they prefer metrics that is helpful for them to make decision to maximize the 7of11
8 Figure 7. The calculated geoelectric field (E x, E y ) from 0100 UT on 22 September to 1000 UT on 23 September The grey line and black line represent geoelectric field calculated from the measured and the simulated geomagnetic variations at SGO, respectively. economic gain. For the complex GIC signal, the basic metric based on data-model correlation can t satisfy the requirement of the users only interested in economic utility. The utility metric (for details, see Weigel et al. [2006] and Pulkkinen et al. [2007]) is more suitable for them. The metric was developed to evaluate a model on the basis of its ability to predict an event. An event is defined as follows: Within a forecast window 0 t T F, the absolute value of GIC exceeds event threshold GIC thres. The windows are moved over the time series in non-overlapping parts, events are recorded for both the measured and the simulated GIC for given T F and GIC thres [Pulkkinen et al., 2007]. The forecast model has utility if the monetary gain U F is greater than zero: U F BN H CN H > 0; ð8þ where N H is the number of correct forecasts, N H is the number of false alarms, C is the cost of taking mitigating action and B is the benefit from having taken mitigating action when an event occurred. A few general assumptions are made to restrict the problem, which are (1) the user takes the same mitigating action following each forecast of an event, (2) an always mitigate strategy yields a net monetary loss for the user, and (3) the user seeks to maximize monetary gain U F. The missed events are not taken into account because here considering utility with respect to a system that was never mitigated, a missed event will cause the same monetary loss for both the reference system and the system taking mitigation actions based on the forecast. [24] In equation (8), B and C are both system dependent and are estimated by the user. Without the detailed Figure 8. Measured (grey) and simulated (black) GIC at the PIR 400 kv transformer from 0100 UT on 22 September to 1000 UT on 23 September of11
9 in Figure 9 are small compared with the values in practical engineering issues. Even so the test approach is considered acceptable, since our primary interest is to test the capability of the MHD-based model in predicting GIC events compared with the persistence models. Practical engineering problems can be taken into account in future work. Figure 9. Forecast ratio R F as a function of the event threshold GIC thres for the simulation result (asterisks) and two different persistence models, model 1 (triangles) and model 2 (crosses). knowledge of these two parameters, Weigel et al. [2006] suggested that the developer of the model can determine whether the model is potentially useful or not by reporting the forecast ratio which is defined as R F N H N H if the inequality R F > C/B is satisfied, then taking action following every forecast will yield U F > 0. The model with a higher R F will provide greater economic utility U F in model comparisons. [25] For the investigated period in this study, we choose the disturbed period from 1700 UT on September 22 to 0700 UT on September 23 to calculate the forecast ratio R F for the simulated GIC as a function of GIC thres. In Figure 9, we show the curve along with the curves for two different persistence models for T F = 60 min as what Pulkkinen et al. [2007] did. The two persistence models are (1) always an event forecast in which during disturbed periods the alarm is switched on, and (2) nearest neighbor forecast in which the forecast for the next forecast window is the same as the observation for current window. As seen from Figure 9, R F generated from the simulation is generally higher than that from the other two persistence models for the PIR 400 kv transformer during the period we examined for T F = 60 min. This means that the application of the MHD-based simulation is capable of producing superior economic benefit than the used persistence models. In the calculation, the GIC thresholds GIC thres are chosen arbitrarily; in reality the threshold should depend on the type of transformer since the sensitivity of different transformer types to the effects of GIC is different. Due to the transformer type used in Finland has a very low GIC sensitivity, large GICs have never caused problems in the Finnish power grid [Elovaara, 2007]. The GIC thres values included ð9þ 5. Discussion and Summary [26] In this paper we propose a new approach to combine the global MHD model with the local 1-D method to model GIC. On one hand it gives a method to extend the framework of space weather modeling which is based on MHD model down to the surface of the Earth; on the other hand it provides a way to predict GIC flowing in technology systems on ground. As an example, we utilizes the new approach to simulate the GIC flowing from the neutral lead of the PIR 400 kv transformer to the Earth in the Finnish high-voltage power grid when an IMC swept over the Earth during the Space Weather Month campaign of September The model captures the main features of the geomagnetic fluctuations which are associated with GIC. The curve of the modeled GIC signal presents similar shape as the measured one, and the timing of major activity is close to perfect. The magnitude of the modeled peak GIC values are somewhat lower than the measured ones, which is expected since we neglect the effect of the telluric currents when calculating the geomagnetic variations on the Earth s surface. Two metrics are applied to estimate the performance of the model, one is the correlation-based metric and the other is user-relevant utility metric. Comparison of the utility metric for the simulation and two different persistence models shows that the global MHD-based modeling outperforms the persistence models. Although using different MHD codes and GIC modeling processes, Pulkkinen et al. [2007] obtained similar result. They found that the MHD-based model used in their work is capable of producing equal or even superior economic benefit to that of the two persistence models during the space weather event they investigated. [27] The main aim of this study is to propose a new MHD-based GIC modeling approach and carry out the approach by a real space weather event. So the specified space weather event and power station are selected because of the data availability of an entire space weather episode including the solar data, solar wind data, geomagnetic data and GIC data. The initial simulation result is encouraging when one power station and one event are considered. It indicates that the MHD-based GIC modeling at least has the potential to forecast the GIC flowing in ground facilities. To further validate the approach, a number of space weather events should be simulated to test its ability in the future work. [28] In spite of the success in reproducing some features of the GIC for the event we investigated, there is still large space to improve the current MHD model for simulating GIC accurately. For numerical consideration, there is a gap 9of11
10 region between the ionosphere model and inner boundary of the global MHD model. A coupling process carried out along the dipole magnetic field lines in that region. The inner boundary is usually set to 2.5R E 3.5R E, which implies that the lowest boundary of the ionosphere does not extend very low in magnetic latitude. This limit the ability of the model to simulate the GIC in technological systems only at region of very high latitude. As mentioned in section 3.2, due to the inner magnetospheric boundary of 3R E set in this study, the model cannot simulate GIC at the RAU 400 kv transformer because its location is beyond the lowest latitudinal boundary (59 in magnetic latitude) of the ionosphere in the MHD-based model. This may be improved by setting the inner magnetospheric boundary as close to the Earth as possible. Additionally, inclusion of model of ring current may promote the capability of the model to simulate the GIC events in middle and low latitude countries. Recently, the number of reports about the GIC-risks in power grids in mid- and low-latitude countries has increased [Gaunt and Coetzee, 2007; Liu et al., 2009]. Modeling GIC events which happen there becomes necessary. Moreover, the inclusion of the ring current model may be important for the development of the region 2 current system in MHD models. [29] To summarize, the focus of this study is to provide a new MHD-based GIC modeling approach. The good performance of the model for the studied space weather event makes us believe this new approach could be a promising way for GIC forecasting in the future. For application purpose, this approach should be tested, validated and improved by more space weather events in the future. [30] Acknowledgments. We thank the institutes who maintain the IMAGE Magnetometer Array for providing the geomagnetic data utilized in this study. We also thank A. Viljanen for providing the GIC data of the PIR 400 kv transformer of the Finnish high-voltage power grid and useful comments on the manuscript. Our acknowledgment also goes to X. J. Wang of Beihang University for helpful discussions. This work was supported by 973 program 2012CB825602, NNSFC grants , , and in part by the Specialized Research Fund for State Key Laboratories of China. References Ahn, B. H., A. D. Richmond, and Y. Kamide (1998), An ionospheric conductance model based on ground magnetic disturbance data, J. Geophys. Res., 103(A7), 14,769 14,780, doi: /97ja Boteler, D. H., and R. J. Pirjola (1998), The complex-image method for calculating the magnetic and electric fields produced at the surface of the Earth by the auroral electrojet, Geophys. J. Int., 132, Boteler, D., L. Trichtchenko, R. Pirjola, J. Parmelee, S. Souksaly, A. Foss, and L. Marti (2007), Real-time simulation of geomagnetically induced currents, in th International Symposium on Elcctromagnetic Compatibility and Electromagnetic Ecology, pp , IEEE Press, New York. Cagniard, L. (1953), Basic theory of the magneto-telluric method of geophysical prospecting, Geophysics, 18(3), Colella, P., and P. R. Woodward (1984), The piecewise parabolic method (PPM) for gas-dynamical simulations, J. Comput. Phys., 54, , doi: / (84) Elovaara, J. (2007), Finnish experiences with grid effects of GIC s, in Space Weather: Research Towards Applications in Europe, Astrophys. Space Sci. Libr. Ser., vol. 344, edited by J. Lilensten, pp , Springer, New York. Friis-Christensen, E., Y. Kamide, A. D. Richmond, and S. Matsushita (1985), Interplanetary magnetic field control of high-latitude electric fields and currents determined from Greenland magnetometer data, J. Geophys. Res., 90(A2), Fukushima, N. (1976), Generalized theorem for no ground magnetic effect of vertical currents connected with Pedersen currents in the uniform-conductivity ionosphere, Rep. Ionos. Space Res. Jpn., 30, Gaunt, C. T., and G. Coetzee (2007), Transformer failures in regions incorrectly considered to have low GIC-risk, in Power Tech, 2007 IEEE Lausanne, pp , IEEE Press, New York. Häkkinen, L., and R. Pirjola (1986), Calculation of electric and magnetic fields due to an electrojet current system above a layered Earth, Geophysica, 22(31 44), Hu, Y. Q., X. C. Guo, and G. Q. Li (2005), Oscillation of quasi-steady Earth s magnetosphere, Chin. Phys. Lett., 22(10), , doi: / x/22/10/073. Hu, Y. Q., X. C. Guo, and C. Wang (2007), On the ionospheric and reconnection potentials of the Earth: Results from global MHD simulations, J. Geophys. Res., 112, A07215, doi: /2006ja Janhunen, P. (1996), GUMICS-3 a global ionosphere-magnetosphere coupling simulation with high ionospheric resolution, in Symposium Proceedings on Environment Modeling for Space-Based Applications, Rep. ESA-SP 392, pp , Eur. Space Agency, Paris. Kamide, Y., A. D. Richimond, and S. Matsushita (1981), Estimation of ionospheric electric fields, ionospheric currents, and field-aligned currents from ground magnetic records, J. Geophys. Res., 86, , doi: /ja086ia02p Lam, H. L., D. H. Boteler, and L. Trichtchenko (2002), Case studies of space weather events from their launching on the Sun to their impacts on power systems on the Earth, Ann. Geophys., 20, , doi: /angeo Lehtinen, M., and R. Pirjola (1985), Currents produced in earthed conductor networks by geomagnetically-induced electric fields, Ann. Geophys., 3(4), Liu, C. M., L. G. Liu, R. Pirjola, and Z. Z. Wang (2009), Calculation of geomagnetically induced currents in mid- to low-latitude power grids based on the plane wave method: A preliminary case study, Space Weather, 7, S04005, doi: /2008sw Lühr, H. (1994), The IMAGE magnetometer network, STEP Int. Newsl., 4(10), 4 6. Lyon, J. G., R. E. Lopez, C. C. Goodrich, M. Wiltberger, and K. Papadopoulos (1998), Simulation of the March 9, 1995, substorm: Auroral brightening and the onset of lobe reconnection, Geophys. Res. Lett., 25(15), , doi: /98gl Moen, J., and A. Brekke (1993), The solar flux influence on quiet time conductances in the auroral ionosohere, Geophys. Res. Lett., 20(10), , doi: /92gl Pirjola, R., and A. Viljanen (1998), Complex image method for calculating electric and magnetic fields produced by an auroral electrojet of finite length, Ann. Geophys., 16(11), , doi: /s Powell, K. G., P. L. Roe, T. J. Linde, T. I. Gombosi, and D. L. De Zeeuw (1999), A solution-adaptive upwind scheme for ideal magnetohydrodynamics, J. Comput. Phys., 154(2), , doi: /jcph Pulkkinen, A., O. Amm, and A. Viljanen (2003), Ionospheric equivalent current distributions determined with the method of spherical elementary current systems, J. Geophys. Res., 108(A2), 1053, doi: /2001ja Pulkkinen, A., A. Viljanen, and R. Pirjola (2006), Estimation of geomagnetically induced current levels from different input data, Space Weather, 4, S08005, doi: /2006sw Pulkkinen, A., M. Hesse, M. Kuznetsova, and L. Rastätter (2007), Firstprinciples modeling of geomagnetically induced electromagnetic fields and currents from upstream solar wind to the surface of the Earth, Ann. Geophys., 25(4), , doi: /angeo Pulkkinen, A., M. Hesse, S. Habib, L. Van der Zel, B. Damsky, F. Policelli, D. Fugate, W. Jacobs, and E. Creamer (2010a), Solar shield: Forecasting and mitigating space weather effects on highvoltage power transmission systems, Nat. Hazards, 53, , doi: /s x. Pulkkinen, A., L. Rastätter, M. Kuznetsova, M. Hesse, A. Ridley, J. Reader, H. Singer, and A. Chulaki (2010b), Systematic evaluation of ground and geostationary magnetic field predictions generated 10 of 11
11 by global magnetohydrodynamic models, J. Geophys. Res., 115, A03206, doi: /2009ja Pulkkinen, A., M. Kuznetsova, A. Ridley, J. Raeder, A. Vapirev, D. Weimer, R. Weigel, M. Wiltberger, G. Millward, and L. Rastätter (2011), Geospace environment modeling challenge: Ground magnetic field perturbations, Space Weather, 9, S02004, doi: /2010sw Raeder, J., Y. L. Wang, T. J. Fuller-Rowell, and H. J. Singer (2001), Global simulation of space weather effects of the Bastille Day storm, Sol. Phys., 204, , doi: /a: Tanskanen, E. I., A. Viljanen, T. I. Pulkkinen, R. Pirjola, L. Häkkinen, A. Pulkkinen, and O. Amm (2001), At substorm onset, 40% of AL comes from underground, J. Geophys. Res., 106(A7), 13,119 13,134, doi: /2000ja Viljanen, A., O. Amm, and R. Pirjola (1999), Modeling geomagnetically induced currents during different ionospheric situations, J. Geophys. Res., 104(A12), 28,059 28,071, doi: /1999ja Viljanen, A., et al. (2004), Fast computation of the geoelectric field using the method of elementary current systems and planar Earth models, Ann. Geophys., 22(1), , doi: /angeo Viljanen, A., A. Pulkkinen, R. Pirjola, K. Pajunpää, P. Posio, and A. Koistinen (2006), Recordings of geomagnetically induced currents and a nowcasting service of the Finnish natural gas pipeline system, Space Weather, 4, S10004, doi: /2006sw Wang, C., J. J. Zhang, B. B. Tang, and S. Y. Fu (2011), Comparison of equivalent current systems for the substorm event of 8 March 2008 derived from the global PPMLR-MHD model and the KRM algorithm, J. Geophys. Res., 116, A07207, doi: /2011ja Weigel, R. S., A. J. Klimas, and D. Vassiliadis (2003), Solar wind coupling to and predictability of ground magnetic fields and their time derivatives, J. Geophys. Res., 108(A7), 1298, doi: /2002ja Weigel, R. S., T. Detman, E. J. Rigler, and D. N. Baker (2006), Decision theory and the analysis of rare event space weather forecasts, Space Weather, 4, S05002, doi: /2005sw Wik, M., A. Viljanen, R. Pirjola, A. Pulkkinen, P. Wintoft, and H. Lundstedt (2008), Calculation of geomagentically induced currents in the 400 kv power grid in southern Sweden, Space Weather, 6, S07005, doi: /2007sw Wintoft, P. (2005), Study of the solar wind coupling to the time difference horizontal geomagnetic field, Ann. Geophys., 23(5), , doi: /angeo Yu, Y. Q., and A. J. Ridley (2008), Validation of the space weather modeling framework using ground-based magnetometers, Space Weather, 6, S05002, doi: /2007sw of 11
Statistics on Geomagnetically-Induced Currents in the Finnish 400 kv Power System Based on Recordings of Geomagnetic Variations
J. Geomag. Geoelectr., 41, 411-420,1989 Statistics on Geomagnetically-Induced Currents in the Finnish 400 kv Power System Based on Recordings of Geomagnetic Variations Ari VILJANEN and Risto PIRJOLA Finnish
More informationSolar wind drivers of large geomagnetically induced currents during the solar cycle 23
Click Here for Full Article SPACE WEATHER, VOL. 6,, doi:10.1029/2007sw000374, 2008 Solar wind drivers of large geomagnetically induced currents during the solar cycle 23 K. E. J. Huttunen, 1 S. P. Kilpua,
More informationRole of IMF B x in the solar wind magnetosphere ionosphere coupling
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010ja015454, 2010 Role of IMF B x in the solar wind magnetosphere ionosphere coupling Z. Peng, 1 C. Wang, 1 and Y. Q. Hu 2 Received 14 March 2010;
More informationRecordings of geomagnetically induced currents and a nowcasting service of the Finnish natural gas pipeline system
SPACE WEATHER, VOL. 4,, doi:10.1029/2006sw000234, 2006 Recordings of geomagnetically induced currents and a nowcasting service of the Finnish natural gas pipeline system A. Viljanen, 1 A. Pulkkinen, 2,3
More informationDependence of magnetic field just inside the magnetopause on subsolar standoff distance: Global MHD results
Article SPECIAL ISSUE Basic Plasma Processes in Solar-Terrestrial Activities April 2012 Vol.57 No.12: 1438 1442 doi: 10.1007/s11434-011-4961-6 SPECIAL TOPICS: Dependence of magnetic field just inside the
More informationThe Solar wind - magnetosphere - ionosphere interaction
The Solar wind - magnetosphere - ionosphere interaction Research seminar on Sun-Earth connections Eija Tanskanen Friday January 27, 2006 12-14 a.m., D115 Outline 1. Basics of the Earth s magnetosphere
More informationGeneration of 100-year geomagnetically induced current scenarios
SPACE WEATHER, VOL. 10,, doi:10.1029/2011sw000750, 2012 Generation of 100-year geomagnetically induced current scenarios A. Pulkkinen, 1,2 E. Bernabeu, 3 J. Eichner, 4 C. Beggan, 5 and A. W. P. Thomson
More informationTime history effects at the magnetopause: Hysteresis in power input and its implications to substorm processes
219 Time history effects at the magnetopause: Hysteresis in power input and its implications to substorm processes M. Palmroth, T. I. Pulkkinen, T. V. Laitinen, H. E. J. Koskinen, and P. Janhunen 1. Introduction
More informationEstimation of geomagnetically induced current levels from different input data
SPACE WEATHER, VOL. 4,, doi:10.1029/2006sw000229, 2006 Estimation of geomagnetically induced current levels from different input data Antti Pulkkinen, 1 Ari Viljanen, 2 and Risto Pirjola 3 Received 21
More informationEffects of neutral point reactors and series capacitors on geomagnetically induced currents in a high voltage electric power transmission system
SPACE WEATHER, VOL. 9,, doi:10.1029/2011sw000715, 2011 Effects of neutral point reactors and series capacitors on geomagnetically induced currents in a high voltage electric power transmission system Esko
More informationA Correlation Study of Steady Magnetospheric Convection in the Northern and Southern Hemispheres
A Correlation Study of Steady Magnetospheric Convection in the Northern and Southern Hemispheres Caroline M. McElhenny 1, Dr. Anna DeJong 1 1 Christopher Newport University Abstract: The purpose of this
More informationRisto Pirjola. Finnish Meteorological Institute, P. O. Box 503, FI Helsinki, Finland
Earth Planets Space, 63, 999 1008, 2011 Eigenvalues and eigenvectors of the transfer matrix involved in the calculation of geomagnetically induced currents in an electric power transmission network Risto
More informationMetrics of model performance for electron fluxes (<200 kev) at geostationary orbit
Metrics of model performance for electron fluxes (
More informationEffect of CME Events of Geomagnetic Field at Indian Station Alibag and Pondicherry
Effect of CME Events of Geomagnetic Field at Indian Station Alibag and Pondicherry Babita Chandel Sri Sai University Palampur, Himachal Pradesh, India Abstract: Space weather activity CMEs, and solar energetic
More informationSimulation of the 23 July 2012 extreme space weather event: What if this extremely rare CME was Earth directed?
SPACE WEATHER, VOL. 11, 671 679, doi:1.12/213sw99, 213 Simulation of the 23 July 212 extreme space weather event: What if this extremely rare CME was Earth directed? Chigomezyo M. Ngwira, 1,2 Antti Pulkkinen,
More informationMagnetospheric modes and solar wind energy coupling efficiency
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009ja014737, 2010 Magnetospheric modes and solar wind energy coupling efficiency T. I. Pulkkinen, 1 M. Palmroth, 1 H.
More informationChapter 8 Geospace 1
Chapter 8 Geospace 1 Previously Sources of the Earth's magnetic field. 2 Content Basic concepts The Sun and solar wind Near-Earth space About other planets 3 Basic concepts 4 Plasma The molecules of an
More informationV r : A new index to represent the variation rate of geomagnetic activity
Earthq Sci (2010)23: 343 348 343 Doi: 10.1007/s11589-010-0731-9 V r : A new index to represent the variation rate of geomagnetic activity Dongmei Yang 1, Yufei He 1 Chuanhua Chen 2 and Jiadong Qian 3 1
More informationINTERPLANETARY ASPECTS OF SPACE WEATHER
INTERPLANETARY ASPECTS OF SPACE WEATHER Richard G. Marsden Research & Scientific Support Dept. of ESA, ESTEC, P.O. Box 299, 2200 AG Noordwijk, NL, Email: Richard.Marsden@esa.int ABSTRACT/RESUME Interplanetary
More informationGlobal Simulation of Magnetospheric Space Weather Effects of the Bastille Day Storm
Global Simulation of Magnetospheric Space Weather Effects of the Bastille Day Storm J. Raeder and Y. L. Wang Institute of Geophysics and Planetary Physics, University of California, Los Angeles T. J. Fuller-Rowell
More informationMitsubishi Electric Power Products Inc. and Duquesne Light Co.
Mitsubishi Electric Power Products Inc. and Duquesne Light Co. Overview on Geomagnetically Induced Current Revision #00 March 2017 Presented by: Elizabeth Cook (DLC) and Wesley Terek (MEPPI) Overview on
More informationThe Dynamic Magnetosphere. Ioannis A. Daglis. National Observatory of Athens, Greece
310/1749-42 ICTP-COST-USNSWP-CAWSES-INAF-INFN International Advanced School on Space Weather 2-19 May 2006 The Dynamic Magnetosphere: Reaction to and Consequences of Solar Wind Variations Yannis DAGLIS
More informationEffect of the altitudinal variation of the gravitational acceleration on the thermosphere simulation
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2008ja013081, 2008 Effect of the altitudinal variation of the gravitational acceleration on the thermosphere simulation Yue Deng, 1 Aaron J. Ridley,
More informationGeomagnetic Disturbances (GMDs) History and Prediction
Geomagnetic Disturbances (GMDs) History and Prediction J. Patrick Donohoe, Ph.D., P.E. Dept. of Electrical and Computer Engineering Mississippi State University Box 9571 Miss. State, MS 39762 donohoe@ece.msstate.edu
More informationDoes the polar cap area saturate?
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L09107, doi:10.1029/2007gl029357, 2007 Does the polar cap area saturate? V. G. Merkin 1 and C. C. Goodrich 1 Received 15 January 2007;
More informationGEOMAGNETICALLY INDUCED CURRENTS IN A BRAZILIAN POWER NETWORK OVER THE SOLAR CYCLES 23 AND 24
GEOMAGNETICALLY INDUCED CURRENTS IN A BRAZILIAN POWER NETWORK OVER THE SOLAR CYCLES 23 AND 24 Cleiton Barbosa 1,2 *, Gelvam A. Hartmann 1, Livia Alves 3, Ramon Caraballo 4, Andrés Papa 1, R.Pirjola 5,6
More informationSupporting Information for The substorm cycle as reproduced by community-available global MHD models
SPACE WEATHER Supporting Information for The substorm cycle as reproduced by community-available global MHD models E. Gordeev 1, V. Sergeev 1, N. Tsyganenko 1, M. Kuznetsova 2, L. Rastätter 2, J. Raeder
More informationExternal Magnetic Field Variations and Aeromagnetic Surveys Experiences, Problems, Potential Solutions
External Magnetic Field Variations and Aeromagnetic Surveys Experiences, Problems, Potential Solutions Jurgen Watermann Hans Gleisner & Thorkild Rasmussen chercheur associé, Le STUDIUM hosted by LPCE/CNRS
More informationVariations of time derivatives of the horizontal geomagnetic field and horizontal geomagnetic field along 96-degree magnetic meridian
Covenant Journal of Physical and Life Sciences (CJPL) Vol. 4 No 1 June, 216 Variations of time derivatives of the horizontal geomagnetic field and horizontal geomagnetic field along 96-degree magnetic
More informationWhat determines when and where reconnection begins
What determines when and where reconnection begins Robert L. McPherron Invited presentation at Unsolved Problems in Magnetospheric Physics, Scarborough, UK, Sept. 6-12. Factors That Might Affect Tail Reconnection
More informationGeomagnetic Disturbance Report Reeve Observatory
Event type: Various geomagnetic disturbances including coronal hole high-speed stream, coronal mass ejection, sudden impulse and reverse shock effects Background: This background section defines the various
More informationObservations of an interplanetary slow shock associated with magnetic cloud boundary layer
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L15107, doi:10.1029/2006gl026419, 2006 Observations of an interplanetary slow shock associated with magnetic cloud boundary layer P. B.
More informationReal-time Earth magnetosphere simulator with three-dimensional magnetohydrodynamic code
SPACE WEATHER, VOL. 4,, doi:10.1029/2004sw000100, 2006 Real-time Earth magnetosphere simulator with three-dimensional magnetohydrodynamic code M. Den, 1 T. Tanaka, 2,3 S. Fujita, 4,5 T. Obara, 5,6 H. Shimazu,
More informationNICT SPACE WEATHER ACTIVITIES
Prepared by NICT Agenda Item: 3 Discussed in Ad-Hoc Meeting on Space Weather NICT SPACE WEATHER ACTIVITIES Tsutomu Nagatsuma Space Weather and Environment Informatics Laboratory, Applied Electromagnetic
More informationDo Countries Under the Equatorial Electrojet Belt Worry About GIC?
Do Countries Under the Equatorial Electrojet Belt Worry About GIC? Endawoke Yizengaw 1 and Brett Carter 2 1 Institute for Scientific Research, Boston College 2 RMIT University, Melbourne, Australia Content
More informationSTUDY ON RELATIONSHIP OF MAGNETOSPHERIC SUBSTORM AND MAGNETIC STORM
Prosiding Seminar Nasional Penelitian, Pendidikan dan Penerapan MIPA Fakultas MIPA, Universitas Negeri Yogyakarta, 16 Mei 2009 STUDY ON RELATIONSHIP OF MAGNETOSPHERIC SUBSTORM AND MAGNETIC STORM L. Muhammad
More informationTwo types of geomagnetic storms and relationship between Dst and AE indexes
Two types of geomagnetic storms and relationship between Dst and AE indexes Lyudmila P. Shadrina 1, * 1 Academy of sciences of Sakha (Yakutia) Republic, Yakutsk, Russia Abstract. The study of the relationship
More informationRelation between substorm characteristics and rapid temporal variations of the ground magnetic field
Relation between substorm characteristics and rapid temporal variations of the ground magnetic field A. Viljanen, E. I. Tanskanen, A. Pulkkinen To cite this version: A. Viljanen, E. I. Tanskanen, A. Pulkkinen.
More informationSubstorm studies with magnetic field measurements from ground-based networks
Substorm studies with magnetic field measurements from ground-based networks Eija I. Tanskanen University of Bergen,Institute of Physics and Technology, Bergen, Norway Finnish Meteorological Institute,
More informationMagnetospheric Model Performance during Conjugate Aurora
18 Magnetospheric Model Performance during Conjugate Aurora William Longley 1, Patricia Reiff 1, Jone Peter Reistad 2, and Nikolai Østgaard 2 Video of Yosemite Talk, URL: http://dx.doi.org/10.15142/t3h01c
More informationSubstorm-associated effects in the variations of low energy electron fluxes in the inner magnetosphere: Does the substorm s strength matter?
Substorm-associated effects in the variations of low energy electron fluxes in the inner magnetosphere: Does the substorm s strength matter? N. Ganushkina (1, 2), S. Dubyagin (1), I. Sillanpää (1), D.
More informationEXTREMELY QUIET 2009 STATE OF GEOMAGNETIC FIELD AS A REFERENCE LEVEL FOR AMPLITUDE OF THE LOCAL GEOMAGNETIC DISTURBANCES
EXTREMELY QUIET 2009 STATE OF GEOMAGNETIC FIELD AS A REFERENCE LEVEL FOR AMPLITUDE OF THE LOCAL GEOMAGNETIC DISTURBANCES A.E. Levitin, L.I. Gromova, S.V. Gromov, L.A. Dremukhina Pushkov Institute of Terrestrial
More informationPC index as a standard of magnetospheric disturbances in the auroral zone
PC index as a standard of magnetospheric disturbances in the auroral zone Oleg Troshichev, Arctic and Antarcrtic Research Institute, St.Petersburg olegtro@aari.ru The Solar-Terrestrial Physics Symposium
More informationA Semi-Empirical Model for Forecasting Relativistic Electrons at Geostationary Orbit
2008 Annual Meeting Theme, 20 24 January 2008, New Orleans, Louisiana Fifth Space Weather Symposium A Semi-Empirical Model for Forecasting Relativistic Electrons at Geostationary Orbit Wladislaw Lyatsky
More information1 Introduction. Cambridge University Press Physics of Space Plasma Activity Karl Schindler Excerpt More information
1 Introduction Space plasma phenomena have attracted particular interest since the beginning of the exploration of space about half a century ago. Already a first set of pioneering observations (e.g.,
More informationComment on Effects of fast and slow solar wind on the correlation between interplanetary medium and geomagnetic activity by P.
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. A10, 1386, doi:10.1029/2002ja009746, 2003 Correction published 20 January 2004 Comment on Effects of fast and slow solar wind on the correlation between interplanetary
More informationSpace Physics. An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres. May-Britt Kallenrode. Springer
May-Britt Kallenrode Space Physics An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres With 170 Figures, 9 Tables, Numerous Exercises and Problems Springer Contents 1. Introduction
More informationInvestigation of the viscous potential using an MHD simulation
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011ja017022, 2012 Investigation of the viscous potential using an MHD simulation R. Bruntz, 1 R. E. Lopez, 1 M. Wiltberger, 2 and J. G. Lyon 3 Received
More informationSpace Weather effects observed on the ground geomagnetic effects
Space Weather effects observed on the ground geomagnetic effects J. Watermann Danish Meteorological Institute, Copenhagen, Denmark my thanks to 1 Ground-based systems affected by space weather ( physics
More information2-2-5 Application of the Magnetospheric Model for Numerical Forecast - Validation Using the AE Index-
2-2-5 Application of the Magnetospheric Model for Numerical Forecast - Validation Using the AE Index- KITAMURA Kentarou, SHIMAZU Hironori, FUJITA Shigeru, WATARI Shinichi, KUNITAKE Manabu, SHINAGAWA Hiroyuki,
More informationInfluence of spatial variations of the geoelectric field on geomagnetically induced currents
J. Space Weather Space Clim. 2017, 7, A22 A. Viljanen and R. Pirjola, Published by EDP Sciences 2017 DOI: 10.1051/swsc/2017024 Available online at: www.swsc-journal.org RESEARCH ARTICLE Influence of spatial
More informationStudy of Geomagnetic Field Variations at Low Latitude of African Equatorial Region
Study of Geomagnetic Field Variations at Low Latitude of African Equatorial Region Agbo G. A 1 ; Azi A. O. 2, Okoro N. O. 3 Industrial Physics Department, Ebonyi State University, P.M.B 053 Abakaliki Abstract:
More informationSaturation of polar cap potential: Nonlinearity in quasi steady solar wind magnetosphere ionosphere coupling
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009ja014389, 2010 Saturation of polar cap potential: Nonlinearity in quasi steady solar wind magnetosphere ionosphere coupling J. R. Kan, 1,2 H.
More informationGeomagnetically induced currents during intense storms driven by coronal mass ejections and corotating interacting regions
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007ja012487, 2008 Geomagnetically induced currents during intense storms driven by coronal mass ejections and corotating interacting regions Ryuho
More informationLow energy electrons in the inner Earth s magnetosphere
Low energy electrons in the inner Earth s magnetosphere Natalia Ganushkina (1, 2) (1) University of Michigan, Ann Arbor MI, USA (2) Finnish Meteorological Institute, Helsinki, Finland The research leading
More informationModeling Interactions between the Magnetosphere, Ionosphere & Thermosphere. M.Wiltberger NCAR/HAO
Modeling Interactions between the Magnetosphere, Ionosphere & Thermosphere M.Wiltberger NCAR/HAO Outline Overview of MIT circuit Modeling Magnetospheric impacts on the Ionosphere Energetic Particle Fluxes
More informationSpecification of electron radiation environment at GEO and MEO for surface charging estimates
Specification of electron radiation environment at GEO and MEO for surface charging estimates Natalia Ganushkina (University of Michigan/FMI) Collaborators: S. Dubyagin (FMI), J.-C. Matéo Vélez, A. Sicard
More informationSpace Weather and Satellite System Interaction
Space Engineering International Course, Kyutech, 4 th Quarter Semester 2017 Space Weather and Satellite System Interaction Lecture 2: Space Weather Concept, Reporting and Forecasting Assoc. Prof. Ir. Dr.
More informationHow changes in the tilt angle of the geomagnetic dipole affect the coupled magnetosphere-ionosphere-thermosphere system
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2012ja018056, 2012 How changes in the tilt angle of the geomagnetic dipole affect the coupled magnetosphere-ionosphere-thermosphere system Ingrid
More informationReal-time forecast service for geomagnetically induced currents
Real-time forecast service for geomagnetically induced currents ESA/ESTEC Contract No. 16953/2/NL/LvH WP 4 Technical Note Forecasting local geomagnetic field fluctuations and GIC from solar wind data P.
More informationTime Series of Images of the Auroral Substorm
ESS 7 Lecture 13 October 27, 2010 Substorms Time Series of Images of the Auroral Substorm This set of images in the ultra-violet from the Polar satellite shows changes that occur during an auroral substorm.
More informationSteady Magnetospheric Convection Selection Criteria: Implications of Global SuperDARN Convection Measurements
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, 1 2 3 Steady Magnetospheric Convection Selection Criteria: Implications of Global SuperDARN Convection Measurements K. A. McWilliams and J. B.
More informationSolar-Wind/Magnetosphere Coupling
Solar-Wind/Magnetosphere Coupling Joe Borovsky Space Science Institute --- University of Michigan 1. Get a feeling for how the coupling works 2. Get an understanding of how reconnection works 3. Look at
More informationReport from Finland. Kirsti Kauristie Finnish Meteorological Institute
Report from Finland Kirsti Kauristie Finnish Meteorological Institute IPT-SWeISS-2, Tokyo, Japan, 21-23 May 2018 Current SWE activities Contributors: Universities Aalto Helsinki (UH) Oulu (UO) Turku (UTU)
More informationEFFECT OF SOLAR AND INTERPLANETARY DISTURBANCES ON SPACE WEATHER
Indian J.Sci.Res.3(2) : 121-125, 2012 EFFECT OF SOLAR AND INTERPLANETARY DISTURBANCES ON SPACE WEATHER a1 b c SHAM SINGH, DIVYA SHRIVASTAVA AND A.P. MISHRA Department of Physics, A.P.S.University, Rewa,M.P.,
More informationResponse of the Earth s magnetosphere and ionosphere to the small-scale magnetic flux rope in solar wind by the MHD simulation
Response of the Earth s magnetosphere and ionosphere to the small-scale magnetic flux rope in solar wind by the MHD simulation Kyung Sun Park 1, Dae-Young Lee 1, Myeong Joon Kim 1, Rok Soon Kim 2, Kyungsuk
More informationWhat can I do with the TIEGCM?
What can I do with the TIEGCM? Astrid Maute and lots of people at HAO, and the community High Altitude Observatory NCAR High Altitude Observatory (HAO) National Center for Atmospheric Research (NCAR) The
More informationChapter 15. The Role and Form of Modelling in Space Weather
Chapter 15 The Role and Form of Modelling in Space Weather Konstantinos Papadopoulos Departments of Physics and Astronomy, University of Maryland College Park, MD 20782, USA Abstract Keywords A critical
More informationReal time Kp predictions from solar wind data using neural networks
MS No.: EGS.3-3 First author: Boberg Real time predictions from solar wind data using neural networks Fredrik Boberg, Peter Wintoft, and Henrik Lundstedt Lund Observatory, Box 3, SE- Lund, Sweden Swedish
More informationEvaluating the applicability of the finite element method for modelling of geoelectric fields
Ann. Geophys.,, 89 98, 0 www.ann-geophys.net//89/0/ doi:0.9/angeo--89-0 Author(s) 0. CC Attribution.0 License. Annales Geophysicae Open Access Evaluating the applicability of the finite element method
More informationSubstorms at Mercury: Old Questions and New Insights. Daniel N. Baker Laboratory for Atmospheric and Space Physics (LASP)
Substorms at Mercury: Old Questions and New Insights Daniel N. Baker Laboratory for Atmospheric and Space Physics (LASP) Outline of Presentation Introduction Substorms in the Earth s Magnetosphere Prior
More informationFirst European Space Weather Week. Space weather - atmospheres, drag, global change future needs. 29 November-3 December 2004
First European Space Weather Week Space weather - atmospheres, drag, global change future needs 29 November-3 December 2004 Timescales of important phenomena Weather Climate No single statement of requirement
More information2 Preliminary Results Achieved by the Meridian Project
Space Science Activities in China cycle peak year ( ), magnetic storm activities increased significantly, the Meridian Project has repeatedly observed the responses of the space environment to solar storms
More informationNumerical Methods Used in the Lyon-Fedder-Mobarry Global Code to Model the Magnetosphere
Advanced Methods for Space Simulations, edited by H. Usui and Y. Omura, pp. 101 109. c TERRAPUB, Tokyo, 2007. Numerical Methods Used in the Lyon-Fedder-Mobarry Global Code to Model the Magnetosphere John
More informationELECTROMAGNETIC SOURCE EQUIVALENCE AND EXTENSION OF THE COMPLEX IMAGE METHOD FOR GEOPHYSICAL APPLICATIONS
Progress In Electromagnetics Research B, Vol. 16, 57 84, 29 ELECTROMAGNETIC SOURCE EQUIVALENCE AND EXTENSION OF THE COMPLEX IMAGE METHOD FOR GEOPHYSICAL APPLICATIONS A. Pulkkinen University of Maryland,
More informationImpact of the altitudinal Joule heating distribution on the thermosphere
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010ja016019, 2011 Impact of the altitudinal Joule heating distribution on the thermosphere Yue Deng, 1 Timothy J. Fuller Rowell, 2,3 Rashid A. Akmaev,
More informationThis document is a preview generated by EVS
INTERNATIONAL STANDARD ISO 16698 First edition 2013-05-01 Space environment (natural and artificial) Methods for estimation of future geomagnetic activity Environnement spatial (naturel et artificiel)
More information[2] Observatorio Geofísico de Aiguá. Ruta 39 km 60 Maldonado Uruguay
in : The : The [1] Laboratorio de Geofísica-Geotectónica, Instituto de Ciencias Geológicas Universidad de la República, Facultad de Ciencias Iguá 4225 CP 11400 Montevideo, Uruguay [2] Observatorio Geofísico
More informationMagnetopause erosion: A global view from MHD simulation
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. A6, 1235, doi:10.1029/2002ja009564, 2003 Magnetopause erosion: A global view from MHD simulation M. Wiltberger High Altitude Observatory, National Center
More informationGlobal modeling of the magnetosphere in terms of paraboloid model of magnetospheric magnetic field
Global modeling of the magnetosphere in terms of paraboloid model of magnetospheric magnetic field I. Alexeev, V. Kalegaev The solar wind influence on the magnetospheric state is sufficiently nonlinear
More informationGeomagnetic storms. Measurement and forecasting
Geomagnetic storms. Measurement and forecasting Anna Gustavsson 17 October 2006 Project of the Space Physics Course 2006 Umeå University 1 Introduction Effects of magnetic storms on technology Geomagnetic
More informationInverse and normal coronal mass ejections: evolution up to 1 AU. E. Chané, B. Van der Holst, C. Jacobs, S. Poedts, and D.
A&A 447, 727 733 (2006) DOI: 10.1051/0004-6361:20053802 c ESO 2006 Astronomy & Astrophysics Inverse and normal coronal mass ejections: evolution up to 1 AU E. Chané, B. Van der Holst, C. Jacobs, S. Poedts,
More informationGeomagnetically induced pipe-to-soil voltages in the Czech oil pipelines during October November 2003
Annales Geophysicae, 3, 389 393, 5 SRef-ID: 43-576/ag/53-389 European Geosciences Union 5 Annales Geophysicae Geomagnetically induced pipe-to-soil voltages in the Czech oil pipelines during October November
More informationAn empirical model of ground-level geomagnetic perturbations
SPACE WEATHER, VOL. 11, 17, doi:1./swe.3, 13 An empirical model of ground-level geomagnetic perturbations D. R. Weimer 1 Received November ; revised 1 January 13; accepted 11 February 13; published March
More informationSpectral Studies of Ionospheric and Solar Wind Electric Fields: Case Studies and Long-Duration Analysis Using Magnetometer Data
Spectral Studies of Ionospheric and Solar Wind Electric Fields: Case Studies and Long-Duration Analysis Using Magnetometer Data M. C. Kelley, M. J. Nicolls, and G. Sachdeva Cornell University, Ithaca,
More informationWhy Study Magnetic Reconnection?
Why Study Magnetic Reconnection? Fundamental Process Sun: Solar flares, Flare loops, CMEs Interplanetary Space Planetary Magnetosphere: solar wind plasma entry, causes Aurora Ultimate goal of the project
More informationEnergy Analysis During the Collision of Two Successive CMEs
Numerical Modeling of Space Plasma Flows: ASTRONUM-2013 ASP Conference Series, Vol. 488 N.V.Pogorelov, E.Audit,and G.P.Zank,eds. c 2014 Astronomical Society of the Pacific Energy Analysis During the Collision
More informationPractical Model Applicable to Investigating the Coast. Effect on the Geoelectric Field in Connection with. Studies of Geomagnetically Induced Currents
Advances in Applied Physics, Vol. 1, 2013, no. 1, 9-28 HIKARI Ltd, www.m-hikari.com Practical Model Applicable to Investigating the Coast Effect on the Geoelectric Field in Connection with Studies of Geomagnetically
More informationLow energy electron radiation environment for extreme events
Low energy electron radiation environment for extreme events Natalia Ganushkina (1, 2) and Stepan Dubyagin (1) Special thanks to Jean-Charles Matéo-Vélez (3) (1) Finnish Meteorological Institute, Helsinki,
More informationSignatures of Geomagnetic Storms and Coronal Mass Ejections on Electron and Ion Temperatures At Low Latitude Upper Ionosphere
International Journal of Physics and Applications. ISSN 0974-3103 Volume 7, Number 1 (2015), pp. 43-48 International Research Publication House http://www.irphouse.com Signatures of Geomagnetic Storms
More informationSpace Weather Effects of Coronal Mass Ejection
J. Astrophys. Astr. (2006) 27, 219 226 Space Weather Effects of Coronal Mass Ejection K. N. Iyer 1,, R. M. Jadav 1, A. K. Jadeja 1, P. K. Manoharan 2, Som Sharma 3 and Hari Om Vats 3 1 Department of Physics,
More informationРС INDEX AS INDICATOR OF THE SOLAR WIND ENERGY ENTERED INTO THE MAGNETOSPHERE: RELATION TO INTERPLANETARY ELECTRIC FIELD AND MAGNETIC DISTURBANCES
РС INDEX AS INDICATOR OF THE SOLAR WIND ENERGY ENTERED INTO THE MAGNETOSPHERE: RELATION TO INTERPLANETARY ELECTRIC FIELD AND MAGNETIC DISTURBANCES О. A.Troshichev Arctic and Antarcrtic Research Institute,
More informationGeomagnetically Induced Currents in Austria. Rachel Bailey, ZAMG Georg Achleitner, APG Thomas Halbedl, TU Graz Roman Leonhardt, ZAMG
Geomagnetically Induced Currents in Austria Rachel Bailey, ZAMG Georg Achleitner, APG Thomas Halbedl, TU Graz Roman Leonhardt, ZAMG Geomagnetically Induced Currents (GIC) Enhanced electrojet and rapid
More informationRationale for a European Space Weather Programme
Rationale for a European Space Weather Programme Hannu Koskinen Finnish Meteorological Institute ESWS Final Presentation ESTEC, 6 December, 2001 Scope WP 300 of ESWS: Establishment of detailed rationale
More informationA global MHD simulation of an event with a quasi-steady northward IMF component
European Geosciences Union 2007 Annales Geophysicae A global MHD simulation of an event with a quasi-steady northward IMF component V. G. Merkin 1, J. G. Lyon 2, B. J. Anderson 3, H. Korth 3, C. C. Goodrich
More informationNatalia Ganushkina (1, 2), Stepan Dubyagin (1), Ilkka Sillanpää (1)
From studying electron motion in the electromagnetic fields in the inner magnetosphere to the operational nowcast model for low energy (< 200 kev) electron fluxes responsible for surface charging Natalia
More informationBenchmark GMD Event. Project (GMD Mitigation) Standard Drafting Team. GMD Task Force In-person meeting March 18-19, 2014
Benchmark GMD Event Project 2013-03 (GMD Mitigation) Standard Drafting Team GMD Task Force In-person meeting March 18-19, 2014 Topics General Description Luis Marti, Hydro One Statistical Considerations
More informationObservations of geomagnetically induced currents in the Australian power network
SPACE WEATHER, VOL. 11, 6 16, doi:10.1029/2012sw000849, 2013 Observations of geomagnetically induced currents in the Australian power network R. A. Marshall, 1 H. Gorniak, 2 T. Van Der Walt, 2 C. L. Waters,
More informationChaos and the limits of predictability for the solar-wind-driven magnetosphere ionosphere system
PHYSICS OF PLASMAS VOLUME 8, NUMBER 6 JUNE 2001 Chaos and the limits of predictability for the solar-wind-driven magnetosphere ionosphere system W. Horton and R. S. Weigel Department of Physics and Institute
More informationJoule heating and nitric oxide in the thermosphere, 2
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010ja015565, 2010 Joule heating and nitric oxide in the thermosphere, 2 Charles A. Barth 1 Received 14 April 2010; revised 24 June 2010; accepted
More information