Search USGS⌕ Search

SEARCH · Search USGS

Results for “Earth in Space”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 487 records · Page 27Linked to original sources

Intensity and impact of the New York Railroad superstorm of May 1921

Analysis is made of low‐latitude ground‐based magnetometer data recording the magnetic superstorm of May 1921. By inference, the storm was driven by a series of interplanetary coronal mass ejections, one of which produced a maximum pressure on the magnetopause of ~64.5 nPa, sufficient to compress the subsolar magnetopause radius to ~5.3 Earth radii. Over the course of the storm, low‐latitude geomagnetic disturbance exhibited extreme local time (longitude) asymmetry that can be attributed to substorm disturbance extending to low latitudes. The storm attained an estimated maximum − Dst on 15 May of 907 ± 132 nT, an intensity comparable to that of the Carrington event of 1859. The May 1921 storm brought spectacular aurorae to the nighttime sky. It also interfered with and damaged telephone and telegraph systems associated with railroad systems in New York City and State. These later effects were due to a combination of three factors: the localized details of geomagnetic vector disturbance, the geographic expression of the Earth's surface impedance tensor, and the configurations and physical parameters of the electrical networks of the day.

New York↗

EAARL topography - George Washington Birthplace National Monument 2008

These remotely sensed, geographically referenced elevation measurements of Lidar-derived bare earth (BE) and first surface (FS) topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL; the National Park Service (NPS), Northeast Coastal and Barrier Network, Kingston, RI; and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of the George Washington Birthplace National Monument in Virginia, acquired on March 26, 2008. The datasets are made available for use as a management tool to research scientists and natural resource managers. An innovative airborne Lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL) was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) Lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive Lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for submeter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a Lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of Lidar data in an interactive or batch mode. Modules for presurvey flight line definition, flight path plotting, Lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is routinely used to create maps that represent submerged or first surface topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations.

Virginia↗

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗

EAARL Coastal Topography and Imagery-Naval Live Oaks Area, Gulf Islands National Seashore, Florida, 2007

These remotely sensed, geographically referenced color-infrared (CIR) imagery and elevation measurements of lidar-derived bare-earth (BE) topography, first-surface (FS) topography, and canopy-height (CH) datasets were produced collaboratively by the U.S. Geological Survey (USGS), St. Petersburg Science Center, St. Petersburg, FL; the National Park Service (NPS), Gulf Coast Network, Lafayette, LA; and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of the Naval Live Oaks Area in Florida's Gulf Islands National Seashore, acquired June 30, 2007. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multispectral CIR camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for presurvey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations. For more information about similar projects, please visit the Decision Support for Coastal Science and Management website.

Data Series↗

Evaluation of candidate Landsat Data Gap Sensors

The capabilities of the currently operational Landsat satellites may be lost before the launch of the follow-on Landsat Data Continuity Mission (LDCM), thus producing a gap in the Landsat data record and the National Satellite Land Remote Sensing Data Archive (NSLRSDA). In anticipation of a gap, the Federal agencies responsible for Landsat program management, the National Aeronautics and Space Administration (NASA) and the Department of Interior (DOI) U. S. Geological Survey (USGS), convened a Landsat Data Gap Study Team (LDGST). The study team assessed the basic characteristics of multiple systems and identified sensors aboard the China-Brazil Earth Resources Satellite (CBERS-2) and the Indian Remote Sensing (IRS-P6) ResourceSat-1 satellite as the most promising sources of Landsat-like data. The sensors include the combination of CBERS-2 Infrared Multi-spectral Scanner (IRMSS) and High Resolution Charged Coupled Device (CCD), as well as the IRS-P6 Advanced Wide Field Sensor (AWiFS) and the IRS-P6 Linear Imaging Self Scanning Sensor (LISS-III). The study team concluded that more robust technical evaluations of data and sensor performance are required before gap mitigation strategies can be fully formulated. A technical report is made available that summarizes the results from those evaluations, including the initial data characterization and science utility evaluation. The report can be accessed at http://calval.cr.usgs.gov/LDGST.php.

