SEARCH · USGS Science

Results for “Article”

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.

979 records · Page 29Linked to original sources

USGS2025

Generalized Bancroft algorithm for locating earthquakes with P- and S-wave arrival times

Because of similarities between locating an earthquake with seismic stations and locating a Global Positioning System (GPS) receiver from satellites, the Bancroft algorithm developed for GPS processing can be used to locate earthquakes. Such an approach to earthquake location differs from the conventional method of choosing an initial or trial solution and then iteratively improving the solution until convergence. The Bancroft algorithm has the advantage of being a direct, noniterative solution but with the disadvantage of only being able to accommodate a homogeneous velocity model. An additional limitation of the standard Bancroft algorithm is that it considers arrival times in a medium with a single propagation velocity. This poses no problem for GPS processing because electromagnetic waves travel at the speed of light; however, for seismic waves it means the algorithm can be applied to collections of either P ‐ or S ‐wave arrival times. Here, I show how the Bancroft algorithm can be generalized to handle both P ‐ and S ‐wave arrival‐time measurements simultaneously. I also show how to accommodate depth‐varying P ‐ and S ‐wave velocity models. I apply the generalized Bancroft algorithm to microearthquakes beneath Tanaga Volcano in Alaska and compare standard locations from the widely used HYPOINVERSE location code to Bancroft locations and to the output of HYPOINVERSE when setting the trial location to the Bancroft location. I find the Bancroft locations outperform the results from the other methods for shallow earthquakes near sea level, where a quantity known as the geometric dilution of precision is large and linearized approaches such as HYPOINVERSE are expected to struggle.

Alaska↗

USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗

USGS2026

The geologic history of the Chehalis Forearc Basin, Washington State, USA

The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.

Washington↗

USGS2026

Egg thiamine concentrations reflect dietary differences between four hatchery strains, including two morphotypes, of Lake Trout in Lake Ontario

Objective Test whether dietary differences among hatchery strains of Lake Trout Salvelinus namaycush were associated with variation in egg thiamine concentrations and assess the potential for using strain-specific traits to mitigate thiamine deficiency complex (TDC) in Lake Ontario restoration efforts. Methods Stomach content analysis and stable isotope mixing models were used to compare diets among four hatchery strains of Lake Trout, representing lean and humper morphotypes. Egg thiamine concentrations were measured to test whether dietary differences correlated with susceptibility to TDC. Results Significant differences in diet and egg thiamine concentrations were observed among Lake Trout strains. A humper-like morphotype strain exhibited elevated thiamine levels, which were associated with reduced consumption of prey known to degrade thiamine levels when consumed, Alewife Alosa pseudoharengus and Rainbow Smelt Osmerus mordax , and increased reliance on Round Goby Neogobius melanostomus or sculpin species (Slimy Sculpin Cottus cognatus and Deepwater Sculpin Myoxocephalus thompsonii ). These findings support the link between the consumption of Alewife and Rainbow Smelt and thiamine deficiency. Conclusions Our findings support a link between strain-specific dietary preferences and susceptibility to TDC in Lake Trout. Preferential stocking of strains less prone to TDC, such as the strain derived from a humper morph from Lake Superior, could enhance restoration success in Lake Ontario. This strategy may be especially effective given the continued dominance of Alewife in the lake's prey base and the current management priorities supporting fisheries for Pacific salmon Oncorhynchus spp.

North American Journal of Fisheries Management↗

USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗

USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗

USGS2024

Advancing sustainable groundwater management with a hydro-economic system model: Investigations in the Harney Basin, Oregon

Groundwater resources frequently trend toward unsustainable levels because, absent effective institutions, individual water users generally act independently without considering the impacts on other users. Hydro-economic models (HEMs) of human-natural systems can play a positive role toward successful groundwater management by yielding valuable knowledge and insight. The current study explores how an HEM that captures essential physical and economic characteristics of a system can shed light on the system's processes and dynamics to benefit stakeholders, managers, and also researchers. These propositions are illustrated using the Harney Basin, Oregon, which has seen large groundwater declines in the past 20 years. The HEM shows that: (a) although current groundwater pumping rates will gradually raise costs and reduce well yields, irrigators gain the highest aggregate economic return by continuing current pumping; (b) lowland areas of the basin are hydrologically connected, which limits the efficacy of remedies focused on regulations only in some portions of the basin; (c) community expectations regarding the efficacy of several proposed solutions are overly optimistic; and (d) the study's scenarios identify interventions that would stabilize the groundwater system and prevent additional adverse impacts on residential and livestock wells and groundwater-dependent ecosystems. These interventions would require limiting groundwater pumping by nearly half and reducing annual profits by $7.5–$9.0M. The HEM also demonstrated its value to researchers: its insights shifted attention toward questions about Oregon's existing groundwater institutions and their inability to adaptively manage the transition from abundant groundwater to scarce groundwater in a timely manner.

Oregon↗

USGS2025

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin↗

USGS2025

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

California↗

USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗

USGS2025

To heal or not to heal?: 1. The effect of pore fluid pressure on the frictional healing behavior of lithologies in Oklahoma

