USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Observer”

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 1,261 records · Page 70Linked to original sources

The Border Ranges fault system in Glacier Bay National Park, Alaska: Evidence for major early Cenozoic dextral strike-slip motion

The Border Ranges fault system of southern Alaska, the fundamental break between the arc basement and the forearc accretionary complex, is the boundary between the Peninsular–Alexander–Wrangellia terrane and the Chugach terrane. The fault system separates crystalline rocks of the Alexander terrane from metamorphic rocks of the Chugach terrane in Glacier Bay National Park. Mylonitic rocks in the zone record abundant evidence for dextral strike-slip motion along north-northwest-striking subvertical surfaces. Geochronologic data together with regional correlations of Chugach terrane rocks involved in the deformation constrain this movement between latest Cretaceous and Early Eocene (~50 Ma). These findings are in agreement with studies to the northwest and southeast along the Border Ranges fault system which show dextral strike-slip motion occurring between 58 and 50 Ma. Correlations between Glacier Bay plutons and rocks of similar ages elsewhere along the Border Ranges fault system suggest that as much as 700 km of dextral motion may have been accommodated by this structure. These observations are consistent with oblique convergence of the Kula plate during early Cenozoic and forearc slivering above an ancient subduction zone following late Mesozoic accretion of the Peninsular–Alexander–Wrangellia terrane to North America.

Alaska↗

Hydroacoustic, meteorologic and seismic observations of the 2016 Nansen Ice Shelf calving event and iceberg formation

On 7 April 2016, the Nansen Ice Shelf (NIS) front calved into two icebergs, the first large-scale calving event in >30 years. Three hydrophone moorings were deployed seaward of the NIS in December 2015 and over the following months recorded hundreds of short duration, broadband (10-400 Hz) cryogenic signals, likely caused by fracturing of the ice-shelf. The majority of these icequakes occur between January and early March 2016, several weeks prior to the calving observed by satellite on 7 April. Barometric pressure and wind speed records show the day the icebergs drifted from the NIS coincided with the largest low-pressure storm system recorded in the previous 7 months. A nearby seismic station also shows an increase in low-frequency energy, harmonic tremor, and microseisms on 7 April. Our interpretation is the northern segment of the NIS leading edge broke free during mid-January to February, producing high acoustic energy, but the icebergs remained stationary until the combination of a strong low-pressure system, with high winds freed the icebergs. As the unpinning of Antarctic ice shelves is not a well-documented process, our observations show that storm systems may play an under-appreciated role in Antarctic ice shelf break-up.

Frontiers in Earth Science↗

Statistical analysis of water-quality data containing multiple detection limits II: S-language software for nonparametric distribution modeling and hypothesis testing

Analysis of low concentrations of trace contaminants in environmental media often results in left-censored data that are below some limit of analytical precision. Interpretation of values becomes complicated when there are multiple detection limits in the data-perhaps as a result of changing analytical precision over time. Parametric and semi-parametric methods, such as maximum likelihood estimation and robust regression on order statistics, can be employed to model distributions of multiply censored data and provide estimates of summary statistics. However, these methods are based on assumptions about the underlying distribution of data. Nonparametric methods provide an alternative that does not require such assumptions. A standard nonparametric method for estimating summary statistics of multiply-censored data is the Kaplan-Meier (K-M) method. This method has seen widespread usage in the medical sciences within a general framework termed "survival analysis" where it is employed with right-censored time-to-failure data. However, K-M methods are equally valid for the left-censored data common in the geosciences. Our S-language software provides an analytical framework based on K-M methods that is tailored to the needs of the earth and environmental sciences community. This includes routines for the generation of empirical cumulative distribution functions, prediction or exceedance probabilities, and related confidence limits computation. Additionally, our software contains K-M-based routines for nonparametric hypothesis testing among an unlimited number of grouping variables. A primary characteristic of K-M methods is that they do not perform extrapolation and interpolation. Thus, these routines cannot be used to model statistics beyond the observed data range or when linear interpolation is desired. For such applications, the aforementioned parametric and semi-parametric methods must be used.

Computers & Geosciences↗

The consequences of climate change for dryland biogeochemistry