Conference Paper↗

Calcite U-Pb dating unravels the age and hydrothermal history of the giant Shuiyindong Carlin-type gold deposit in the Golden Triangle, south China

The ages of Carlin-type gold deposits in the Golden Triangle of South China have long been questioned due to the general lack of minerals unequivocally linked to gold deposition that can be precisely dated using conventional radiogenic isotope techniques. Recent advances in U-Pb methods show that calcite can be used to constrain the ages of hydrothermal processes, but few studies have been applied to ore deposits. Herein, we show that this approach can be used to constrain the timing of hydrothermal activity that generated and overprinted the giant Shuiyindong Carlin-type gold deposit in the Golden Triangle. Three stages of calcite (Cal-1, Cal-2, and Cal-3) have been recognized in this deposit based on crosscutting relationships, cathodoluminescence colors, and chemical (U, Pb, and rare earth element [REE]) and isotope (C, O, Sr) compositions. Cal-1 is texturally associated with ore-stage jasperoid and disseminated Au-bearing arsenian pyrite in hydrothermally altered carbonate rocks, which suggests it is synmineralization. Cal-2 fills open spaces and has a distinct orange cathodoluminescence, suggesting that it precipitated during a second fluid pulse. Cal-1 and Cal-2 have similar carbonate rock-buffered chemical and isotopic compositions. Cal-3 occurs in veins that often contain realgar and/or orpiment and are chemically (low U, Pb, and REE) and isotopically (higher δ 13 C, lower δ 18 O and Sr i values) distinct from Cal-1 and Cal-2, suggesting that it formed from a third fluid. U-Pb isotope analyses, by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) for U-rich Cal-1 and Cal-2 and by LA-multicollector (MC)-ICP-MS for U-poor Cal-3, yield well-defined age constraints of 204.3 to 202.6, 191.9, and 139.3 to 137.1 Ma for Cal-1, Cal-2, and Cal-3, respectively. These new ages suggest that the Shuiyindong gold deposit formed in the late Triassic and was overprinted by hydrothermal events in the early Jurassic and early Cretaceous. Given the association of Cal-3 with orpiment and realgar, and previous geochronologic studies of several other major gold deposits in the Golden Triangle, we infer that the latest stage of calcite may be associated with an early Cretaceous regional gold metallogenic event. Combined with existing isotopic ages in the region, these new ages lead us to propose that Carlin-type gold deposits in the Golden Triangle formed during two metallogenic episodes in extensional settings, associated with the late Triassic Indochina orogeny and early Cretaceous paleo-Pacific plate subduction. This study shows that the calcite U-Pb method can be used to constrain the timing of Carlin-type gold deposits and successive hydrothermal events.

Golden Triangle, south China↗

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation↗

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗

Topographic lidar survey of Dauphin Island, Alabama and Chandeleur, Stake, Grand Gosier and Breton Islands, Louisiana, July 12-14, 2013

This Data Series Report contains lidar elevation data collected on July 12 and 14, 2013, for Dauphin Island, Alabama, and Chandeleur, Stake, Grand Gosier and Breton Islands, Louisiana. Classified point cloud data—data points described in three dimensions—in lidar data exchange format (LAS) and bare earth digital elevation models (DEMs) in ERDAS Imagine raster format (IMG) are available as downloadable files. Photo Science, Inc., was contracted by the U.S. Geological Survey (USGS) to collect and process these data. The lidar data were acquired at a horizontal spacing (or nominal pulse spacing) of 1 meter (m) or less. The USGS surveyed points within the project area from July 14–23, 2013, for use in ground control and accuracy assessment. Photo Science, Inc., calculated a vertical root mean square error (RMSEz) of 0.012 m by comparing 10 surveyed points to an interpolated elevation surface of unclassified lidar data. The USGS also checked the data using 80 surveyed points and unclassified lidar point elevation data and found an RMSEz of 0.073 m. The project specified an RMSEz of 0.0925 m or less. The lidar survey was acquired to document the short- and long-term changes of several different barrier island systems. Specifically, this survey supports detailed studies of Chandeleur and Dauphin Islands that resolve annual changes in beaches, berms and dunes associated with processes driven by storms, sea-level rise, and even human restoration activities. These lidar data are available to Federal, State and local governments, emergency-response officials, resource managers, and the general public.

