Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

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 865 records · Page 48Linked to original sources
USGS2012

Spatial distribution of damage around faults in the Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah: A mechanical analog for faulting in pyroclastic deposits on Mars

Volcanic ash is thought to comprise a large fraction of the Martian equatorial layered deposits and much new insight into the process of faulting and related fluid flow in these deposits can be gained through the study of analogous terrestrial tuffs. This study identifies a set of fault-related processes that are pertinent to understanding the evolution of fault systems in fine-grained, poorly indurated volcanic ash by investigating exposures of faults in the Miocene-aged Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah. The porosity and granularity of the host rock are found to control the style of localized strain that occurs prior to and contemporaneous with faulting. Deformation bands occur in tuff that was porous and granular at the time of deformation, while fractures formed where the tuff lost its porous and granular nature due to silicic alteration. Non-localized deformation of the host rock is also prominent and occurs through compaction of void space, including crushing of pumice clasts. Significant off-fault damage of the host rock, resembling fault pulverization, is recognized adjacent to one analog fault and may reflect the strain rate dependence of the resulting fault zone architecture. These findings provide important new guidelines for future structural analyses and numerical modeling of faulting and subsurface fluid flow through volcanic ash deposits on Mars.

Utah↗
USGS2012

Soil genesis on the island of Bermuda in the Quaternary: the importance of African dust transport and deposition

The origin of terra rossa, red or reddish-brown, clay-rich soils overlying high-purity carbonate substrates, has intrigued geologists and pedologists for decades. Terra rossa soils can form from accumulation of insoluble residues during dissolution of the host limestones, addition of volcanic ash, or addition of externally derived, long-range-transported (LRT) aeolian particles. We studied soils and paleosols on high-purity, carbonate aeolianites of Quaternary age on Bermuda, where terra rossa origins have been debated for more than a century. Potential soil parent materials on this island include sand-sized fragments of local volcanic bedrock, the LRT, fine-grained (<20 &mu;m) component of distal loess from the lower Mississippi River Valley, and LRT dust from Africa. These parent materials can be characterized geochemically using trace elements that are immobile in the soil-forming environment. Results indicate that local volcanic bedrock on Bermuda has Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that can be distinguished from African dust and lower Mississippi River valley loess. Bermuda soils have Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that indicate derivation from a combination of LRT dust from Africa and local volcanic bedrock. Our results indicate that soils on islands in a very broad latitudinal belt of the western Atlantic margin have been influenced by African LRT dust inputs over much of the past &ndash;500 ka.

Journal of Geophysical Research↗
USGS2012

Calving seismicity from iceberg-sea surface interactions

Iceberg calving is known to release substantial seismic energy, but little is known about the specific mechanisms that produce calving icequakes. At Yahtse Glacier, a tidewater glacier on the Gulf of Alaska, we draw upon a local network of seismometers and focus on 80 hours of concurrent, direct observation of the terminus to show that calving is the dominant source of seismicity. To elucidate seismogenic mechanisms, we synchronized video and seismograms to reveal that the majority of seismic energy is produced during iceberg interactions with the sea surface. Icequake peak amplitudes coincide with the emergence of high velocity jets of water and ice from the fjord after the complete submergence of falling icebergs below sea level. These icequakes have dominant frequencies between 1 and 3 Hz. Detachment of an iceberg from the terminus produces comparatively weak seismic waves at frequencies between 5 and 20 Hz. Our observations allow us to suggest that the most powerful sources of calving icequakes at Yahtse Glacier include iceberg-sea surface impact, deceleration under the influence of drag and buoyancy, and cavitation. Numerical simulations of seismogenesis during iceberg-sea surface interactions support our observational evidence. Our new understanding of iceberg-sea surface interactions allows us to reattribute the sources of calving seismicity identified in earlier studies and offer guidance for the future use of seismology in monitoring iceberg calving.

Journal of Geophysical Research F: Earth Surface↗
USGS2012

Nitrate removal in deep sediments of a nitrogen-rich river network: A test of a conceptual model

