Geology topics

Research about California

Source-linked reports with geographic coverage including California.

At least 613 records · Page 34Linked to original sources

USGS2023

Changes in habitat suitability for wintering dabbling ducks during dry conditions in the Central Valley of California

In arid and Mediterranean regions, landscape-scale wetland conservation requires understanding how wildlife responds to dynamic freshwater availability and conservation actions to enhance wetland habitat. Taking advantage of Landsat satellite data and structured and community science bird survey data, we built species distribution models to describe how three duck species, the Northern Pintail ( Anas acuta ), Green-winged Teal ( Anas crecca ), and Northern Shoveler ( Anas clypeata ), respond to freshwater supply and food resources on different flooded land cover types in the Central Valley of California. Specifically, our models compared duck habitat suitability between the wettest and driest conditions in each month from September through April. Using abundance-weighted boosted regression trees, we created three sets of species occurrence models based on different covariates: (1) near real-time (hereafter “real-time”) covariates in which duck observations were matched to the water availability within the 16-day window of a Landsat observation, (2) a combination of real-time covariates and waterfowl food resource covariates describing annual corn and rice biomass and managed wetland moist soil seed yield estimates derived from Landsat data, and (3) long-term average covariates—the most common approach to species distribution modeling—in which long-term average surface water availability was used. We modeled the monthly occurrence of three duck species as a function of surface water availability, land cover type, road density, temperature, and bird data source. We found that dry conditions result in reduced habitat suitability, with the biggest reductions in November through January and in agricultural fields; in contrast, suitability of flooded wetland habitat was relatively robust to surface water availability. When models of habitat suitability based on long-term average climate conditions were compared to models based on real-time conditions, the highest long-term suitability values occurred in areas where suitability was high regardless of whether it was a wet or a dry year. While all models performed well, the inclusion of crop and wetland plant yield covariates resulted in slightly higher model performance. Overall, species distribution models created using data on the environmental conditions present at the time of bird observations can aid conservation efforts under extreme conditions over large spatial scales.

California↗

USGS2023

Hydrodynamics and habitat interact to structure fish communities within terminal channels of a tidal freshwater delta

Terminal channels were historically a common feature of tidal delta ecosystems but have become increasingly rare as landscapes have been modified. Tidal hydrodynamics are a defining feature in tidal terminal channel ecosystems from which native aquatic communities have evolved. However, few studies have explored the relationship between fish community structure and hydrodynamics in these tidal terminal channel ecosystems. We sampled fish communities throughout a network of terminal channels within the northeasternmost region of the San Francisco Estuary to determine the relationship between fish community structure and hydrodynamics within these environments. We collected two years (2017 and 2018) of fish community samples using gill nets and analyzed data using multivariate community analyses and count models. We found metrics of fish diversity and counts of native fishes to be greatest upstream (farthest from tidal influence) of the tidal excursion within terminal channels. Counts of non-native fishes were less affected by this hydrodynamic feature of terminal channels and more tightly correlated to local habitat conditions (e.g., water temperature, depth). Our results suggest that channel hydrodynamics plays a role in structuring fish communities within terminal channels, particularly native fishes. These results indicate that hydrodynamics in tidal delta ecosystems may be able to be altered in ways that benefit native fishes without the cost of water pumping.

California↗

USGS2023

Quality of groundwater used for domestic drinking-water supply in the Coachella Valley, 2020

Groundwater is the primary source of drinking water in the Coachella Valley in the desert region of southern California. Although most people in Coachella Valley are served by public drinking-water systems, about 20,000 people rely on private domestic or small-system wells (referred to herein as domestic wells). Recently, the U.S. Geological Survey (USGS) found that 39 percent of the groundwater resources used by domestic wells in Coachella Valley contained arsenic, fluoride, or both constituents at concentrations greater than the maximum contaminant levels established for public drinking-water systems. Uranium, chromium, nitrate, and perchlorate were detected at moderate concentrations below maximum contaminant levels. Elevated (above background) perchlorate concentrations in some areas indicate that domestic wells may receive recharge from Colorado River water used for irrigation or aquifer replenishment. Moderate total dissolved solids (TDS) concentrations throughout the study area and the co-occurrence of high concentrations of TDS and perchlorate indicates that Colorado River water is a source of recharge in the southeastern Indio groundwater subbasin. Four volatile organic compounds were detected at low concentrations, and pesticides and per- and polyfluoroalkyl substances were not detected.

