Search

Study’s Novelty/Excerpt • This study investigates the prevalence of Enterobacteriaceae in clinical samples from the Federal Teaching Hospital Gombe, providing a comprehensive analysis of 420 non-duplicate isolates over a four-month period. • The research uniquely identifies the distribution of various Enterobacteriaceae species, with Escherichia coli being the most prevalent, and correlates these findings with the types of clinical specimens, revealing urine as the most common source. • This work underscores the necessity for further research into the antimicrobial resistance profiles of these isolates, aiming to inform more effective treatment strategies for bacterial infections in the region. Full Abstract Enterobacteriaceae is a family of Gram-negative, oxidase-negative, and catalase-positive bacteria mostly found in Humans and animals' intestines. Some of these organisms are enteric opportunistic pathogens associated with urinary tract infections, respiratory tract infections, and wound infections, whereas others are regularly pathogenic for humans. This study aimed to determine the prevalence of Enterobacteriaceae in clinical samples from Federal Teaching Hospital Gombe between August 2022 and November 2022. A total of 420 non-duplicate isolates from various clinical samples were analyzed in the study. The isolates were identified based on cultural characteristics, Gram staining, and standard biochemical tests. Out of the 420 isolates identified, Escherichia coli was the most prevalent with 163(38.8%) isolates, followed by Klebsiella pneumoniae with 69(16.4%), Klebsiella oxytoca with 61(14.5%), Proteus spp. 28(6.7%), Citrobacter spp. 27(6.4%), Enterobacter spp. 22(5.2%), Serratia marcescens 13(3.1%), Providencia spp. 12(2.9%), Yersinia enterocolitica 11(2.6%), Morganella morganii 3(0.7%) and Salmonella spp. with 2(0.5%). Based on clinical specimens, urine had the highest percentage of isolates with 53.4%, followed by wound swab (19.1%), stool (10.1%), High vaginal swab (6.8%), Endocervical swab (3.6%), sputum (3.1%), blood (1.4%), Cerebrospinal fluid (1.0%) and semen (1.0%) and then pleural fluid (0.5%). In conclusion, Enterobacteriaceae clinical isolates were highly prevalent in Federal Teaching Hospital Gombe. Further research to assess the antimicrobial resistance profile of these clinical bacterial organisms in the study area is recommended for effective treatment options for bacterial infections.

The seismometer deployed by the InSight lander measured the seismic velocity of the Martian crust. We use a rock physics model to interpret those velocities and constrain hydrogeological properties. The seismic velocity of the upper ∼10 km is too low to be ice-saturated. Hence there is no cryosphere that confines deeper aquifers and possibly no aquifers locally. An increase in seismic velocity at depths of ∼10 km could be explained by a few volume percent of mineral cement (1%–5%) in pore space and may document the past depth of aquifers.

Tidal interactions between planets or stars and the bodies that orbit them dissipate energy in their interiors. The dissipated energy heats the interior and a fraction of that energy will be released as seismic energy. Here we formalize a model to describe the tidally-driven seismic activity on planetary bodies based on tidal dissipation. To constrain the parameters of our model we use the seismic activity of the Moon, driven by tidal dissipation from the Earth-Moon interactions. We then apply this model to predict the amount of seismic energy release and largest seismic events on other moons in our Solar System and exoplanetary bodies. We find that many moons in the Solar System should be more seismically active than the Earth's Moon and many exoplanets should exhibit more seismic activity than the Earth. Finally, we examine how temporal-spatial variations in tidal dissipation manifest as variations in the locations and timing of seismic events on these bodies.

Active eruptions from the south polar region of Saturn’s ~500-km-diameter moon Enceladus are concentrated along a series of lineaments known as the ‘tiger stripes’1,2, thought to be partially open fissures that connect to the liquid water ocean beneath the ice shell3,4. To date, no study simultaneously explains why the tiger stripes should be located only at the south pole, why there are multiple approximately parallel and regularly spaced fractures, what accounts for their spacing of about 35 km, and why similarly active fissures have not been observed on other icy bodies. Here we propose that secular cooling, which leads to a thickening of the ice shell and building of global tensile stresses5,6, causes the first fracture to form at one of the poles, where the ice shell is thinnest owing to tidal heating7. The tensile stresses are thereby relieved, preventing a similar failure at the opposite pole. The steadily erupting water ice loads the flanks of the open fissure, causing bending in the surrounding elastic plate and further tensile failure in bands parallel to the first fracture—a process that may be unique to Enceladus, where the gravity is too weak for compressive stresses to prevent fracture propagation through the thin ice shell. The sequence of fissures then cascades outwards until the loading becomes too weak or the background shell thickness becomes too great to permit through-going fractures.

Tidal interactions between planets or stars and the bodies that orbit them dissipate energy in their interiors. The dissipated energy heats the interior and a fraction of that energy will be released as seismic energy. Here we formalize a model to describe the tidally-driven seismic activity on planetary bodies based on tidal dissipation. To constrain the parameters of our model we use the seismic activity of the Moon, driven by tidal dissipation from the Earth-Moon interactions. We then apply this model to predict the amount of seismic energy release and largest seismic events on other moons in our Solar System and exoplanetary bodies. We find that many moons in the Solar System should be more seismically active than the Earth's Moon and many exoplanets should exhibit more seismic activity than the Earth. Finally, we examine how temporal-spatial variations in tidal dissipation manifest as variations in the locations and timing of seismic events on these bodies.