Drylands, which cover more than 40% of Earth’s terrestrial surface, are dominant drivers of global biogeochemical cycling and home to more than one third of the human population. Climate projections predict warming, drought frequency and severity, and evaporative demand will increase in drylands at faster rates than global means. Due to extreme temperatures and high biological dependency on limited water availability, drylands are predicted to be exceptionally sensitive to climate change and, indeed, significant climate impacts are already being observed. Yet our understanding and ability to forecast climate change effects on dryland biogeochemistry and ecosystem functions lag behind many mesic systems. To improve our capacity to forecast ecosystem change, we propose focusing on the controls and consequences of two key characteristics affecting dryland biogeochemistry: i) high spatial and temporal heterogeneity in environmental conditions and ii) generalized resource scarcity. In addition to climate change, drylands are experiencing accelerating land use change. Building our understanding of dryland biogeochemistry in both intact and disturbed systems will better equip us to address the interacting effects of climate change and landscape degradation. Responding to these challenges will require a diverse, globally distributed, and interdisciplinary community of dryland experts united towards better understanding these vast and important ecosystems.

New Phytologist Foundation↗

Evolution of Ore Deposits and Technology Transfer Project: Isotope and Chemical Methods in Support of the U.S. Geological Survey Science Strategy, 2003-2008

Principal functions of the U.S. Geological Survey (USGS) Mineral Resources Program are providing assessments of the location, quantity, and quality of undiscovered mineral deposits, and predicting the environmental impacts of exploration and mine development. The mineral and environmental assessments of domestic deposits are used by planners and decisionmakers to improve the stewardship of public lands and public resources. Assessments of undiscovered mineral deposits on a global scale reveal the potential availability of minerals to the United States and other countries that manufacture goods imported to the United States. These resources are of fundamental relevance to national and international economic and security policy in our globalized world economy. Performing mineral and environmental assessments requires that predictions be made of the likelihood of undiscovered deposits. The predictions are based on geologic and geoenvironmental models that are constructed for the diverse types of mineral deposits from detailed descriptions of actual deposits and detailed understanding of the processes that formed them. Over the past three decades the understanding of ore-forming processes has benefited greatly from the integration of laboratory-based geochemical tools with field observations and other data sources. Under the aegis of the Evolution of Ore Deposits and Technology Transfer Project (referred to hereinafter as the Project), a 5-year effort that terminated in 2008, the Mineral Resources Program provided state-of-the-art analytical capabilities to support applications of several related geochemical tools to ore-deposit-related studies. The analytical capabilities and scientific approaches developed within the Project have wide applicability within Earth-system science. For this reason the Project Laboratories represent a valuable catalyst for interdisciplinary collaborations of the type that should be formed in the coming years for the United States to meet its natural-resources and natural-science needs. This circular presents an overview of the Project. Descriptions of the Project laboratories are given first including descriptions of the types of chemical or isotopic analyses that are made and the utility of the measurements. This is followed by summaries of select measurements that were carried out by the Project scientists. The studies are grouped by science direction. Virtually all of them were collaborations with USGS colleagues or with scientists from other governmental agencies, academia, or the private sector.

Circular↗

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth↗

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

Episodic fluid flow in the Nankai accretionary complex: Timescale, geochemistry, flow rates, and fluid budget

Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 m M , but measured chlorinity is as low as 447 m M . We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3–4 orders of magnitude along a flow conduit (in this case the décollement), using steady state results as initial conditions. Transient results with an increase in décollement permeability from 10 −16 m 2 to 10 −13 m 2 and 20% smectite reproduce the observed chloride profile at site 808 after 80–160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the décollement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2–3 orders of magnitude smaller than observed. Fluid budget calculations indicate that nearly 71% of the incoming water in the sediments leaves the accretionary wedge via diffuse flow out the seafloor, 0–5% escapes by focused flow along the décollement, and roughly 1% is subducted.

Journal of Geophysical Research B: Solid Earth↗

Trace element abundances of high-MgO glasses from Kilauea, Mauna Loa and Haleakala volcanoes, Hawaii