California↗

USGS2023

Fracture-mesh faulting in the swarm-like 2020 Maacama sequence revealed by high-precision earthquake detection, location, and focal mechanisms

In August of 2020, an earthquake sequence initiated within the Maacama fault zone in northern California, raising questions about its relationship with the larger-scale fault. To investigate the faulting geometry and its implications for physical processes driving seismicity, we applied an integrated, multi-faceted seismic analysis including waveform-correlation-based event detection, relative relocation, and high-precision focal mechanisms. To determine mechanisms for a large population of very small earthquakes (predominantly M < 1), we combined correlation-derived polarity analysis with a recently developed technique to determine S/P amplitude ratios from single seismic components. Finally, we applied an iterative stress inversion to distinguish between the likely fault and auxiliary planes. This analysis reveals that although the sequence initiated on a right-lateral fault, plausibly a strand of the Maacama Fault, it also activated numerous en echelon left-lateral conjugate faults. Together, these interlocking faults form a fracture mesh, consistent with structures associated elsewhere with fluid-induced faulting and earthquake swarms.

California↗

USGS2023

Ages of the granitic basement of Long Valley Caldera, California, USA, and siting of the Quaternary granite-rhyolite pluton

The leucogranitic crystal-mush pluton beneath the iconic Long Valley Caldera, California, USA, released >820 km 3 of crystal-poor Pleistocene rhyolite, which was hosted by numerous Mesozoic granitic plutons, only a few of which had been dated until now. Reported here are U-Pb zircon ages, determined by sensitive high-resolution ion microprobe−reverse geometry (SHRIMP-RG), for 11 circumcaldera granitoids, all of them either Triassic or Cretaceous. Growth of the 35-km-wide Quaternary rhyolite-leucogranite plutonic reservoir was fostered by collocation of (1) a dense swath of late Pliocene basaltic vents, (2) a left-stepping extensional reentrant in the rangefront fault zone of the Sierra Nevada batholith, and (3) a sharp offset of the Proterozoic continental margin as represented by the Sr-isotope 0.706 line. We further consider whether the basement architecture of as many as 26 separate Triassic and Cretaceous plutons and intervening septa and pendants of Paleozoic metasedimentary rocks influenced siting of the Quaternary pluton and whether the ragged margin of Proterozoic lithosphere helped to focus asthenospheric edge upwelling that intensified crustal melting and intrusion in both the Triassic and the Quaternary.

California↗

USGS2023

Simulating post-dam removal effects of hatchery operations and disease on juvenile Chinook salmon (Oncorhynchus tshawytscha) production in the Lower Klamath River, California

