Geochemistry of waters in the Valley of Ten Thousand Smokes region, Alaska

Meteoric waters from cold springs and streams outside of the 1912 eruptive deposits filling the Valley of Ten Thousand Smokes (VTTS) and in the upper parts of the two major rivers draining the 1912 deposits have similar chemical trends. Thermal springs issue in the mid-valley area along a 300-m lateral section of ash-flow tuff, and range in temperature from 21 to 29.8°C in early summer and from 15 to 17°C in mid-summer. Concentrations of major and minor chemical constituents in the thermal waters are nearly identical regardless of temperature. Waters in the downvalley parts of the rivers draining the 1912 deposits are mainly mixtures of cold meteoric waters and thermal waters of which the mid-valley thermal spring waters are representative. The weathering reactions of cold waters with the 1912 deposits appear to have stabilized and add only subordinate amounts of chemical constituents to the rivers relative to those contributed by the thermal waters. Isotopic data indicate that the mid-valley thermal spring waters are meteoric, but data is inconclusive regarding the heat source. The thermal waters could be either from a shallow part of a hydrothermal system beneath the 1912 vent region or from an incompletely cooled, welded tuff lens deep in the 1912 ash-flow sheet of the upper River Lethe area.Bicarbonate-sulfate waters resulting from interaction of near-surface waters and the cooling 1953–1968 southwest Trident plug issue from thermal springs south of Katmai Pass and near Mageik Creek, although the Mageik Creek spring waters are from a well-established, more deeply circulating hydrothermal system. Katmai caldera lake waters are a result of acid gases from vigorous drowned fumaroles dissolving in lake waters composed of snowmelt and precipitation.     

Helium isotope geochemistry of mid-ocean ridge basalts from the South Atlantic

We report new helium isotope results for 49 basalt glass samples from the Mid-Atlantic Ridge between 1°N and 47°S.³He/⁴He in South Atlantic mid-ocean ridge basalts (MORB) varies between 6.5 and 9.0 R_A (R_A is the atmospheric ratio of1.39 × 10⁻⁶), encompassing the range of previously reported values for MORB erupted away from high³He/⁴He hotspots such as Iceland. He, Sr and Pb isotopes show systematic relationships along the ridge axis. The ridge axis is segmented with respect to geochemical variations, and local spike-like anomalies in³He/⁴He, Pb and Sr isotopes, and trace element ratios such as(La/Sm)_N are prevalent at the latitudes of the islands of St. Helena, Tristan da Cunha and Gough to the east of the ridge. The isotope systematics are consistent with injection beneath the ridge of mantle “blobs” enriched in radiogenic He, Pb and Sr, derived from off-axis hotspot sources. The variability in³He/⁴He along the ridge can be used to refine the hotspot source-migrating-ridge sink model.

MORB from the 2–7°S segment are systematically the least radiogenic samples found along the mid-ocean ridge system to date. Here the depleted mantle source is characterized by⁸⁷Sr/⁸⁶Sr of 0.7022, Pb isotopes close to the geochron and with²⁰⁶Pb/²⁰⁴Pb of 17.7, and³He/⁴He of 8.6–8.9 R_A. The “background contamination” of the subridge mantle, by radiogenic helium derived from off-ridge hotspots, displays a maximum between 20 and 24°S. The HePb and HeSr isotope relations along the ridge indicate that the³He/⁴He ratios are lower for the hotspot sources of St. Helena, Tristan da Cunha and Gough than for the MORB source, consistent with direct measurements of³He/⁴He ratios in the island lavas. Details of the HeSrPb isotope systematics between 12 and 22°S are consistent with early, widespread dispersion of the St. Helena plume into the asthenosphere, probably during flattening of the plume head beneath the thick lithosphere prior to continental breakup. The geographical variation in theHe/Pbratio deduced from the isotope systematics suggests only minor degassing of the plume during this stage. Subsequently, it appears that the plume component reaching the mid-Atlantic ridge was partially outgassed of He during off-ridge hotspot volcanism and related melting activity.

Overall, the similar behavior of He and Pb isotopes along the ridge indicates that the respective mantle sources have evolved under conditions which produced related He and Pb isotope variations.     

Fracture in Westerly Granite under AE Feedback and Constant Strain Rate Loading: Nucleation, Quasi-static Propagation, and the Transition to Unstable Fracture Propagation