We performed an ion-microprobe study of eleven high-MgO (6.7-14.8 wt%) tholeiite glasses from the Hawaiian volcanoes Kilauea, Mauna Loa and Haleakala. We determined the rare earth (RE), high field strength, and other selected trace element abundances of these glasses, and used the data to establish their relationship to typical Hawaiian shield tholeiite and to infer characteristics of their source. The glasses have trace element abundance characteristics generally similar to those of typical shield tholeiites, e.g. L(light)REE/H(heavy)REE(C1) > 1. The Kilauea and Mauna Loa glasses, however, display trace and major element characteristics that cross geochemical discriminants observed between Kilauea and Mauna Loa shield lavas. The glasses contain a blend of these discriminating chemical characteristics, and are not exactly like the typical shield lavas from either volcano. The production of these hybrid magmas likely requires a complexly zoned source, rather than two unique sources. When corrected for olivine fractionation, the glass data show correlations between CaO concentration and incompatible trace element abundances, indicating that CaO may behave incompatibly during melting of the tholeiite source. Furthermore, the tholeiite source must contain residual garnet and clinopyroxene to account for the variation in trace element abundances of the Kilauea glasses. Inversion modeling indicates that the Kilauea source is flat relative to C1 chondrites, and has a higher bulk distribution coefficient for the HREE than the LREE.

Hawaii↗

Nucleation and triggering of earthquake slip: Effect of periodic stresses

Results of stability analyses for spring and slider systems, with state variable constitutive properties, are applied to slip on embedded fault patches. Unstable slip may nucleate only if the slipping patch exceeds some minimum size. Subsequent to the onset of instability the earthquake slip may propagate well beyond the patch. It is proposed that the seismicity of a volume of the earth's crust is determined by the distribution of initial conditions on the population of fault patches that nucleate earthquake slip, and the loading history acting upon the volume. Patches with constitutive properties inferred from laboratory experiments are characterized by an interval of self-driven accelerating slip prior to instability, if initial stress exceeds a minimum threshold. This delayed instability of the patches provides an explanation for the occurrence of aftershocks and foreshocks including decay of earthquake rates by time −1 . A population of patches subjected to loading with a periodic component results in periodic variation of the rate of occurrence of instabilities. The change of the rate of seismicity for a sinusoidal load is proportional to the amplitude of the periodic stress component and inversely proportional to both the normal stress acting on the fault patches and the constitutive parameter, A 1 , that controls the direct velocity dependence of fault slip. Values of A 1 representative of laboratory experiments indicate that in a homogeneous crust, correlation of earthquake rates with earth tides should not be detectable at normal stresses in excess of about 8 MPa. Correlation of earthquakes with tides at higher normal stresses can be explained if there exist inhomogeneities that locally amplify the magnitude of the tidal stresses. Such amplification might occur near magma chambers or other soft inclusions in the crust and possibly near the ends of creeping fault segments if the creep or afterslip rates vary in response to tides. Observations of seismicity rate variations associated with seasonal fluctuations of reservoir levels appear to be consistent with the model.

Tectonophysics↗

Ash3d: A finite-volume, conservative numerical model for ash transport and tephra deposition

We develop a transient, 3-D Eulerian model (Ash3d) to predict airborne volcanic ash concentration and tephra deposition during volcanic eruptions. This model simulates downwind advection, turbulent diffusion, and settling of ash injected into the atmosphere by a volcanic eruption column. Ash advection is calculated using time-varying pre-existing wind data and a robust, high-order, finite-volume method. Our routine is mass-conservative and uses the coordinate system of the wind data, either a Cartesian system local to the volcano or a global spherical system for the Earth. Volcanic ash is specified with an arbitrary number of grain sizes, which affects the fall velocity, distribution and duration of transport. Above the source volcano, the vertical mass distribution with elevation is calculated using a Suzuki distribution for a given plume height, eruptive volume, and eruption duration. Multiple eruptions separated in time may be included in a single simulation. We test the model using analytical solutions for transport. Comparisons of the predicted and observed ash distributions for the 18 August 1992 eruption of Mt. Spurr in Alaska demonstrate to the efficacy and efficiency of the routine.

Journal of Geophysical Research↗

Three-dimensional electrical resistivity characterization of Mountain Pass, California and surrounding region

