SEARCH · USGS Science

Results for “Geophysics”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,549 records · Page 86Linked to original sources

USGS2014

Report of committee on relation of inch and meter

Those of you who attended the annual meeting of this Section on April 27, 1938, heard a paper entitled “A method for introducing a new standard of length” that was presented by Professor Philip Klssara, calling attention to the Bill then in Congress proposing to redefine the length of the inch. The paper has been published in the “Transactions of the American Geophysical Union“ [19th Annual Meeting, 1938, pp. 94–96]. Its closing paragraph is as follows: “As the American Geophysical Union has as one of its functions the duty of offering advice to the Congress on technical matters, I move, Mr. Chairman, that the Chairman be empowered to appoint a committee to study the matter and to report what recommendations should be made to Congress by the Union.”

Eos, Transactions, American Geophysical Union↗

USGS2014

Discussion of “The relation of geology to dry weather stream flow in Ohio”

The concepts presented in this paper are of great value to the ground‐water hydrologist. They indicate one way to analyze and conveniently use the vast accumulation of stream‐flow records collected by governmental agencies as a tool for geophysical reconnaisance. To be usable as a method of geophysical prospecting for ground water, stream‐flow records must first be transformed into flow‐duration curves. Then, by means of the 90 per cent index, chosen by Cross, or other, more suitable indices, the ground‐water storage characteristics of the surficial deposits of drainage basins can readily be determined. Where there is good indication of extensive groundwater storage, it is probable that large supplies of ground water can be developed for use.

Ohio↗

USGS2025

