USGS Science⌕ Search

SEARCH · USGS Science

Results for “Annual Review of Environment and Resources”

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.

1,696 records · Page 48Linked to original sources
USGS2024

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California↗
USGS2025

Noble and base metal distribution and processes affecting ore tenors in the disrupted lower stratigraphy of the Stillwater Complex, USA

Exploration continues for contact-style Ni-Cu sulfide and chromitite-associated PGE mineralization in ultramafic rocks of the Stillwater Complex. At the Iron and Chrome Mountain areas, massive sulfides occur along the complex’s footwall contact and anomalous concentrations of PGE+Au are associated with the three lowermost chromitite seams. Southeast of Chrome Mountain, magmatic layering is highly disrupted by the presence of faults, magmatic breccias, serpentinized discordant dunites, pyroxenite pegmatoids, and disaggregated chromitite seams. The bulk rock chemistry, sulfide chemistry, and noble metal mineralogy of samples from this area were examined to determine the deportment of PGE and processes that led to enrichments in PGE, Au, Cu, Co and Ni. Results show that a sulfide liquid was the principal collector of PGE. If sulfide liquid was initially deposited with chromite, it was disaggregated or redistributed by subsequent melt or fluid infiltration, which may have resulted in the offset of peak PGE(+Cu, Ni) from peak Cr 2 O 3 concentrations, and upgraded PGE tenors. Upon cooling, PGE exsolved from sulfides to form discrete bismuth tellurides, arsenides, arsenic sulfides, antimonides, and alloys, commonly along the margins of sulfide globules. Calculated metal tenors are highest in the disseminated sulfides southwest of Chrome Mountain, whereas massive and net-textured sulfides near the Iron Mountain-Camp zone represent monosulfide solid solution cumulates. At progressively shallower levels, higher metal tenors combined with lower S/Se ratios are consistent with increasing R-factors from 100 to 100,000. Serpentinization and talc-tremolite alteration resulted in S loss through partial replacement of sulfides by secondary silicate+carbonate+magnetite+sulfide assemblages, further upgrading Ni-Cu-PGE tenors. The present work shows that processes responsible for the disruption of magmatic layering and post-magmatic fluid alteration along the intrusion’s lower contact led to noble and base metal enrichments.

Montana↗
USGS2026

Divisions of geologic time—Major chronostratigraphic and geochronologic units

Introduction The reports and maps of our Nation’s geological surveys inform and benefit the public, private industry, government officials, and scientists. The use of clear and consistent nomenclature and classifications can improve communication of data and interpretations. Since 1899, the U.S. Geological Survey (USGS) Geologic Names Committee (GNC) has been responsible for defining standards that promote uniform geologic nomenclature and classifications among geoscientists. The GNC periodically publishes a geologic time scale, the “Divisions of Geologic Time,” that serves as the national standard for USGS publications (for example, refer to Orndorff and others, 2023). Authors may use other published geologic time scales, such as those of the Geological Society of America (GSA) or the International Commission on Stratigraphy (ICS), provided that they are clearly specified and referenced. Access to the USGS, GSA, and ICS geologic time scales is also available from the U.S. Geologic Names Lexicon (Geolex) website https://ngmdb.usgs.gov/Geolex/stratres/timescales. The geologic time scale serves a dual purpose by enabling authors to distinguish earth material units by position (chronostratigraphic) and time (geochronologic), as outlined in order of decreasing rank.

Fact Sheet↗
USGS2025

Bayesian ETAS modeling for the Pacific Northwest: Uncovering effects of tectonic regimes, regional differences, and swarms on aftershock parameters