Crystals and bubbles nucleate and grow in a magma that experiences a range of temperatures, pressures and strain-rates. We have a good conceptual and sometimes quantitative understanding of how crystallization and bubble nucleation are controlled by decompression and cooling. Here we explore the effect of strain-rate on the crystallization kinetics of magmas. In order to understand the interaction between deformation and crystallization, samples of basalt were deformed during their crystallization. We made measurements at subliquidus conditions (1160°C) and deformed samples in compression at strain-rates varying from 0 to 2 × 10–4 s–1 for a total strain of 0.31. We simultaneously imaged the samples using X-ray micro-tomography. Without deformation, no crystallization was observed over the course of a 260 min experiment. Once deformation was applied, crystallization initiated. Deformation increased the nucleation rate, increased crystal growth rates, and decreased the incubation time. Increasing the strain-rate, however, does not show a discernable effect of crystallization kinetics. We hypothesize that deformation may have an effect on the parameters that govern the crystallization kinetics of magmas, such as activation energy and diffusion by changing chemical potentials.

Faecal Sludge (FS) contains high concentrations of pathogenic microorganisms that are 10-100 times higher than those in domestic wastewater. Proper and sustainable treatment is required to inactivate these pathogens if FS is to be recycled in agriculture, so as to minimise public health and environmental risks. Composting is one of the common low-cost technologies of sanitising FS in Urban Africa; however, it is associated with longer pathogen inactivation periods that make it commercially uneconomical. This study investigated the effect of different organic wastes types and their mixing ratios with FS on the inactivation efficiency of viable Ascaris eggs (suum and lumbricoides) during composting. Dewatered FS was mixed with market waste (MW), chicken feathers (CF) and sawdust (SD) in different ratios. Compost piles of FS:MW:SD and FS:CF:SD both in volumetric ratios of 1:2:1 and 1:3:1 were set-up in duplicate (3 m3 each), composted and monitored weekly for viable Ascaris eggs presence for a period of 15 weeks. The results suggest that the organic waste types have a significant effect on the temperature evolution and pathogen inactivation efficiency while their mixing ratios do not. Piles containing CF achieved the shortest pathogen survival period of 4 weeks compared with 6-8 weeks for those with MW. The temperature-time factor was found to be the most important variable responsible for viable Ascaris eggs inactivation. However, other mechanisms such as microbial antagonism or antibiotic action induced by indigenous microorganisms and toxic by-products such as free ammonia were found to have also played an important role in Ascaris eggs inactivation. All piles attained 100% Ascaris eggs inactivation from FS, and therefore, the compost was safe for use in agriculture. The study findings suggest that composting of FS with CF can reduce Ascaris eggs inactivation periods by 42%, which may thus reduce the operational costs of FS treatment facilities.

Volcanic eruptions are energetic events driven by the imbalance of magmatic forces. The magnitudes of these forces remain poorly constrained because they operate in regions that are inaccessible, either underground or dangerous to approach. New techniques are needed to quantify the processes that drive eruptions and to probe magma storage conditions. Here we present X-ray microdiffraction measurements of volcanic stress imparted as lattice distortions to the crystal cargo of magma from Yellowstone and Long Valley eruptions. Elevated residual stresses between 100 and 300 MPa are preserved in erupted quartz. Multiple volcanic forces could be culpable for the deformation so we analyzed crystals from pyroclastic falls, pyroclastic density currents, and effusive lavas. Stresses are preserved in all quartz but cannot be attributed to differences in eruption style. Instead, lattice deformation likely preserves an in situ measurement of the deviatoric stresses required for the brittle failure of viscous, crystal-bearing glass during ascent.

Rocks are heterogeneous multiscale porous media: two rock samples with identical bulk properties can vary widely in microstructure. The advent of digital rock technology and modern 3-D printing provides new opportunities to replicate rocks. However, the inherent trade-off between imaging resolution and sample size limits the scales over which microstructure and macrostructure can be identified and related to each other. Here, we develop a multiscale digital rock construction strategy by combining X-ray computed microtomography and focused-ion beam (FIB)-scanning electron microscope (SEM) images, and we apply the technique to a tight sandstone. The computed tomography (CT) scanning images characterize macroscale pore structures, while the FIB-SEM images capture microscale pore textures. The FIB-SEM images are then coupled to CT images via a template-matching algorithm and superposition. Bulk properties, including porosity and pore and throat size distribution, can be recovered with this approach. Permeability prediction with a pore network model for the largest connected pore network are 3 orders and 1 order of magnitude greater than the bulk rock measured value using the CT-only and the SEM-CT coupled images, respectively.

Hollow reentrants in quartz phenocrysts from Yellowstone (western United States) caldera's Lava Creek Tuff are preserved vestiges of bubbles in the supereruption's pre-eruptive magma reservoir. We characterized the reentrants using a combination of petrographic techniques, synchrotron X-ray microtomography, and cathodoluminescence imagery. One or more reentrants occur in ~20% of quartz, and up to ~90% of those reentrants are hollow. The earliest-erupted parts of the Lava Creek Tuff have the most empty reentrants. The hollow reentrants provide direct, physical evidence for volatile saturation, exsolution, and retention in a magma reservoir. Quartz-melt surface tension permits bubbles to attach to quartz only when bubbles have been able to nucleate and grow in the melt. Prior to eruption, the Lava Creek Tuff existed as a bubbly, volatile-saturated magma reservoir. The exsolved volatiles increased magma compressibility, helping to prevent the ever-accumulating magma from reaching a critical, eruptive overpressure until it reached a tremendous volume.

