Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial mapping”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

Mapping Soil pH Buffering Capacity of Selected Fields

Soil pH buffering capacity, since it varies spatially within crop production fields, may be used to define sampling zones to assess lime requirement, or for modeling changes in soil pH when acid forming fertilizers or manures are added to a field. Our objective was to develop a procedure to map this soil property. One hundred thirty six soil samples (0 to 15 cm depth) from three Georgia Coastal Plain fields were titrated with calcium hydroxide to characterize differences in pH buffering capacity of the soils. Since the relationship between soil pH and added calcium hydroxide was approximately linear for all samples up to pH 6.5, the slope values of these linear relationships for all soils were regressed on the organic C and clay contents of the 136 soil samples using multiple linear regression. The equation that fit the data best was b (slope of pH vs. lime added) = 0.00029 - 0.00003 * % clay + 0.00135 * % O/C, r(exp 2) = 0.68. This equation was applied within geographic information system (GIS) software to create maps of soil pH buffering capacity for the three fields. When the mapped values of the pH buffering capacity were compared with measured values for a total of 18 locations in the three fields, there was good general agreement. A regression of directly measured pH buffering capacities on mapped pH buffering capacities at the field locations for these samples gave an r(exp 2) of 0.88 with a slope of 1.04 for a group of soils that varied approximately tenfold in their pH buffering capacities.

Weaver, A. R.↗

Features of Point Clouds Synthesized from Multi-View ALOS/PRISM Data and Comparisons with LiDAR Data in Forested Areas

LiDAR waveform data from airborne LiDAR scanners (ALS) e.g. the Land Vegetation and Ice Sensor (LVIS) havebeen successfully used for estimation of forest height and biomass at local scales and have become the preferredremote sensing dataset. However, regional and global applications are limited by the cost of the airborne LiDARdata acquisition and there are no available spaceborne LiDAR systems. Some researchers have demonstrated thepotential for mapping forest height using aerial or spaceborne stereo imagery with very high spatial resolutions.For stereo imageswith global coverage but coarse resolution newanalysis methods need to be used. Unlike mostresearch based on digital surface models, this study concentrated on analyzing the features of point cloud datagenerated from stereo imagery. The synthesizing of point cloud data from multi-view stereo imagery increasedthe point density of the data. The point cloud data over forested areas were analyzed and compared to small footprintLiDAR data and large-footprint LiDAR waveform data. The results showed that the synthesized point clouddata from ALOSPRISM triplets produce vertical distributions similar to LiDAR data and detected the verticalstructure of sparse and non-closed forests at 30mresolution. For dense forest canopies, the canopy could be capturedbut the ground surface could not be seen, so surface elevations from other sourceswould be needed to calculatethe height of the canopy. A canopy height map with 30 m pixels was produced by subtracting nationalelevation dataset (NED) fromthe averaged elevation of synthesized point clouds,which exhibited spatial featuresof roads, forest edges and patches. The linear regression showed that the canopy height map had a good correlationwith RH50 of LVIS data with a slope of 1.04 and R2 of 0.74 indicating that the canopy height derived fromPRISM triplets can be used to estimate forest biomass at 30 m resolution.

LiDARD↗

Ice Roughness and Thickness Evolution on a Swept NACA 0012 Airfoil

Several recent studies have been performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center focusing on the evolution, spatial variations, and proper scaling of ice roughness on airfoils without sweep exposed to icing conditions employed in classical roughness studies. For this study, experiments were performed in the IRT to investigate the ice roughness and thickness evolution on a 91.44-cm (36-in.) chord NACA 0012 airfoil, swept at 30-deg with 0deg angle of attack, and exposed to both Appendix C and Appendix O (SLD) icing conditions. The ice accretion event times used in the study were less than the time required to form substantially three-dimensional structures, such as scallops, on the airfoil surface. Following each ice accretion event, the iced airfoils were scanned using a ROMER Absolute Arm laser-scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger to determine the spatial roughness variations along the surfaces of the iced airfoils. The resulting measurements demonstrate linearly increasing roughness and thickness parameters with ice accretion time. Further, when compared to dimensionless or scaled results from unswept airfoil investigations, the results of this investigation indicate that the mechanisms for early stage roughness and thickness formation on swept wings are similar to those for unswept wings.

glaze icing↗

Complex Spatial Structure in a Population of Didymopanax pittieri, A Tree of Wind-Exposed Lower Montane Rain Forest