Executive Summary The Federal Energy Regulatory Commission has been considering the approval to breach four dams on lower Klamath River in southern Oregon and northern California. Approval of this application would allow for Strikeouts indicate text deletion hereafter. decommissioning and dam removal, beginning as early as 2023. This action would affect Klamath River salmon ( Oncorhynchus ssp.) populations, a critical food source for federally endangered Southern Resident Killer Whales ( Orcinus orca ). In the long run, reintroduction of salmon populations to the upper Klamath River Basin may increase salmon abundance available to Southern Resident Killer Whales, but in the near term, it is uncertain how changes in hatchery management and disease-caused mortality by the myxosporean parasite Ceratonova shasta will influence abundance of salmon populations entering the ocean. To assess this uncertainty, we used the Stream Salmonid Simulator (S3) to simulate population dynamics of juvenile Chinook salmon ( Oncorhynchus tshawytscha ) for nine different population sources that rear and migrate through the Klamath River. S3 is a spatially explicit population model that runs on a daily time-step and simulates daily growth, survival, and movement of juvenile Chinook salmon from the time of spawning through ocean entry. The key features of this model relevant to this report include (1) a C. shasta disease submodel; (2) a temperature-dependent bioenergetics model that calculates daily growth rates; (3) size-dependent movement; (4) density-dependent dynamics that are influenced by the effect of flow on suitable habitat area; and (5) habitat, river flow, and water temperature specific to each scenario. We constructed and ran four scenarios: two scenarios for dams in place (Dams In) and dams removed (Dams Out), and given these dam-removal conditions, a low- and high-spore scenario for C. shasta . Each scenario was run for nine water years representing a range of conditions from dry to wet. Previously published daily river flows and water temperatures for Dams In and Dams Out provided physical inputs for each scenario. Daily spore concentrations were simulated using a three-part mechanistic model that used river discharge, water temperature, and the prevalence of infection (POI) of hatchery-origin Chinook salmon juveniles with C. shasta in the previous year 3 . We constructed two spore scenarios for each Dams In and Dams Out scenario, a “Low Spore” scenario and a “High Spore” scenario resulting in four scenarios for comparison. Spore scenarios were established by setting the prior-year POI of hatchery fish to 0.15 and 0.75 in the estimation of spore concentrations. Hatchery releases under Dams Out differed from those under the current Dams In scenario. Hatchery releases under the Dams Out scenario were modified to emulate changes in hatchery production that would occur under Dams Out conditions. This included moving hatchery production and releases from Iron Gate Dam to a proposed hatchery at Fall Creek, which would be located about 11 kilometers (km) upstream of Iron Gate Dam. It is anticipated that the Fall Creek hatchery would produce fewer fish at smaller and larger sizes at different release timings. For salmon inputs, we used observed historical abundance of main-stem spawners from brood year 2009 and juvenile salmon entering from tributaries in water year 2010, which represented an average return year for the 2005–18 period. Main-stem spawning was allowed to shift upstream from Iron Gate Dam under the Dams Out scenario. We also included hatchery-origin fish as natural spawners that would have otherwise returned to Iron Gate Hatchery in the first 3 years following dam removal. The S3 model simulated considerably higher total abundance for Dams Out relative to the respective Dams In scenarios, and higher abundance for the Low Spore scenario relative to the High Spore scenario. The difference in abundance between the four combinations of the dam-removal and spore scenarios varied among population groups. For main-stem natural production, juvenile abundance at ocean entry was 2–3 times higher for Dams Out scenarios than for Dams In scenarios, and juvenile abundance for High Spore scenarios was lower than that for the Dams Out Low Spores scenario. For hatchery releases, abundance at ocean entry was similar between Dams In and Dams Out scenarios for most water years, despite lower release sizes from Fall Creek Hatchery under Dams Out. For tributary populations, abundance for the High Spore scenarios was consistently lower than for the Low Spore scenarios, but differences between dam-removal scenarios varied among water years, with Dams Out scenarios having similar or higher abundance than Dams In scenarios, and dry water years having the largest difference between Dams In and Dams Out scenarios. We determined that different factors affected the response of each population group. For main-stem natural production, survival from fry emergence to ocean entry was higher under Dams Out scenarios compared to Dams In scenarios because juveniles emerged later and tended to arrive at the ocean sooner and at larger sizes, causing the population to have less time-dependent in-river mortality. Owing to their late release timing, hatchery populations had high disease-caused mortality in Dams In and Dams Out High Spore scenarios. Furthermore, a high proportion of infected fish (those that would be expected to die at some future point) survived to the ocean. Iron Gate Hatchery fish had lower survival rates than releases from Fall Creek Hatchery because the last mid-June release group from the 2010 Iron Gate Hatchery release incurred nearly total mortality in most water years owing to water temperatures exceeding 24 degrees Celsius. Our analysis shows how the S3 model was able to track different populations and provide insights on how the differential response of each population combined to influence the simulated number of juvenile Chinook salmon arriving at the Pacific Ocean where they become available as a food source for Southern Resident Killer Whales.

California↗

USGS2023

Quality of groundwater used for domestic supply in the Modesto, Turlock, and Merced Subbasins of the San Joaquin Valley, California