Alabama;Louisiana↗

Current activity on the Martian surface: A key subject for future exploration

One of the fundamental discoveries in Mars science in the last decade has been the extent and importance of current surface activity. Recent results have shifted our view of Mars from a world where the most interesting geologic events were in the distant past (similar to the Moon) to a world that undergoes active evolution and one where understanding the present is key to deciphering the planet’s history. When input was requested for the last Planetary Science Decadal Survey, some observations of surface changes had been published, but the number of detections was small and their significance not fully appreciated. Since that time, detections have proliferated, driven primarily by the long-term operation of the Mars Reconnaissance Orbiter (MRO) and the High Resolution Imaging Science Experiment (HiRISE) as well as landed observations of aeolian activity. In addition to observed changes, theory suggests that additional important surface processes are likely active but not yet observed because orbital data are limited in space and time and landed studies are rare. Understanding current Martian surface processes is a fundamental science question in itself, as it provides a test for physical and terrestrial analog-based models of specific geological processes acting under non-Earth planetary and environmental conditions. It is also an essential step for reading Mars’ geologic history and providing input to climate models: without understanding current dynamic processes, we cannot understand how they have varied during recent climate cycles, nor how they are reflected in ancient rock or modern ice records. Understanding the rates and types of current surface activity is also highly relevant to selecting geological samples, setting Planetary Protection rules, and understanding the hazards and environment that would be experienced by future human explorers.

Bulletin of the AAS↗

Characteristics of the Landsat Multispectral Data System

Landsat satellites were launched into orbit in 1972 and 1975. Additional Landsat satellites are planned for launch in 1978 and 1981. The satellites orbit the Earth at an altitude of approximately 900 km and each can obtain repetitive coverage of cloud-free areas every 18 days. A sun-synchronous orbit is used to insure repeatable illumination conditions. Repetitive satellite coverage allows optimal cover conditions for geologic applications to be identified. Seasonal variations in solar illumination must be analyzed to select the best Landsat data for geologic applications. Landsat data may be viewed in stereo where there is sufficient sidelap and sufficient topographic relief. Landsat-1 ceased operation on January 10, 1978. Landsat-2 detects, only solar radiation that is reflected from the Earth's surface in visible and near-visible wavelengths. The third Landsat will also detect emitted thermal radiation. The multispectral scanner (MSS) was the only sensing instrument used on the first two satellites. The MSS on Landsats-1 and -2 detect radiation which is reflected from a 79 m by 79 m area, and the data are formatted as if the measurement was made from a 56 m by 79 m area. The MSS integrates spectral response from all cover types within the 79 m by 79 m area. The integrated spectral signature often does not resemble the spectral signature from individual cover types, and the integrated signature is also modified by the atmosphere. Landsat-1 and -2 data are converted to 70 mm film and computer compatible tapes (CCT's) at Goddard Space Flight Center (GSFC); these are shipped to the EROS Data Center (EDC) for duplication and distribution to users. Landsat-C data will be converted to 241 mm-wide film and CCT's at EDC. Landsat-D data will be relayed from the satellite directly to geosynchronous satellites and then to the United States from any location on Earth.

Open-File Report↗

Ten ways remote sensing can contribute to conservation