The response of unconfined aquifers to Earth's solid tides (Earth tides) has been used to study a number of hydrogeological, environmental, and ecological processes. The traditional interpretation of such responses has been based on a model that assumes free flow of groundwater to the water table, even though capillary effects have long been recognized to influence the response of coastal aquifers to ocean tides. Here we study the capillary effect on the response of unconfined aquifers to Earth tides, using a numerical procedure that couples groundwater flow and poroelastic strain in both the saturated and the unsaturated zones. We show that when the water table rises to the surface and the capillary zone disappears, the predicted tidal response is identical to that predicted by the traditional model. But when the capillary zone is present when the water table falls below the surface, the predicted amplitude ratio and phase shift depart notably from those predicted by the traditional model. We apply the model to a field case where large seasonal changes of tidal response have been documented. We show that the large seasonal changes in the tidal response may be due to the change of the water table from a capillary-affected state to a capillary-free state when the water table rises annually to the ground surface. Thus, the omission of the capillary effects in the interpretation of the response of unconfined aquifers to Earth tides may lead to incorrect inference of subsurface processes and properties.

Groundwater-level changes after earthquakes provide insight into changes in hydrogeological properties such as permeability and pore pressure. Quantifying such changes, both their location and magnitude, is usually hindered by limited data. Using extensive high-resolution water-level monitoring records, we provide direct evidence of significant groundwater drawdown (4.74-m maximum) over a 160-km2 area along crustal ruptures after the Mw 7.0, 2016, Kumamoto earthquake. Approximately 106 m3 of water disappeared within 35 min after the main shock. The loss of water was not caused by static-strain driven pore-pressure decrease nor by releasing of water through structural pathways, but most likely by water transfer downwards through open cracks. Such changes may impact the security of water resources, the safety of underground waste repositories, and contaminant transport in seismically active areas.

Meter-scale vesicular blocks, termed “giant pumice,” are characteristic primary products of many subaqueous silicic eruptions. The size of giant pumices allows us to describe meter-scale variations in textures and geochemistry with implications for shearing processes, ascent dynamics, and thermal histories within submarine conduits prior to eruption. The submarine eruption of Havre volcano, Kermadec Arc, in 2012, produced at least 0.1 km3 of rhyolitic giant pumice from a single 900-m-deep vent, with blocks up to 10 m in size transported to at least 6 km from source. We sampled and analyzed 29 giant pumices from the 2012 Havre eruption. Geochemical analyses of whole rock and matrix glass show no evidence for geochemical heterogeneities in parental magma; any textural variations can be attributed to crystallization of phenocrysts and microlites, and degassing. Extensive growth of microlites occurred near conduit walls where magma was then mingled with ascending microlite-poor, low viscosity rhyolite. Meter- to micron-scale textural analyses of giant pumices identify diversity throughout an individual block and between the exteriors of individual blocks. We identify evidence for post-disruption vesicle growth during pumice ascent in the water column above the submarine vent. A 2D cumulative strain model with a flared, shallow conduit may explain observed vesicularity contrasts (elongate tube vesicles vs spherical vesicles). Low vesicle number densities in these pumices from this high-intensity silicic eruption demonstrate the effect of hydrostatic pressure above a deep submarine vent in suppressing rapid late-stage bubble nucleation and inhibiting explosive fragmentation in the shallow conduit.

Abstract: The maximum solid phase concentration in geophysical turbulent gas‐particle mixtures is essential for understanding the flow dynamics but is poorly known. We present laboratory experiments on turbulent mixtures of air and ceramic particles generated in a vertical pipe. The mixtures had maximum bulk concentrations Cmax = 0.7–2.4 vol. % set by the onset of clustering and that increased with the degree of turbulence. Comparison with results of similar experiments with less dense (glass) particles reveals that Cmax increases with the particle Reynolds number according to Cmax = 0.78 × Rep0.17. Published results of experiments at specific Rep in different configurations are consistent with our data, suggesting that the model for Cmax may be generally applicable to turbulent mixtures. From our empirical laws we infer that natural turbulent gas‐particle flows with Rep ~ 102–105 have maximum solid concentrations of ~2–5 vol. %.

Pore pressure in aquifers confined below a cryosphere will increase as Mars cools and the cryosphere thickens. Increased pore pressure decreases the effective stress and hence promotes seismicity. We calculate the rate of pore pressure change from cooling of Mars's interior and the modulation of pore pressure from solar and Phobos tides and barometric loading. Using the time-varying pressure and tidal stresses, we compute Coulomb stress changes and the expected seismicity rate from a rate-and-state friction model. Seismicity rate will vary by several tens of percent to 2 orders of magnitude if the mean pore pressure is within 0.2 and 0.01 MPa of lithostatic, respectively. Seismic events promoted by high pore pressure may be tremor-like. Documenting (or not) tidally modulated shallow seismicity would provide evidence for (or against) water-filled confined aquifers, that pore pressure is high, and that the state of stress is close to failure—with implications for processes that can deliver water to the Martian surface.

Wastewater injection over the past decade has increased seismicity in the central USA, in some cases accompanied by detectable surface uplift. Here, we use this uplift to constrain subsurface properties and pore pressure evolution. We apply an advanced multitemporal interferometric algorithm to 35 synthetic aperture radar images acquired by ALOS satellite over four years before the 2012 earthquake sequence in east Texas, where large volumes of wastewater are disposed at depths of ∼800 m and ∼1800 m. To solve for the hydraulic diffusivity of the injection layers, we jointly inverted the injected volume and uplift data, considering a poroelastic layered half space. We find diffusivity values of 0.3±0.1 m 2 /s and 0.7±0.15 m 2 /s for shallow and deep injection layers, respectively, which combined with seismically-derived bulk moduli yields permeability values of 5.5±2.6×10 −14 m 2 and 1.9±0.25×10 −13 m 2 for these layers, consistent with permeability range reported for Rodessa formation and well test values. Hydraulic conductivity determines the evolution of pore pressure and thus the origin and location of induced seismicity. This study highlights the value of geodetic observations to constrain key hydrogeological properties of injection layers and to monitor the evolution of the subsurface pressure change.