Summary More than 2 million Californians rely on groundwater from privately owned domestic wells for drinking-water supply. This report summarizes a water-quality survey of domestic and small-system drinking-water supply wells in the Modesto, Turlock, and Merced subbasins of the San Joaquin Valley where more than 78,000 residents are estimated to use privately owned domestic wells. Results indicate that inorganic and organic constituents in groundwater were respectively present above regulatory (maximum contaminant level, MCL) benchmarks for public drinking-water quality in 37 percent and 9 percent of the aquifer area used for domestic drinking-water supplies (herein, “domestic groundwater resources”). The most prevalent inorganic constituents exceeding regulatory benchmarks were nitrate, uranium, and arsenic. The only organic constituents exceeding regulatory benchmarks were the fumigant constituents 1,2,3-trichloropropane (1,2,3-TCP) and 1,2-dibromo-3-chloropropane (DBCP), but the herbicides atrazine and simazine were detected at low concentrations below one-tenth of regulatory benchmarks in 30 percent of domestic groundwater resources. Total dissolved solids (TDS) and manganese exceeded aesthetic-based (secondary maximum contaminant level [SMCL]) benchmarks for drinking water in 3 percent and 13 percent of domestic groundwater resources, respectively. Per- and polyfluoroalkyl substances (PFAS) were detected in 23 percent of domestic groundwater resources, with 4 percent exceeding California state notification or response levels for specific compounds. Total coliform bacteria were detected in 20 percent of domestic groundwater resources. Elevated concentrations of nitrate, uranium, TDS, and pesticides (fumigant constituents and herbicides) are related to agricultural land use and were typically present at shallow depths up to 75 meters below land surface. Agriculturally derived constituents were detected in wells screened below the Corcoran Clay Member of the Tulare Formation (herein, “Corcoran Clay”) in the southeastern part of the study area, where the Corcoran Clay tends to be shallower and thinner than in areas to the northwest. Nitrate, uranium, and TDS were most prevalent in the northwest part of the study area proximal to the valley trough where soils are poorly drained and agricultural land uses are predominantly grain, alfalfa, and dairy farms. Pesticides tended to occur in groundwater below coarse-grained surficial deposits and within a northwest to southeast trending band along the eastern extent of the Corcoran Clay that typically demarcates the western extent of well-drained soils associated with perennial orchard crops. Elevated concentrations of arsenic tended to occur west of this band in reducing groundwater but also sometimes co-occurred with elevated nitrate in oxic groundwater, most likely because of geochemical conditions in agriculturally affected groundwater that can enhance the mobility of arsenic from aquifer sediments.

California↗

USGS2023

Magmatic record of changing Cordilleran plate-boundary conditions—Insights from Lu-Hf isotopes in the Mojave Desert

Belts of Cordilleran arc plutons in the eastern part of the Mojave crustal province, inboard from the southwestern North American plate boundary, record major magmatic pulses at ca. 180–160 and 75 Ma and smaller pulses at ca. 100 and 20 Ma. This cyclic magmatism likely reflects evolving plate-margin processes. Zircon Lu-Hf isotopic characteristics and inherited zircons for different-age plutons may relate magma sources to evolving tectonics. Sources similar in age to the bulk of the exposed Mojave crust (1.6–1.8 Ga) dominated the magmas. Rare zircons having ε Hf (t) values as low as −52 indicate that Cretaceous melt sources also included more ancient crustal components, such as Archean-derived detritus in supracrustal gneisses of the Vishnu basin. Some rocks signal contributions from mantle lithosphere (in the Miocene) or asthenosphere (middle Cretaceous). Temporal shifts in isotopic pattern in this sample of the Cordillera relate to cyclic pulses of magmatic flux. Hf-isotopic pull-downs suggestive of dominantly crustal sources characterize the Jurassic and Late Cretaceous flare-ups. The Late Cretaceous flare-up, occurring near the onset of flat-slab subduction, produced abundant Proterozoic xenocrystic zircon and Hf isotopes implicating derivation largely from heterogeneous deep Mojave crust. Isotopic pull-ups characterize the lower-flux middle Cretaceous and Miocene magmatic episodes. The middle Cretaceous pulse ca. 105–95 Ma produced Mojave crust signals but also the isotopically most juvenile magmatic zircons, ranging upward to barely positive ε Hf values and suspected to signal an asthenosphere contribution. This may point toward transtension or slab retreat causing 105–95 Ma backarc extension in the Mojave hinterland of the Cordillera. That possibility of backarc extension raises questions about the tectonic environment of the contemporaneous main Sierra Nevada high-flux arc closer to the continental margin.