In an effort to increase conservation effectiveness through the use of Earth observation technologies, a group of remote sensing scientists affiliated with government and academic institutions and conservation organizations identified 10 questions in conservation for which the potential to be answered would be greatly increased by use of remotely sensed data and analyses of those data. Our goals were to increase conservation practitioners’ use of remote sensing to support their work, increase collaboration between the conservation science and remote sensing communities, identify and develop new and innovative uses of remote sensing for advancing conservation science, provide guidance to space agencies on how future satellite missions can support conservation science, and generate support from the public and private sector in the use of remote sensing data to address the 10 conservation questions. We identified a broad initial list of questions on the basis of an email chain-referral survey. We then used a workshop-based iterative and collaborative approach to whittle the list down to these final questions (which represent 10 major themes in conservation): How can global Earth observation data be used to model species distributions and abundances? How can remote sensing improve the understanding of animal movements? How can remotely sensed ecosystem variables be used to understand, monitor, and predict ecosystem response and resilience to multiple stressors? How can remote sensing be used to monitor the effects of climate on ecosystems? How can near real-time ecosystem monitoring catalyze threat reduction, governance and regulation compliance, and resource management decisions? How can remote sensing inform configuration of protected area networks at spatial extents relevant to populations of target species and ecosystem services? How can remote sensing-derived products be used to value and monitor changes in ecosystem services? How can remote sensing be used to monitor and evaluate the effectiveness of conservation efforts? How does the expansion and intensification of agriculture and aquaculture alter ecosystems and the services they provide? How can remote sensing be used to determine the degree to which ecosystems are being disturbed or degraded and the effects of these changes on species and ecosystem functions?

Conservation Biology↗

Mapping of the Moon by Clementine

The "faster, cheaper, better" Clementine spacecraft mission mapped the Moon from February 19 to May 3, 1994. Global coverage was acquired in 11 spectral bandpasses from 415 to 2792 nm and at resolutions of 80-330 m/pixel; a thermal-infrared camera sampled ???20% of the surface; a high-resolution camera sampled selected areas (especially the polar regions); and a lidar altimeter mapped the large-scale topography up to latitudes of ??75??. The spacecraft was in a polar, elliptical orbit, 400-450 km periselene altitude. Periselene latitude was -28.5?? for the first month of mapping, then moved to +28.5??. NASA is supporting the archiving, systematic processing, and analysis of the ???1.8 million lunar images and other datasets. A new global positional network has been constructed from 43,000 images and ???0.5 million match points; new digital maps will facilitate future lunar exploration. In-flight calibrations now enable photometry to a high level of precision for the uv-visible CCD camera. Early science results include: (1) global models of topography, gravity, and crustal thicknesses; (2) new information on the topography and structure of multiring impact basins; (3) evidence suggestive of water ice in large permanent shadows near the south pole; (4) global mapping of iron abundances; and (5) new constraints on the Phanerozoic cratering rate of the Earth. Many additional results are expected following completion of calibration and systematic processing efforts. ?? 1997 COSPAR. Published by Elsevier Science Ltd.

Advances in Space Research↗

A stochastic fault model. 2. Time-dependent case

A random model of fault motion in an earthquake is formulated by assuming that the slip velocity is a random function of position and time truncated at zero, so that it does not have negative values. This random function is chosen to be self-affine; that is, on change of length scale, the function is multiplied by a scale factor but is otherwise unchanged statistically. A snapshot of slip velocity at a given time resembles a cluster of islands with rough topography; the final slip function is a smoother island or cluster of islands. In the Fourier transform domain, shear traction on the fault equals the slip velocity times an impedance function. The fact that this impedance function has a pole at zero frequency implies that traction and slip velocity cannot have the same spectral dependence in space and time. To describe stress fluctuations of the order of 100 bars when smoothed over a length of kilometers and of the order of kilobars at the grain size, shear traction must have a one-dimensional power spectrum is space proportional to the reciprocal wave number. Then the one-dimensional power spectrum for the slip velocity is proportional to the reciprocal wave number squared and for slip to its cube. If slip velocity has the same power law spectrum in time as in space, then the spectrum of ground acceleration will be flat (white noise) both on the fault and in the far field.

Journal of Geophysical Research Solid Earth↗

Ultraviolet and visible remote sensing of volcanic gases

