SEARCH · USGS Science

Results for “JGR - Planets”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗

USGS2024

Seasonality of retreat rate of a wave-exposed marsh edge

Wave-driven erosion of marsh boundaries is a major cause of marsh loss, but little research has captured the effect of seasonal differences on marsh-edge retreat rates to illuminate temporal patterns of when the majority of this erosion is occurring. Using five surface models captured over a study year of a marsh with a steep escarped boundary in South San Francisco Bay, we find a pronounced seasonal signal, where rapid marsh retreat in the spring and summer is driven by a strong sea breeze but little change is found in the marsh-edge position in the fall and winter. We found accretion in the mudflat transition region close to the marsh boundary in the calmer seasons however, suggesting intertwined morphodynamics of mudflats and the eroding marsh-scarp. We observed large spatial heterogeneity in retreat rates within seasons, but less on longer (annual and decadal) timescales. The relationship between marsh-edge retreat rates and properties of the wave field nearby is explored and contextualized against extant relationships, but our results speak to the difficulty in addressing spatial erosion/accretion variability on short (seasonal) timescales with simple models.

JGR Earth Surface↗

USGS2024

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface↗

USGS2023

Calibration and validation for the Surface Biology and Geology (SBG) mission concept: Recommendations for a multi-sensor system for imaging spectroscopy and thermal imagery

The primary objective of the National Aeronautics and Space Administration (NASA) Surface Biology and Geology (SBG) mission is to measure biological, physical, chemical, and mineralogical features of the Earth's surface, realizing a key conceptual component of the envisioned NASA Earth System Observatory (ESO). SBG is planned to launch as a two-platform mission in the late 2020s, the first of the ESO satellites. Targeted science and applications objectives based on observations of the Earth's SBG helped to define the mission architecture and instrument capabilities for the SBG mission concept. These objectives further drove the need for enabling change detection and trending of surface biological and geological features. These needs implied fundamental calibration goals to achieve the necessary science data quality characteristics. To meet those goals, calibration and validation pre-launch and on-orbit methods formed a basis of the calibration and validation concept, including the combined use of on-board references, vicarious techniques, and routine lunar imaging. International collaboration with space agencies in other countries, an important feature of the recommended SBG mission architecture, uncovered and emphasized the need for inter-calibration techniques that underscored the importance of collaborative instrument characterization data sharing and the use of common calibration references that are International System of Units (SI) traceable in pre-launch and post-launch on orbit calibration mission phases. International collaboration through the use of terrestrial and aquatic networks on six continents for vicarious calibration and validation activities will further assure necessary science data quality while in orbit.

JGR Biogeosciences↗

USGS2023

Biophysical drivers of coastal treeline elevation

Sea level rise is leading to the rapid migration of marshes into coastal forests and other terrestrial ecosystems. Although complex biophysical interactions likely govern these ecosystem transitions, projections of sea level driven land conversion commonly rely on a simplified “threshold elevation” that represents the elevation of the marsh-upland boundary based on tidal datums alone. To determine the influence of biophysical drivers on threshold elevations, and their implication for land conversion, we examined almost 100,000 high-resolution marsh-forest boundary elevation points, determined independently from tidal datums, alongside hydrologic, ecologic, and geomorphic data in the Chesapeake Bay, the largest estuary in the U.S. located along the mid-Atlantic coast. We find five-fold variations in threshold elevation across the entire estuary, driven not only by tidal range, but also salinity and slope. However, more than half of the variability is unexplained by these variables, which we attribute largely to uncaptured local factors including groundwater discharge, microtopography, and anthropogenic impacts. In the Chesapeake Bay, observed threshold elevations deviate from predicted elevations used to determine sea level driven land conversion by as much as the amount of projected regional sea level rise by 2050. These results suggest that local drivers strongly mediate coastal ecosystem transitions, and that predictions based on elevation and tidal datums alone may misrepresent future land conversion.

Delaware, Maryland, Virginia↗

USGS2023

Submersed macrophyte density regulates aquatic greenhouse gas emissions