California↗

USGS2023

“Aftershock Faults” and what they could mean for seismic hazard assessment

We study stress‐loading mechanisms for the California faults used in rupture forecasts. Stress accumulation drives earthquakes, and that accumulation mechanism governs recurrence. Most moment release in California occurs because of relative motion between the Pacific plate and the Sierra Nevada block; we calculate relative motion directions at fault centers and compare with fault displacement directions. Dot products between these vectors reveal that some displacement directions are poorly aligned with plate motions. We displace a 3D finite‐element model according to relative motions and resolve stress tensors onto defined fault surfaces, which reveal that poorly aligned faults receive no tectonic loading. Because these faults are known to be active, we search for other loading mechanisms. We find that nearly all faults with no tectonic loading show increase in stress caused by slip on the San Andreas fault, according to an elastic dislocation model. Globally, faults that receive a sudden stress change respond with triggered earthquakes that obey an Omori law rate decay with time. We therefore term this class of faults as “aftershock faults.” These faults release ∼4% of the moment release in California, have ∼0.1%–5% probability of M 6.7 earthquakes in 30 yr, and have a 0.001%–1% 30 yr M 7.7 probability range.

California↗

USGS2023

Simulating the migration dynamics of juvenile salmonids through rivers and estuaries using a hydrodynamically driven enhanced particle tracking model

Juvenile salmonids migrate hundreds of kilometers from their natal streams to mature in the ocean. Throughout this migration, they respond to environmental cues such as local water velocities and other stimuli to direct and modulate their movements, often through heavily modified riverine and estuarine habitats. Management strategies in an uncertain future of climate change and altered land use regimes depend heavily on being able to reliably predict their ocean entry timings, route use, and survival rates through rivers and estuaries. We developed a spatially-explicit agent-based model of fish movement in response to hydrodynamic flows that uses movement dynamics gleaned from multi-dimensional tracking datasets of acoustically tagged juveniles moving through an urbanized, branched tidal estuary. We demonstrate how such models can be calibrated, and we apply it to the Sacramento-San Joaquin Delta in Central California. The quality of the out-of-sample validation of the model to predict juvenile salmon survival and route selection indicates that the model is versatile and flexible enough to be used in novel hydroclimatological conditions.

California↗

USGS2023

Flushing time variability in a short, low-inflow estuary

Flushing time, the time scale for exchange and mixing between embayed and oceanic waters in an estuary, plays an integral role in determining water quality and aquatic ecosystem health. Here, we investigated the spatiotemporal variability of flushing times throughout Morro Bay, a short, low-inflow estuary (LIE) on the California coast, using a calibrated and validated hydrodynamic model (Delft3D). Morro Bay has historically supported an extensive eelgrass ( Zostera marina ) habitat, which declined substantially from 139 to 5.4 ha during 2007–2017. Eelgrass decline motivated the current research into the role of changing bed roughness and oceanic drivers (i.e., tide and sea-level rise) on estuarine hydrodynamics and flushing times. We found that tidal variability exerts the strongest control on flushing times compared to other effects, i.e., bed roughness or sea-level rise. Additionally, we found that increasing sea level and decreasing bed roughness (associated with declining seagrass coverage) yielded higher rates of mixing (lower flushing times). We detected a strong correspondence between areas having shorter flushing times (e.g., near the estuary mouth) and areas occupied by resilient eelgrass populations in Morro Bay. Our findings further indicated that flushing times in short LIEs are particularly sensitive to several factors (e.g., bed roughness, sea level) that are susceptible to anthropogenic disturbance and future climate change.

California↗

USGS2023

A decision-support framework for dam removal planning and its application in northern California