Broadband seismic data were recorded on the ground surface around an exceptionally regular eruptive system, geyser El Jefe, in the El Tatio geyser field, Chile. We identify two stages in the eruption, recharge and discharge, characterized by a radial expansion and contraction, respectively, of the surface around the geyser. We model the deformation with spherical sources that vary in size, location, and pressure, constrained by pressure observations inside the conduit that are highly correlated with deformation signals. We find that in order to fit the data, the subsurface pressure sources must be laterally offset from the geyser vent during the recharge phase and that they must migrate upward toward the vent during the eruption phase. This pattern is consistent with models in which ascending fluids accumulate and then are released from a bubble trap that is horizontally offset from the shallow conduit of the geyser.

Hydrologic responses to earthquakes such as streamflow increases, water-level changes, and changes in geyser eruption frequency often reflect changes in permeability caused by seismic waves. The dynamic nature of permeability, as revealed by coseismic hydrologic phenomena, holds implications for groundwater systems, geothermal resources, mineral resources, and geologic hazards. Analysis of water-level responses to solid Earth tides and changes in atmospheric pressure provides a passive way to continuously monitor changes in permeability and storage properties in tectonically active regions.

Transitions in eruptive style are common at volcanoes. Understanding how and why these transitions occur remain open questions. The 2012 eruption of the submarine Havre volcano in the Kermadec arc (southwest Pacific Ocean) produced a raft of floating pumice followed by a pair of domes from the same vent. Here, we used measurements on erupted magmas and constraints on the eruption rate, combined with a model for magma ascent, to identify the dominant controls on the transition in eruption style. During the raft-forming stage, magma ascent was fast enough that little gas was lost. Magma reached the seafloor with great enough vesicularity to be buoyant and produce clasts that could float. As the eruption waned, the eruption rate decreased, and the conduit narrowed. Sufficient gas was then lost to the surrounding country rocks during ascent such that the erupted magma was no longer buoyant relative to seawater. Most of the original dissolved water in the magma was lost to the crust surrounding the conduit during the dome-forming stage.

In the submarine setting pyroclast cooling and saturation, by controlling the buoyancy of individual clasts and surrounding fluid, exert fundamental controls on the dispersal and fate of volcanic material. In this study we use laboratory experiments to understand and quantify how hot and air-filled pyroclasts cool and saturate in liquid water. By measuring internal temperature and clast submerged weight in over thirty experiments, we find that pumice cools and saturates in two stages. We observe that the first stage is marked by much higher rates of cooling and saturation than the second stage and that the transition between stages is sharp. We propose conceptual and quantitative models for each stage of pyroclast cooling and saturation. Specifically we argue that the first stage is defined by the presence of internal steam. That is, hot pyroclasts ingest liquid water and generate steam. Once a clast's internal pore space is filled with steam, heat loss to the surroundings causes condensation of internal vapor and clast saturation. Stage 1 ends when all internal steam has condensed. We test this idea with a Stefan Model with advection and conclude that heat loss, rather than permeability, controls the saturation of pyroclasts above the boiling temperature. We find that permeability matters for air-filled clasts with initial temperatures below the boiling point. Stage 2 begins after all internal vapor has condensed and we propose that stage 2 cooling and saturation are controlled by heat conduction and thermal contraction of trapped gas, respectively. From our understanding of the processes that govern the stage 1 to 2 transition, we derive an empirical average pyroclast cooling rate of q = 7.5 ± 0.5 W cm−2. Finally, we develop quantitative models for cooling-controlled porous clast saturation and for buoyant clast rise to the ocean surface. We find that meter size clasts stay buoyant for more than 10 min - long enough to reach the surface from a depth of 1000 m. The models developed here test our understanding of clast-scale processes and demonstrate that - in hot and porous clasts - the details of internal texture do not matter because clast saturation, and thus buoyancy and fate in the water column, are governed by heat transfer.

The cause of the Lusi mud eruption remains controversial. The review by Miller and Mazzini (2017) firmly dismisses a role of drilling operations at the adjacent Banjarpanji-1 well and argues that the eruption was triggered by the M6.3 Yogyakarta earthquake some 254 km away. We disagree with these conclusions. We review drilling data and the daily drilling reports, which clearly confirm that the wellbore was not intact and that there was a subsurface blowout. Downhole pressure data from Lusi directly witness the birth of Lusi at the surface on the 29th of May 2006, indicating a direct connection between the well and the eruption. Furthermore, the daily drilling reports specifically state that Lusi activity was visibly altered on three separate occasions by attempts to kill the eruption by pumping dense fluid down the BJP-1 well, providing further evidence of a connection between the wellbore and Lusi. By comparison with other examples of newly initiated mud eruptions elsewhere by other earthquakes, the Yogyakarta earthquake was far away given its magnitude. The seismic energy density of the Yogyakarta earthquake was only 0.0043 J/m3, which is less than a quarter of the minimum 0.019 J/m3 seismic energy density that has ever been inferred to trigger other mud eruptions. We show that the Lusi area had previously experienced other shallow earthquakes with similar frequencies and stronger ground shaking that did not trigger an eruption. Finally, the data from the BJP-1 well indicates that there was no prior hydrodynamic connection between deep overpressured hydrothermal fluids and the shallow Kalibeng clays, and that there was no evidence of any liquefaction or remobilization of the Kalibeng clays induced by the earthquake. We thus strongly favor initiation by drilling and not an earthquake.

