Search USGS⌕ Search

Geology topics

E. Joshua Rigler

Publications and source records attributed to E. Joshua Rigler.

At least 19 recordsLinked to original sources

Near-Earth solar wind speed of fast coronal mass ejections

Using solar wind and ground magnetometer data for solar cycles 23–25, extreme-value power-law models are developed that quantify relationships between near-Earth solar wind speed and Sun-to-Earth transit times of interplanetary coronal mass ejections (ICMEs). The models of near-Earth solar wind speed, conditional on ICME transit time, are (a) approximately stable for transit times as short as 15 hr, (b) sublinear, consistent with hydrodynamic drag acting on ICMEs during transit, with linear (ballistic) models not supported by the data, (c) representative of most of the information content of the data, and (d) insensitive to the solar cycle and interplanetary preconditioning. Applying the models to the September 1859 Carrington event ( 𝑇 =17.6 hours), we estimate an ICME-average speed of 𝑉 𝑎 =102⁢1 1164 896 km/s and a 1-hr ICME-maximum speed of 𝑉 𝑚 =157⁢5 1795 1382 km/s. These estimates indicate that Carrington-class magnetic storms do not require exceptionally extreme solar wind speeds, with values lower than some previous estimates.

JGR Space Physics↗

Interplanetary electric fields for extreme magnetic storms

Using a list of sudden-commencement storms, the ring-current index, and 1-h near-Earth solar-wind measurements from solar cycles 20–25, we develop extreme-value statistical models relating storm intensity 𝐷 =max⁡{−𝐷⁢𝑠⁢𝑡} to the storm main-phase maximum duskward interplanetary electric field 𝐸 . The conditional relationship 𝐷|𝐸 is demonstrably sublinear—linear models are confidently rejected—indicating saturation of magnetospheric response under extreme solar-wind forcing. An event like that of July 2012 ( 𝐸 =69.6 mV/m), if Earth-directed, would be associated with a median storm intensity of 𝐷 =49⁢5 648 378 nT. Storms comparable to March 1989 ( 𝐷 =594 nT) correspond to electric fields of 𝐸 =5⁢4 76 39 mV/m, while Carrington-class storms ( 𝐷 =964 nT) correspond to 𝐸 =9⁢5 133 68 mV/m—substantially lower than several previous estimates. These results indicate that solar-wind conditions capable of driving extremely intense magnetic storms are less exceptional, and potentially more frequent, than previously thought.

Geophysical Research Letters↗

Magnetic storms and geoelectric hazards

Magnetic storms induce geoelectric fields at Earth's surface that can interfere with grounded long-line systems. The September 1859 storm disrupted global telegraph operations, the March 1989 storm caused a blackout in Canada and interfered with electric-power-transmission systems in the United States, and other storms have had related impacts. The geographic and temporal dependence of geoelectric fields are functions of both geomagnetic variation and local surface impedance, which differ considerably across different geological regions. These dependencies can be mapped across the contiguous United States by combining magnetotelluric impedance tensors with ground magnetometer time series. This review illustrates such mapping for the 1989 storm and shows that power-system interference was experienced where surface impedance is high, and when and where geoelectric fields were intense. Statistical analyses indicate that storms comparable to that of March 1989 occur roughly once every four solar cycles. Ongoing developments in numerical modeling and real-time monitoring are anticipated to enable prediction of geoelectric hazards. ▪ Magnetic storms can induced electric fields in the solid Earth that interfere with electric-power-transmission systems. ▪ Geoelectric hazards depend on the storm-time geomagnetic disturbance and the electrical conductivity structure of Earth. ▪ Historically, impacts on telecommunication and power-transmission systems in the United States have been concentrated in the East and Midwest. ▪ The future occurrence of a magnetic superstorm could cause widespread disruption of electric-power-transmission systems.

Annual Review of Earth and Planetary Sciences↗

U.S. Geological Survey geomagnetic variometer data: Capitalizing on seismic infrastructure

The U.S. Geological Survey’s Geomagnetism Program is collaborating with the Earthquake Hazards Program and Global Seismographic Network Program to densify magnetic field observations. This collaboration focuses on the installation of magnetometers, or magnetic variometers, at existing seismic stations. Along with improving the density of space weather observations for hazard monitoring, these data can be used to correct colocated magnetic field induced noise in seismic data. Such corrections are especially useful during time periods of large magnetic storms where the magnetic field‐induced instrument noise can be of similar amplitude to earthquake ground‐motion records.