Dam removals are occurring more frequently with the rising cost of maintaining aging infrastructure, public safety concerns, and growing interest in river restoration. So far, most dam-removals have been unsystematic in their approach. Given the several thousand dam removals expected over the coming decades, a systematic approach to plan future dam removals holds potential for aligning and delivering multiple benefits. Despite multi-sector factors driving decision-making, most existing prioritization frameworks tend to operate within single or related disciplines. Here we present a hierarchical, multi-disciplinary decision-support framework to prioritize dam removals based on opportunistic factors (Tier 1), hydro-ecological variables (Tier 2), and socio-cultural considerations (Tier 3). This framework integrates multiple decision criteria under data availability constraints, incorporates value-driven weights, and can be applied to a portfolio of dams at various spatial scales. The final output facilitates the identification of dam removal projects that align opportunistic, environmental, and social benefits. We recommend the application of this framework as a critical first step to identifying high-priority candidates for removal, recognizing that removal decisions will ultimately require detailed feasibility studies and stakeholder engagement. To illustrate its utility, we apply this framework to California's North Coast region and identify a small number of “good” candidates to be considered for removal. We conclude with recommendations for filling critical knowledge gaps and advancing systematic dam removal planning in the United States and beyond.

California↗

USGS2023

Deep root activity overprints weathering of petrogenic organic carbon in shale