Many estimates of nitrogen removal in streams and watersheds do not include or account for nitrate removal in deep sediments, particularly in gaining streams. We developed and tested a conceptual model for nitrate removal in deep sediments in a nitrogen-rich river network. The model predicts that oxic, nitrate-rich groundwater will become depleted in nitrate as groundwater upwelling through sediments encounters a zone that contains buried particulate organic carbon, which promotes redox conditions favorable for nitrate removal. We tested the model at eight sites in upwelling reaches of lotic ecosystems in the Waupaca River Watershed that varied by three orders of magnitude in groundwater nitrate concentration. We measured denitrification potential in sediment core sections to 30 cm and developed vertical nitrate profiles to a depth of about 1 m with peepers and piezometer nests. Denitrification potential was higher, on average, in shallower core sections. However, core sections deeper than 5 cm accounted for 70%, on average, of the depth-integrated denitrification potential. Denitrification potential increased linearly with groundwater nitrate concentration up to 2 mg NO 3 -N/L but the relationship broke down at higher concentrations (> 5 mg NO 3 -N/L), a pattern that suggests nitrate saturation. At most sites groundwater nitrate declined from high concentrations at depth to much lower concentrations prior to discharge into the surface water. The profiles suggested that nitrate removal occurred at sediment depths between 20 and 40 cm. Dissolved oxygen concentrations were much higher in deep sediments than in pore water at 5 cm sediment depth at most locations. The substantial denitrification potential in deep sediments coupled with the declines in nitrate and dissolved oxygen concentrations in upwelling groundwater suggest that our conceptual model for nitrate removal in deep sediments is applicable to this river network. Our results suggest that nitrate removal rates can be high in deep sediments of upwelling stream reaches, which may have implications for efforts to understand and quantify nitrogen transport and removal at larger scales.

Journal of Geophysical Research: Biogeosciences↗
USGS2012

Climatic limits on foliar growth during major droughts in the Southwestern U.S.A.

Pronounced droughts during the 1950s and 2000s in the Southwestern U.S.A. (SW) provide an opportunity to compare mesoscale ecosystem responses to anomalously dry conditions before and during the regional warming that started in the late 1970s. This year-round warming has produced fewer cool season freezes, losses in regional snowpack, an 8-10 day advance in spring onset, and hotter summers, all of which should affect vegetation differently across seasons and elevations. Here, we examine indices that represent climatic limits on foliar growth for both drought periods, and evaluate these indices for areas that experienced tree mortality during the 2000s drought. Relative to the 1950s drought, warmer conditions during the 2000s drought decreased the occurrence of temperatures too low for foliar growth at lower elevations in winter and higher elevations in summer. Higher vapor pressure deficits (VPDs) largely driven by warmer temperatures in the more recent drought were more limiting to foliar growth from spring through summer at lower and middle elevations. At many locations where tree mortality occurred during the 2000s drought, low-temperature constraints on foliar growth were extremely unlimiting, whereas VPD constraints were extremely limiting from early spring through late autumn. Our analysis shows that in physiographically complex regions like the SW, seasonality and elevational gradients are important for understanding vegetative responses to warming. It also suggests that continued warming will increase the degree to which VPD limits foliar growth during future droughts, and expand its reach to higher elevations and other seasons.

Journal of Geophysical Research G: Biogeosciences↗
USGS2012

Quantifying soil surface change in degraded drylands: shrub encroachment and effects of fire and vegetation removal in a desert grassland

Woody plant encroachment, a worldwide phenomenon, is a major driver of land degradation in desert grasslands. Woody plant encroachment by shrub functional types ultimately leads to the formation of a patchy landscape with fertile shrub patches interspaced with nutrient-depleted bare soil patches. This is considered to be an irreversible process of land and soil degradation. Recent studies have indicated that in the early stages of shrub encroachment, when there is sufficient herbaceous connectivity, fires (prescribed or natural) might provide some reversibility to the shrub encroachment process by negatively affecting shrub demography and homogenizing soil resources across patches within weeks to months after burning. A comprehensive understanding of longer term changes in microtopography and spatial patterning of soil properties following fire in shrub-encroached grasslands is desirable. Here, we investigate the changes in microtopography with LiDAR (light detection and ranging), vegetation recovery, and spatial pattering of soil properties in replicated burned, clipped, and control areas in a shrub-grass transition zone in the northern Chihuahuan Desert four years after prescribed fire or clipping. Results indicate a greater homogeneity in soil, microtopography, and vegetation patterning on burned relative to clipped and control treatments. Findings provide further evidence that disturbance by prescribed fire may allow for reversal of the shrub encroachment process, if the event occurs in the early stages of the vegetation shift. Improved understanding of longer-term effects of fire and associated changes in soil patterning can inform the use and role of fire in the context of changing disturbance regimes and climate.