The Pacific Northwest (PNW) of North America has high seismic hazard due to numerous earthquake sources under populated areas. It hosts several tectonic regimes and subregional seismic zones that are hypothesized to have different patterns of earthquake and aftershock occurrence. It is also predisposed to earthquake swarms, which can complicate the statistical modeling of these patterns. We present the first statistical seismicity model of the PNW catalog using the epidemic‐type aftershock sequence (ETAS) framework. We develop a Bayesian inference procedure that provides a stable estimation of both ETAS parameters and their uncertainties for different sets of PNW earthquakes, even those with very sparse catalogs. The Bayesian approach allows us to investigate how parameter estimates change between the intraslab and crustal tectonic regimes, the northern and southern PNW, and when swarms are included and excluded from the catalog. We also utilize our Bayesian framework to calculate parameter estimates under different prior beliefs about PNW seismicity, as well as to propagate catalog measurement errors into ETAS parameter estimates. We discuss the implications of parameter differences across the region for aftershock forecasting for the PNW.

Pacific Northwest↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2024

Collision structures of the Prince William terrane and Chugach terrane docking along the Shumagin and Unimak convergent margins, Alaska, USA

Western Alaska’s convergent margins are composed of tectonostratigraphic terranes. On land, terrane assembly is recognized along boundaries or sutures between neighboring geologic elements with distinctly different origins. In marine areas where rock outcrops are covered by sediment, recognizing terrane sutures is problematic. A fault in seismic dip line 5 of the ALEUT project has been interpreted as a terrane suture. It is imaged intermittently down to the 30+-km-deep plate interface. Processing of ALEUT strike line 7 revealed the suture at ~18 km depths extending 300 km along the margin. Upper structures in line 5 are like the structures of adjacent seismic transects where imaging is only 8−10 km deep. They were previously not recognized as the upper reaches of terrane sutures and show structural details obscured at greater depths. The composite data are the basis for a simple tectonic model of terrane docking.

Alaska↗
USGS2026

Earthquake-hazard exposure of residents with potential access and functional needs in the United States

Earthquake response plans and earthquake early warning (EEW) systems designed for general populations may not consider potential access and functional needs (AFN) of individuals with physical, sensory, cognitive, or social limitations. Previous efforts to map the distribution of these populations have focused on social-vulnerability indices that ignore or oversimply these limitations. The descriptive and exploratory analysis summarized in this United States (U.S.) case study addresses this gap by identifying and integrating spatially explicit data for AFN-related residential populations, earthquake hazards, and county and county equivalents for the conterminous U.S., Alaska, Puerto Rico, and Hawaii. We focus on 13 AFN-related attributes that relate to an individual's ability to access information contained in an EEW alert, to understand and process earthquake information or observed ground shaking, and to take self-protective actions based on this information and physical cues of an earthquake. Depending on the demographic attribute, there are millions to tens of millions of U.S. residents with AFN-related attributes in areas considered to have varying likelihoods (2%, 10%, and 50%) of exceedance of a damaging earthquake in the next 50 years. Although these amounts represent low percentages at the national level, the percentage of individuals with AFN-related attributes in many counties and county equivalents substantially exceeds national percentages. No one county, county equivalent, U.S. state, or U.S. territory has the highest percentage of individuals in all AFN-related attributes; therefore, future efforts to increase individual resilience to earthquakes may benefit from understanding the local context of individuals with potential access and functional needs.

International Journal of Disaster Risk Reduction↗
USGS2026

Non-destructive sampling of larval amphibians shows tail tissues reflect whole-body methylmercury concentrations and trophic-related differences in bioaccumulation

Mercury (Hg), and specifically methylmercury (MeHg), is a contaminant of global concern to humans and wildlife, but there is still limited understanding of its effects on many taxa. Amphibians are closely associated with water, where Hg is converted to MeHg, and high concentrations of MeHg can reduce survival of amphibians in the wild. However, limited non-lethal proxies exist for estimating MeHg bioaccumulation in amphibian larvae, which could be used to reduce destructive sampling and allow sampling of imperiled species. We evaluated whether tail clips were indicative of whole-body MeHg concentrations of larval amphibians, whether dragonflies were an effective bioindicator of Hg concentrations of larval amphibians, and compared MeHg bioaccumulation in predatory caudate (salamanders, newts) and non-predatory anuran (frogs, toads) larvae. Tail-clip and whole-body MeHg concentrations were strongly correlated, especially for caudates. There was high variability among sites and species, but caudates had higher MeHg than anurans, and MeHg concentrations of caudates and anurans were strongly correlated. Dragonflies were a weak indicator of MeHg for anurans but strong indicators for larval caudates, likely because caudates and dragonflies are carnivorous and occupy similar trophic positions. Our study revealed that sampling tails can be an effective non-destructive index of whole-body MeHg for larval amphibians and demonstrates differences in MeHg bioaccumulation related to trophic position in larval amphibians and dragonflies. All post-embryo life stages for a broad suite of amphibian species can now be sampled non-lethally for MeHg, allowing for a broader examination of population-level effects and focus on imperiled species that cannot be sampled destructively.