Shallow freshwater ecosystems emit large amounts of greenhouse gases (GHGs), such as carbon dioxide (CO 2 ) and methane (CH 4 ), yet emissions are highly variable. The role that aquatic macrophytes play in regulating aquatic GHG emissions is uncertain despite their ability to dominate shallow waterbodies. Here, we studied the effects of submersed macrophyte ( Ceratophyllum demersum ) density on CO 2 and CH 4 concentrations and fluxes. We conducted a 61-days experiment using mesocosms containing one of the following C . demersum density treatments: 0, 10, 20, or 30 individual shoots ( n = 3). We found that high density C . demersum had the highest CO 2 and CH 4 surface water concentrations and emissions while there was no significant difference in CH 4 in the low and medium densities and no plant control. The high density treatment lost biomass over the course of the experiment, indicating die-off and additions of organic matter to the sediment. High organic matter loading and low dissolved oxygen likely stimulated GHG production in the high density treatment. Our results emphasize that submersed macrophyte density and periods of growth and dieback are important in regulating GHG emissions, which may help explain why shallow waterbodies are high yet variable sources of GHGs to the atmosphere.

JGR Biogeosciences↗

USGS2024

The effect of drying boreal lakes on plants, soils, and microbial communities in lake margin habitats

Decadal scale lake drying in interior Alaska results in lake margin colonization by willow shrub and graminoid vegetation, but the effects of these changes on plant production, biodiversity, soil properties, and soil microbial communities are not well known. We studied changes in soil organic carbon (SOC) and nitrogen (N) storage, plant and microbial community composition, and soil microbial activities in drying and non-drying lakes in the Yukon Flats National Wildlife Refuge. Historic changes in lake area were determined using Landsat imagery. Results showed that SOC storage in drying lake margins declined by 0.13 kg C m −2 yr −1 over 30 years of exposure of lake sediments, with no significant change in soil N. Lake drying resulted in an increase in graminoid and shrub aboveground net primary production (ANPP, +3% yr −1 ) with little change in plant functional composition. Increases in ANPP were similar in magnitude (but opposite in sign) to losses in SOC over a 30-year drying trend. Potential decomposition rates and soil enzyme activities were lower in drying lake margins compared to stable lake margins, possibly due to high salinities in drying lake margin soils. Microbial communities shifted in response to changing plant communities, although they still retained a legacy of the previous plant community. Understanding how changing lake hydrology impacts the ecology and biogeochemistry of lake margin terrestrial ecosystems is an underexamined phenomenon with large impacts to landscape processes.

Alaska↗

USGS2024

Vegetation loss following vertical drowning of Mississippi River deltaic wetlands leads to faster microbial decomposition and decreases in soil carbon

Wetland ecosystems hold nearly a third of the global soil carbon pool, but as wetlands rapidly disappear the fate of this stored soil carbon is unclear. The aim of this study was to quantify and then link potential rates of microbial decomposition after vertical drowning of vegetated tidal marshes in coastal Louisiana to known drivers of anaerobic decomposition altered by vegetation loss. Profiles of potential CH 4 and CO 2 production (surface to 60 cm deep) were measured during anaerobic incubations, organic matter chemistry was assessed with infrared spectroscopy, and soil porewater nutrients and redox potentials were measured in the field along a chronosequence of wetland loss. After vertical drowning, pond soils had lower redox potentials, higher pH values, lower soil carbon and nitrogen concentrations, lower lignin: polysaccharide ratios, more NH 4 + and PO 4 3− , and higher rates of potential CO 2 release than vegetated marsh soils. Potential CH 4 production was similar in vegetated marshes and open water ponds, with depth-dependent decreases in CH 4 production as soil carbon concentrations increased. In these anoxic soils, vegetation loss exerts a primary control on decomposition rates because flooding drives sustained increases in porewater nutrient availability (NH 4 + and PO 4 3 , dissolved organic carbon) and decreases in redox potential (from −150 to −500 mV) that lead to higher potential CO 2 fluxes within a few years. Without new carbon inputs following wetland loss, the sustained decomposition in open water ponds may lead to losses of stored soil carbon and could influence global carbon budgets.

Louisiana↗

USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗

USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗

USGS2025

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

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

Oklahoma↗

USGS2024

Large eddy simulation of cross-shore hydrodynamics under random waves in the inner surf and swash zones