Chihuahuan Desert↗
USGS2012

Hydrogeomorphology of the hyporheic zone: stream solute and fine particle interactions with a dynamic streambed

Hyporheic flow in streams has typically been studied separately from geomorphic processes. We investigated interactions between bed mobility and dynamic hyporheic storage of solutes and fine particles in a sand-bed stream before, during, and after a flood. A conservatively transported solute tracer (bromide) and a fine particles tracer (5 μm latex particles), a surrogate for fine particulate organic matter, were co-injected during base flow. The tracers were differentially stored, with fine particles penetrating more shallowly in hyporheic flow and retained more efficiently due to the high rate of particle filtration in bed sediment compared to solute. Tracer injections lasted 3.5 h after which we released a small flood from an upstream dam one hour later. Due to shallower storage in the bed, fine particles were rapidly entrained during the rising limb of the flood hydrograph. Rather than being flushed by the flood, we observed that solutes were stored longer due to expansion of hyporheic flow paths beneath the temporarily enlarged bedforms. Three important timescales determined the fate of solutes and fine particles: (1) flood duration, (2) relaxation time of flood-enlarged bedforms back to base flow dimensions, and (3) resulting adjustments and lag times of hyporheic flow. Recurrent transitions between these timescales explain why we observed a peak accumulation of natural particulate organic matter between 2 and 4 cm deep in the bed, i.e., below the scour layer of mobile bedforms but above the maximum depth of particle filtration in hyporheic flow paths. Thus, physical interactions between bed mobility and hyporheic transport influence how organic matter is stored in the bed and how long it is retained, which affects decomposition rate and metabolism of this southeastern Coastal Plain stream. In summary we found that dynamic interactions between hyporheic flow, bed mobility, and flow variation had strong but differential influences on base flow retention and flood mobilization of solutes and fine particulates. These hydrogeomorphic relationships have implications for microbial respiration of organic matter, carbon and nutrient cycling, and fate of contaminants in streams.

Journal of Geophysical Research: Biogeosciences↗
USGS2018

Three-dimensional modeling of fine sediment transport by waves and currents in a shallow estuary

A suspended sediment transport model is implemented in the unstructured‐grid SUNTANS model and applied to study fine‐grained sediment transport in South San Francisco Bay. The model enables calculation of suspension of bottom sediment based on combined forcing of tidal currents and wind waves. We show that accurate results can be obtained by employing two‐size classes which are representative of microflocs and macroflocs in the Bay. A key finding of the paper is that the critical calibration parameter is the ratio of the erosion of the microflocs to macroflocs from the bed. Different values of this erosion ratio are needed on the shallow shoals and deeper channels because of the different nature of the sediment dynamics in these regions. Application of a spatially variable erosion ratio and critical shear stress for erosion is shown to accurately reproduce observed suspended sediment concentration at four‐field sites located along a cross‐channel transect. The results reveal a stark contrast between the behavior of the suspended sediment concentration on the shoals and in the deep channel. Waves are shown to resuspend sediments on the shoals, although tidal and wind‐generated currents are needed to mix the thin wave‐driven suspensions into the water column. The contribution to the suspended sediment concentration in the channel by transport from the shoals is similar in magnitude to that due to local resuspension. However, the local contribution is in phase with strong bottom currents which resuspend the sediments, while the contribution from the shoals peaks during low‐water slack tide.

San Francisco Bay↗
USGS2018

