SEARCH · USGS Science
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The unstable nature of many chemical and physical constituents in ground water requires special collection procedures and field analysis immediately after collection. This report describes the techniques and equipment commonly used in the collection and field analysis of samples for pH, temperature, carbonate, bicarbonate, specific conductance, Eh, and dissolved oxygen.
Groundwater depletion for agricultural irrigation poses significant environmental and economic challenges. This study introduces a proof-of-concept that combines hydro-economic modeling, scenario-based modeling, and multi-objective optimization to manage pumping curtailment in an over-allocated basin in the western United States. Three optimization scenarios were evaluated, each offering different degrees of management flexibility. Results reveal that scenarios with finer spatial resolution achieved greater environmental benefits per unit profit loss. Additionally, strategies allowing fractional reductions in curtailed wells–rather than complete shutdowns based on water rights seniority–substantially improved efficiency, highlighting the value of increased decision-making flexibility. Although scenario testing can aid stakeholder engagement and strategy exploration, multi-objective optimization provides a systematic framework to quantify tradeoffs between competing objectives. This combined approach demonstrates promise for building consensus and supporting the design of sustainable water management strategies that balance agricultural livelihoods with ecosystem preservation.
Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.
Foreward This document contains design specifications of a navigational standard for the Moon, including a Lunar Transverse Mercator system, a Lunar Polar Stereographic system, a Lunar Grid Reference System, and a unique coordinate structure, Artemis Condensed Coordinates, for Artemis mission navigation and lunar surface science. The National Aeronautics and Space Administration (NASA) Artemis campaign seeks to place humans on the Moon for the first time since the Apollo missions. Early Artemis missions are heavily focused on the lunar south pole, which promises to return valuable data on the Moon’s geologic record, amongst other mission objectives. Coordinate systems in use today for the lunar south pole provides crew members on the surface neither an efficient nor intuitive means to communicate their position and orientation. A novel grid coordinate system, the Lunar Grid Reference System, is proposed to address these concerns for use in real-time extravehicular activity operations on the lunar surface. The many stakeholders involved in the Artemis missions will need a common system to communicate position and orientation while astronauts are operating on the lunar surface. To that end, Artemis crew members will need that system to be efficient and intuitive to promote efficient extravehicular activity timelines and reduce confusion. In the context of this document, these characteristics are addressed on the design of lunar coordinate systems: Efficient.—The number of characters required to communicate a location within a desired precision level in both local and global contexts, and how many steps are required for a recipient or sender to interpret a location. Intuitive.—How well the system aligns with human perceptual abilities, and whether the system yields distances that have the same relationship to actual lunar surface distance in all directions from the point where a person is located. Technological systems are currently being investigated to supplement the crew members’ ability to locate and orient themselves and other assets on the lunar surface; however, it is unlikely that those systems will be fully operational for the first few landed missions. Even with future positional aids, crew members will still need an efficient and intuitive means to communicate position and orientation. In addition, if technological systems fail, the crews will require land navigation skills and have maps available, thus providing further motivation for a crew-centric coordinate system. The contents of this U.S. Geological Survey (USGS) document detail a comprehensive framework for standardizing lunar crewed surface navigation within NASA and outlines the protocols, methods, and designs necessary for achieving consistency and interoperability across relevant space mission teams and lunar surface navigators. Key components of this document include designs of map projections, projected coordinate reference systems (Lunar Transverse Mercator and Lunar Polar Stereographic systems), and a grid system (Lunar Grid Reference System and Artemis Condensed Coordinates) for the Moon. The work proposed in this document seeks to accomplish something similar to the National Geospatial-Intelligence Agency (NGA) document SIG 0012 (NGA, 2014a), but for using grid systems for the Moon. This report incorporates initial feedback and input from NASA’s Artemis Geospatial Data Team, NASA’s Flight Operations Directorate, National Geodetic Survey, USGS Astrogeology Science Center, and NGA and is intended to serve as a resource for all involved with the Artemis missions, as well as for engineers designing and operating lunar infrastructure.
The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.
We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.
While deterministically predicting the time and location of earthquakes remains impossible, earthquake forecasting models can provide estimates of the probabilities of earthquakes occurring within some region over time. To enable informed decision-making of civil protection, governmental agencies, or the public, Operational Earthquake Forecasting (OEF) systems aim to provide authoritative earthquake forecasts based on current earthquake activity in near-real time. Establishing OEF systems involves several nontrivial choices. This review captures the current state of OEF worldwide and analyzes expert recommendations on the development, testing, and communication of earthquake forecasts. An introductory summary of OEF-related research is followed by a description of OEF systems in Italy, New Zealand, and the United States. Combined, these two parts provide an informative and transparent snapshot of today's OEF landscape. In Section 4, we analyze the results of an expert elicitation that was conducted to seek guidance for the establishment of OEF systems. The elicitation identifies consensus and dissent on OEF issues among a non-representative group of 20 international earthquake forecasting experts. While the experts agree that communication products should be developed in collaboration with the forecast user groups, they disagree on whether forecasting models and testing methods should be user-dependent. No recommendations of strict model requirements could be elicited, but benchmark comparisons, prospective testing, reproducibility, and transparency are encouraged. Section 5 gives an outlook on the future of OEF. Besides covering recent research on earthquake forecasting model development and testing, upcoming OEF initiatives are described in the context of the expert elicitation findings.
Continental remobilization is a crucial driver of metallogenesis and the formation of ore deposits. Some of the world’s largest mineral deposits of the economically valuable elements tin (Sn), tungsten (W), gold (Au), copper (Cu), lead (Pb), and zinc (Zn) formed during the Mesozoic Era. Additionally, the chemistry and distribution of the elements Sn and W have been investigated in previous studies to understand planetary formation and differentiation processes. These two elements are largely co-located during certain South China Mesozoic metallogenic events but are not co-located during other time periods in the same regions. Here, we investigated the mineral chemistry network similarities and dissimilarities of Sn and W to understand their mineral formation and distribution during the Mesozoic Era and throughout Earth history. Mineral chemistry network community detection analysis and electronegativity associations among mineral constituent elements of Sn minerals and W minerals indicate that the elements have similar chemistry among their oxide minerals. However, Sn forms a much wider range of minerals that also contain S compared to W, which occurs in a limited number of S-containing minerals. The divergent constituent element interactions among S-containing Sn minerals and W minerals reflect the redox sensitivity and importance of oxygen (O) fugacity in Sn mineral formation. Conversely, extensive W mineral deposits are known to form at both high and low O fugacities. The similarities and differences between the mineral chemistry networks of Sn and W reflect the mineral distribution of the two elements in the Sn-W mineralization event from 160 to 139 Ma vs. the Sn–uranium (U) mineralization event from 125 to 98 million years ago (Ma). The mineral chemistry and distribution of Mesozoic Sn and W deposits illustrate the contrasting importance of redox and O fugacity on the mineral formation of different elements, and the dynamic crustal evolution that took place during this period of Earth history.
Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.
Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.
We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.
The City of Boca Raton in southeastern Palm Beach County, Florida, is an urban residential area that has sustained a constant population growth with subsequent increase in water use. The Boca Raton network of canals is controlled to provide for drainage of excess water, to maintain proper coastal ground-water levels to prevent saltwater intrusion, and to recharge the surficial aquifer system from which the city withdraws potable water. Most of the water supplied to the Boca Raton canal system and the surficial aquifer system, other than rainfall and runoff, is pumped from the Hillsboro Canal. The Biscayne aquifer, principal hydrogeologic unit of the surficial aquifer system, is highly permeable and there is a close relation between water levels in the canals and the aquifer. The amount of water supplied by seepage from the conservation areas is unknown. Because the Hillsboro Canal flows from Lake Okeechobee and Water Conservation Areas 1 and 2, which are places of more highly mineralized ground water and surface water, the canal is a possible source of contamination. Water samples were collected at 10 canal sites during wet and dry seasons and analyzed for major inorganic ions and related characteristics, nutrients, and trace elements. All concentrations were generally within or less than the drinking-water standards established by the Florida Department of Environmental Protection. The high concentrations of sodium and chloride that were detected in samples from the Boca Raton canal system are probably from the more mineralized water of the Hillsboro Canal. Other water-quality data, gathered from various sources from 1982 through 1991, did not indicate any significant changes nor trends. The data include pesticide and metal analyses of water samples and bottom sediments collected at four canal sites in the Boca Raton study area by the U.S. Geological Survey during 1982-84. The effects of the Hillsboro Canal on the water quality of the Boca Raton canal system are indicated by increased concentrations of sodium, chloride, dissolved solids, and total organic carbon. Concentrations of the constituents in the canal water generally decrease with distance from the Hillsboro Canal pumping station and are the result of dilution by receiving canal waters.
Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.
The Puget Sound basin encompasses the 13,700-square-mile area that drains to the Puget Sound and the adjacent marine waters of Washington State. Well more than 4 million people live within the basin, with numbers continuing to increase, who rely on the basin’s natural resources including groundwater. The Puget Sound Partnership was created by a Washington State statute to implement a science-based recovery of the Puget Sound to help address impacts to these resources. As part of the recovery, the partnership developed the Puget Sound Vital Signs as measures of ecosystem health that guide the assessment of progress toward Puget Sound recovery goals. The Puget Sound Partnership Leadership Council adopted a Drinking Water Vital Sign associated with human health and quality of life, recognizing certain indicators as integral to the sustainability of Puget Sound recovery efforts. One such Vital Sign indicator was the vulnerability of groundwater throughout the aquifers of the Puget Sound basin to elevated nitrate concentrations as defined by the probability of exceeding 2 milligrams/liter (mg/L) at a specific location and well depth. The U.S. Geological Survey (USGS) led the effort to characterize groundwater vulnerability. For this study, groundwater vulnerability refers to a probability with which a contaminant applied at or near the land surface can migrate to the aquifer of interest for a given set of land-use practices. Nitrate concentration data were selected for evaluation because elevated nitrate concentrations are typically caused by anthropogenic activities and have been associated with deleterious impacts on human health. To identify groundwater vulnerability to elevated nitrate concentrations, logistic regression was used to relate anthropogenic (human associated) and natural variables to the occurrence of elevated nitrate concentrations in untreated groundwater from large public water supply system wells found within the Washington State Department of Health Sentry database. Variables that were analyzed included well depth, soil hydraulic conductivity, precipitation, population density, fertilizer application amounts, and land-use types. Statistically significant models that predicted the probabilities of groundwater nitrate concentrations greater than 2 mg/L based on the predictor variables were created for the time periods 2000–04, 2005–09, 2010–14, and 2015–19. For all time periods, well depth and a measure of the abundance of urban and agricultural land over or near the well consistently helped explain the vulnerability of the well to elevated nitrate concentrations defined as a probability of exceeding 2 mg/L of nitrate. Precipitation and (or) soil hydraulic conductivity were also important predictor variables in the models. The models for each time period were used to create maps of groundwater vulnerability at 150- and 300-foot depths throughout the Puget Sound basin. As expected, the most vulnerable locations were associated with shallower well depths and increased agriculture and urban land cover. Across all four time periods, groundwater vulnerability throughout the Puget Sound was low, with probabilities of exceeding 2 mg/L concentrations of nitrate at depths at 150 and 300 feet typically less than 50 percent. Results also found a slight decrease in probabilities of elevated nitrate concentrations throughout the basin over time. More specifically, additional statistical tests found that groundwater with probabilities of less than about 60 percent declined from 2000 to 2019 and represented more than 75 percent of the modeled Puget Sound basin aquifer. Wells with greater than 60 percent probability increased over the same time period but represented only about 25 percent of the aquifer. The maps and statistical analysis presented in the study provide valuable and informative evaluation of the vulnerability of groundwater in the Puget Sound basin to elevated nitrate concentrations. The probability maps do not represent measured nitrate concentrations in groundwater, but rather they present the probability that nitrate concentrations exceed 2 mg/L. The models and predictions from this study are a viable indicator for the Puget Sound Partnership’s Healthy Human Population—Drinking Water Vital Sign. The logistic regression modeling approach presented here benefits water managers by allowing them to assess temporal trends in a range of probabilities, explore vulnerability changes as new regional land cover and anthropogenic data are generated, and distinguish vulnerabilities at different depths within the aquifer.
The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions
Background seismic noise fundamentally sets a lower bound on our ability to record signals arising from earthquakes. The background noise spectrum at a station is a combination of cultural noise, ocean-generated microseism noise, intrinsic instrument self-noise, and the sensitivity of the instrument to nonseismic noise sources. The USArray-Transportable Array Alaska deployed 195 stations across Alaska and parts of Canada (Yukon, British Columbia, and Northwest Territories). These stations were all installed using similar techniques and made use of instruments with similar self-noise levels. As such, this network provides an opportunity to look at how geographic location influences seismic background. Using these broadband stations, we report background noise levels from 0.2 to 75 s period in six discrete bands. By constructing “noise maps,” we depict both spatial and temporal changes in the background noise field. Using these maps, combined with targeted analysis, we infer sources and contributing factors to noise levels in these different period bands. These include cultural noise, the formation of sea ice, seasonal changes in permafrost and wave activity in the Gulf of Alaska, and magnetic field variability. We use this study as an opportunity to review several previous studies examining seismic noise in Arctic regions.
We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.
Anthropogenic nitrogen (N) deposition is unequally distributed across space and time, with inputs to terrestrial ecosystems impacted by industry regulations and variations in human activity. Soil carbon (C) content normally controls the fraction of mineralized N that is nitrified ( ƒ nitrified ), affecting N bioavailability for plants and microbes. However, it is unknown whether N deposition has modified the relationships among soil C, net N mineralization, and net nitrification. To test whether N deposition alters the relationship between soil C and net N transformations, we collected soils from coniferous and deciduous forests, grasslands, and residential yards in 14 regions across the contiguous United States that vary in N deposition rates. We quantified rates of net nitrification and N mineralization, soil chemistry (soil C, N, and pH), and microbial biomass and function (as beta-glucosidase (BG) and N -acetylglucosaminidase (NAG) activity) across these regions. Following expectations, soil C was a driver of ƒ nitrified across regions, whereby increasing soil C resulted in a decline in net nitrification and ƒ nitrified . The ƒ nitrified value increased with lower microbial enzymatic investment in N acquisition (increasing BG:NAG ratio) and lower active microbial biomass, providing some evidence that heterotrophic microbial N demand controls the ammonium pool for nitrifiers. However, higher total N deposition increased ƒ nitrified , including for high soil C sites predicted to have low ƒ nitrified , which decreased the role of soil C as a predictor of ƒ nitrified . Notably, the drop in contemporary atmospheric N deposition rates during the 2020 COVID-19 pandemic did not weaken the effect of N deposition on relationships between soil C and ƒ nitrified . Our results suggest that N deposition can disrupt the relationship between soil C and net N transformations, with this change potentially explained by weaker microbial competition for N. Therefore, past N inputs and soil C should be used together to predict N dynamics across terrestrial ecosystems.