The Sulphide Queen carbonatite deposit at Mountain Pass in southeast California is a world class rare earth element (REE) resource. This study images electrical resistivity structure of the REE deposit and surrounding area to characterize resources under cover. An east-west elongated grid (35 × 15 km) of 65 wideband magnetotelluric stations spanning from eastern Shadow Valley to eastern Ivanpah Valley were collected and modeled in three-dimensions (3-D). Gravity, aeromagnetic, and geologic data are used to inform interpretation of structures in the resistivity model, including the following observations. Shadow Valley is filled with conductive sediment that locally dips southward to a depth of 1 km. The Kingston Range-Halloran Hills detachment fault dips westward at ∼15 degrees. The REE deposit is a moderate low resistivity zone dipping southwest to a possible depth of ∼1 km, and is bounded by the North and South faults and bisected by the Middle fault. Ivanpah Dry Lake is underlain by a north striking southward dipping sedimentary basin. Two possible zones of mineralization are observed in Ivanpah Valley, one along the western edge of Ivanpah Dry Lake and one on the western edge of valley along a new inferred fault. The brittle-ductile transition is imaged at ∼10 km below mean sea level. No deep electrically conductive structures are imaged to be related to the REE deposit likely due to the complex geologic history of the Mojave terrane. Future studies should regional target Proterozoic rocks and search within for geophysical signatures similar to Mountain Pass.

California↗

Initial pulse of Siberian Traps sills as the trigger of the end-Permian mass extinction

Mass extinction events are short-lived and characterized by catastrophic biosphere collapse and subsequent reorganization. Their abrupt nature necessitates a similarly short-lived trigger, and large igneous province magmatism is often implicated. However, large igneous provinces are long-lived compared to mass extinctions. Therefore, if large igneous provinces are an effective trigger, a subinterval of magmatism must be responsible for driving deleterious environmental effects. The onset of Earth’s most severe extinction, the end-Permian, coincided with an abrupt change in the emplacement style of the contemporaneous Siberian Traps large igneous province, from dominantly flood lavas to sill intrusions. Here we identify the initial emplacement pulse of laterally extensive sills as the critical deadly interval. Heat from these sills exposed untapped volatile-fertile sediments to contact metamorphism, likely liberating the massive greenhouse gas volumes needed to drive extinction. These observations suggest that large igneous provinces characterized by sill complexes are more likely to trigger catastrophic global environmental change than their flood basalt- and/or dike-dominated counterparts.

Siberia, Siberian Traps↗

Quaternarymagmatism in the Cascades - Geologic perspectives

Foreward The Cascade magmatic arc is a belt of Quaternary volcanoes that extends 1,250 km from Lassen Peak in northern California to Meager Mountain in Canada, above the subduction zone where the Juan de Fuca Plate plunges beneath the North American Plate. This Professional Paper presents a synthesis of the entire volcanic arc, addressing all 2,300 known Quaternary volcanoes, not just the 30 or so visually prominent peaks that comprise the volcanic skyline. Study of Cascade volcanoes goes back to the geological explorers of the late 19th century and the seminal investigations of Howel Williams in the 1920s and 1930s. However, major progress and application of modern scientific methods and instrumentation began only in the 1970s with the advent of systematic geological, geophysical, and geochemical studies of the entire arc. Initial stimulus from the USGS Geothermal Research Program was enhanced by the USGS Volcano Hazards Program following the 1980 eruption of Mount St. Helens. Together, these two USGS Programs have provided more than three decades of stable funding, staffing, and analytical support. This Professional Paper summarizes the resultant USGS data sets and integrates them with the parallel contributions of other investigators. The product is based upon an all-encompassing and definitive geological database, including chemical and isotopic analyses to characterize the rocks and geochronology to provide the critical time constraints. Until now, this massive amount of data has not been summarized, and a systematic and uniform interpretation firmly grounded in geological fact has been lacking. Herein lies the primary utility of this Cascade volume. It not only will be the mandatory starting point for new workers, but also will provide essential geological context to broaden the perspectives of current investigators of specific Cascade volcanoes. Wes Hildreth's insightful understanding of volcanic processes and his uncompromising scientific integrity make him uniquely qualified to present this synthesis. During more than three decades of volcanological studies, he has carried out comprehensive investigations of Mount Adams, Mount Baker, the Three Sisters, and the Simcoe Mountains Volcanic Field. He also brings a broad experience in other volcanic arcs, having conducted integrated field and laboratory investigations at several major volcanic centers in the Andes and the Aleutian arcs. His expertise and perspective have been further enhanced by in-depth petrologic studies of caldera environments, primarily in Long Valley, California, and Yellowstone. On the basis of all these field and laboratory investigations and exhaustive literature searches, he has published three definitive petrologic syntheses addressing the passage and transformation of basaltic magmas from their mantle sources through the crust to form the many types of volcanic manifestations at the Earth's surface. A major strength of this Professional Paper is that it adheres to data first and foremost, and only then correlates these data with relevant theories. Petrological and geophysical interpretation is left to the later sections of the volume, and even there is never allowed to stray from the pertinent databases. Hildreth's interpretations are not just idle speculations, but are carefully reasoned inferences firmly based on his thorough evaluation of the observational geological data. Professional Paper 1744 should not be skimmed lightly, in the hope that the salient points will quickly rub off. Instead, every section, indeed every paragraph, presents scholarly observations and insightful interpretations that demand careful and thoughtful study. This volume will influence and guide the course of Cascade investigations for decades to come.