Didymopanax pittieri is a common shade-intolerant tree colonizing treefall gaps in the elfin forests on windswept ridgecrests in the lower montane rain forests of the Cordillera de Tilarain, Costa Rica. All D. pittieri taller than > 0.5 m in a 5.2-ha elfin forested portion of a gridded study watershed in the Monteverde Cloud Forest Preserve were located, mapped, and measured. This local population of D. pittieri is spatially inhomogeneous, in that density increases with increasing wind exposure; D. pittieri are more abundant near ridge crests than lower on windward slopes. The important and ubiquitous phenomenon of spatial inhomogeneity in population density is addressed and corrected for in spatial analyses by the application of the inhomogeneous version of Ripley's K. The spatial patterns of four size classes of D. pittieri (<5 cm dbh, 5-10 cm dbh, 10-20 cm dbh, and> 20 cm dbh) were investigated. Within the large-scale trend in density driven by wind exposure, D. pittieri saplings are clumped at the scale of treefall gaps and at the scale of patches of aggregated gaps. D. pittieri 5-10 cm dbh are randomly distributed, apparently due to competitive thinning of sapling clumps during the early stages of gap-phase regeneration. D. pittieri larger than 10 cm dbh are overdispersed at a scale larger than that of patches of gaps. Natural disturbance can influence the distribution of shade intolerant tree populations at several different spatial scales, and can have discordant effects at different life history stages.

Lawton, Robert M.↗

Observations of deformation and mixing of the total ozone field in the Antarctic polar vortex

Total Ozone Mapping Spectrometer (TOMS) images of the springtime Southern Hemisphere commonly show concentric layers in the total ozone field outside the Antarctic polar vortex. The layering appears to result from horizontal folding and stretching of regions on the equatorward flank of the polar vortex near the midlatitude ozone maximum. This folding and stretching interleaves low and high ozone air from the subtropics and midlatitudes, respectively. Occaisional large amplitude wave events can extract very low ozone air from the interior of the polar vortex (the Antarctic ozone hole), but the folding and stretching occurs in midlatitiudes even when wave amplitudes are not exceptionally large. The folding and stretching results in relatively rapid horizontal mixing of the atmosphere on the equatorward flank of the jet. This type of Lagrangian behavior may be common in the atmosphere, but is visible when local tracer gradients are large and observations with high spatial resolution are available. Also, experimentation has shown that gray-scale images of TOMS data show the details of the spatial distribution of ozone much more clearly than contour maps or false-color images.

Bowman, Kenneth P.↗

'Nano' Morphology and Element Signatures of Early Life on Earth: A New Tool for Assessing Biogenicity

The relatively young technology of NanoSIMS is unlocking an exciting new level of information from organic matter in ancient sediments. We are using this technique to characterize Proterozoic organic material that is clearly biogenic as a guide for interpreting controversial organic structures in either terrestrial or extraterrestrial samples. NanoSIMS is secondary ion mass spectrometry for trace element and isotope analysis at sub-micron resolution. In 2005, Robert et al. [1] combined NanoSIMS element maps with optical microscopic imagery in an effort to develop a new method for assessing biogenicity of Precambrian structures. The ability of NanoSIMS to map simultaneously the distribution of organic elements with a 50 nm spatial resolution provides new biologic markers that could help define the timing of life s development on Earth. The current study corroborates the work of Robert et al. and builds on their study by using NanoSIMS to map C, N (as CN), S, Si and O of both excellently preserved microfossils and less well preserved, non-descript organics in Proterozoic chert from the ca. 0.8 Ga Bitter Springs Formation of Australia.

Oehler, D. Z.↗

Multielement mapping of alpha-SiC by scanning Auger microscopy

Fine second-phase particles, numerous in sintered alpha-SiC, were analyzed by scanning Auger microscopy and conventional techniques. The Auger analysis utilized computer-controlled data acquisition, multielement correlation diagrams, and a high spatial resolution of 100 nm. This procedure enabled construction of false color maps and the detection of fine compositional details within these particles. Carbon, silicon oxide, and boron-rich particles (qualitatively as BN or B4C) predominated. The BN particles, sometimes having a carbon core, are believed to result from reaction between B4C additives and nitrogen sintering atmospheres.

Browning, Ray↗

Mapping compositional and particle size variations across Silver Lake Playa: Relevance to analyses of Mars TIR data

