Geology topics

Anna Kelbert

Publications and source records attributed to Anna Kelbert.

34 records · Page 2Linked to original sources

USGS2020

Simultaneous observations of geoelectric and geomagnetic ﬁelds 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↗

USGS2020

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↗

USGS2020

Modified GIC estimation using 3-D Earth conductivity

Geomagnetically induced currents (GICs) are quasi-direct current (DC) electric currents that flow in technological conductors during geomagnetic storms. Extreme GICs are hazardous to man-made infrastructure. GICs enter and exit the technological systems, such as the electric power grid, at grounding points, and their magnitudes depend on the currents that flow underground. They are, therefore, a function of the Earth's electrical conductivity, represented at ground level as Earth impedances, as well as the resistance parameters of the power network. Traditional GIC estimation practices are based on Earth impedances obtained from laterally homogeneous or piecewise layered-Earth models. We refer to these methods, collectively, as the 1-D approximation. However, GIC hazard mitigation can be improved with more accurate GIC modeling that takes the spatially heterogeneous Earth's conductivity into account. Here, we propose a modified approximation for GIC estimation that is very similar to the 1-D approximation but is instead derived from empirical 3-D Earth impedances. Our formulation sets up the computation of static, frequency-dependent power line telluric response functions, which, once computed, may be considered part of the power grid system model. These response functions may then be used for historical scenario analysis of GIC hazards and for simplified real-time, albeit approximate, GIC estimation in a power grid. This modest modification to the simpler local field formulation approach avoids real-time integration of geoelectric fields along power lines while taking the realistic 3-D Earth into account in a rigorous manner. Once implemented, the method provides a power grid operator with the benefits of convenience and computational speed for a first look real-time operational GIC hazard assessment. We estimate that the proposed modified 3-D GIC modeling approach produces GIC values that are well within 50% of those obtained with the full-scale power line integration of spatially variable geoelectric fields, for storms comparable in scale to the 2003 Halloween storm, all geological structures, and power lines located in the contiguous United States and other low- to middle-latitude regions.

Space Weather↗

USGS2020

A 100-year geoelectric hazard analysis for the U.S. high-voltage power grid

A once-per-century geoelectric hazard map is created for the United States high-voltage power grid. A statistical extrapolation from 31 years of magnetic field measurements is made by identifying 84 geomagnetic storms with the Kp and Dst indices. Data from 24 geomagnetic observatories, 1079 magnetotelluric survey sites, and 17,258 transmission lines are utilized to perform a geoelectric hazard analysis with the most comprehensive data publicly available. With this data we estimate once-per-century geoelectric fields at the magnetotelluric survey sites and calculate the theoretical voltages within transmission lines in the United States power grid. Once-per-century geoelectric field strengths span more than three orders of magnitude from a minimum of 0.02 V/km at a site in Idaho to a maximum of 26.8 V/km at a site in Maine, with nearly 30% of the surveyed land area ex- ceeding 1 V/km. We show the influence that geoelectric field polarization has on geoelec- tric hazards when viewed on a power transmission network. The calculated transmission line voltages can exceed 1000 V in some transmission lines. Four regions in the United States with particularly notable geoelectric hazards are identified and discussed: the East coast, Pacific Northwest, upper Midwest, and the Denver metropolitan area.

Space Weather↗

USGS2019

Empirical estimation of natural geoelectric hazards

Geoelectric field time series can be estimated by convolving estimates of Earth‐surface impedance, such as those obtained from magnetotelluric survey measurements, with historical records of geomagnetic variation obtained at magnetic observatories. This straightforward procedure permits the mapping of geoelectric field variation during magnetic storms. Statistical analysis of the time series allows extrapolation to extreme‐value amplitudes, such as might be realized during an intense magnetic storm in the future. The development of these products is illustrated for the Mid‐Atlantic United States, using impedances obtained from EarthScope survey data and geomagnetic variation records obtained at the Fredericksburg observatory operated by the U.S. Geological Survey. For this region, 100‐year geoelectric exceedance amplitudes have a range of almost three orders of magnitude (from 0.04 V/km at a site in southern Pennsylvania to 24.29 V/km at a site in central Virginia), and they have significant geographic granularity, which is due to site‐to‐site differences in surface impedance (and subsurface electrical conductivity structure). Maps of 100‐year exceedance amplitudes resemble those of geoelectric amplitudes for the March 1989 magnetic storm, and, in that sense, the March 1989 storm resembles what might be loosely called a “100‐year” event.