The impressive number of stream gauges in Chile, combined with a suite of past and recent large earthquakes, makes Chile a unique natural laboratory to study several streams that recorded responses to multiple seismic events. We document changes in discharge in eight streams in Chile following two or more large earthquakes. In all cases, discharge increases. Changes in discharge occur for peak ground velocities greater than about 7–11 cm/s. Above that threshold, the magnitude of both the increase in discharge and the total excess water do not increase with increasing peak ground velocities. While these observations are consistent with previous work in California, they conflict with lab experiments that show that the magnitude of permeability changes increases with increasing amplitude of ground motion. Instead, our study suggests that streamflow responses are binary.

El Tatio in northern Chile is one of the best-studied geothermal fields in South America. However, there remain open questions about the mass and energy budgets, water recharge rates and residence time in the subsurface, origin of dissolved solutes, and processes affecting the phase and chemical composition of groundwater and surface water. We measured and sampled surface manifestations of the geothermal system (geysers, perpetual spouters, mud pools/volcanoes, and non-eruptive hot springs) and meteoric water. From the isotopic composition we infer that the thermal water has a meteoric origin that is different from the composition of local meteoric water. The absence of detectable tritium in thermal waters indicates that most of the recharge occurred pre-1950. Boiling and steam separation from the deep reservoir appear to be the main subsurface processes affecting the thermal fluids. A large amount of heat is lost to the atmosphere by evaporation from surface water and by steam emitted from erupting geysers. Using the chloride inventory method, we estimate thermal water discharge to be 218 to 234 L/s, and the advective heat flow to be 120 to 170 MW.

We image seismic intrinsic (Q−1i) and scattering (Q−1s) attenuations in Long Valley Caldera, California, by analyzing more than 1,700 vertical component waveforms from local earthquakes. Observed energy envelopes are fit to the diffusion model and seismic attenuation images are produced using two-dimensional space weighting functions. Low intrinsic and low scattering attenuation anomalies in the center south of the caldera correspond to the location of an earthquake swarm in 2014. We identify high intrinsic and high scattering attenuation anomalies in the fluid-rich western and eastern areas of the caldera. From a comparison with other geophysical images (magnetotellurics and seismic tomography) we attribute these anomalies to a hydrothermal system (high attenuation). Average to high attenuation values are also observed at the adjacent Mammoth Mountain (southwest of the caldera) and may also have a hydrothermal origin. High intrinsic attenuation at low frequencies to the west of the Hartley Springs Fault may be produced by the magmatic system that produced the Inyo Craters. Seismic attenuation imaging provides insights into subsurface structures that are complementary to velocity and conductivity images.

Clinopyroxene and plagioclase (andesine) microlites in an obsidian flow from Glass Mountain (NE California, USA) display strong alignment. Synchrotron X-ray diffraction, coupled with Rietveld analysis, was used to quantify crystallographic-preferred orientation (CPO). Clinopyroxene, with a rod-shaped morphology, shows a strong alignment of [001] in the flow direction and (010) aligned parallel to the inferred flow plane. Andesine, with a platy morphology, displays an alignment of (010) platelets in the flow plane. Some pole densities exceed 90 multiples of random distribution. Applying a model of rigid ellipsoidal inclusions in a viscous matrix, the local pure shear strains are between 2 and 3.

It has been noted for centuries that earthquakes appear to trigger the eruption of volcanoes. For example, analyses of global volcanic and seismic records since 1500 AD have shown that explosive eruptions with Volcanic Explosivity Index (VEI) values ≥2 are preceded within days by nearby major earthquakes (magnitude M8 or larger) about 4 times more often than expected due to coincidence, suggesting that large earthquakes can trigger eruptions. We expand the definition of a triggered eruption to include the possibility of M6 or greater earthquakes within 5 days and 800 km of a VEI 2 or greater eruption. Removing pre-1964 records, to ensure complete and accurate catalogs, we find 30 volcanoes that at some point experienced a potentially triggered eruption and define these volcanoes as “sensitive” volcanoes. Within this group of sensitive volcanoes, normalized distributions of volcano-centric factors such as tectonic setting, dominant rock type, and type of volcano are practically indistinguishable from those of sensitive volcanoes in which the time of eruption is randomized. Comparisons of sensitive volcanoes and insensitive volcanoes (i.e., volcanoes that have never experienced a triggered eruption) reveal that sensitive volcanoes are simply more active than insensitive volcanoes: They erupt more frequently, are located solely in subduction zones, and erupt primarily andesites and basaltic-andesites. The potentially triggered eruptions do not show the magnitude-distance relationship expected for seismically induced responses (e.g., hydrologic responses), and eruptions that do occur within days of nearby earthquakes do so within rates expected by random chance. There is, however, a 5–12% increase in the number of explosive eruptions in the 2 months to 2 years following major earthquakes. We conclude that short-term seismically triggered explosive eruptions occur less frequently than previously inferred, an important conclusion when considering volcanic hazards and for understanding the nature of earthquake-volcano interactions.