The high spectral and spatial resolution thermal infrared (TIR) data to be acquired from the upcoming Mars Observer-Thermal Emission Spectra (TES) mission will map the composition and texture of the Martian sediments. To prepare for these data, portions of two remote sensing experiments were conducted to test procedures for extracting surface property information from TIR data. Reported here is the continuing analysis of Thermal Infrared Multispectral Scanner (TIMS) data, field emission spectra, laboratory Fourier Transform Infrared (FTIR) reflectance spectra, and field observations with respect to the physical characteristics (composition, emissivity, etc.) of Silver Lake playa in southern California.

Petroy, S. B.↗

Acousto-Optic Imaging Spectrometers for Mars Surface Science

NASA's long term plan for Mars sample collection and return requires a highly streamlined approach for spectrally characterizing a landing site, documenting the mineralogical make-up of the site and guiding the collections of samples which represent the diversity of the site. Ideally, image data should be acquired at hundreds of VIS and IR wavelengths, in order to separately distinguish numerous anticipated species, using principal component analysis and linear unmixing. Cameras with bore-sighted point spectrometers can acquire spectra of isolated scene elements, but it requires 10(exp 2) to 10(exp 2) successive motions and precise relative pointing knowledge in order to create a single data cube which qualifies as a spectral map. These and other competing science objectives have to be accomplished within very short lander/rover operational lifetime (a few sols). True, 2-D imaging spectroscopy greatly speeds up the data acquisition process, since the spectra of all pixels in the scene are collected at once. This task can be accomplished with cameras that use electronically tunable acousto-optic tunable filters (AOTFs) as the optical tuning element. AOTFs made from TeO2 are now a mature technology, and operate at wavelengths from near-UV to about 5 microns. Because of incremental improvements in the last few years, present generation devices are rugged, radiation-hard and operate at temperatures down to at least 150K so they can be safely integrated into the ambient temperature optics of in-situ instruments such as planetary or small-body landers. They have been used for ground-based astronomy, and were also baselined for the ST-4 Champollion IR comet lander experiment (CIRCLE), prior to cancellation of the ST-4 mission last year. AIMS (for Acousto-optic Imaging spectrometer), is a prototype lander instrument which is being built at GSFC with support by the NASA OSS Advanced Technologies and Mission Studies, Mars Instrument Definition and Development Program (MIDP). AIMS is capable of tunable spectroscopic imaging of surface mineralogy, ices and dust between 0.5 and 2.4 microns, at a resolving power (lambda/delta lambda) which is typically several hundred. The design spatial resolution, similar to IMP and SSI, will allow mapping at scales down to about 1 cm.

Glenar, D. A.↗

Comparison of multi-temporal NOAA-AVHRR and SPOT-XS satellite data for mapping land-cover dynamics in the West African Sahel

Multi-resolution and multi-temporal remote sensing data (SPOT-XS and AVHRR) were evaluated for mapping local land-cover dynamics in the Sahel of West Africa. The aim of this research was to evaluate the agricultural information that could be derived from both high and low spatial resolution data in areas where there is very often limited ground information. A combination of raster-based image processing and vector-based geographical information system mapping was found to be effective for understanding both spatial and spectral land-cover dynamics. The SPOT data proved useful for mapping local land-cover classes in a dominantly recessive agricultural region. The AVHRR-LAC data could be used to map the dynamics of riparian vegetation, but not the changes associated with recession agriculture. In areas where there was a complex mixture of recession and irrigated agriculture, as well as riparian vegetation, the AVHRR data did not provide an accurate temporal assessment of vegetation dynamics.

Marsh, S. E.↗

Spatial Analysis of Great Lakes Regional Icing Cloud Liquid Water Content

Abstract Clustering of cloud microphysical conditions, such as liquid water content (LWC) and drop size, can affect the rate and shape of ice accretion and the airworthiness of aircraft. Clustering may also degrade the accuracy of cloud LWC measurements from radars and microwave radiometers being developed by the government for remotely mapping icing conditions ahead of aircraft in flight. This paper evaluates spatial clustering of LWC in icing clouds using measurements collected during NASA research flights in the Great Lakes region. We used graphical and analytical approaches to describe clustering. The analytical approach involves determining the average size of clusters and computing a clustering intensity parameter. We analyzed flight data composed of 1-s-frequency LWC measurements for 12 periods ranging from 17.4 minutes (73 km) to 45.3 minutes (190 km) in duration. Graphically some flight segments showed evidence of consistency with regard to clustering patterns. Cluster intensity varied from 0.06, indicating little clustering, to a high of 2.42. Cluster lengths ranged from 0.1 minutes (0.6 km) to 4.1 minutes (17.3 km). Additional analyses will allow us to determine if clustering climatologies can be developed to characterize cluster conditions by region, time period, or weather condition. Introduction