The thermophysical properties of the Bagnold Dunes, Mars: Ground truthing orbital data

We compare the thermophysical properties and particle sizes derived from the Mars Science Laboratory rover's Ground Temperature Sensor of the Bagnold dunes, specifically Namib dune, to those derived orbitally from Thermal Emission Imaging System, ultimately linking these measurements to ground truth particle sizes determined from Mars Hand Lens Imager images. In general, we find that all three datasets report consistent particle sizes for the Bagnold dunes (~110–350 μm and are within measurement and model uncertainties), indicating that particle sizes of homogeneous materials inferred from temperature measurements and thermophysical models are reliable. Furthermore, we examine the effects of two physical characteristics that could influence the modeled thermal inertia and particle sizes, including (1) fine‐scale (centimeter to meter scale) ripples and (2) thin layering of indurated/armored materials. To first order, we find that small‐scale ripples and thin (approximately centimeter scale) layers do not significantly affect the determination of bulk thermal inertia from orbital thermal data using a single nighttime temperature. Modeling of a layer of coarse or indurated material reveals that a thin layer (< ~5 mm; similar to what was observed by the Curiosity rover) would not significantly change the observed thermal properties of the surface and would be dominated by the properties of the underlying material. Thermal inertia and particle sizes of relatively homogeneous materials derived from nighttime orbital data should be considered as reliable, as long as there are no significant subpixel anisothermality effects (e.g., lateral mixing of multiple thermophysically distinct materials).

Journal of Geophysical Research: Planets↗
USGS2018

Consequences of abrading bed load on vertical and lateral bedrock erosion in a curved experimental channel

In this study, we conducted multiple physical experiments to estimate the efficacy and spatial pattern of erosion by abrading sediment moving through a simple U‐shaped channel bend with erodible bed and banks. The experiments showed that in the bend, lateral abrasion followed a monotonically increasing linear relationship with sediment feed rate. However, vertical incision had a more complex relation with the sediment feed rate, with an initial increase in abrasion as the feed rate increased followed by a decrease in abrasion of the bed as cover effects became dominant at higher feed rates. Bank erosion was large in places where the width and the lateral slope of the point bar were relatively large. On the other hand, in places where the width of the point bar was smaller, the bedrock bed was eroded primarily along the boundary of the point bar, resulting in a bedrock bench near the outer bank.

Shimanto River↗
USGS2018

A model ensemble for projecting multi‐decadal coastal cliff retreat during the 21st century

Sea cliff retreat rates are expected to accelerate with rising sea levels during the 21 st century. Here we develop an approach for a multi‐model ensemble that efficiently projects time‐averaged sea cliff retreat over multi‐decadal time scales and large (>50 km) spatial scales. The ensemble consists of five simple 1‐D models adapted from the literature that relate sea cliff retreat to wave impacts, sea level rise (SLR), historical cliff behavior, and cross‐shore profile geometry. Ensemble predictions are based on Monte Carlo simulations of each individual model, which account for the uncertainty of model parameters. The consensus of the individual models also weights uncertainty, such that uncertainty is greater when predictions from different models do not agree. A calibrated, but unvalidated, ensemble was applied to the 475 km‐long coastline of Southern California (USA), with 4 SLR scenarios of 0.5, 0.93, 1.5, and 2 m by 2100. Results suggest that future retreat rates could increase relative to mean historical rates by more than two‐fold for the higher SLR scenarios, causing an average total land loss of 19 – 41 m by 2100. However, model uncertainty ranges from +/‐ 5 – 15 m, reflecting the inherent difficulties of projecting cliff retreat over multiple decades. To enhance ensemble performance, future work could include weighting each model by its skill in matching observations in different morphological settings

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Enhancement of a parsimonious water balance model to simulate surface hydrology in a glacierized watershed