Book chapter↗

USGS2019

The first 3D conductivity model of the contiguous US: Reflections on geologic structure and application to induction hazards

Estimation of ground level geoelectric fields has been identified by the National Space Weather Action Plan as a key component of assessment and mitigation of space weather impacts on critical infrastructure. Estimates of spatially and temporally variable electric fields are used to generate statistically based hazard maps and show promise toward monitoring and responding to geomagnetic disturbances in near real‐time. One approach to geoelectric field estimation is to employ three‐dimensional (3D) Earth conductivity models. These data‐constrained conductivity models are the results of regional magnetotelluric inversions based primarily on NSF’s Earthscope USArray impedances, which to date cover ~60% of the contiguous United States. Here, we present the first‐ever composite conductivity model of the contiguous United States and describe its compilation from 3D regional conductivity models, a global mantle conductivity model, offshore bathymetry, and sediment thickness data. We discuss structures within the conductivity model and how they relate to the complex geologic tapestry of the continent. Finally, we discuss the utility of this synthesis model for estimation and mitigation of geomagnetically induced currents.

Book chapter↗

USGS2019

Data sharing in magnetotellurics

Here, we introduce the first openly available comprehensive database of magnetotelluric (MT) and related electromagnetic data that we developed and matured over the past decade, explain how to access the data, and describe the challenges that had to be overcome to make MT data sharing possible. The database is a helpful tool for MT scientists, and is widely used by the international scientific community in teaching, planning future surveys, and then to leverage the value of new data at the inversion and interpretation stage. It has also turned into a primary resource for researchers in space weather, and for power grid industry. MT impedance, a transfer function that relates geoelectric fields to geomagnetic fields at the Earth’s surface, provides a critical missing component for estimation of geomagnetically induced currents: quasi-direct currents caused by geomagnetic storms that can destroy valuable equipment and cause extended blackouts. We envision that this database will also be useful for researchers in other areas of geophysics, such as seismology, who would like to use real MT data for hands-on learning and joint inversion analysis.

Eos, Earth and Space Science News↗

USGS2018

Geoelectric hazard maps for the Pacific Northwest

Maps of extreme value, horizontal component geoelectric field amplitude are constructed for the Pacific Northwest United States (and parts of neighboring Canada). Multidecade long geoelectric field time series are calculated by convolving Earth surface impedance tensors from 71 discrete magnetotelluric survey sites across the region with historical 1‐min (2‐min Nyquist) geomagnetic variation time series obtained from two nearby observatories. After fitting statistical models to 1‐min geoelectric amplitudes realized during magnetic storms, extrapolations are made to estimate threshold amplitudes that are only exceeded, on average, once per century. One hundred‐year geoelectric exceedance amplitudes range from 0.06 V/km at a survey site in western Washington State to 9.47 V/km at a site in southeast British Columbia; 100‐year geoelectric exceedance amplitudes equal 7.10 V/km at a site north of Seattle and 2.28 V/km at a site north of Portland. Systematic and random errors are estimated to be less than 20%, much less than site‐to‐site differences in geoelectric amplitude that arise from site‐to‐site differences in surface impedance. Maps of 100‐year exceedance amplitudes are compared with the peak geoelectric amplitudes realized during the March 1989 magnetic superstorm; it is noted that some storms of relatively modest intensity can generate localized geoelectric fields of relatively high amplitude. The geography of geoelectric hazard across the Pacific Northwest is closely related to known geologic and tectonic structures.

Pacific Northwest↗

USGS2018

On the feasibility of real-time mapping of the geoelectric field across North America

A review is given of the present feasibility for accurately mapping geoelectric fields across North America in near-realtime by modeling geomagnetic monitoring and magnetotelluric survey data. Should this capability be successfully developed, it could inform utility companies of magnetic-storm interference on electric-power-grid systems. That real-time mapping of geoelectric fields is a challenge is reflective of (1) the spatiotemporal complexity of geomagnetic variation, especially during magnetic storms, (2) the sparse distribution of ground-based geomagnetic monitoring stations that report data in realtime, (3) the spatial complexity of three-dimensional solid-Earth impedance, and (4) the geographically incomplete state of continental-scale magnetotelluric surveys.