The oxidation of organic carbon in sedimentary bedrock (petrogenic OC, OC p e t r o "> petro ) is increasingly recognized as a potential source of CO 2 to the atmosphere. Recent studies provide evidence for the mobilization and oxidation of OC p e t r o "> petro in sedimentary bedrock during rock weathering. However, the mechanisms and rates remain uncertain, particularly where overlying soils and vegetation drive contemporaneous oxidation of recently fixed organic carbon. Here, we quantify OC p e t r o "> petro weathering across a 16 m shale depth profile in a steep, rapidly eroding forested hillslope in the Northern California Coast Ranges. We report solid and gas phase radiocarbon and stable isotope analyses of samples extracted from specialized in-situ samplers, and a supporting laboratory incubation experiment of the shale regolith. OC p e t r o "> petro is removed from the weathered bedrock at a rate of approximately 0.12 gC/m 3 yr, which is orders of magnitude lower than the rate of OC p e t r o "> petro oxidation we achieved in the laboratory with crushed samples (557.1 gC/m 3 /yr). This disparity occurs despite high O 2 ( g ) "> 2(g) content across the depth profile, indicating that physical accessibility of OC p e t r o "> petro can regulate oxidative weathering. There is no direct radiocarbon evidence of OC p e t r o "> petro oxidation in CO 2 ( g ) "> 2(g) across the upper 13 m of the weathering profile during both wet and dry seasons. Instead, vadose zone CO 2 ( g ) "> 2(g) production at the site is dominated by respiration of recently fixed carbon associated with deep rooting. OC p e t r o "> petro is clearly mobilized across the vadose zone during weathering in this rapidly eroding, oxygen-rich, biologically dynamic hillslope, but at rates far below what can be measured given the contribution of root-derived CO 2 ( g ) "> 2(g) .

California↗

USGS2023

A life cycle model for evaluating estuary residency and restoration potential in Chinook salmon

Understanding the spatial and temporal habitat use of a population is a necessary step for recovery planning. For Chinook salmon ( Oncorhynchus tshawytscha ), variation in their migration and habitat use complicate predicting how restoring habitats could impact total recruitment. To evaluate how juvenile life history variation affects a population’s response to potential restoration, we developed a stage-structured model for a Chinook salmon population in a northern California river with a seasonally closed estuary. We modeled the timing of juvenile migration and estuarine use as a function of freshwater conditions and fish abundance. We used the model to evaluate the sensitivity of the population to different estuary and freshwater restoration scenarios that could affect population parameters at different life stages. The population’s run size increased most in response to freshwater restoration that enhanced spawning productivity (egg and fry survival), followed by spawner capacity. In contrast, estuary restoration scenarios affected only a subset of Chinook salmon (average 15%), and as a result, did not have a large impact on the total recruitment of a cohort. Under current condition, estuary rearing fish were over six times less likely to survive than fish that migrate to the ocean in the spring or early summer before estuary closure. Because estuary residents experienced low survival in the estuary and in the ocean, improvements to both estuary survival and growth would be needed to increase their total survival. When life cycle monitoring data is available, life cycle models such as ours generate predictions at scales relevant to conservation and are an advantageous approach to managing and conserving anadromous salmon that use multiple habitats throughout their life cycle.

California↗

USGS2023

Interactions among rainfall, fire, forbs and non-native grasses predict occupancy dynamics for the endangered Pacific pocket mouse (Perognathus longimembris pacificus) in a Mediterranean-type ecosystem

It is important to understand species-habitat relationships to implement effective adaptive management for rare species. However, it can be challenging to assess habitat associations and their relationships to abiotic stressors in dynamic habitats without the insights that can be gained from long-term monitoring. We report results from the first six years of extensive track tube monitoring of the largest two of three remaining extant populations of federally endangered Pacific pocket mouse ( Perognathus longimembris pacificus ) in a coastal Mediterranean-type ecosystem on Marine Corps Base, Camp Pendleton in southern California, USA. We used dynamic occupancy and structural equation modeling to assess potential drivers of population trends that included habitat, fire history, rainfall, disturbance, and the presence of other small mammals. We found that the variables that best predicted mouse occupancy were moderate to high forb and perennial herb cover (40–80%), and moderate to high open ground (20–70%) and low non-native grass cover (<20%), Non-native grass cover (>20%) was also a strong predictor of lower PPM colonization and increased extinction probabilities, with the extent of non-native grass cover being strongly influenced by annual rainfall and recency of fire. Our study adds to the growing literature on effects of invasive annual grasses on native species in Mediterranean-type ecosystems. We suggest that habitat management could be based upon promotion of open forb and perennial herb dominated habitats with reduction of non-native grasses by prescribed fire and other methods. These types of spatial and temporal monitoring programs can support land managers by creating a monitoring and management feedback loop. They can reveal landscape and environmental variables associated with species persistence, inform habitat management goals, and help managers to assess the success of management actions on populations of conservation concern.

California↗

USGS2023

Water-soluble organic carbon release from mineral soils and sediments in an irrigated agricultural system

Water interactions with soil and vegetation are greatly altered in agricultural watersheds compared to natural landscapes, which impacts sources and fates of organic carbon (OC). While mineral soil horizons in natural ecosystems primarily act as filters for dissolved organic carbon (DOC) leached from organic surface horizons, tilled soils largely lack an organic horizon and their mineral horizons therefore act as a source for both DOC and sediment to surface waters. Irrigated watersheds highlight this difference, as DOC and total suspended sediment (TSS) concentrations simultaneously increase during the low-discharge irrigation season, suggesting that sediment-associated OC may constitute a significant source of DOC. While water-soluble OC (WSOC) from sediments and soils has been found to be compositionally similar to stream DOC, these contributions remain poorly quantified in agricultural streams. To address this, we conducted abiotic solubilization experiments using sediments (suspended and bed) and soils from an irrigated agricultural watershed in northern California, USA. Sediments (R 2 > 0.99) and soils (0.74 < R 2 < 0.89) displayed linear solubilization behaviors over the range of concentrations tested. Suspended sediment from the irrigation season exhibited the largest solubilization efficiency (10.9 ± 1.6% TOC sediment solubilized) and potential (1.79 ± 0.26 mg WSOC g −1 dry sediment), followed by suspended sediment from a winter storm, then bed sediment and soils. Successive solubilization experiments increased the total release of WSOC by ∼50%, but most (88–97%) of the solid-phase OC remained insoluble in water. Using these solubilization potential estimates and measured TSS concentrations, we estimated that WSOC from suspended sediment in streams represented 4–7% of the annual DOC export from the watershed. However, field sediment export is much higher than what is represented by suspended sediment in the water column, therefore field-scale contributions from sediments could be much higher than estimated.

California↗

USGS2023

Modeling the dynamic penetration depth of post-1950s water in unconfined aquifers using environmental tracers: Central Valley, California

The penetration depth of post-1950s recharge (D-1950) in aquifers is a marker that is frequently used to identify groundwater that is susceptible to anthropogenic contamination. Here, we compute D-1950 values at wells, interpolate them in space, and project them across time to map the moving front of modern recharge in four dimensions in the Central Valley aquifer system, California, USA. Tracers of groundwater age (tritium, carbon-14, noble gases, sulfur hexafluoride, and chlorofluorocarbons) were collected at 650 wells spatially distributed throughout the Central Valley and were fit to a lumped-parameter model that assumes a logarithmic age-depth profile in the aquifer. For samples where tritium was present (>0.3 tritium units), the model was used to predict D-1950 at wells screened above or across the modern-premodern interface (n = 484). Wells with samples where tritium was absent (≤0.3 tritium units) were used to define the depth beyond which groundwater is completely premodern (n = 166). Predicted D-1950 values were below the depth of screen bottoms for wells where groundwater is completely modern, and above the depth of screen tops for wells where groundwater is completely premodern. The interpolated surface of D-1950 is dynamic, less prone to extreme values, and produces maps with lower interpolation errors due to a higher spatial density of wells than maps based on the depth of premodern groundwater. Between 2005 and 2025, D-1950 is expected to deepen by 11 and 12 m in the northern and southern parts of the Central Valley, respectively. Areas where D-1950 increases rapidly are likely to see increases in nitrate and other anthropogenic contaminants associated with the downward moving front of modern water.

California↗

USGS2023

Precaldera mafic magmatism at Long Valley, California: Magma-tectonic siting and incubation of the Great Rhyolite System

The iconic volcanic center at Long Valley has released ∼820 km 3 of rhyolite in at least 110 eruptions. From 2.2 Ma until 0.23 Ma, products were exclusively rhyolitic, and ∼ 700 km 3 were high-silica rhyolite severely depleted in Sr, Ba, and Eu. The rhyolitic interval was preceded by an interval from 3.9 to 2.6 Ma with numerous basalt-andesite-dacite eruptions accompanied by no rhyolite at all. We have now mapped the circumcaldera products of this interval, defined 107 eruptive units, characterized them all chemically and petrographically, and dated many by 40 Ar/ 39 Ar. Here we display and describe them by sector around the caldera, interpret the nature of the transcrustal magma system that eventuated in the 35-km-wide Long Valley granite-rhyolite pluton, and analyze regional tectonic factors that did or did not contribute to siting the system. Nine Miocene (12–6 Ma) eruptive units close to Long Valley were followed by a Pliocene flare-up that released >300 mafic eruptions in a SW–NE swath 170 km long, centered across the later site of Long Valley. The basalts and their fractionates are intraplate alkalic products dominated by a continental lithosphere that had long been fluxed by Mesozoic subduction. Tertiary arc volcanism had not impinged on the area of Long Valley. Volumes estimated for the 107 Neogene precaldera eruptive units (only 40 of which exceeded 0.1 km 3 ) total ∼ 27 km 3 ± 50%—only ∼3% of the subsequent volume of rhyolite erupted. Such a volume of high-silica rhyolite with ultra-low Sr and Eu is not a product of partial melting but requires as proximate parent a leucogranitic crystal mush that is itself the upper level of a long-lived plutonic reservoir that extends to the lower crust. The 27 km 3 of Neogene magma that erupted was a small contingent of the mantle-derived basaltic flux needed to energize (and contribute its fractionated melt to) a 30-km-deep compositionally graded crustal column, which culminated in ∼10,000 km 3 of granitoid mush from which 820 km 3 of Quaternary high-silica rhyolitic melt escaped and erupted. Pliocene basaltic eruptions ceased at ∼2.6 Ma, probably because the basaltic flux intensified sufficiently to render the mushy upper crust impenetrable. The 2.6–2.2 Ma quiescent interval represented culmination of thermal activation of the plutonic column and refinement of its leucogranitic mushy upper layer, from which extreme melts escaped for the next 2 Myr. The Pliocene mafic swath crosses the Sierran rangefront fault zone coincident with a left-stepping extensional reentrant that also began developing at ∼3 Ma. Moreover, Long Valley overlies a dextral offset in the initial Sr-isotope 0.706 line, which may reflect the rifted or attenuated edge of Proterozoic crust and mantle lithosphere. Concatenation of these three influences may account for siting of intensified edge-focused magmatism that produced the great Quaternary pluton.

California↗