The U.S. Geological Survey monthly water balance model (MWBM) was enhanced with the capability to simulate glaciers in order to make it more suitable for simulating cold region hydrology. The new model, MWBMglacier, is demonstrated in the heavily glacierized and ecologically important Copper River watershed in Southcentral Alaska. Simulated water budget components compared well to satellite‐based observations and ground measurements of streamflow, evapotranspiration, snow extent, and total water storage, with differences ranging from 0.2% to 7% of the precipitation flux. Nash Sutcliffe efficiency for simulated and observed streamflow was greater than 0.8 for six of eight stream gages. Snow extent matched satellite‐based observations with Nash Sutcliffe efficiency values of greater than 0.89 in the four Copper River ecoregions represented. During the simulation period 1949 to 2009, glacier ice melt contributed 25% of total runoff, ranging from 12% to 45% in different tributaries, and glacierized area was reduced by 6%. Statistically significant ( p < 0.05) decreasing and increasing trends in annual glacier mass balance occurred during the multidecade cool and warm phases of the Pacific Decadal Oscillation, respectively, reinforcing the link between climate perturbations and glacier mass balance change. The simulations of glaciers and total runoff for a large, remote region of Alaska provide useful data to evaluate hydrologic, cryospheric, ecologic, and climatic trends. MWBM glacier is a valuable tool to understand when, and to what extent, streamflow may increase or decrease as glaciers respond to a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2018

A global empirical model for near real-time assessment of seismically induced landslides

Earthquake-triggered landslides are a significant hazard in seismically active regions, but our ability to assess the hazard they pose in near real-time is limited. In this study, we present a new globally applicable model for seismically induced landslides based on the most comprehensive global dataset available; we use 23 landslide inventories that span a range of earthquake magnitudes and climatic and tectonic settings. We use logistic regression to relate the presence and distribution of earthquake-triggered landslides with spatially distributed estimates of ground shaking, topographic slope, lithology, land-cover type, and a topographic index designed to estimate variability in soil wetness to provide an empirical model of landslide distribution. We tested over 100 combinations of independent predictor variables to find the best-fitting model, using a diverse set of statistical tests. Blind validation tests show the model accurately estimates the distribution of available landslide inventories. The results indicate that the model is reliable and stable, with high “balanced accuracy” (correctly vs. incorrectly classified pixels) for the majority of test events. A cross validation analysis shows high balanced accuracy for a majority of events as well. By combining near-real time estimates of ground shaking with globally available landslide susceptibility data, this model provides a tool to estimate the distribution of co-seismic landslide hazard within minutes of the occurrence of any earthquake worldwide for which a USGS ShakeMap is available.

Journal of Geophysical Research: Earth Surface↗
USGS2018

Long-term evolution of sand transport through a river network: Relative influences of a dam versus natural changes in grain size from sand waves

Temporal and spatial nonuniformity in supplies of water and sand in a river network leads to sand transport that is in local disequilibrium with the upstream sand supply. In such river networks, sand is transported downstream as elongating waves in which coupled changes in grain size and transport occur. Depending on the magnitude of each sand‐supplying event and the interval between such events, changes in bed‐sand grain size associated with sand‐wave passage may more strongly regulate sand transport than do changes in water discharge. When sand transport is controlled more by episodic resupply of sand than by discharge, upstream dam construction may exacerbate or mitigate sand‐transport disequilibria, thus leading to complicated and difficult‐to‐predict patterns of deposition and erosion. We analyzed all historical sediment‐transport data and embarked on a 4‐year program of continuous sediment‐transport measurements to describe disequilibrium sand transport in a river network. Results indicate that sand transport in long river segments can evolve over ≥50‐year timescales following rare large sand‐supplying events. These natural changes in sand transport in distal downstream river segments can be larger than those caused by an upstream dam. Because there is no way to know a priori whether sand transport in a river has changed in response to changes in the upstream sand supply, contemporary continuous measurements of sand transport are required for accurate sand loads and budgeting. Analysis of only historical sediment‐transport measurements, as is common in the literature, may lead to incorrect conclusions with respect to current or future sediment‐transport conditions.

Journal of Geophysical Research: Earth Surface↗
USGS2018

Fire, flood, and drought: Extreme climate events alter flow paths and stream chemistry