Open-File Report↗

USGS2018

Calculation of voltages in electric power transmission lines during historic geomagnetic storms: An investigation using realistic earth impedances

Commonly, one-dimensional (1-D) Earth impedances have been used to calculate the voltages induced across electric power transmission lines during geomagnetic storms under the assumption that much of the three-dimensional structure of the Earth gets smoothed when integrating along power transmission lines. We calculate the voltage across power transmission lines in the mid-Atlantic region with both regional 1-D impedances and 64 empirical 3-D impedances obtained from a magnetotelluric survey. The use of 3-D impedances produces substantially more spatial variance in the calculated voltages, with the voltages being more than an order of magnitude different, both higher and lower, than the voltages calculated utilizing regional 1-D impedances. During the March 1989 geomagnetic storm 62 transmission lines exceed 100 V when utilizing empirical 3-D impedances, whereas 16 transmission lines exceed 100 V when utilizing regional 1-D impedances. This demonstrates the importance of using realistic impedances to understand and quantify the impact that a geomagnetic storm has on power grids.

Space Weather↗

USGS2018

Geoelectric hazard maps for the Mid-Atlantic United States: 100 year extreme values and the 1989 magnetic storm

Maps of extreme value geoelectric field amplitude are constructed for the Mid‐Atlantic United States, a region with high population density and critically important power grid infrastructure. Geoelectric field time series for the years 1983–2014 are estimated by convolving Earth surface impedances obtained from 61 magnetotelluric survey sites across the Mid‐Atlantic with historical 1 min (2 min Nyquist) measurements of geomagnetic variation obtained from a nearby observatory. Statistical models are fitted to the maximum geoelectric amplitudes occurring during magnetic storms, and extrapolations made to estimate threshold amplitudes only exceeded, on average, once per century. For the Mid‐Atlantic region, 100 year geoelectric exceedance amplitudes have a range of almost 3 orders of magnitude (from 0.04 V/km at a site in southern Pennsylvania to 24.29 V/km at a site in central Virginia), and they have significant geographic granularity, all of which is due to site‐to‐site differences in magnetotelluric impedance. Maps of these 100 year exceedance amplitudes resemble those of the estimated geoelectric amplitudes attained during the March 1989 magnetic storm, and, in that sense, the March 1989 storm resembles what might be loosely called a “100 year” event. The geoelectric hazard maps reported here stand in stark contrast with the 100 year geoelectric benchmarks developed for the North American Electric Reliability Corporation.

Geophysical Research Letters↗

USGS2018

Extreme‐value geoelectric amplitude and polarization across the northeast United States

Maps are presented of extreme‐value geoelectric field amplitude and horizontal polarization for the Northeast United States. These maps are derived from geoelectric time series calculated for sites across the Northeast by frequency‐domain multiplication (time‐domain convolution) of 172 magnetotelluric impedance tensors, acquired during a survey, with decades‐long, 1‐min resolution time series of geomagnetic variation, acquired at three magnetic observatories. The maps show that, during intense magnetic storms, high geoelectric amplitude hazards are realized across electrically resistive, igneous and metamorphic rock of the Appalachian Mountains and the New England Highlands, while low geoelectric hazards are realized across electrically conductive, sedimentary rock of the Appalachian Plateau and the Mid‐Atlantic Coastal Plain. From statistical extrapolation, once‐per‐century (100‐year) geoelectric amplitudes are highest at a site in Virginia at 25.44 V/km (followed by a site in Maine at 21.75 V/km and a site in Connecticut at 19.39 V/km); 100‐year geoelectric amplitude exceeds 10 V/km at 12 sites across the northeast; geoelectric amplitude is lowest at a site in Virginia at 0.05 V/km. Average errors for these values are estimated to be about 38%, or much less than the more than 2 orders of magnitude range seen in geoelectric amplitudes from one survey site to another across the northeast. It is noteworthy that geoelectric fields tend to be most (least) polarized at locations with high (low) geoelectric hazard. Furthermore, geoelectric fields over the Appalachians tend to be polarized southeast‐to‐northwest, or generally in a direction orthogonal to the southwest‐to‐northeast geological strike. Results reported here inform utility companies in projects for evaluating and managing the response of power grid systems to the deleterious effects of geomagnetic disturbance.