The competition between fault healing (i.e., re-strengthening) and fault loading determines the timing and magnitude of fault failure within the seismic cycle. Repeating earthquakes can give observational estimates of fault healing rates, however, it is difficult to link laboratory studies of frictional healing and observed healing rates from repeating earthquakes in part because of uncertainty in lithology at depth. Due to well-constrained and relatively simple geology, earthquakes in Oklahoma can be linked to the granitic basement rock and to the Arbuckle Group, which is primarily composed of dolomite at earthquake depths. Here, we conduct friction experiments to measure healing rates of the two earthquake-bearing lithologies at confining pressures representative of earthquake depths and pore pressures ranging from 0% to 80% of the confining pressure. We measure frictional healing by executing slide-hold-slide tests with hold times ranging from 3 s to 3000 s. The friction experiments on the Troy Granite indicate that pore fluid pressure does not greatly affect healing rate. On the other hand, the dolomite of the Arbuckle Group exhibits decreased healing with increased pore fluid pressure, with weakening at the highest pore pressure. We hypothesize that this is due to an increase in dissolution of dolomite at high pore pressures/low effective normal stress. These healing rates are used in the companion paper to understand the moment-recurrence time behavior of repeating earthquakes in Prague, Oklahoma. This work has implications for possible enhanced dissolution and weakening behavior of the Arbuckle Group during wastewater injection activities.

Oklahoma↗

USGS2026

Potential warming-induced changes in stream water sources between proglacial and non-glacial streams, south-central Alaska

Snow and ice-melt are essential for global water resources. Rising air temperatures are causing vegetation encroachment and shrubification of previously snow- and ice-covered landscapes, and precipitation regimes are changing. Collectively, these changes will likely affect the hydrology of mountain environments, although the exact ways in which the hydrology will change are poorly understood. Changing source waters (i.e., the proportion of ice-melt, snowmelt, rain and groundwater) affect the timing and magnitude of streamflow and affect stream temperature, sediment, solute and nutrient fluxes. We investigate spatial and seasonal source water contributions across sub-watersheds of the Nellie Juan River watershed on Alaska's Kenai Peninsula, which differ in glacier cover, elevation and land cover, using a space-for-time approach. To do so, we use a three-endmember mixing model to quantify the contribution of rain, snow/ice-melt and deep groundwater. We show that snow/ice-melt peak later in proglacial streams compared to non-glacial streams. Endmember contributions in non-glacial streams generally vary with elevation and season. We also found that groundwater is a major contributor to streamflow across the study area, contributing between 12% and 47% seasonally, and may become increasingly important as snowpack and ice-melt decrease. This research characterises the shift from a glacierised to a deglaciating landscape through the lens of source water contribution. Broadly, these findings could help improve understanding of how water resources in glacierised watersheds are affected by warming air temperatures.

Alaska↗

USGS2024

Prioritizing US Geological Survey science on salinization and salinity in candidate and selected priority river basins

The US Geological Survey (USGS) is selecting and prioritizing basins, known as Integrated Water Science basins, for monitoring and intensive study. Previous efforts to aid in this selection process include a scientifically defensible and quantitative assessment of basins facing human-caused water resource challenges (Van Metre et al. in Environmental Monitoring and Assessment , 192(7), 458 2020). In the present work, we explore this ranking process based on water quality considerations, specifically salinity and salinization. We selected top candidate basins to study salinity and salinization issues in 18 hydrologic regions that include 163 candidate basins. Our prioritization is based on quantitative assessment of sources of salinity, drivers of change, and receptors that must respond to those sources and drivers. Source terms represented in the prioritization include geology, depth to brackish groundwater, stream conductivity, chloride in precipitation, urban and agricultural land use, application of road salt as a deicer, and irrigation. Drivers represented in prioritization include changes in chemical weathering as a result of changes in rainwater chemistry. Receptors include measures of water stress, measurements of stream ecological health, and socioeconomic factors. In addition, we present research activities for the USGS on salinity and salinization that can be pursued in these basins including assessment of sources, pathways, and loadings; predicting and understanding changes in sources, peaks, and trends; understanding the components of salinity and mobilization of contaminants; understanding the relationship between salinization and changing ecosystems; and developing knowledge on the causes and distribution of groundwater salinity, brackish water resources, and challenges related to desalination.

Enviornmental Monitoring and Assessment↗

USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗

USGS2025

Characterizing sedimentary organic carbon in a hydrothermal spreading center, the Escanaba Trough

Sediments in critical marine mineral environments are of wide importance due to their preservation of both marine minerals and organic carbon (OC) stocks. However, OC storage and cycling is often overlooked in mineral system studies. This work characterizes sedimentary OC within the Escanaba Trough, a hydrothermal sulfide system off the coast of northern California. By utilizing ROV-based push coring, we collected sediments near and far from high temperature, low temperature, and inactive vents. We applied a multipronged organic geochemical approach, measuring bulk sediment, OC fractions of varying labilities, and biomarkers to tease apart the storage, source, and cycling of OC within this complex system. In contrast to past work indicating a primarily terrestrial source to deeper, Pleistocene sediments, our results suggest a primarily marine source in surface sediments at least 50 m away from active venting. Near active venting, we see evidence of locally produced chemosynthetic OC in addition to the background marine source. This chemosynthetic OC appears to be rapidly remineralized and supports more complex deep-sea food webs compared to hemipelagic sites. Still, the greatest contribution of labile OC was observed at inactive vent sites, which we suggest is due to hydrothermal fluid alteration of more labile OC at actively venting sites. Depending on fluid temperature, hydrothermally altered OC is either preserved in the sediments as petroleum or migrated with fluid flow. Considering inactive sites are those most likely to be targeted by potential mining, we suggest additional studies could verify if these results apply at other sedimented seafloor massive sulfide systems.

California↗

USGS2025

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗

USGS2026

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences↗