contiguous United States↗

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗

Multi-scale intense geoelectric and geomagnetic field perturbations observed after an interplanetary magnetic field turning

Intense geoelectric fields during geomagnetic storms generate geomagnetically induced currents in power grids and other infrastructure, necessitating an understanding of their causes, for example, through coordinated space and ground observations. This study investigates localized intense geoelectric ( E ) and geomagnetic ( B ) field perturbations following an Interplanetary Magnetic Field (IMF) turning during a geomagnetic storm on 25 October 2011. Observations from EarthScope magnetotelluric sites in the upper Midwest United States revealed shorter period ( ∼ 1 min) ultra-low-frequency (ULF) waves superimposed on longer period ( ∼ 10 min) perturbations in both E and B fields. These sites, located at ∼ 19 hr magnetic local time and 56 −57 ° magnetic latitude, recorded large amplitude E and B perturbations. Ground-based all-sky imagers showed auroral brightening with sunward and poleward propagation, while upstream spacecraft linked the perturbations to an IMF turning and solar wind dynamic pressure impulse. The longer-period E and B field perturbations likely stem from localized ionospheric currents tied to substorm auroral activity post-IMF turning. The combination of ionospheric currents, ULF waves, and the Earth's varying conductivity produces intense geoelectric fields of ≥ 2 V/km in the upper Midwest. A comparison using input data and software compatible with the NOAA/USGS geoelectric field nowcast model revealed its limitations in capturing such events due to the temporal and spatial resolution of the underlying data. Using 1-s geomagnetic field data can improve geoelectric field models by capturing short-period and large spatial scale waves, although localized magnetic perturbations remain underestimated due to insufficient ground magnetometer density.

Space Weather↗

The solar cycle, geology, and geoelectric hazards for power grids

When sunspots are large and numerous, intense magnetic storms are likely to occur on the Earth. Magnetic storms can generate electric fields in the Earth, and these fields can, in turn, interfere with electric power transmission grids that are grounded at the Earth’s surface. Across the contiguous United States, geoelectric hazards are highest in the Upper Midwest and in the East. These regions correspond to geological structures that are electrically resistive, and they have, historically, experienced the most interference to electric power systems.

Fact Sheet↗

On the uncertain intensity estimate of the 1859 Carrington storm

A study is made of the intensity of the Carrington magnetic storm of September 1859 as inferred from visual measurements of horizontal-component geomagnetic disturbance made at the Colaba observatory in India. Using data from modern observatories, a lognormal statistical model of storm intensity is developed, to characterize the maximum-negative value of the storm-time disturbance index (maximum – Dst) versus geomagnetic disturbance recorded at low-latitude observatories during magnetic storms. With this model and a recently published presentation of the Colaba data, the most likely maximum – Dst of the Carrington storm and its credibility interval are estimated. A related model is used to examine individual Colaba disturbance values reported for the Carrington storm. Results indicate that only about one in a million storms with maximum – Dst like the Carrington storm would result in local disturbance greater than that reported from Colaba. This indicates that either the Colaba data were affected by magnetospheric-ionospheric current systems in addition to the ring current, or there might be something wrong with the Colaba data. If the most extreme Colaba disturbance value is included in the analysis, then, of all hypothetical storms generating the hourly average disturbance recorded at Colaba during the Carrington storm, the median maximum – Dst = 964 nT, with a 68% credibility interval of [855,1087] nT. If the most extreme Colaba disturbance value is excluded from the analysis, then the median maximum – Dst = 866 nT, with a 68% credibility interval of [768,977] nT. The widths of these intervals indicate that estimates of the occurrence frequency of Carrington-class storms are very uncertain, as are related estimates of risk for modern technological systems.

Journal of Space Weather and Space Climate↗

Correction to “The March 1940 superstorm: Geoelectromagnetic hazards and impacts on American communication and power systems”