Infrasonic directivity of monopole, dipole, and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local ( ⁠< 5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source ( ⁠< 5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Geophysical Journal International↗

USGS2010

Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity andinduced-polarization data

Electrical geophysical methods have found wide use in the growing discipline of hydrogeophysics for characterizing the electrical properties of the subsurface and for monitoring subsurface processes in terms of the spatiotemporal changes in subsurface conductivity, chargeability, and source currents they govern. Presently, multichannel and multielectrode data collections systems can collect large data sets in relatively short periods of time. Practitioners, however, often are unable to fully utilize these large data sets and the information they contain because of standard desktop-computer processing limitations. These limitations can be addressed by utilizing the storage and processing capabilities of parallel computing environments. We have developed a parallel distributed-memory forward and inverse modeling algorithm for analyzing resistivity and time-domain induced polar-ization (IP) data. The primary components of the parallel computations include distributed computation of the pole solutions in forward mode, distributed storage and computation of the Jacobian matrix in inverse mode, and parallel execution of the inverse equation solver. We have tested the corresponding parallel code in three efforts: (1) resistivity characterization of the Hanford 300 Area Integrated Field Research Challenge site in Hanford, Washington, U.S.A., (2) resistivity characterization of a volcanic island in the southern Tyrrhenian Sea in Italy, and (3) resistivity and IP monitoring of biostimulation at a Superfund site in Brandywine, Maryland, U.S.A. Inverse analysis of each of these data sets would be limited or impossible in a standard serial computing environment, which underscores the need for parallel high-performance computing to fully utilize the potential of electrical geophysical methods in hydrogeophysical applications.

Geophysics↗

USGS2011

Quantifying solute transport processes: Are chemically "conservative" tracers electrically conservative?

The concept of a nonreactive or conservative tracer, commonly invoked in investigations of solute transport, requires additional study in the context of electrical geophysical monitoring. Tracers that are commonly considered conservative may undergo reactive processes, such as ion exchange, thus changing the aqueous composition of the system. As a result, the measured electrical conductivity may reflect not only solute transport but also reactive processes. We have evaluated the impacts of ion exchange reactions, rate-limited mass transfer, and surface conduction on quantifying tracer mass, mean arrival time, and temporal variance in laboratory-scale column experiments. Numerical examples showed that (1) ion exchange can lead to resistivity-estimated tracer mass, velocity, and dispersivity that may be inaccurate; (2) mass transfer leads to an overestimate in the mobile tracer mass and an underestimate in velocity when using electrical methods; and (3) surface conductance does not notably affect estimated moments when high-concentration tracers are used, although this phenomenon may be important at low concentrations or in sediments with high and/or spatially variable cation-exchange capacity. In all cases, colocated groundwater concentration measurements are of high importance for interpreting geophysical data with respect to the controlling transport processes of interest.

Geophysics↗

USGS2013

Application and evaluation of electromagnetic methods for imaging saltwater intrusion in coastal aquifers: Seaside Groundwater Basin, California

Developing effective resource management strategies to limit or prevent saltwater intrusion as a result of increasing demands on coastal groundwater resources requires reliable information about the geologic structure and hydrologic state of an aquifer system. A common strategy for acquiring such information is to drill sentinel wells near the coast to monitor changes in water salinity with time. However, installation and operation of sentinel wells is costly and provides limited spatial coverage. We studied the use of noninvasive electromagnetic (EM) geophysical methods as an alternative to installation of monitoring wells for characterizing coastal aquifers. We tested the feasibility of using EM methods at a field site in northern California to identify the potential for and/or presence of hydraulic communication between an unconfined saline aquifer and a confined freshwater aquifer. One-dimensional soundings were acquired using the time-domain electromagnetic (TDEM) and audiomagnetotelluric (AMT) methods. We compared inverted resistivity models of TDEM and AMT data obtained from several inversion algorithms. We found that multiple interpretations of inverted models can be supported by the same data set, but that there were consistencies between all data sets and inversion algorithms. Results from all collected data sets suggested that EM methods are capable of reliably identifying a saltwater-saturated zone in the unconfined aquifer. Geophysical data indicated that the impermeable clay between aquifers may be more continuous than is supported by current models.

California↗

USGS2011

Hydrogeophysical investigations at Hidden Dam, Raymond, California

Self-potential and direct current resistivity surveys are carried out at the Hidden Dam site in Raymond, California to assess present-day seepage patterns and better understand the hydrogeologic mechanisms that likely influence seepage. Numerical modeling is utilized in conjunction with the geophysical measurements to predict variably-saturated flow through typical two-dimensional dam cross-sections as a function of reservoir elevation. Several different flow scenarios are investigated based on the known hydrogeology, as well as information about typical subsurface structures gained from the resistivity survey. The flow models are also used to simulate the bulk electrical resistivity in the subsurface under varying saturation conditions, as well as the self-potential response using petrophysical relationships and electrokinetic coupling equations. The self-potential survey consists of 512 measurements on the downstream area of the dam, and corroborates known seepage areas on the northwest side of the dam. Two direct-current resistivity profiles, each approximately 2,500 ft (762 m) long, indicate a broad sediment channel under the northwest side of the dam, which may be a significant seepage pathway through the foundation. A focusing of seepage in low-topography areas downstream of the dam is confirmed from the numerical flow simulations, which is also consistent with past observations. Little evidence of seepage is identified from the self-potential data on the southeast side of the dam, also consistent with historical records, though one possible area of focused seepage is identified near the outlet works. Integration of the geophysical surveys, numerical modeling, and observation well data provides a framework for better understanding seepage at the site through a combined hydrogeophysical approach.

California↗

USGS2021

Ground‐penetrating radar, electromagnetic induction, terrain, and vegetation observations coupled with machine learning to map permafrost distribution at Twelvemile Lake, Alaska

We collected ground‐penetrating radar (GPR) and frequency‐domain electromagnetic induction (FDEM) profiles in 2011 and 2012 to identify the extent of permafrost relative to surface biomass and solar insolation around Twelvemile Lake near Fort Yukon, Alaska. We compared a Landsat‐derived biomass estimate and modeled solar insolation from a digital elevation model to the geophysical measurements. We show correspondence between vegetation type and biomass relative to permafrost extent and seasonal freeze–thaw. Thicker permafrost (≥25 m) was covered by greater biomass, and seasonal thaw depths in these regions were minimal (1 m). Shallow (1–3 m depth) and thin (20–50 cm) newly forming permafrost or frozen layers from the previous winter occurred below northward oriented slopes with thin biomass cover. South‐facing slopes exhibited permafrost when there was enough biomass to shield incoming solar energy. We developed an artificial neural network to predict permafrost extent across the broader region by mapping GPR‐observed instances of permafrost to FDEM, biomass, and terrain observations with 90.2% accuracy. We identified a strong linear correlation ( r = −0.77) between permafrost probability and seasonal thaw depth, indicating that our models may also be used to explore thaw patterns and variability in active layer thickness. This study highlights the combined influence of biomass and terrain on the presence of permafrost and the value of evaluating such parameters via remote sensing to predict permafrost spatial or temporal variability. Incorporating diverse geophysical datasets with in‐situ validation into machine learning models demonstrates a useful approach to upscale estimated permafrost extent across large Arctic expanses.

Alaska↗

USGS2008

Information fusion in regularized inversion of tomographic pumping tests

In this chapter we investigate a simple approach to incorporating geophysical information into the analysis of tomographic pumping tests for characterization of the hydraulic conductivity (K) field in an aquifer. A number of authors have suggested a tomographic approach to the analysis of hydraulic tests in aquifers - essentially simultaneous analysis of multiple tests or stresses on the flow system - in order to improve the resolution of the estimated parameter fields. However, even with a large amount of hydraulic data in hand, the inverse problem is still plagued by non-uniqueness and ill-conditioning and the parameter space for the inversion needs to be constrained in some sensible fashion in order to obtain plausible estimates of aquifer properties. For seismic and radar tomography problems, the parameter space is often constrained through the application of regularization terms that impose penalties on deviations of the estimated parameters from a prior or background model, with the tradeoff between data fit and model norm explored through systematic analysis of results for different levels of weighting on the regularization terms. In this study we apply systematic regularized inversion to analysis of tomographic pumping tests in an alluvial aquifer, taking advantage of the steady-shape flow regime exhibited in these tests to expedite the inversion process. In addition, we explore the possibility of incorporating geophysical information into the inversion through a regularization term relating the estimated K distribution to ground penetrating radar velocity and attenuation distributions through a smoothing spline model. ?? 2008 Springer-Verlag Berlin Heidelberg.

Studies in Computational Intelligence↗

USGS1993

Estimation of undiscovered deposits in quantitative mineral resource assessments-examples from Venezuela and Puerto Rico

Quantitative mineral resource assessments used by the United States Geological Survey are based on deposit models. These assessments consist of three parts: (1) selecting appropriate deposit models and delineating on maps areas permissive for each type of deposit; (2) constructing a grade-tonnage model for each deposit model; and (3) estimating the number of undiscovered deposits of each type. In this article, I focus on the estimation of undiscovered deposits using two methods: the deposit density method and the target counting method. In the deposit density method, estimates are made by analogy with well-explored areas that are geologically similar to the study area and that contain a known density of deposits per unit area. The deposit density method is useful for regions where there is little or no data. This method was used to estimate undiscovered low-sulfide gold-quartz vein deposits in Venezuela. Estimates can also be made by counting targets such as mineral occurrences, geophysical or geochemical anomalies, or exploration "plays" and by assigning to each target a probability that it represents an undiscovered deposit that is a member of the grade-tonnage distribution. This method is useful in areas where detailed geological, geophysical, geochemical, and mineral occurrence data exist. Using this method, porphyry copper-gold deposits were estimated in Puerto Rico. ?? 1993 Oxford University Press.

Nonrenewable Resources↗

USGS2017

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, 1 / f α "> 1 / f α 1/fα with frequency, f . With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013 . doi: 10.1007/s00190-012-0605-0 ) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.

Journal of Geodesy↗

USGS2016

End of the chain? Rugosity and fine-scale bathymetry from existing underwater digital imagery using structure-from-motion (SfM) technology

The rugosity or complexity of the seafloor has been shown to be an important ecological parameter for fish, algae, and corals. Historically, rugosity has been measured either using simple and subjective manual methods such as ‘chain-and-tape’ or complicated and expensive geophysical methods. Here, we demonstrate the application of structure-from-motion (SfM) photogrammetry to generate high-resolution, three-dimensional bathymetric models of a fringing reef from existing underwater video collected to characterize the seafloor. SfM techniques are capable of achieving spatial resolution that can be orders of magnitude greater than large-scale lidar and sonar mapping of coral reef ecosystems. The resulting data provide finer-scale measurements of bathymetry and rugosity that are more applicable to ecological studies of coral reefs than provided by the more expensive and time-consuming geophysical methods. Utilizing SfM techniques for characterizing the benthic habitat proved to be more effective and quantitatively powerful than conventional methods and thus might portend the end of the ‘chain-and-tape’ method for measuring benthic complexity.

Coral Reefs↗

USGS2009

Holocene evolution of Apalachicola Bay, Florida

A program of geophysical mapping and vibracoring was conducted to better understand the geologic evolution of Apalachicola Bay. Analyses of the geophysical data and sediment cores along with age control provided by 34 AMS 14 C dates on marine shells and wood reveal the following history. As sea level rose in the early Holocene, fluvial deposits filled the Apalachicola River paleochannel, which extended southward under the central part of the bay and seaward across the continental shelf. Sediments to either side of the paleochannel contain abundant wood fragments, with dates documenting that those areas were forested at 8,000 14 C years B.P . As sea level continued to rise, spits formed of headland prodelta deposits. Between ∼6,400 and ∼2,500 14 C years B.P ., an Apalachicola prodelta prograded and receded several times across the inner shelf that underlies the western part of the bay. An eastern deltaic lobe was active for a shorter time, between ∼5,800 and 5,100 14 C years B.P . Estuarine benthic foraminiferal assemblages occurred in the western bay as early as 6,400 14 C years B.P ., and indicate that there was some physical barrier to open-ocean circulation and shelf species established by that time. It is considered that shoals formed in the region of the present barrier islands as the rising sea flooded an interstream divide. Estuarine conditions were established very early in the post-glacial flooding of the bay.

Florida↗

USGS2020

Constraints on eruption processes and event masses for the 2016–2017 eruption of Bogoslof volcano, Alaska, through evaluation of IASI satellite SO2 masses and complementary datasets

Bogoslof volcano, Alaska, experienced at least 70 explosive eruptions between 12 December 2016 and 31 August 2017. Due to its remote location and limited local monitoring network, this eruption was monitored and characterized primarily using remote geophysical and satellite techniques. SO 2 emissions from Bogoslof were persistently detected by the Infrared Atmospheric Sounding Interferometer (IASI) satellite sensors. Of Bogoslof’s 70 explosive events, 50% produced measurable SO 2 masses ranging from 0.1 to 21.5 kt, with a median and standard deviation of 0.7 ± 4.0 kt SO 2 , respectively. Here, we compare IASI-derived SO 2 masses from Bogoslof events to complementary geophysical datasets to provide insights into eruption source processes, namely the degree of seawater scrubbing of water-soluble SO 2 and variations in magma flux. Correlations with the number of lightning strokes and infrasound energy are expected to indicate magma-flux as a controlling process, while correlations with infrasound frequency index are expected to indicate variations in vent-water content as a controlling factor. These comparisons suggest that the measured SO 2 masses are primarily a function of eruption magnitude (degassed magma mass) and that scrubbing of SO 2 emissions by vent seawater may have exerted a minor effect on the observed SO 2 masses. SO 2 masses were combined with petrologic constraints on melt inclusion and matrix glass S concentrations to calculate degassed magma masses and volumes. The cumulative SO 2 -derived degassed magma mass and estimated volume (dense-rock equivalent) for the full Bogoslof eruption were found to be 2.8 × 10 10 kg and 9.3 × 10 6 m 3 , respectively. When individual event masses are compared against event masses calculated using an empirical plume-height method, a strong correlation is found ( R 2 = 0.83), with better than order-of-magnitude agreement in most cases. These estimates of eruption masses provide useful information on the magnitude, behavior, and associated hazards of the 2016–2017 eruption, and potentially future unrest at Bogoslof volcano.

Alaska↗

USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Hawaii↗

USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part I: Lavas and tephras

There are few petrological constraints on magma storage depths at Mauna Loa, Hawai‘i. Yet understanding the geometry of the magmatic plumbing system is critical for interpreting geophysical signals of unrest at this very high-threat volcano. We address this gap by examining CO₂-rich fluid inclusions (FI) in lava and tephra from seven eruptions (8870 ± 56 14C yr BP, 1852, 1855, 1868, 1949, 1950, and 1984), supplemented with published data from 2022. Raman spectroscopy was used to determine FI densities, from which entrapment pressures were calculated using a CO₂-H₂O equation of state. Most FI record pressures of ~ 0.25–1.25 kbar (~ 2–5 km depth below the summit), consistent with geophysical estimates from the past 40 years. In summit eruptions, FI hosted in more evolved olivine and orthopyroxene clots (Fo and Mg# < 84) record slightly shallower pressures than those in more primitive olivines (Fo > 84) from rift zone eruptions, suggesting a crystal-poor evolved cap near the top of the reservoir (~ 2 km). The similarity in storage depths across all eight eruptions indicates that Mauna Loa’s magmas have tapped a quasi-stable reservoir over the past two centuries, and possibly over 10 kyr. Electron backscatter diffraction reveals deformations to the crystal lattice in Fo82-83 olivines, likely due to deformation during storage in mush piles. The intensity of deformation is comparable to that seen at Kīlauea, implying that mush pile stress may be decoupled from edifice size or longevity. Finally, SO₂ contents in FI increase from ~ 2 mol% at 2 kbar to ~ 15 mol% at 0.5 kbar, suggesting sulphur degassing begins far deeper than the 0.2–0.3 kbar commonly assumed for Hawaiian systems. This validates the newest generation of S degassing models (e.g., Sulfur_X), and explains precursory SO2 emissions in the ~3 hours prior to the onset of the 2022 eruption (Esse et al. 2025).

Hawaii↗

USGS2025

Utilizing crystal records to determine the underlying architecture and timescales of magmatic processes during the 2022 Mauna Loa eruption

Details of Mauna Loa’s plumbing system and magma storage dynamics are relatively poorly constrained, relative to nearby Kīlauea volcano. Recent studies have begun to address this gap in knowledge using geophysical imaging or olivine-focused petrology. Here, we use the crystal cargo from the 2022 summit and Northeast Rift Zone eruption to advance our understanding of Mauna Loa’s storage regions and pre˗eruptive processes. We identify textural and geochemical populations of olivine, orthopyroxene, clinopyroxene, and plagioclase macrocrysts as well as evolved glomerocrysts consisting of orthopyroxene, clinopyroxene, and plagioclase. Summit versus rift zone trends can be distinguished between specific crystal populations and are used to infer magma transfer pathways. Melt inclusion barometry from three orthopyroxene populations indicate crystallization depths of ~2–3 km which coincide with previous olivine˗hosted fluid inclusion barometry across Mauna Loa’s recent 200˗yr eruptive history. We combine orthopyroxene˗hosted melt inclusion barometry with diffusion chronometry on the same crystals using Fe˗Mg interdiffusion to estimate timescales of mixing and recharge prior to the eruption. Summit orthopyroxene yield timescales of weeks to months, aligning with geophysically recognized pre˗eruptive intrusion events and previously determined olivine timescales. Northeast Rift Zone orthopyroxene crystals have resorbed cores and yield short pre˗eruptive mixing timescales on the order of hours to days. Chemistry and timescales from the orthopyroxene populations suggest some of the 2022 magmas may have bypassed the shallowest portions of the summit reservoir system and were transferred directly into the rift zone syn˗eruptively.

Hawaii↗

USGS2019

Hydrogeology of an alpine talus aquifer: Cordillera Blanca, Peru

The dramatic loss of glacial mass in low latitudes is causing shifts in downstream water availability and use during the driest months of the year. The world’s largest concentration of tropical glaciers lies in the Cordillera Blanca range of Peru, where glacial runoff is declining and regional stresses are emerging over water resources. Throughout the Cordillera Blanca, groundwater inputs from alpine meadow–talus complexes, locally known as pampas, supply proglacial streams with up to 80% of their flow during the region’s dry season. Structural knowledge of the pampa aquifers is needed to estimate their drainable groundwater storage capacity and residence time, to elucidate the role and importance of alpine groundwater storage in the regional water budget of the Cordillera Blanca. To understand the structure of these proglacial aquifers, multiple near-surface geophysical methods were implemented in a proglacial valley near dense networks of spring-fed tributaries. Geophysical results and borehole logs suggest groundwater is stored in a confined aquifer composed of buried talus deposits overlain by lacustrine clay, while deeper portions of the unit, 10–15 m in depth, are relatively clay-free and more hydraulically conductive. Based on these findings and assumptions of aquifer porosity, the pampas of the Callejon de Huaylas may store from 0.006 to 0.02 km 3 of groundwater. Furthermore, these findings suggest that the talus aquifers of the Cordillera Blanca were formed in proglacial lakes, followed by infilling with fine lacustrine sediments that confine lower units and allow for groundwater discharge to springs via macropores and preferential flow.

Hydrogeology Journal↗