The Martian crustal dichotomy is the stark ∼5 km difference in surface elevation and ∼26 km difference in crustal thickness between the northern lowlands and southern highlands that originated within 100s of Myr of Mars’ formation. The origin of the dichotomy has broad implications for the geodynamic history of Mars, but purely exogenic or endogenic theories so far cannot explain all of the large scale geophysical observations associated with dichotomy formation. A giant impact can produce the shape and slope of the dichotomy boundary, but struggles to explain Mars’ remanent crustal magnetic signatures and the ultimate formation of Tharsis. Degree-1 mantle convection can relate the crustal dichotomy to the formation of Tharsis, but does not explain the elliptical dichotomy shape and must be initiated by a large pre-existing viscosity jump in the mantle. We propose a hybrid model of dichotomy formation in which a giant impact induces degree-1 convection with an upwelling antipodal to the impact site. In this scenario, a giant impact in the northern hemisphere excavates crust, creating an initial difference in crustal thickness and possibly composition between the two hemispheres. Over 10s to 100s of Myr, the dominant upwelling(s) would migrate to be under the thicker, insulating crust in the southern hemisphere, generating melt that further thickens the southern crust. We examine this process using 3-D mantle convection simulations, and find that a hemispherical difference in crustal thickness and composition caused by a giant impact can induce degree-1 convection with the upwelling(s) antipodal to the impact site in

Dilute gas-particle mixtures in which the particles are carried by the turbulent fluid are found in various geophysical contexts, from cold snow avalanches to hot pyroclastic density currents. Though previous studies suggest that such mixtures have maximum particle concentrations of a few volume percent, the dependence of this maximum concentration on the Reynolds number is unclear. We addressed this issue through laboratory experiments in a vertical pipe, where dilute gas-particle mixtures were created by injecting a turbulent air flow from below. Nearly monodisperse mixtures of glass beads of different grain sizes (77 to 1,550 μm) were used with varying bulk concentrations from 0.025 to 8 vol. %. To create quasi-static mixtures, the mean air velocity matched the terminal settling velocity for the grain sizes investigated. The maximum Reynolds numbers of the mixtures were ~104–106. The air pressure indicated full support of the particle weight at concentrations down to 0.025 vol. %. Above a critical particle concentration, at which all the particles were suspended, subsequent additional particles were not maintained in the mixture and led to the formation of clusters that settled downward in the pipe to form a dense fluidized bed. Maximum mean particle concentrations of the dilute mixtures increased from ~1 to ~2.8 vol. % and reached a plateau at increasing mixture Reynolds number. These results give insights into the maximum particle concentrations of geophysical turbulent gas-particle mixtures and may serve to constrain observations as well as the input and output data of models.

A long-standing conceptual model for deep submarine eruptions is that high hydrostatic pressure hinders degassing and acceleration, and suppresses magma fragmentation. The 2012 submarine rhyolite eruption of Havre volcano in the Kermadec arc provided constraints on critical parameters to quantitatively test these concepts. This eruption produced a >1 km3 raft of floating pumice and a 0.1 km3 field of giant (>1 m) pumice clasts distributed down-current from the vent. We address the mechanism of creating these clasts using a model for magma ascent in a conduit. We use water ingestion experiments to address why some clasts float and others sink. We show that at the eruption depth of 900 m, the melt retained enough dissolved water, and hence had a low enough viscosity, that strain-rates were too low to cause brittle fragmentation in the conduit, despite mass discharge rates similar to Plinian eruptions on land. There was still, however, enough exsolved vapor at the vent depth to make the magma buoyant relative to seawater. Buoyant magma was thus extruded into the ocean where it rose, quenched, and fragmented to produce clasts up to several meters in diameter. We show that these large clasts would have floated to the sea surface within minutes, where air could enter pore space, and the fate of clasts is then controlled by the ability to trap gas within their pore space. We show that clasts from the raft retain enough gas to remain afloat whereas fragments from giant pumice collected from the seafloor ingest more water and sink. The pumice raft and the giant pumice seafloor deposit were thus produced during a clast-generating effusive submarine eruption, where fragmentation occurred above the vent, and the subsequent fate of clasts was controlled by their ability to ingest water.

On 24 August 2016 a seismic event (Mw 6.0) was the first of the long Central Italy sequence (ongoing at the end of 2017) of medium-to-high magnitude earthquakes, with nine Mw ≥5 up to October 2017, and with about 74.000 seismic events registered after 1 year. The largest was the Mw 6.5 30 October 2016 event near Norcia. After the major seismic events, 17 mud volcanoes erupted around Monteleone di Fermo village (Marche region). Mud volcano eruptions generally occurred a few hours to a few days after the main earthquakes, suggesting a seismic triggering. We analyzed the peak ground velocities and dynamic stresses during the three largest earthquakes. We also evaluated the static stress changes in order to assess the potential influence of normal stress changes on the feeder system of the activated mud volcanoes. We find a correlation with dynamic stresses, whereas static stress changes are negligible or negative (with values reaching −0.44 bar, clamping feeder dykes). We conclude that seismic shaking (up to ~3.9 bar during Norcia earthquake) is the dominant driver for these eruptions. Finally, we evaluated the response ratio as a function of the dynamic stress. It increases exponentially with peak dynamic stress varying from 0.3 bar to >50% for peak dynamic stress >2 bar, indicating a link between earthquake shaking and mud volcano activity.