In Love et al. (2023), 5 estimated geopotentials on long lines for the March 1989 storm are incorrect. This was the result of incorrect calculation of line integrals. This error affects the two paragraphs of Section 11 and Figure 10 of that paper. The most significant impact of this error is that we overestimated the average ratio of the geopotentials measured the long lines during the March 1940 storm to the estimated geopotentials for the March 1989 storm. Instead of being a factor of 9 higher, the average ratio is only 9% higher. This means the inductive geoeffectiveness of the March 1940 storm was comparable to, or just slightly greater than, the geoeffectiveness of the March 1989 storm. Corrected paragraphs and a corrected figure are provided here.

Space Weather↗

Determining ULF wave contributions to geomagnetically induced currents: The important role of sampling rate

Past studies found that large-amplitude geomagnetically induced current (GIC) related to magnetospheric Ultra Low Frequency (ULF) waves tend to be associated with periods >120 s at magnetic latitudes >60°, with comparatively (a) smaller GIC amplitudes at lower latitudes and shorter wave periods and (b) fewer reports of waves associated with GIC at lower latitudes. ULF wave periods generally decrease with decreasing latitude; thus, we examine whether these trends might be due, in part, to the undersampling of ULF wave fields in commonly available measurements with 60 s sampling intervals. We use geomagnetic field ( B ), geoelectric field ( E ), and GIC measurements with 0.5–10 s sampling intervals during the 29–31 October 2003 geomagnetic storm to show that waves with periods <∼120 s were present during times with the largest amplitude E and GIC variations. These waves contributed to roughly half the maximum E and GIC values, including during times with the maximum GIC values reported over a 14-year monitoring interval in New Zealand. The undersampling of wave periods <120 s in 60 s measurements can preclude identification of the cause of the GIC during some time intervals. These results indicate (a) ULF waves with periods ≤120 s are an important contributor to large amplitude GIC variations, (b) the use of 0.1–1.0 Hz sampling rates reveals their contributions to B , E , and GIC, and (c) these waves' contributions are likely strongest at magnetic latitudes <60° where ULF waves often have periods <120 s.

Space Weather↗

Mapping a magnetic superstorm: March 1989 geoelectric hazards and impacts on United States power systems

A study is made of the relationships between geomagnetic and geoelectric field variation, Earth-surface impedance, and operational interference (anomalies) experienced on electric-power systems across the contiguous United States during the March 13-14, 1989 magnetic storm. For this, a 1-minute-resolution sequence of geomagnetic field maps is constructed from magnetometer time series acquired at ground-based observatories. Induced geoelectric field maps are calculated by convolving the geomagnetic maps with magnetotelluric impedance tensors. During the storm, anomalies were concentrated where the lithosphere is electrically resistive, and when and where geoelectric field amplitudes were high. This was particularly true in the Mid-Atlantic, Northeast, and the upper Midwest. Few anomalies were experienced in other parts of the Midwest and across much of the West, where the lithosphere is more conductive, and when and where geoelectric field amplitudes were low. Peak 1-minute-resolution geoelectric field amplitude ranged from 21.66 V/km in Maine and 19.02 V/km in Virginia to < 0.02 V/km in Idaho. Latitude-dependent organization of geoelectric hazards by auroral-zone electrojet currents is detectable, but it is much weaker than geographic organization due to surface impedance. Hazardous geoelectric fields were induced during different storm phases, at different local times, and, by inference, by a variety of ionospheric currents. Compared to geoelectric field amplitudes realized across the United States during March 1989, hazard maps used by utility companies to estimate systems exposure have much less geographic detail and a much smaller maximum-to-minimum range in geoelectric field amplitude. Future research will benet from denser geomagnetic monitoring, additional magnetotelluric surveying, and access to power-system impact data.

Space Weather↗

Synthesizing ground magnetic disturbance using dipole-aligned loop elementary currents and Biot-Savart relationship

This report presents a method for constructing a simplified numerical description of the electric current distributions in the ionosphere and gap region based on dipole-aligned loop elementary currents (DALECs). A theoretical basis for DALECs is presented, along with a prototypical algorithm for constructing an elementary numerical DALEC. The algorithm is verified and validated by combining DALECs with an efficient Biot-Savart solver in order to estimate magnetic disturbance on the Earth’s surface. We examine (1) simple scenarios with known solutions and (2) hemispherical magnetic disturbance fields obtained from a state-of-the-art global geospace circulation model.