British Columbia, California, Oregon, Washington↗

Magmatic (silicates/saline/sulfur-rich/CO2) immiscibility and zirconium and rare-earth element enrichment from alkaline magma chamber margins : Evidence from Ponza Island, Pontine Archipelago, Italy

Fluid inclusions were measured from a feldspathoid-bearing syenite xenolith entrained in trachyte from Ponza, one of the islands of the Pontine Archipelago, located in the Gulf of Gaeta, Italy. The feldspathoid-bearing syenite consists mainly of potassium feldspar, clinopyroxene, amphibole, biotite, titanite, manganoan magnetite, apatite with minor nosean, Na-rich feldspar, pyrrhotite, and rare cheralite. Baddeleyite and zirkelite occur associated with manganoan magnetite. Detailed electron-microprobe analysis reveals enrichments in REE, Y, Nb, U, Th as well as Cl and F in appropriate phases. Fluid inclusions observed in potassium feldspar are either silicate-melt or aqueous inclusions. The aqueous inclusions can be further classified as. (1) one-phase vapor, (2) two-phase (V + L) inclusions, vapor-rich inclusions with a small amount of CO2 in most cases; homogenization of the inclusions always occurred in the vapor phase between 359 and 424??C, salinities vary from 2.9 to 8.5 wt. % NaCl equivalent; and. (3) three-phase and multiphase inclusions (hypersaline/sulfur-rich aqueous inclusions sometimes with up to 8 or more solid phases). Daughter minerals dissolve on heating before vapor/liquid homogenization. Standardless quantitative scanning electron microscope X-ray fluorescence analysis has tentatively identified the following chloride and sulfate daughter crystals; halite, sylvite, glauberite. arcanite, anhydrite, and thenardite. Melting of the daughter crystals occurs between 459 and 536??C (54 to 65 wt. % NaCI equivalent) whereas total homogenization is between 640 and 755??C. The occurrence of silicate-melt inclusions and high-temperature, solute-rich aqueous inclusions suggests that the druse or miarolitic texture of the xenolith is late-stage magmatic. The xenolith from Ponza represents a portion of the peripheral magma chamber wall that has recorded the magmatic/hydrothermal transition and the passage of high solute fluids enriched in chlorides, sulfur, and incompatible elements.

European Journal of Mineralogy↗

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): A new observational database to support explosive eruptive column model validation and development

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Journal of Volcanology and Geothermal Research↗

A synthesis of the basal thermal state of the Greenland Ice Sheet

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Journal of Geophysical Research F: Earth Surface↗

Removing sun glint from optical remote sensing images of shallow rivers

Sun glint is the specular reflection of light from the water surface, which often causes unusually bright pixel values that can dominate fluvial remote sensing imagery and obscure the water-leaving radiance signal of interest for mapping bathymetry, bottom type, or water column optical characteristics. Although sun glint is ubiquitous in fluvial remote sensing imagery, river-specific methods for removing sun glint are not yet available. We show that existing sun glint-removal methods developed for multispectral images of marine shallow water environments over-correct shallow portions of fluvial remote sensing imagery resulting in regions of unreliable data along channel margins. We build on existing marine glint-removal methods to develop a river-specific technique that removes sun glint from shallow areas of the channel without overcorrection by accounting for non-negligible water-leaving near-infrared radiance. This new sun glint-removal method can improve the accuracy of spectrally-based depth retrieval in cases where sun glint dominates the at-sensor radiance. For an example image of the gravel-bed Snake River, Wyoming, USA, observed-vs.-predicted R 2 values for depth retrieval improved from 0.66 to 0.76 following sun glint removal. The methodology presented here is straightforward to implement and could be incorporated into image processing workflows for multispectral images that include a near-infrared band.

Earth Surface Processes and Landforms↗