Ryerson, Charles C.↗

4 Vesta in Color: High Resolution Mapping from Dawn Framing Camera Images

Rotational surface variations on asteroid 4 Vesta have been known from ground-based and HST observations, and they have been interpreted as evidence of compositional diversity. NASA s Dawn mission entered orbit around Vesta on July 16, 2011 for a year-long global characterization. The framing cameras (FC) onboard the Dawn spacecraft will image the asteroid in one clear (broad) and seven narrow band filters covering the wavelength range between 0.4-1.0 microns. We present color mapping results from the Dawn FC observations of Vesta obtained during Survey orbit (approx.3000 km) and High-Altitude Mapping Orbit (HAMO) (approx.950 km). Our aim is to create global color maps of Vesta using multi spectral FC images to identify the spatial extent of compositional units and link them with other available data sets to extract the basic mineralogy. While the VIR spectrometer onboard Dawn has higher spectral resolution (864 channels) allowing precise mineralogical assessment of Vesta s surface, the FC has three times higher spatial resolution in any given orbital phase. In an effort to extract maximum information from FC data we have developed algorithms using laboratory spectra of pyroxenes and HED meteorites to derive parameters associated with the 1-micron absorption band wing. These parameters will help map the global distribution of compositionally related units on Vesta s surface. Interpretation of these units will involve the integration of FC and VIR data.

Reddy, V.↗

Trade-Offs of Simplified Versus Comprehensive Representation of Mineralogy When Studying Dust Impacts on Earth’s Climate Systems

The intensity and direction of dust impacts on Earth’s climate systems depend on mineral composition. For example, the presence or absence of a few percent of iron oxides in dust will determine if dust is warming or cooling the atmosphere. Similarly, feldspar will enhance ice cloud formation, while acid gases in the atmosphere will react on the surface of dust calcite limiting acid rain. Still, most climate models use a simplified representation of dust mineralogy. They assume a fixed composition at emission which stays invariant during transport and removal. Such simplification assumes spatially and temporally constant physical and chemical properties of dust, and appears to provide satisfactory results when comparing some properties with observations. The trade-off is their lack of spatial gradients, which will fail to induce circulation, cloud and precipitation changes. The two reasons to omit mineral variations are the uncertainty of current atlases of soil mineral composition in arid regions, and, more practically, an improved runtime efficiency. The former reason is losing ground with the recent launch (July 2022) of a dedicated mission (NASA/JPL EMIT) to retrieve global soil mineralogy of dust sources at high spatial resolution. While the EMIT science team is finalizing a satisfactory global map of mineral composition of dust sources, we analyzed the interaction of dust mineralogy on radiation and its impact on the fast temperature response using different representations of mineral composition from detailed and spatially varying to simplified and globally uniform, assuming different hematite contents and methods to calculate optical properties. Our results show that resolving dust mineralogy reduces dust absorption, and results in improved agreement with observation-based single scattering albedo (SSA), radiative fluxes from CERES (the Clouds and the Earth’s Radiant Energy System), and land surface temperature from CRU (Climatic Research Unit), compared to the baseline bulk dust model version. It also results in distinct radiative impacts on Earth’s climate over North Africa. From our 19-year simulation, we will show that it leads to a reduction of over 50% in net downward radiation at top of atmosphere (TOA) across the Sahara and an approximately 20% reduction over the Sahel. We will explain how the surface temperature response affects the monsoon flow from the Gulf of Guinea. Interestingly, we find similar results by simply fixing the hematite content of dust to a globally uniform value of 0.9% by volume. We will discuss the underlying reasons for such results and show that they may be unrelated to the distribution of soil mineralogy. Still, an accurate representation of soil mineralogy is necessary to better understand dust impacts on the Earth’s climate systems.

dust impacts↗

Geometric processing of digital images of the planets

New procedures and software have been developed for geometric transformation of images to support digital cartography of the planets. The procedures involve the correction of spacecraft camera orientation of each image with the use of ground control and the transformation of each image to a Sinusoidal Equal-Area map projection with an algorithm which allows the number of transformation calculations to vary as the distortion varies within the image. When the distortion is low in an area of an image, few transformation computations are required, and most pixels can be interpolated. When distortion is extreme, the location of each pixel is computed. Mosaics are made of these images and stored as digital databases. Completed Sinusoidal databases may be used for digital analysis and registration with other spatial data. They may also be reproduced as published image maps by digitally transforming them to appropriate map projections.