Scientific Investigations Report↗

Down to Earth with nuclear electromagnetic pulse: Realistic surface impedance affects mapping of the E3 geoelectric hazard

An analysis is made of Earth-surface geoelectric fields and voltages on electricity transmission power-grids induced by a late-phase E3 nuclear electromagnetic pulse (EMP). A hypothetical scenario is considered of an explosion of several hundred kilotons set several hundred kilometers above the eastern-midcontinental United States. Ground-level E3 geoelectric fields are estimated by convolving a standard parameterization of E3 geomagnetic field variation with magnetotelluric Earth-surface impedance tensors derived from wideband measurements acquired across the study region during a recent survey. These impedance tensors are a function of subsurface three-dimensional electrical conductivity structure. Results, presented as a movie-map, demonstrate that localized differences in surface impedance strongly distort the amplitude, polarization, and variational phase of induced E3 geoelectric fields. Locations with a high degree of E3 geoelectric polarization tend to have high geoelectric amplitude. Uniform half-space models and one-dimensional, depth-dependent models of Earth-surface impedance, such as those widely used in government and industry reports informing power-grid vulnerability assessment projects, do not provide accurate estimates of the E3 geoelectric hazard in complex geological settings. In particular, for the Eastern-Midcontinent, half-space models can lead to (order-one) overestimates/underestimates of EMP-induced geovoltages on parts of the power grid by as much as 1,000 volts (a range of 2,000 volts)—comparable to the amplitudes of the geovoltages themselves.

Earth and Space Sciences↗

Magnetotelluric sampling and geoelectric hazard estimation: Are national-scale surveys sufficient?

At present, the most reliable information for inferring storm-time ground electric fields along electrical transmission lines comes from coarsely sampled, national-scale magnetotelluric (MT) data sets, such as that provided by the EarthScope USArray program. An underlying assumption in the use of such data is that they adequately sample the spatial heterogeneity of the surface relationship between geomagnetic and geoelectric fields. Here, we assess the degree to which the density of MT data sampling affects geoelectric hazard assessments. For electrical transmission networks in each of four focus regions across the contiguous United States, we perform two parallel band-limited (10 1 –10 3 s) hazard analyses: one using only USArray-style (∼70-km station spacing) MT data, and one incorporating denser (≪70-km station spacing) MT data. We find that the use of USArray-style MT sampling alone provides a useful first-order estimate of integrated geoelectric fields along electrical transmission lines. However, we also find that the use of higher density MT data can in some areas lead to order-of-magnitude differences in line-averaged electric field estimates at the level of individual transmission lines and can also yield significant differences in subregional hazard patterns. As we demonstrate using variogram plots, these differences reflect short-spatial-scale variability in Earth conductivity, which in turn reflects regional lithotectonic structure and history. We also provide a cautionary example in the use of electrical conductivity models to predict dense MT data; although valuable for hazard applications, models may only be able to reproduce surface geoelectric fields as captured by the MT data from which they were derived.

AGU Space Weather↗

Numerical simulations of the geospace response to the arrival of an idealized perfect interplanetary coronal mass ejection