Connecticut, Delaware, Maine, Maryland, Massachuse↗

USGS2017

Methodology for time-domain estimation of storm time geoelectric fields using the 3-D magnetotelluric response tensors

Geoelectric fields at the Earth's surface caused by magnetic storms constitute a hazard to the operation of electric power grids and related infrastructure. The ability to estimate these geoelectric fields in close to real time and provide local predictions would better equip the industry to mitigate negative impacts on their operations. Here we report progress toward this goal: development of robust algorithms that convolve a magnetic storm time series with a frequency domain impedance for a realistic three-dimensional (3-D) Earth, to estimate the local, storm time geoelectric field. Both frequency domain and time domain approaches are presented and validated against storm time geoelectric field data measured in Japan. The methods are then compared in the context of a real-time application.

Space Weather↗

USGS2017

Some results from ModEM3DMT, the freely available OSU 3D MT inversion code

At the 3DEM-5 workshop in 2013, we presented a paper entitled "ModEM: developing 3D EM inversion for the masses", outlining our then recent development of a modular system for inversion of EM geophysical data, called ModEM. As promised in that presentation, we made a version of the code that is suitable for 3D modeling and inversion of magnetotelluric data freely available for academic use shortly thereafter. There are now over 250 registered users, of ModEM3DMT from around the globe. To date at least 50 publications cite use of ModEM for 3D inversion of real MT datasets to address diverse problems in applied and basic Earth Science research at a range of scales. Here we present an overview of some of these results, focusing on studies that the authors have been involved in, and are thus most familiar to us.

Conference Paper↗

USGS2016

Geoelectric hazard maps for the continental United States

In support of a multiagency project for assessing induction hazards, we present maps of extreme-value geoelectric amplitudes over about half of the continental United States. These maps are constructed using a parameterization of induction: estimates of Earth surface impedance, obtained at discrete geographic sites from magnetotelluric survey data, are convolved with latitude-dependent statistical maps of extreme-value geomagnetic activity, obtained from decades of magnetic observatory data. Geoelectric amplitudes are estimated for geomagnetic waveforms having 240 s sinusoidal period and amplitudes over 10 min that exceed a once-per-century threshold. As a result of the combination of geographic differences in geomagnetic activity and Earth surface impedance, once-per-century geoelectric amplitudes span more than 2 orders of magnitude and are an intricate function of location. For north-south induction, once-per-century geoelectric amplitudes across large parts of the United States have a median value of 0.26 V/km; for east-west geomagnetic variation the median value is 0.23 V/km. At some locations, once-per-century geoelectric amplitudes exceed 3 V/km.

Geophysical Research Letters↗

USGS2015

Ionospheric current source modeling and global geomagnetic induction using ground geomagnetic observatory data

Long-period global-scale electromagnetic induction studies of deep Earth conductivity are based almost exclusively on magnetovariational methods and require accurate models of external source spatial structure. We describe approaches to inverting for both the external sources and three-dimensional (3-D) conductivity variations and apply these methods to long-period ( T ≥1.2 days) geomagnetic observatory data. Our scheme involves three steps: (1) Observatory data from 60 years (only partly overlapping and with many large gaps) are reduced and merged into dominant spatial modes using a scheme based on frequency domain principal components. (2) Resulting modes are inverted for corresponding external source spatial structure, using a simplified conductivity model with radial variations overlain by a two-dimensional thin sheet. The source inversion is regularized using a physically based source covariance, generated through superposition of correlated tilted zonal (quasi-dipole) current loops, representing ionospheric source complexity smoothed by Earth rotation. Free parameters in the source covariance model are tuned by a leave-one-out cross-validation scheme. (3) The estimated data modes are inverted for 3-D Earth conductivity, assuming the source excitation estimated in step 2. Together, these developments constitute key components in a practical scheme for simultaneous inversion of the catalogue of historical and modern observatory data for external source spatial structure and 3-D Earth conductivity.

Journal of Geophysical Research↗