Edwards, Kathleen↗

Flood Detection with Synthetic Aperture Radar: A Case Study of Houston, Texas Following Hurricane Harvey (2017) using C- and X-Band Observations

Hurricane Harvey produced record-breaking rainfall of up to 60 inches resulting in extensive flooding in Houston, Texas, in late August and early September of 2017. The slow forward motion of the storm following landfall left much of the area unobservable to optical remote sensing instruments for several days due to cloud cover. The active nature of Synthetic Aperture Radar (SAR) instruments allows for observations through clouds, which can supplement efforts to estimate hurricane-induced flood impacts. A growing fleet of SAR constellations has helped lower the latency of imagery following a hurricane, allowing for more timely detections of flooding to help support emergency response efforts. In this study, we leverage publicly available C-band SAR observations from the European Space Agency’s Sentinel-1B (S1B) satellite, collected on 30 August, and X-band SAR imagery collected on 1 September by the Airbus TanDEM-X (TDX) satellite made available through the NASA Commercial Smallsat Data Acquisition (CSDA) program. For each dataset, one co-polarized, StripMap, Radiometric Terrain Corrected (RTC) image was used to create a binary water/no water classification map by referencing permanent water in the Cropland Data Layer (CDL) dataset to determine thresholds. A validation dataset of randomly distributed “ground truth” points was generated using optical imagery from PlanetScope on 31 August, where the domain had relatively little cloud cover. We found that the S1B- and TDX-derived open water maps achieved overall accuracies of 94.17% and 91.57%, respectively. The variation in performance is attributed to both the penetrative abilities of C- and X-band SAR wavelengths in vegetated areas and the increased spatial resolution of the commercial SAR (~3 m) over Sentinel-1 (30 m). While publicly available Sentinel-1 observations are commonly relied upon in response and recovery efforts, these results suggest that including commercial X-band SAR imagery could be beneficial by providing both increased spatial resolution and more frequent revisits when deriving post-event flood mapping products.

Alexander M Melancon↗

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (httl://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Until recently, the forecasters created each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent was to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗

The Temporal and Spatial Development of dB/dt for Substorms

Ground induced currents (GICs) due to space weather are a threat to high voltage power transmission systems. However, knowledge of ground conductivity is the largest source of errors in the determination of GICs. A good proxy for GICs is dB/dt obtained from the Bx and By components of the magnetic field fluctuations. It is known that dB/dt values associated with magnetic storms can reach dangerous levels for power transmission systems. On the other hand, it is not uncommon for dB/dt values associated with substorms to exceed prior Pulkkinen and Molinski critical thresholds of 1.5 nT/s and 5 nT/s, respectively, and the temporal and spatial changes of the dB/dt associated with substorms, unlike storms, are not well understood. Using two dimensional maps of dB/dt over North America and Greenland derived from the spherical elementary currents, we investigate the temporal and spatial change of dB/dt for both a single substorm event and a two dimensional superposed epoch analysis of many substorms. Both the single event and the statistical analysis shows a sudden increase of dB/dt at substorm onset followed by an expansion poleward, westward, and eastward after the onset during the expansion phase. The area of dB/dt values exceeding the two critical thresholds from the initial onset dB/dt values showed little to no expansion equatorward. The temporal and spatial development of the dB/dt resembles the temporal and spatial change of the auroral emissions. Substorm values of dB/dt peak shortly after the auroral onset time and in at least one event exceeded 35 nT/s for a non-storm time substorm. In many of our 81 cases the area that exceeds the threshold of 1.5 nT/s is over several million square kilometers and after about 30 minutes the dB/dt values fall below the threshold level. These results address one of goals of the Space Weather Action Plan, which are to establish benchmarks for space weather events and improve modeling and prediction of their impacts on infrastructure. Plain language: The change in the ground magnetic field with respect to time (dB/dt) associated with magnetic storms (a large disturbance of the magnetic field of the earth) can reach dangerous levels for power transmission systems. On the other hand, substorms, which are a smaller more localized disturbance of the Earth’s magnetic field, are more common. It is not uncommon for substorm dB/dt values to also exceed dangerous levels and the temporal and spatial changes of the dB/dt associated with substorms, unlike storms, are not well understood. Our analysis shows a sudden increase of dB/dt at substorm onset, which peaks shortly after the start of the substorm, followed shortly after by an expansion northward, westward, and eastward after the onset.

ionospheric currents↗