A 3D large eddy simulation coupled with a free surface tracking scheme was used to simulate cross-shore hydrodynamics as observed in a large wave flume experiment. The primary objective was to enhance the understanding of wave-backwash interactions and the implications for observed morphodynamics. Two simulation cases were carried out to elucidate key processes of wave-backwash interactions across two distinct stages: berm erosion and sandbar formation, during the early portion of a modeled storm. The major difference between the two cases was the bathymetry: one featuring a berm without a sandbar (Case I), and the other, featuring a sandbar without a berm (Case II) at similar water depth. Good agreement (overall Willmott's index of agreement greater than 0.8) between simulations and measured data in free surface elevation, wave spectrum, and flow velocities validated the model skill. The findings indicated that the bottom shear stress, represented by the Shields parameter, was significant in both cases, potentially contributing substantial sediment transport. Notably, the occurrence of intense wave-backwash interactions were more frequent in the absence of a sandbar. These intense wave-backwash interactions resulted in a pronounced horizontal pressure gradient, quantified by high Sleath parameters, exceeding the criteria for momentary bed failure. Additionally, a more vigorous turbulence-bed interaction, characterized by near-bed turbulent kinetic energy, was observed in the case lacking a sandbar, potentially augmenting sediment suspension. These insights are pivotal in understanding the mechanisms underlying berm erosion and how sandbar formation serves to protect further beach erosion.

JGR Oceans↗

USGS2024

Microplastic and associated black particles from road-tire wear: Implications for radiative effects across the cryosphere and in the atmosphere

The environmental effects of airborne micro- and nano-size plastic particles are poorly understood. Microscopy and chemical analyses of atmospherically deposited particles on snow surfaces at high elevation (2,865–3,690 m) in the Upper Colorado River basin (UCRB; Colorado Rocky Mountains) revealed the presence of black substances intimately associated with microplastic fibers, particles interpreted to have originated as tire matter. Identical and similar particles occur in shredded tires and road-surface samples. The substance responsible for the black color of all tires is carbon black, a graphitic light-absorbing tire additive produced by hydrocarbon combustion that homogeneously permeates the mixture of tire polymers and other additives. Such black tire matter may thus exert radiative effects closely similar to those of black carbon. The presence in snow of many organic compound types common to tires, measured by two-dimensional gas chromatography, suggests that atmospherically deposited black road-tire-wear matter is among the light-absorbing particulates that advance the onset and rate of snow melt in the UCRB. The mass of road-tire-wear particles shed from vehicles may be estimated by multiplying measured amounts of eroded tire-per-distance traveled by vehicular distances. Under a combination of measurements and assumptions about the amounts and radiative properties of atmospheric tire-wear particles, the radiative effects of these particles might add about 10%–30% to those effects from black carbon, an estimate ripe for revision. On regional and global scales, the amounts and effects of emitted and deposited tire-wear matter likely vary by factors of geographic source, transport pathway, and depositional setting.

JGR - Atmospheres↗

USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗

USGS2025

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2025

Shoreline seasonality of California’s beaches

We report on remote sensing techniques developed to characterize seasonal shoreline cycles from satellite-derived shoreline measurements. These techniques are applied to 22-yr of shoreline measurements for over 777 km of beach along California's 1,700-km coast, for which the general understanding is that shorelines exhibit winter-narrow and summer-recovery seasonality. We find that approximately 90% of beach transects exhibit significant and recurring seasonal cycles in the shoreline position. Seasonal shoreline excursions are twice as large in northern and central California (17.5–32.2 m) than southern California (7.3–15.9 m; interquartile ranges). Clustering analyses were effective at characterizing the temporal patterns of the seasonality, revealing that ∼459 km of beach (59%) exhibit winter-narrow conditions, whereas ∼189 km (24%) and ∼50 km (6.4%) exhibit spring-narrow and summer-narrow conditions, respectively. These spring- and summer-narrow conditions are most common in southern California, where they represent over half of the total length of beach shoreline. Multivariate analyses reveal that wave climate and geomorphic setting are significantly related to the magnitude and timing of shoreline seasonal cycles. Combinations of these variables explain 44% of the seasonality variance of the complete data set and 85% of the variance for a subset of 93 long (>1 km) continuous beaches. We conclude that diversity in waves and geomorphic setting along California cause a broad range of seasonal patterns in the shoreline. Combined, this indicates that the overly generalized “winter-narrow/summer-recovery” conventions for California beaches are not expressed universally and that shoreline seasonality is far more diverse than these simple canonical rules.

California↗