A few square kilometers of the Bishop Tuff in eastern California (USA) have evenly spaced columns that are more resistant to erosion than the surrounding tuff owing to the precipitation of mordenite, a low-temperature (100-130 °C) zeolite. We hypothesize that the columns are a result of instabilities at the liquid water and steam interface as cold water seeped into the still-cooling Bishop Tuff. We use two methods to quantitatively assess this hypothesis. First, scaling shows which hydrodynamic instabilities exist in the system. Second, to account for the effects of multiphase flow, latent heat, and the finite amplitude and temporal evolution of these instabilities we use two-dimensional numerical models of liquid water infiltrating hot tuff. These tests highlight several features of boiling hydrothermal systems. (1) The geometry of at least some convection appears to be broadly captured by linear stability theory that neglects reactive transport, heterogeneity of the host rock, and the finite amplitude of instabilities. (2) Slopes > 10% set the wavelength of convection, meaning that these columns formed somewhere with relatively gentle topography. (3) For permeabilities of > 10-13 m2, the wavelength of the instability changes through time, slowing infiltration, while for permeabilities < 10-15 m2, cooling is dominated by conduction. The spacing and stability of columns increase with higher vertical permeability and decrease with higher horizontal permeability. These columns are a rare window into hydrothermal processes that may be widespread.

Geysers episodically erupt liquid and vapor. Despite two centuries of scientific study, basic questions persist-why do geysers exist? What determines eruption intervals, durations, and heights? What initiates eruptions? Through monitoring eruption intervals, analyzing geophysical data, taking measurements within geyser conduits, performing numerical simulations, and constructing laboratory models, some of these questions have been addressed. Geysers are uncommon because they require a combination of abundant water recharge, magmatism, and rhyolite flows to supply heat and silica, and large fractures and cavities overlain by low-permeability materials to trap rising multiphase and multicomponent fluids. Eruptions are driven by the conversion of thermal to kinetic energy during decompression. Larger and deeper cavities permit larger eruptions and promote regularity by isolating water from weather variations. The ejection velocity may be limited by the speed of sound of the liquid + vapor mixture.

Documentation and analysis of groundwater responses to induced earthquakes are important to better understand their influence on shallow groundwater systems and hydrogeological properties and processes. Here we show that induced seismicity in Oklahoma can cause changes of groundwater level over distances >150 km from the epicenter. We test existing models for the cause of the observed responses and find that the model most consistent with observations is enhanced crustal permeability produced by seismic waves, changing aquifer recharge. Simulation suggests that the sources of this recharge are close to the responding wells and have lateral dimensions of ∼100 m. Continuous monitoring of pressure and temperature in wells, installing clustered wells to monitor multiple water levels near injection sites, and isotopic and chemical analysis of groundwater near injection sites are required to better understand and quantify the recharging sources.

Groundwater temperature changes will lag surface temperature changes from a changing climate. Steady state solutions of the heat-transport equations are used to identify key processes that control the long-term thermal response of springs and other groundwater discharge to climate change, in particular changes in (1) groundwater recharge rate and temperature and (2) land-surface temperature transmitted through the vadose zone. Transient solutions are developed to estimate the time required for new thermal signals to arrive at ecosystems. The solution is applied to the volcanic Medicine Lake highlands, California, USA, and associated springs complexes that host groundwater-dependent ecosystems. In this system, upper basin groundwater temperatures are strongly affected only by recharge conditions. However, as the vadose zone thins away from the highlands, changes in the average annual land-surface temperature also influence groundwater temperatures. Transient response to temperature change depends on both the conductive time scale and the rate at which recharge delivers heat. Most of the thermal response of groundwater at high elevations will occur within 20 years of a shift in recharge temperatures, but the large lower elevation springs will respond more slowly, with about half of the conductive response occurring within the first 20 years and about half of the advective response to higher recharge temperatures occurring in approximately 60 years.

Europa's surface contains numerous quasi-elliptical features called pits, domes, spots and small chaos. We propose that these features, collectively referred to as lenticulae, are the surface expression of saucer-shaped sills of liquid water in Europa's ice shell. In particular, the inclined sheets of water that surround a horizontal inner sill limit the lateral extent of intrusion, setting the lateral dimension of lenticulae. Furthermore, the inclined sheets disrupt the ice above the intrusion allowing the inner sill to thicken to produce the observed relief of lenticulae and to fracture the crust to form small chaos. Scaling relationships between sill depth and lateral extent imply that the hypothesized intrusions are, or were, 1–5 km below the surface. Liquid water is predicted to exist presently under pits and for a finite time under chaos and domes.

Pyroclastic density currents have been observed to both enter the sea, and to travel over water for tens of kilometers. Here, we identified a 1.2-m-thick, stratified pumice lapilli-ash cored at Site U1396 offshore Montserrat (Integrated Ocean Drilling Program [IODP] Expedition 340) as being the first deposit to provide evidence that it was formed by submarine deposition from pumice-rich pyroclastic density currents that traveled above the water surface. The age of the submarine deposit is ca. 4 Ma, and its magma source is similar to those for much younger Soufrière Hills deposits, indicating that the island experienced large-magnitude, subaer-ial caldera-forming explosive eruptions much earlier than recorded in land deposits. The deposit's combined sedimentological characteristics are incompatible with deposition from a submarine eruption, pyroclastic fall over water, or a submarine seafloor-hugging turbidity current derived from a subaerial pyroclastic density current that entered water at the shoreline. The stratified pumice lapilli-ash unit can be subdivided into at least three depositional units, with the lowermost one being clast supported. The unit contains grains in five separate size modes and has a >12 phi range. Particles are chiefly sub-rounded pumice clasts, lithic clasts, crystal fragments, and glass shards. Pumice clasts are very poorly segregated from other particle types, and lithic clasts occur throughout the deposit; fine particles are weakly density graded. We interpret the unit to record multiple closely spaced (