New observations of fracture nucleation are presented from three triaxial compression experiments on intact samples of Westerly granite, using Acoustic Emission (AE) monitoring. By conducting the tests under different loading conditions, the fracture process is demonstrated for quasi-static fracture (under AE Feedback load), a slowly developing unstable fracture (loaded at a `slow' constant strain rate of 2.5 × 10<sup>−6</sup> /s) and an unstable fracture that develops near instantaneously (loaded at a `fast' constant strain rate of 5 × 10⁻⁵ /s). By recording a continuous ultrasonic waveform during the critical period of fracture, the entire AE catalogue can be captured and the exact time of fracture defined. Under constant strain loading, three stages are observed: (1) An initial nucleation or stable growth phase at a rate of ~ 1.3 mm/s, (2) a sudden increase to a constant or slowly accelerating propagation speed of ~ 18 mm/s, and (3) unstable, accelerating propagation. In the ~ 100 ms before rupture, the high level of AE activity (as seen on the continuous record) prevented the location of discrete AE events. A lower bound estimate of the average propagation velocity (using the time-to-rupture and the existing fracture length) suggests values of a few m/s. However from a low gain acoustic record, we infer that in the final few ms, the fracture propagation speed increased to 175 m/s. These results demonstrate similarities between fracture nucleation in intact rock and the nucleation of dynamic instabilities in stick slip experiments. It is suggested that the ability to constrain the size of an evolving fracture provides a crucial tool in further understanding the controls on fracture nucleation.     

Noble gases in the earth's mantle

The abundances and isotopic compositions of noble gases in two samples from ultramafic xenoliths in alkali basalt, a young kaersutitic amphibole separated from a peridotite xenolith from Dish Hill, California and an ancient whole-rock lherzolite xenolith from Baja California, are reported and compared with the results of analyses on other mantle samples. In addition to previously recognized excesses of ³He and ¹²⁹Xe, our results indicate that ambient gases in the mantle show a general enrichment of the lighter-mass nonradiogenic isotopes of Ar, Kr and Xe, and Ar with ⁴⁰Ar/³⁶Ar = 3 · 10².     

Under the surface: Pressure-induced planetary-scale waves,volcanic lightning,and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha'apai volcano

We present a narrative of the eruptive events culminating in the cataclysmic January 15, 2022 eruption of Hunga Tonga-Hunga Ha'apai Volcano by synthesizing diverse preliminary seismic, volcanological, sound wave, and lightning data available within the first few weeks after the eruption occurred. The first hour of eruptive activity produced fast-propagating tsunami waves, long-period seismic waves, loud audible sound waves, infrasonic waves, exceptionally intense volcanic lightning and an unsteady volcanic plume that transiently reached—at 58 ?km—the Earth's mesosphere. Energetic seismic signals were recorded worldwide and the globally stacked seismogram showed episodic seismic events within the most intense periods of phreatoplinian activity, and they correlated well with the infrasound pressure waveform recorded in Fiji. Gravity wave signals were strong enough to be observed over the entire planet in just the first few hours, with some circling the Earth multiple times subsequently. These large-amplitude, long-wavelength atmospheric disturbances come from the Earth's atmosphere being forced by the magmatic mixture of tephra, melt and gasses emitted by the unsteady but quasi-continuous eruption from 0402±1–1800 UTC on January 15, 2022. Atmospheric forcing lasted much longer than rupturing from large earthquakes recorded on modern instruments, producing a type of shock wave that originated from the interaction between compressed air and ambient (wavy) sea surface. This scenario differs from conventional ideas of earthquake slip, landslides, or caldera collapse-generated tsunami waves because of the enormous (～1000x) volumetric change due to the supercritical nature of volatiles associated with the hot, volatile-rich phreatoplinian plume. The time series of plume altitude can be translated to volumetric discharge and mass flow rate. For an eruption duration of ～12 ?h, the eruptive volume and mass are estimated at 1.9 ?km³ and ～2 900 ?Tg, respectively, corresponding to a VEI of 5–6 for this event. The high frequency and intensity of lightning was enhanced by the production of fine ash due to magma—seawater interaction with concomitant high charge per unit mass and the high pre-eruptive concentration of dissolved volatiles. Analysis of lightning flash frequencies provides a rapid metric for plume activity and eruption magnitude. Many aspects of this eruption await further investigation by multidisciplinary teams. It represents a unique opportunity for fundamental research regarding the complex, non-linear behavior of high energetic volcanic eruptions and attendant phenomena, with critical implications for hazard mitigation, volcano forecasting, and first-response efforts in future disasters.     

δ18O record and temperature change over the past 100 years in ice cores on the Tibetan Plateau

The 213 m ice core from the Puruogangri Ice Field on the Tibetan Plateau facilitates the study of the regional temperature changes with its δ18O record of the past 100 years. Here we combine information from this core with that from the Dasuopu ice core (from the southern Tibetan Plateau), the Guliya ice core (from the northwestern Plateau) and the Dunde ice core (from the northeastern Plateau) to learn about the regional differences in temperature change across the Tibetan Plateau. The δ18O changes vary with region on the Plateau, the variations being especially large between South and North and between East and West. Moreover, these four ice cores present increasing δ18O trends, indicating warming on the Tibetan Plateau over the past 100 years. A comparative study of Northern Hemisphere (NH) temperature changes, the δ18O-reflected temperature changes on the Plateau, and available meteorological records show consistent trends in overall warming during the past 100 years.     

Altitude variation of the plasmapause signature in the main ionospheric trough

The projection of the plasmapause magnetic-field lines to low altitudes, where the light-ion chemistry is dominated by O⁺, tends to occur near the minimum electron density in the main (midlatitude) electron density trough at night. With increasing altitude in the trough, where H⁺ emerges as the dominant ion on the low-latitude boundary, we have found cases where the plasmapause field lines are located on the sharp low-latitude side of the trough as expected if this topside ionosphere H⁺ distribution varies in step with the plasmapause gradient in the distant plasmasphere. These conclusions are based on near-equatorial crossings of the plasmapause (corresponding to the steep gradient in the dominant species H⁺) by the Explorer-45 satellite as determined from electric-field measurements by Maynard and Cauffman in the early 1970s and ISIS-2 ionospheric topside-sounder measurements. The former data have now been converted to digital form and made available at http://nssdcftp.gsfc.nasa.gov. The latter provide samples of nearly coincident observations of ionospheric main trough crossings near the same magnetic-field lines of the Explorer 45-determined equatorial plasmapause. The ISIS-2 vertical electron density profiles are used to infer where the F-region transitions from an O⁺ to a H⁺ dominated plasma through the main trough boundaries.     

An evaluation of the health status of the Lavaca Bay, Texas ecosystem using Crassostrea virginica as the sentinel species

Locational risks for compromised ecosystem health for the eastern oysters (Crassostrea virginica) harvested from Lavaca Bay, Texas were estimated. Flow cytometric evaluation of variations in DNA content and the lysosomal destabilization assay were used for evaluation of genotoxicity and stress, respectively. Bayesian geo-statistical methods were utilized to estimate and evaluate spatial effects. For models with spatial risks, continuous surface maps of predicted parameter values were created to evaluate risk location. Lysosomal destabilization assay results were spatially oriented whereas flow cytometry results were fit best with the random effects model. While not spatially oriented, the highest levels of variations in DNA content were also present near industrial facilities. Locational risks of increased biomarkers of genotoxicity and stress in the eastern oyster (Crassostrea virginica) were increased with proximity to industrial facilities     

Improving at-site flood frequency analysis with additional spatial information: a probabilistic regional envelope curve approach

Extreme flood events have detrimental effects on society, the economy and the environment. Widespread flooding across South East Queensland in 2011 and 2013 resulted in the loss of lives and significant cost to the economy. In this region, flood risk planning and the use of traditional flood frequency analysis (FFA) to estimate both the magnitude and frequency of the 1-in-100 year flood is severely limited by short gauging station records. On average, these records are 42 years in Eastern Australia and many have a poor representation of extreme flood events. The major aim of this study is to test the application of an alternative method to estimate flood frequency in the form of the Probabilistic Regional Envelope Curve (PREC) approach which integrates additional spatial information of extreme flood events. In order to better define and constrain a working definition of an extreme flood, an Australian Envelope Curve is also produced from available gauging station data. Results indicate that the PREC method shows significant changes to the larger recurrence intervals (≥100 years) in gauges with either too few, or too many, extreme flood events. A decision making process is provided to ascertain when this method is preferable for FFA.     

Urban stormwater thermal gene expression models for protection of sensitive receiving streams

Thermal impact of typical high‐density residential, industrial, and commercial land uses is a major concern for the health of aquatic life in urban watersheds, especially in smaller, cold, and cool‐water streams. This is the first study of its kind that provides simple easy‐to‐use equations, developed using gene expression programming (GEP) that can guide the assessment and the design of urban stormwater management systems to protect thermally sensitive receiving streams. We developed 3 GEP models using data collected during 3 years (2009–2011) from 4 urban catchments; the first GEP model predicts event mean temperature at the inlet of the pond; the second model predicts the stormwater temperature at the outlet of the pond; and the third model predicts the temperature of the stormwater after flowing through a cooling trench and before discharging to the receiving stream. The new models have high correlation coefficients of 0.90–0.94 and low prediction uncertainty of less than 4% of the median value of the predicted runoff temperatures. Sensitivity analysis shows that climatic factors have the highest influence on the thermal enrichment followed by the catchment characteristics and the key design variables of the stormwater pond and the cooling trench. The general method presented here is easily transferable to other regions of the world (but not necessarily the exact equations developed here); also through sensitivity and parametric analysis, we gained insight on the key factors and their relative importance in modelling thermal enrichment of urban stromwater runoff.