Extreme climate events—such as hurricanes, droughts, extreme precipitation, and wildfires—have the potential to alter watershed processes and stream response. Yet due to the destructive and hazardous nature and unpredictability of such events, capturing their hydrochemical signal is challenging. A 5‐year postwildfire study of stream chemistry in the Fourmile Creek watershed, Colorado Front Range, USA, focused on high‐frequency storm sampling. During the study, the watershed was impacted by three additional extreme climate events—drought and two periods of extreme rainfall totals. These events altered concentration‐discharge relationships in ways that elucidate how hydrologic flow paths and source material availability affect stream water chemistry. Reduced infiltration after wildfire led to overland flow during thunderstorms, which conveyed ash and soil into streams. This resulted in elevated stream concentrations of constituents elevated in ash—Ca, K, Mg, alkalinity, and dissolved organic carbon—along with sediment and nitrate. Subsurface flow paths were bypassed, leading to low concentrations of Na and SiO 2 , which are bedrock derived and not elevated in ash. During drought conditions, when stream discharge was <20% of average, concentrations of sediment, dissolved organic carbon, and Ca fell below average concentrations, but SiO 2 did not. Extreme rainfall totals saturated the subsurface and led to prolonged elevated stream discharge. Concentration‐discharge relationships for bedrock‐derived constituents, such as Ca and SiO 2 , were altered in that time period, while those for dissolved organic carbon were not. Previous disturbances, including historical mining, also affect stream chemistry, and water‐quality impairment can be exacerbated by extreme climate events.

Colorado↗
USGS2018

Impact of pore fluid chemistry on fine-grained sediment fabric and compressibility

Fines, defined here as grains or particles, less than 75 μm in diameter, exist nearly ubiquitously in natural sediment, even those classified as coarse. Macroscopic sediment properties, such as compressibility, which relates applied effective stress to the resulting sediment deformation, depend on the fabric of fines. Unlike coarse grains, fines have sizes and masses small enough to be more strongly influenced by electrical interparticle forces than by gravity. These electrical forces acting through pore fluids are influenced by pore fluid chemistry changes. Macroscopic property dependence on pore fluid chemistry must be accounted for in sediment studies involving subsurface flow and sediment stability analyses, as well as in engineered flow situations such as groundwater pollutant remediation, hydrocarbon migration, or other energy resource extraction applications. This study demonstrates how the liquid limit‐based electrical sensitivity index can be used to predict sediment compressibility changes due to pore fluid chemistry changes. Laboratory tests of electrical sensitivity, sedimentation, and compressibility illustrate mechanisms linking microscale and macroscale processes for selected pure, end‐member fines. A specific application considered here is methane extraction via depressurization of gas hydrate‐bearing sediment, which causes a dramatic pore water salinity drop concurrent with sediment being compressed by the imposed effective stress increase.

Journal of Geophysical Research↗
USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Thermal, deformation, and degassing remote sensing time-series (A.D. 2000-2017) at the 47 most active volcanoes in Latin America: Implications for volcanic systems

Volcanoes are hazardous to local and global populations, but only a fraction are continuously monitored by ground-based sensors. For example, in Latin America, more than 60% of Holocene volcanoes are unmonitored, meaning long-term multi-parameter datasets of volcanic activity are rare and sparse. We use satellite observations of degassing, thermal anomalies, and surface deformation spanning 17 years at 47 of the most active volcanoes in Latin America, and compare these datasets to ground-based observations archived by the Global Volcanism Program (GVP). This first comparison of multi-satellite time-series on a regional scale provides information regarding volcanic behavior during, non-, pre-, syn- and post-eruptive periods. For example, at Copahue volcano, deviations from background activity in all three types of satellite measurements were manifested months to years in advance of renewed eruptive activity in 2012. By quantifying the amount of degassing, thermal output, and deformation measured at each of these volcanoes, we test the classification of these volcanoes as open or closed volcanic systems. We find that ~28% of the volcanoes do not fall into either classification and the rest show elements of both, demonstrating a dynamic range of behavior that can change over time. Finally, we recommend how volcano monitoring could be improved through better coordination of available satellite-based capabilities and new instruments.

Journal of Geophysical Research↗