Pumice can float on water for months to years – long enough for pumice to travel across oceans and facilitate the spread of species. Long-lived pumice floatation is unexpected, however, because pumice pores are highly connected and water wets volcanic glass. As a result, observations of long floating times have not been reconciled with predictions of rapid sinking. We propose a mechanism to resolve this paradox – the trapping of gas bubbles by water within the pumice. Gas trapping refers to the isolation of gas by water within pore throats such that the gas becomes disconnected from the atmosphere and unable to escape. We use X-ray microtomography to image partially saturated pumice and demonstrate that non-condensable gas trapping occurs in both ambient temperature and hot (500 °C) pumice. Furthermore, we show that the size distribution of trapped gas clusters matches predictions of percolation theory. Finally, we propose that diffusion of trapped gas determines pumice floatation time. Experimental measurements of pumice floatation support a diffusion control on pumice buoyancy and we find that floatation time τ scales as τ∝L2Dθ2 where L is the characteristic length of pumice, D is the gas–water diffusion coefficient, and θ is pumice water saturation. A mechanistic understanding of pumice floatation is a step towards understanding how pumice is partitioned into floating and sinking components and provides an estimate for the lifetime of pumice rafts in the ocean.

Moderate to large earthquakes can increase the amount of water feeding stream flows, mobilizing excess water from deep groundwater, shallow groundwater, or the vadose zone. Here we examine the regional pattern of streamflow response to the Maule M8.8 earthquake across Chile's diverse topographic and hydro-climatic gradients. We combine streamflow analyses with groundwater flow modeling and a random forest classifier, and find that, after the earthquake, at least 85 streams had a change in flow. Discharge mostly increased (n=78) shortly after the earthquake, liberating an excess water volume of >1.1 km3, which is the largest ever reported following an earthquake. Several catchments had increased discharge of >50 mm, locally exceeding seasonal streamflow discharge under undisturbed conditions. Our modeling results favor enhanced vertical permeability induced by dynamic strain as the most probable process explaining the observed changes at the regional scale. Supporting this interpretation, our random forest classification identifies peak ground velocity and elevation extremes as most important for predicting streamflow response. Given the mean recurrence interval of ∼25 yr for >M8.0 earthquakes along the Peru–Chile Trench, our observations highlight the role of earthquakes in the regional water cycle, especially in arid environments.

Individual flood basalt lavas often exceed 103 km3 in volume, and many such lavas erupt during emplacement of flood basalt provinces. The large volume of individual flood basalt lavas implies correspondingly large magma reservoirs within or at the base of the crust. To erupt, some fraction of this magma must become buoyant and overpressure must be sufficient to encourage failure and dike propagation. The overpressure associated with a new injection of magma is inversely proportional to the total reservoir volume, and as a large magma body heats the surrounding rocks thermally activated creep will relax isotropic overpressure more rapidly. Here, we examine the viability of buoyancy overpressure as a trigger for continental flood basalt eruptions. We employ a new one-dimensional model that combines volatile exsolution, bubble growth and rise, assimilation, and permeable fluid escape from Moho-depth and crustal chambers. We investigate the temporal evolution of degassing and the eruptibility of magmas using the Siberian Traps flood basalts as a test case. We suggest that the volatile inventory set during mantle melting and redistributed via bubble motion controls ascent of magma into and through the crust, thereby regulating the tempo of flood basalt magmatism. Volatile-rich melts from low degrees of partial melting of the mantle are buoyant and erupt to the surface with little staging or crustal interaction. Melts with moderate volatile budgets accumulate in large, mostly molten magma chambers at the Moho or in the lower crust. These large magma bodies may remain buoyant and poised to erupt—triggered by volatile-rich recharge or external stresses—for ∼106 yr. If and when such chambers fail, enormous volumes of magma can ascend into the upper crust, staging at shallow levels and initiating substantial assimilation that contributes to pulses of large-volume flood basalt eruption. Our model further predicts that the Siberian Traps may have released 1019–1020 g of CO2 during a number of brief (∼104 yr) pulses, providing a plausible trigger for warming and ocean acidification during the end-Permian mass extinction. The assimilation of carbon-rich crustal rocks strongly enhances both flood basalt eruptibility and CO2 release, and the tempo of eruptions influences the environmental effects of CO2, SO2, and halogen degassing. The eruptive dynamics of flood basalts are thus inextricably linked with their environmental consequences.

The effect of sand filter media thickness on the performance of faecal sludge (FS) drying beds was determined in terms of: dewatering time, contaminant load removal efficiency, solids generation rate, nutrient content and helminth eggs viability in the dried sludge. A mixture of ventilated improved pit latrine sludge and septage in the ratio 1:2 was dewatered using three pilot-scale sludge drying beds with sand media thicknesses of 150, 250 and 350 mm. Five dewatering cycles were conducted and monitored for each drying bed. Although the 150 mm filter had the shortest average dewatering time of 3.65 days followed by 250 mm and 350 mm filters with 3.83 and 4.02 days, respectively, there was no significant difference (p > 0.05) attributable to filter media thickness configurations. However, there was a significant difference for the percolate contaminant loads in the removal and recovery efficiency of suspended solids, total solids, total volatile solids, nitrogen species, total phosphorus, chemical oxygen demand, dissolved chemical oxygen demand and biochemical oxygen demand, with the highest removal efficiency for each parameter achieved by the 350 mm filter. There were also significant differences in the nutrient content (NPK) and helminth eggs viability of the solids generated by the tested filters. Filtering media configurations similar to 350 mm have the greatest potential for optimising nutrient recovery from FS.