Previously, Tsurutani and Lakhina (2014, https://doi.org/10.1002/2013GL058825 ) created estimates for a “perfect” interplanetary coronal mass ejection and performed simple calculations for the response of geospace, including . In this study, these estimates are used to drive a coupled magnetohydrodynamic-ring current-ionosphere model of geospace to obtain more physically accurate estimates of the geospace response to such an event. The sudden impulse phase is examined and compared to the estimations of Tsurutani and Lakhina (2014, https://doi.org/10.1002/2013GL058825 ). The physics-based simulation yields similar estimates for Dst rise, magnetopause compression, and equatorial values as the previous study. However, results diverge away from the equator. values in excess of 30 nT/s are found as low as magnetic latitude. Under southward interplanetary magnetic field conditions, magnetopause erosion combines with strong region one Birkeland currents to intensify the response. Values obtained here surpass those found in historically recorded events and set the upper threshold of extreme geomagnetically induced current activity at Earth.

Space Weather↗

Simultaneous observations of geoelectric and geomagnetic fields produced by magnetospheric ULF waves

Geomagnetic perturbations ( B G E O ) related to magnetospheric ultralow frequency (ULF) waves induce electric fields within the conductive Earth—geoelectric fields ( E G E O )—that in turn drive geomagnetically induced currents. Though numerous past studies have examined ULF wave B G E O from a space weather perspective, few studies have linked ULF waves with E G E O . Using recently available magnetotelluric impedance and E G E O measurements in the contiguous United States, we explore the relationship between ULF waves and E G E O . We use satellite, ground‐based radar, B G E O , and E G E O measurements in a case study of a plasmaspheric virtual resonance (PVR), demonstrating that the PVR E G E O has significant spatial variation in contrast to a relatively uniform B G E O , consistent with spatially varying Earth conductivity. We further show ULF wave E G E O measurements during two moderate storms of ∼ 1 V/km. We use both results to highlight the need for more research characterizing ULF wave E G E O .

Geophysical Research Letters↗

Geomagnetism Program research plan, 2020–2024

The Geomagnetism Program of the U.S. Geological Survey (USGS) monitors geomagnetic field variation through operation of a network of observatories across the United States and its territories, and it pursues scientific research needed to estimate and assess geomagnetic and geoelectric hazards. Over the next five years (2020–2024 inclusive) and in support of national and agency priorities, Geomagnetism Program research scientists plan to pursue an integrated set of research projects broadly encompassing empirical estimation and mapping of geomagnetic disturbance, modeling of solid-Earth conductivity structure and surface impedance, and mapping of magnetic-storm-induced geoelectric fields. Analyses are empirically based, relying on measured time series as well as statistical and numerical modeling of geomagnetic-monitoring data from ground-based observatories and surface-impedance tensors acquired during magnetotelluric surveys. The plan describes augmentation and development of the Geomagnetism Program's existing research portfolio, assuming present funding levels and staffing numbers. Because the projects are interdependent, they cannot be straightforwardly prioritized. They will all be pursued as resources and time permit; additional funding and staffing would enable the projects to be broadened and more rapidly completed. Where appropriate and subject to budgetary constraints and staffing numbers, research on specific projects might be accelerated or even judiciously expanded—some opportunities for expansion are discussed in this plan. Results will provide realistic illumination of the nature of the ground-level expression of space-weather disturbance, a subject of particular importance for projects focused on evaluating the vulnerability of electric-power-grid systems. This plan does not cover Geomagnetism Program operations, which are primarily concerned with the operation of magnetic observatories and, now, magnetotelluric surveys, although the context of such observatories and surveys is discussed. The research element of the program provides guidance for the expansion of program operations and research projects. In addition to the research projects summarized here, program scientists continue to provide leadership to the national and international geomagnetic, magnetotelluric, and space-weather communities.

Circular↗

Adjusted geomagnetic data—Theoretical basis and validation

Adjusted geomagnetic data are magnetometer measurements with provisional correction factors applied such that vector quantities are oriented in a local Cartesian frame in which the X axis points north, the Y axis points east, and the Z axis points down. These correction factors are determined from so-called absolute measurements, which are “ground truth” observations made in the field using specialized magnetometers and survey equipment that are (nearly) colocated with the automated and continuously running magnetic measurement instrumentation. Correction factors can be substantial, up to hundreds of nanoTeslas, depending on the geologic and geomagnetic characteristics of the observatory site. They also tend to evolve over time because of instrument response instability and changing site characteristics. Historically, correction factors were determined offline, up to 1 year or more post-measurement, and applied to raw measurements to produce “Definitive” data for scientific analysis. Growing demand for corrected real-time geomagnetic data to better support space weather operations motivated development of an “Adjusted” geomagnetic data product. Modern computational tools, and some notable practical concerns, dictated a transition to affine transformations in lieu of more traditional baseline corrections, as well as a calibration parameter estimation algorithm that is more robust and statistically optimal, and therefore better suited for automated and unsupervised execution. A theoretical basis for this algorithm is presented, along with a demonstration and validation based on a comparison of results obtained with traditional techniques. Discrepancies between Definitive corrected data and near real-time Adjusted data obtained using affine transformations are minimal, generally much less than 5 nanoTeslas per vector component, and less than 1 nanoTesla for the total field magnitude, which satisfies International Real-Time Magnetic Observatory Network (INTERMAGNET) standards.

Open-File Report↗