Archive of sidescan-sonar data collected off the North Carolina coast during USGS cruise ATSV99045, 7-27 October, 1999
No abstract available.
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.
No abstract available.
This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.
The long-term producing Coso Geothermal Field (CGF) in California operates over 100 wells tapping into a reservoir characterized by an extensive fracture network, complex fluid pathways, and regular seismic activity. Understanding the interaction between seismicity and injection can shed important light on the hydrothermal characteristics of the field. Here, we analyze 15 years of local seismic and daily operational data from the CGF, identifying a strong correlation between short-term increase in seismicity rate and seasonal volumetric and temperature variations in the reinjected fluid. Furthermore, the seismic footprint during peak injection of colder fluids reveals a near-instantaneous response up to 2 km away from the injection well, too rapid for pore pressure diffusion alone. This short-term and distant response is observed to have directional preference, indicating structural or permeability anisotropy within the reservoir. Additionally, the seismic response correlates with the initial volumetric increase of colder fluids, but also with temperature decrease during stable injection periods, suggesting thermal effects alone can play an important role in triggering distant seismicity.
Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.
Biological community samples were collected at 15 sites in the Upper Colorado River Basin (UCOL) in Colorado as part of the National Water-Quality Assessment (NAWQA) Program. Sites sampled in two physiographic provinces, the Southern Rocky Mountains and the Colorado Plateau, represented agriculture, mining, urban and recreation, and mixed land uses and background conditions. Nine measures of water quality, which include information on nutrients, specific conductance (a surrogate for salinity), trace elements in streambed sediment, pesticides in fish tissue, fish communities, and macroinvertebrate richness and composition and stream habitat were used for comparisons among sites within the two physiographic provinces. Sampling sites from three other NAWQA study units—the Rio Grande Valley, the South Platte River Basin, and the Upper Snake River Basin study units—were categorized on the basis of land use and stream size in order to develop a larger data set for comparison to sites in the UCOL. Three categories of land use—forested (includes mining, urban and recreation, and background), agriculture, and mixed—were used for comparison to the UCOL fixed sites. Results indicated that all sites other than the Colorado River below Baker Gulch (a background site) showed some water-quality characteristics to be significantly affected. Results indicated that the concentrations of cadmium and zinc in streambed sediment at mining land-use sites in the Southern Rocky Mountains physiographic province generally were orders of magnitude higher than streambed-sediment concentrations at the background site. Streambed-sediment concentrations at mining land-use sites in the UCOL were greater than the 75 th percentile of concentrations from sites in the three other NAWQA study units. Fish communities and habitat conditions were degraded at mining land-use sites compared to the background site. Ephemeroptera, Plecoptera, and Trichoptera (EPT) richness and the percentage of EPT were lower at mining land-use sites than at the background site and were less than the 50 th percentile of those for sites from the three other NAWQA study units. Nutrient concentrations at urban and recreation sites in the Southern Rocky Mountain physiographic province generally were greater than concentrations at the background site and generally were between the 25 th and 90 th percentile of concentrations for sites from the three other NAWQA study units. Habitat conditions and fish communities at urban and recreation sites were slightly degraded compared to the background site. EPT richness and the percentage of EPT were lower at urban and recreation sites than at the background site and were between the 25 th and 75 th percentile of those for sites from the three other NAWQA study units. The percentage of Chironomidae, which may be indicative of pollutant-tolerant organisms, was higher at urban and recreation sites than at the background site. Mixed land-use sites in the Southern Rocky Mountains physiographic province had similar nutrient concentrations and similar cadmium and zinc streambed-sediment concentrations. Fish-community degradation index values were very different among the three mixed land-use sites in the Southern Rocky Mountains physiographic province. Larger percentages of omnivores and anomalies such as lesions and deformities at two mixed land-use sites resulted in higher degradation values of the fish community. Agriculture land-use sites had higher concentrations of nutrients and selenium than the background site in the Colorado Plateau physiographic province. Concentrations of p,p '-DDE in fish tissue at agriculture sites were higher than the 75 th percentile of concentrations for sites from the three other NAWQA study units. Fish communities had degradation values near the 75 th percentile for agriculture sites. The percentage of EPT was low at agriculture sites when compared to the background site. Two mixed land-use sites in the Colorado Plateau physiographic province had similar concentrations of nutrients, selenium, and p,p '-DDE, and similar EPT richness and composition. These two sites were located downstream from agricultural and urban activities. Some water-quality measures at these two sites indicated degradation compared to a mixed land-use site upstream from most of the agriculture and urban activities in the Colorado Plateau physiographic province.
An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.
Temporal and spatial distributions of ground-water inflow to, and leakage from Lake Five-O, a softwater, seepage lake in northwestern Florida, were evaluated using hydrologic data and simulation models of the shallow ground-water system adjacent to the lake. The simulation models indicate that ground-water inflow to the lake and leakage from the lake to the ground-water system are the dominant components in the total inflow (precipitation plus ground-water inflow) and total outflow (evaporation plus leakage) budgets of Lake Five-O. Simlulated ground-water inflow and leakage were approximately 4 and 5 times larger than precipitation inputs and evaporative losses, respectively, during calendar years 1989-90. Exchanges of water between Lake Five-O and the ground-water system were consistently larger than atmospheric-lake exchanges. A consistent pattern of shallow ground-water inflow and deep leakage was also evident throughout the study period. The mean time of travel from ground-water that discharges at Lake Five-O (time from recharge at the water table to discharge at the lake) was estimated to be within a range of 3 to 6 years. Flow-path evaluations indicated that the intermediate confining unit probably has a negligible influence on the geochemistry of ground-water inflow to Lake Five-O. The hydrologic budgets and flow-path evaluations provide critical information for developing geochemical budgets for Lake Five-O and for improving the understanding of the relative importance of various processes that regulate the acid-neutralizing capacity of softwater seepage lakes in Florida.
A total of 229 ground-water samples were collected from 215 sites as part of the Ozark Plateaus study unit of the National Water-Quality Assessment Program. These samples were collected from 1993 through 1995 using a network of springs and wells with three scale-dependent components. The first component, the study-unit survey, consisted of 99 randomly selected springs and domestic wells in the Springfield Plateau and Ozark aquifers. The second component, two land-use studies, consisted of 42 springs and domestic wells in a poultry-dominated agricultural area and 40 springs and domestic wells in a cattle-dominated agricultural area overlying the Springfield Plateau aquifer. The third component, the small-watershed study, consisted of 4 springs, 18 domestic wells, and 11 monitoring wells in a small basin within the poultry land-use study area. Samples were analyzed for major ions, nutrients, dissolved organic carbon, methylene blue active substances, tritium, and 88 pesticides and metabolites.The water-quality data from these samples were analyzed with descriptive and statistical methods. Nitrite plus nitrate, which was detected more often and in greater concentrations than any of the other nutrients, ranged from less than 0.05 to 25 milligrams per liter as nitrogen. Nitrite plus nitrate concentrations positively correlated to percent agricultural land use around each site. Median nitrite plus nitrate concentrations generally were greater in samples from springs than in samples from wells. Concentrations of nitrite, ammonia, and ammonia plus organic nitrogen were also affected by land use and also by concentrations of dissolved oxygen in the ground water. Concentrations of phosphorus and orthophosphate probably were affected by land use and also by phosphorus solubility. Pesticides were detected in 80 of 229 samples from 73 of 215 sites. A total of 20 pesticides were detected with a maximum of 5 pesticides detected in any 1 sample. The most commonly detected pesticides were tebuthiuron, atrazine, prometon, desethylatrazine, and simazine. Maximum concentrations ranged from 0.003 to 1.0 microgram per liter. The occurrence and distribution of pesticides were related to land use. Percent agricultural land use was greater for samples with pesticides detected than for samples with no pesticides detected. Pesticides were detected more often in samples from springs than in samples from wells. The occurrence of pesticides also was related to seasonality and chemical characteristics, such as solubility and persistence, of the compounds.
Offshore of the Pacific Northwest of the United States is the Cascadia Subduction Zone, a 1,000-kilometer-long tectonic boundary defined by a large fault, called a megathrust, that extends from the Mendocino Junction off northern California to the Nootka Fracture Zone off Vancouver Island, Canada (U.S. Geological Survey, 2023). The Juan de Fuca and Gorda oceanic plates to the west of this boundary subduct under the North America continental plate to the east. Several other smaller faults that cut through the North America plate crust also affect the region. Although their effects upon Ozette Lake are uncertain, geological evidence for past earthquakes, such as underwater landslides, may be found in Pacific Northwest lakes. Underwater landslides caused by past earthquakes should be well preserved in these relatively undisturbed lake environments. The floor of Ozette Lake, Washington, located along the Pacific coast of the United States, west of the Puget Sound region and about 140 kilometers east of the megathrust was mapped by the U.S. Geological Survey in July of 2019 to search for evidence of past earthquakes. Mapping was completed using a SWATHplus-M 234-kHz interferometric side scan sonar system pole-mounted on the U.S. Geological Survey research vessel San Lorenzo. The system collected full-coverage bathymetric and acoustic backscatter data that were processed to 2-meter spatial resolution (Dartnell and others, 2024). This two-map series displays the results of this mapping. A colored shaded relief bathymetry map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities, respectively, that can be analyzed for evidence of past earthquakes.
We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .
The Geochemical Geodatabase of Sedimentary Strata (Coal, Coal-adjacent Rocks, Tuffaceous Oil Shale, Phosphate-rich Rocks) and Produced Water in the Uinta Region, Utah and Colorado, consists of compiled datasets acquired as part of the Carbon Ore, Rare Earth, and Critical Mineral (CORE-CM) Uinta Region assessment funded by the U.S. Department of Energy (DEFE0032046, 2021–2024; Birgenheier et al., 2024). The CORE-CM assessment focused on providing comprehensive geological and geochemical characterization of current and prospective sedimentary-hosted resources including coal, oil shale, phosphatic limestone, and produced water from oil and gas targets present in eastern Utah and northwestern Colorado (Figure 1). This Data Series includes a geodatabase that consists of analytical geochemical data collected September 2021 through December 2024 via portable X-ray fluorescence (pXRF), and laboratory measured analyses produced by inductively coupled plasma mass spectrometry (ICP-MS) and inductively coupled plasma optical emission spectroscopy (ICP-OES). The coal-related geochemical data are derived primarily from the Cretaceous Blackhawk Formation and Ferron Sandstone of Utah, and the Mesaverde Group of Colorado. Additional non-coal resources assessed include oil shale-bearing strata of the Eocene upper Green River Formation (Utah and Colorado), phosphate-rich limestone of the Permian Park City Formation (Utah) and produced water from oil and gas-bearing strata of the Eocene Green River and Wasatch Formations (Uinta Basin) and the Pennsylvanian Paradox Formation (Paradox Basin) (Table 1). The CORE-CM assessment included a wide range of lithologies present in the coal, oil shale, and phosphate geologic resource systems whether or not the specific lithology has current economic value. Geochemical analyses of produced water from oil and gas wells focused on current and emerging hydrocarbon targets in the central Uinta Basin and northern Paradox Basin. A total of 13,092 geochemical analyses from these geologic systems is provided in the included geodatabase. A series of coal quality data (e.g., composition and maceral analyses) is also included in the database and was digitized from archived coal samples from the Utah Geological Survey (Appendix A).
Human population growth and changing settlement patterns fuel the development of urban fringe lands worldwide, with implications for biodiversity. We conducted a 12-year study of birds in the fast-developing urban fringe lands of the central Puget Sound region, Washington, USA, to examine the effect of development configuration on birds. We hypothesized that lower-intensity conservation developments, compared to higher-intensity planned community developments, would benefit the overall bird community, as well as native forest birds (avoiders of human development) and avian generalist species (adapters to human development), but that higher-intensity planned community developments would benefit synanthropic species (exploiters of human development). We fit single-species and multi-species occupancy models to test these hypotheses. Consistent with our hypotheses, we found that a greater proportion of the overall community, avoiders, and adapters occupied lower-intensity conservation developments compared to higher-intensity planned community developments. However, we did not detect an effect of development type on the exploiter guild, and we found that species in the exploiter guild are variable in their response to the configuration of suburban developments. We also hypothesized that human tolerance guilds would be a useful predictor of individual species responses to development type. This hypothesis was somewhat supported: we found that, for avoiders, 87% of species in the guild had the same response to development type as the overall guild; for adapters, 63% had the same response as the overall guild, and for exploiters, only 44% had the same response as the overall guild. Our results indicate that the configuration of suburban developments can have a meaningful impact on bird communities, particularly on those species that are most sensitive to any level of development.
Lithium (Li), classified as a critical mineral and key in energy storage applications because of its high energy density, faces unprecedented demand. This has driven interest in unconventional sources of Li, especially brines from sedimentary basins. Accordingly, the spatial distribution of basinal brinehosted Li resources across North America is becoming better defined, including the Williston Basin underlying Saskatchewan, Canada. This has led to growing investment and exploration by industry throughout the 2020s; however, despite this interest, the Li potential of basinal brines of the Williston Basin in North Dakota remain underexplored. This work presents the results from the 2025 joint brine sampling campaign between the North Dakota Geological Survey, the United States Geological Survey (USGS), and The University of Texas at Austin. Twenty-seven fluid samples were collected from producing oil and water source wells from eight stratigraphic intervals across western North Dakota and analyzed for major cation and anion concentrations. Similar to southeastern Saskatchewan, the highest Li concentrations were identified in the Devonian Duperow Formation, with concentrations up to 179 mg/L being observed in a well commingled between the Duperow and Red River Formations. Some samples from the overlying Birdbear Formation exceeded 75 mg/L, whereas Li concentrations were lower in the Madison Group (33-44 mg/L) and the Three Forks (50 mg/L), Dawson Bay (43 mg/L), Winnipegosis (23 mg/L), and Red River (40.0 mg/L) Formations. This work indicates Li resource potential exists across the Williston Basin and advances our understanding of the Li distribution in basinal brines which could provide a new domestic source of this critical mineral.
We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.
Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.
The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.
Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.
Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.