As magma rises in volcanic systems, volatile species exsolve from the silicate melt and are emitted as gases into the atmosphere. Measuring the magnitude and composition of gas emissions from volcanoes provides insights into processes occurring deep within the Earth and helps constrain the impact of volcanic degassing on atmospheric chemistry. Optical remote sensing techniques allow volcanic gas emissions to be characterized without the need to access hazardous areas near active volcanic vents. This paper reviews the state of the art in ultraviolet and visible volcanic gas remote sensing from the ground, air, and space. Special attention is given to discussing the physics of atmospheric radiative transfer on which these techniques are based. The functionality and limitations of different remote sensing instruments are examined, making clear that the ideal choice of instrumentation will depend on the volcanic system to which it is applied and the sought measurement parameters. Common algorithms for determining trace gas column densities, gas burdens, and volcanic emission rates from measurements of spectral radiance are outlined and compared, showing how some algorithms attempt to model the physics of the measurement while others maximize sensitivity. Several examples demonstrate how remote sensing measurements continue to advance our understanding of volcanic systems and their impact on the atmosphere. Finally, a few promising directions of inquiry are suggested that could lead to improvements in remote sensing instrumentation and analysis techniques. By combining spectroscopic and imaging techniques, improving our understanding of atmospheric radiative transfer, expanding the suite of target gases, and increasing the coverage and frequency of observations, we stand to significantly improve our ability to detect and quantify volcanic gas emissions and gain new insights into important Earth-system processes.

Journal of Volcanology and Geothermal Research↗

Evaluation of a new model of aeolian transport in the presence of vegetation

Aeolian transport is an important characteristic of many arid and semiarid regions worldwide that affects dust emission and ecosystem processes. The purpose of this paper is to evaluate a recent model of aeolian transport in the presence of vegetation. This approach differs from previous models by accounting for how vegetation affects the distribution of shear velocity on the surface rather than merely calculating the average effect of vegetation on surface shear velocity or simply using empirical relationships. Vegetation, soil, and meteorological data at 65 field sites with measurements of horizontal aeolian flux were collected from the Western United States. Measured fluxes were tested against modeled values to evaluate model performance, to obtain a set of optimum model parameters, and to estimate the uncertainty in these parameters. The same field data were used to model horizontal aeolian flux using three other schemes. Our results show that the model can predict horizontal aeolian flux with an approximate relative error of 2.1 and that further empirical corrections can reduce the approximate relative error to 1.0. The level of error is within what would be expected given uncertainties in threshold shear velocity and wind speed at our sites. The model outperforms the alternative schemes both in terms of approximate relative error and the number of sites at which threshold shear velocity was exceeded. These results lend support to an understanding of the physics of aeolian transport in which (1) vegetation's impact on transport is dependent upon the distribution of vegetation rather than merely its average lateral cover and (2) vegetation impacts surface shear stress locally by depressing it in the immediate lee of plants rather than by changing the bulk surface's threshold shear velocity. Our results also suggest that threshold shear velocity is exceeded more than might be estimated by single measurements of threshold shear stress and roughness length commonly associated with vegetated surfaces, highlighting the variation of threshold shear velocity with space and time in real landscapes.

Journal of Geophysical Research F: Earth Surface↗

Insights into shallow magmatic processes at Kīlauea Volcano, Hawaiʻi, from a multiyear continuous gravity time series

Continuous gravity data collected near the summit eruptive vent at Kīlauea Volcano, Hawaiʻi, during 2011–2015 show a strong correlation with summit-area surface deformation and the level of the lava lake within the vent over periods of days to weeks, suggesting that changes in gravity reflect variations in volcanic activity. Joint analysis of gravity and lava level time series data indicates that over the entire time period studied, the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low—approximately 1000–1500 kg/m 3 . The ratio of gravity change (adjusted for Earth tides and instrumental drift) to lava level change measured over 15 day windows rose gradually over the course of 2011–2015, probably reflecting either (1) a small increase in the density of lava within the eruptive vent or (2) an increase in the volume of lava within the vent due to gradual vent enlargement. Superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of Kīlauea's summit and coincident changes in lava level. The unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation, perhaps indicating magma accumulation within, and withdrawal from, void space—a process that might not otherwise be apparent from lava level and deformation data alone. Continuous gravity data thus provide unique insights into magmatic processes, arguing for continued application of the method at other frequently active volcanoes.

Hawaii↗