conterminous United States↗
USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗
USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗
USGS2025

Sources and risk factors for nitrate, pathogens, and fecal contamination of private wells in rural southwestern Wisconsin, USA

Household well water can be degraded by contaminants from the land's surface, but private well owners lack means to protect the source water from neighboring disturbances. Rural residents of southwestern Wisconsin, USA, rely on private well water, and the combination of land use and fractured carbonate bedrock makes groundwater vulnerable to contamination. To identify the extent, sources, and risk factors of private well contamination, randomly selected wells sampled during two-day periods in fall (n = 301) and spring (n = 529) were analyzed for nitrate and indicator bacteria, and a subset (n = 138) was sampled across four seasonal events for analysis of pathogens and microbial source tracking markers by quantitative polymerase chain reaction. Risk factors representing land use, hydrology, geology, and well construction were analyzed for associations with contamination in multivariable models. The importance of risk factors varied by contaminant, illustrating the multifaceted nature of rural groundwater quality. Nitrate contamination was associated with agricultural land use, and wells with casings that extended below a shale aquitard accessed less contaminated water than those drawing water from above it. Human fecal microbes were detected in 64 wells (46%), and rainfall was the key risk factor for contamination, indicating that wastewater from septic systems was available to contaminate wells when transport conditions were favorable. Manure microbes from cattle/ruminants and pigs were detected in 33 and 13 wells, respectively, and concentrations increased with the hectarage of cultivated land near wells. Pathogen genes for viruses, bacteria, and protozoa were detected in 66 wells (48%), including more detections of zoonotic than human-specific pathogens, and human Bacteroides , an indicator of wastewater, was an equivocal predictor of pathogen presence in private wells. Characterizing important elements of the setting, like geology, and identifying sources and risk factors for contaminants can inform landscape-level policies to protect groundwater quality.

Wisconsin↗
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

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.

Mandalay↗
USGS2024

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Geological Society, London, Special Publications↗
USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗
USGS2023

Preliminary geologic map of the Southern Santa Rosa Mountains and Borrego Badlands, San Diego County, Southern California

This investigation delineates the geologic framework of an area of 75 square kilometers (km 2 ) located west of the Salton Sea in southern California (fig. 1, on sheet 1). The study area encompasses the south flank of the Santa Rosa Mountains and the eastern part of the Borrego Badlands (sheet 1). In this study area, regionally important stratigraphic and structural elements collectively inform the late Cenozoic geologic evolution of the Anza-Borrego sector of the Salton Trough province. Critical stratigraphic and structural elements in the map area include the following: The well exposed sequence of late Cenozoic, nonmarine sedimentary rocks that filled the Anza-Borrego subbasin (fig. 1) of the Salton Trough; A tectonic boundary that—in the southern Santa Rosa Mountains—separates the sedimentary strata from underlying crystalline rocks of Peninsular Ranges type. This tectonic boundary, named the West Salton Detachment Fault System by Axen and Fletcher (1998), is projected to underlie all late Cenozoic sedimentary strata in the Anza-Borrego subbasin of the Salton Trough; A variety of transpressional, transtensional, and strike-slip structures that have deformed the late Cenozoic sedimentary strata and collectively guided syntectonic and posttectonic depositional events within the Anza-Borrego subbasin of the Salton Trough; and The southeasternmost surface expression of the Clark Fault, a major strand of the dextral San Jacinto Fault Zone. Geologic mapping and analysis for this investigation focused on clarifying geologic relations among these four stratigraphic and structural aspects in the map area.

California↗