Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ocean Modelling”

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.

At least 487 records · Page 27Linked to original sources
USGS2013

How old is the Isthmus of Panama?

The Standard Model of the formation of the Isthmus of Panama proposes that final closure occurred at 4–3 Ma. The model is based on evidence from studies of marine stratigraphy, fossil sequences, divergent molecular phylogenies, the timing of the Great American Biological Interchange (GABI), and proxies for marine paleosalinity, paleobathymetry, productivity, and paleotemperatures. The new model uses cooling of magmas in the Cretaceous to Early Miocene Central American Volcanic arc to propose Eocene emergence of the discrete structural blocks of the arc and then U/Pb dating, paleomagnetic pole rotations, and atlantic sea-floor anomalies to reconfigure the blocks for different time slices back to 25 Ma. Closure is proposed at 15 Ma, because by this time the alignment of the blocks leaves no space for trans-isthmian marine passages. We propose that the Indonesian Australian Archipelago (IAA) provides a model for the Central American arc between 15 and 3 Ma because it accounts for the extensive marine interchange between the Pacific and Indian oceans through few and narrow passages while maintaining a complete separation of the terrestrial faunas (Wallace and Lydekker lines) of the two continental platforms of Sunda and Sahul. Using the Indonesian Volcanic arc as a model, we can now accommodate the general tectonic configuration and much of the emergence of the new model, while accounting for the marine fossil record, the molecular evidence of rapid recent speciation, and the delayed Great American Biotic Interchange that the New Model fails to explain.

Bulletin of Marine Science↗
USGS2022

Puerto Rico’s state of the climate 2014-2021

The climate of Puerto Rico is influenced by the changing global climate. The following chapters present the current knowledge of the geophysical and chemical drivers and signals of global climate change as they affect the climate of Puerto Rico and influence the climate-dependent services, risks, and vulnerabilities that govern human well-being. These include sustainable economic development, delivery of ecosystem services, the conservation of natural and cultural resources, resiliency in built and natural systems, and food security. The chapters draw on global expertise of land, atmosphere, and ocean geophysical interactions associated with increasing greenhouse gases that drive global warming and on local scientific expertise, data, observations, and modeled projections. They present the global warming scenario (section 1), the contribution of Puerto Rico to global climate change as GHG emissions and aerosols (section 2), the context of natural climate variability (section 3), observed and projected trends in temperature (section 4), rainfall (section 5), sea level rise (section 6), ocean acidification and sea surface warming (section 7), and the expected implications of warming climate on tropical cyclones affecting Puerto Rico (section 8).

Puerto Rico↗
USGS2019

Heat flow in the Western Arctic Ocean (Amerasian Basin)

From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T-3). The resulting measurements, which are unevenly distributed on Alpha-Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T-3 heat flow is ~54.7 ± 11.3 mW m-2, and Nautilus Basin, including Mendeleev Plain, is the only well-surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m-2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near-linear, implying that conductive heat transport dominates and that near-seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A thermal model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T-3 dataset. .

Journal of Geophysical Research B: Solid Earth↗
USGS2024

Hurricane wave energy dissipation and wave-driven currents over a fringing reef

In 2018, two successive tropical cyclones, Hurricane Hector and Hurricane Lane, generated waves that impacted the Hawaiian Islands. This study investigates wave breaking over a broad fringing reef and aims to quantify the magnitudes and length scales of the corresponding wave-driven circulation using detailed field observations and numerical models corresponding to these wave events. Detailed wave and current measurements were collected across a 1200-m wide cross-reef transect off the coral reef-lined south coast of Moloka’i, Hawai’i. High-resolution numerical model grids were developed to resolve reef features and the coupled Delft3D-SWAN modeling system was applied to simulate spectral wave transformation and wave-driven currents for these two energetic ocean wave events generated by distant passing hurricanes. The results indicate that the wave-driven circulation is generally weak, with current speeds typically less than 0.15 m/s for the wave conditions generated by Hurricane Lane, with significant wave heights up to 1.9 m. Higher energy dissipation rates from larger waves up to 2.5 m breaking during Hurricane Hector resulted in stronger observed currents up to approximately 0.3 m/s. However, the model results show that these currents are confined to the wave-breaking region over the upper fore reef in a narrow (100–300 m) region near the reef crest, and weaker flows of less than 0.1 m/s are generated over the shallow and wide reef flat. Wave heights across the reef flat are less than 0.5 m and are controlled by the tidal water levels. The coral reef structures therefore provide significant protection for the coastline even during large wave and variable sea level conditions. Climate change is likely to increase sea level and storm intensity, the combination of which will influence wave transmission over fringing reefs that may degrade habitat, fueling the need for further research on changing conditions on coral reef-lined coasts.

Hawaii↗
USGS2007

Coastal processes study at Ocean Beach, San Francisco, CA: Summary of data collection 2004-2006

Ocean Beach in San Francisco, California, contains a persistent erosional section in the shadow of the San Francisco ebb tidal delta and south of Sloat Boulevard that threatens valuable public infrastructure as well as the safe recreational use of the beach. Coastal managers have been discussing potential mediation measures for over a decade, with little scientific research available to aid in decision making. The United States Geological Survey (USGS) initiated the Ocean Beach Coastal Processes Study in April 2004 to provide the scientific knowledge necessary for coastal managers to make informed management decisions. This study integrates a wide range of field data collection and numerical modeling techniques to document nearshore sediment transport processes at the mouth of San Francisco Bay, with emphasis on how these processes relate to erosion at Ocean Beach. The Ocean Beach Coastal Processes Study is the first comprehensive study of coastal processes at the mouth of San Francisco Bay.

California↗
USGS2025

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗
USGS2019

Ross Ice Shelf response to climate driven by the tectonic imprint on seafloor bathymetry

Ocean melting has thinned Antarctica's ice shelves at an increasing rate over the past two decades, leading to loss of grounded ice. The Ross Ice Shelf is currently close to steady state but geological records indicate that it can disintegrate rapidly, which would accelerate grounded ice loss from catchments equivalent to 11.6 m of global sea level rise. Here, we use data from the ROSETTA-Ice airborne survey and new ocean simulations, to identify the principal threats to Ross Ice Shelf stability. We locate the tectonic boundary between East and West Antarctica from magnetic anomalies and use gravity data to generate a new high-resolution map of sub-ice-shelf bathymetry. The tectonic imprint on bathymetry constrains sub-ice-shelf ocean circulation, protecting the ice shelf grounding line from moderate changes in global ocean heat content. In contrast, local, seasonal production of warm upper-ocean water near the ice front drives rapid ice shelf melting east of Ross Island, where thinning would lead to faster grounded ice loss from both East and West Antarctic ice sheets. We confirm high modelled melt rates in this region using ROSETTA-Ice radar data. Our findings highlight the significance of both the tectonic framework and local ocean-atmosphere exchange processes near the ice front in determining the future of the Antarctic Ice Sheet.

Nature Geoscience↗
USGS2023

Monitoring and modeling dispersal of a submerged nearshore berm at the mouth of the Columbia River, USA

A submerged, low-relief nearshore berm was constructed in the Pacific Ocean near the mouth of the Columbia River, USA, using 216,000 m 3 of sediment dredged from the adjacent navigation channel. The material dredged from the navigation channel was placed on the northern flank of the ebb-tidal delta in water depths between 12 and 15 m and created a distinct feature that could be tracked over time. Field measurements and numerical modeling were used to evaluate the transport pathways, time scales, and physical processes responsible for dispersal of the berm and evaluate the suitability of the location for operational placement of dredged material to enhance the sediment supply to eroding beaches onshore of the placement site. Repeated multibeam bathymetric surveys characterized the initial berm morphology and dispersion of the berm between September 22, 2020, and March 10, 2021. During this time, the volume of sediment within the berm decreased by about 40%to 127,000 m 3 , the maximum height decreased by almost 60%, and the center of the deposit shifted onshore over 200 m. Observations of berm morphology were compared with predictions from a three-dimensional hydrodynamic and sediment transport model application to refine poorly constrained model input parameters including sediment transport coefficients, bed schematization, and grain size. The calibrated sediment transport model was used to predict the amount, timing, and direction of transport outside of the observed survey area. Model simulations predicted that tidal currents were weak in the vicinity of the berm and wave processes including enhanced bottom stresses and asymmetric bottom orbital velocities resulted in dominant onshore movement of sediment from the berm toward the coastline. Roughly 50% of the berm volume was predicted to disperse away from the initial placement site during the 169 day hindcast. Between 9 and 17% of the initial volume of the berm was predicted to accumulate along the shoreface of a shoreline reach experiencing chronic erosion directly onshore of the placement site. Scenarios exploring alternate placement locations suggested that the berm was relatively effective in enhancing the sediment supply along the eroding coastline north of the inlet. The transferable monitoring and modeling framework developed in this study can be used to inform implementation of strategic nearshore placements and regional sediment management in complex, high-energy coastal environments elsewhere.

Oregon, Washington↗
USGS2012

Implications of regional gravity for state of stress in the earth's crust and upper mantle

Topography is maintained by stress differences within the earth. Depending on the distribution of the stress we classify the support as either local or regional compensation. In general, the stresses implied in a regional compensation scheme are an order of magnitude larger than those corresponding to local isostasy. Gravity anomalies, a measure of the earth's departure from hydrostatic equilibrium, can be used to distinguish between the two compensation mechanisms and thus to estimate the magnitude of deviatoric stress in the crust and upper mantle. Topography created at an oceanic ridge crest or in a major continental orogenic zone appears to be locally compensated. Such features were formed on weak crust incapable of maintaining stress differences much greater than the stress from the applied load. Oceanic volcanoes formed on an already cooled, thickened lithosphere are regionally supported with elastic stresses. Similarly, the broad topographic rise seaward of subduction zones is elastically supported as the lithosphere is bent near the plate margin. Although the implied stress is to some degree dependent on the rheological model assumed, the gravity anomalies and surface deformation produced by these features demonstrate that the upper 30–40 km of the oceanic lithosphere is capable of regionally supporting stress differences in the 100-MPa range. Given certain conditions of load emplacement, continental crust can also support loads regionally over 100-m.y. time scales, but the effects of erosion only allow an estimate of a lower bound on stress. Data from space probes indicate that the upper layers of other terrestrial planets also support topographic-induced stress differences in excess of 100 MPa.

Journal of Geophysical Research Solid Earth↗
USGS2022

Development of projected depth-duration frequency curves (2050–89) for south Florida

Planning stormwater projects requires estimates of current and future extreme precipitation depths for events with specified return periods and durations. In this study, precipitation data from four downscaled climate datasets are used to determine changes in precipitation depth-duration-frequency curves from the period 1966–2005 to the period 2050–89 primarily on the basis of Representative Concentration Pathways 4.5 and 8.5 emission scenarios from the Coupled Model Intercomparison Project Phase 5. The four downscaled climate datasets are (1) the Coordinated Regional Downscaling Experiment (CORDEX) dataset, (2) the Localized Constructed Analogs (LOCA) dataset, (3) the Multivariate Adaptive Constructed Analogs (MACA) dataset, and (4) the Jupiter Intelligence Weather Research and Forecasting Model (JupiterWRF) dataset. Change factors—multiplicative changes in expected extreme precipitation magnitude from current to future period—were computed for grid cells from the downscaled climate datasets containing National Oceanic and Atmospheric Administration Atlas 14 stations in central and south Florida. Change factors for specific durations and return periods may be used to scale the National Oceanic and Atmospheric Administration Atlas 14 historical depth-duration-frequency values to the period 2050–89 on the basis of changes in extreme precipitation derived from downscaled climate datasets. Model culling was implemented to select downscaled climate models that best captured observed historical patterns of precipitation extremes in central and south Florida. Overall, a large variation in change factors across downscaled climate datasets was found, with change factors generally greater than one and increasing with return period. In general, median change factors were higher for the south-central Florida climate region (1.05–1.55 depending on downscaled climate dataset, duration, and return period) than for the south Florida climate region (1–1.4 depending on downscaled climate dataset, duration, and return period) when considering best performing models for both areas, indicating a projected overall increase in future extreme precipitation events.

Florida↗
USGS2019

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Proceedings of the Royal Society B: Biological Sci↗
USGS2017

Investigation of input reduction techniques for morphodynamic modeling of complex inlets with baroclinic forcing

The Mouth of the Columbia River (MCR) is a complex estuary inlet system characterized by a buoyant plume created by high freshwater flows from the Columbia River into the Pacific Ocean. Data obtained during two major field campaigns have resulted in a comprehensive dataset of hydrodynamics and sediment transport under high (2013) and low (2005) river flow conditions. Through the analysis of this data and model simulations obtained with the Delft3D (MCR) model application we explored the importance and effect of stratification on sand-sized sediment for short- and long-term sediment transport simulations. Stratification influences the sediment transport through much of the estuary, and significantly reduces sediment export at the MCR. A correlation analysis reveals that a similar representative tide that best approximates the spring-neap averaged transport can be selected for both stratified and non-stratified flow. This correspondence implies that standard morphodynamic tide schematizations (e.g. Lesser, 2009) may also be valid in the stratified conditions found at MCR and other highly stratified estuaries.

Columbia River↗
USGS2022

Silicate volcanism on Europa’s seafloor and implications for habitability

Habitable ocean environments on Europa require an influx of reactants to maintain chemical disequilibrium. One possible source of reactants is seafloor volcanism. Modeling has shown that dissipation of tidal energy in Europa's asthenosphere can generate melt, but melt formation cannot be equated with volcanism. Melt must also be transported through Europa's cold lithosphere to erupt at the seafloor. Here, we use two models of dike propagation to show that dikes can only traverse the lithosphere if either the fracture toughness of the lithosphere or the flux into the dike is large (>500 MPa m 1/2 or ∼1 m 2 s −1 , respectively). We conclude that cyclic volcanic episodes might provide reactants to Europa's ocean if magma accumulates at the base of the lithosphere for several thousand years. However, if dikes form too frequently, or are too numerous, the magma flux into each will be insufficient, and volcanism cannot support a habitable ocean environment.

Geophysical Research Letters↗
USGS2003

Sediment transport and deposition processes near ocean outfalls in southern California

An urbanized coastal ocean that has complex topography and large-scale atmospheric and oceanographic forcing can contain a variety of sediment and pollutant distribution patterns. For example, the central southern California Bight has two large embayments, Santa Monica and San Pedro Bays, that are connected by a short, very narrow shelf off the Palos Verdes peninsula. The complex topography causes quite different oceanographic and sediment distribution patterns in this fairly small region of the coastal ocean. In addition, three sewage outfalls discharge material over the outer shelf. A large suite of sediment cores was obtained and analyzed for contaminants, physical properties, accumulation rates, and grain sizes. Arrays of instrumented moorings that monitor currents, waves, water clarity, water density and collect resuspended materials were deployed. The data and models developed for the Palos Verdes margin suggest that a large reservoir of DDT and its byproducts exists in the coastal ocean sediment and will continue to be exhumed and transported along the shelf for a long time. On the Santa Monica shelf, very large internal waves, or bores, are generated at the shelf break. The near-bottom currents associated with these waves sweep sediment and the associated contaminants from the shelf onto the continental slope. On the San Pedro margin an initial examination of recent data collected in the coastal ocean does not suggest that bacterial contamination on local beaches is primarily caused by transport of material from the adjacent ocean outfall.

California↗
USGS2011

Contribution of global groundwater depletion since 1900 to sea-level rise

Removal of water from terrestrial subsurface storage is a natural consequence of groundwater withdrawals, but global depletion is not well characterized. Cumulative groundwater depletion represents a transfer of mass from land to the oceans that contributes to sea-level rise. Depletion is directly calculated using calibrated groundwater models, analytical approaches, or volumetric budget analyses for multiple aquifer systems. Estimated global groundwater depletion during 1900–2008 totals ~4,500 km 3 , equivalent to a sea-level rise of 12.6 mm (>6% of the total). Furthermore, the rate of groundwater depletion has increased markedly since about 1950, with maximum rates occurring during the most recent period (2000–2008), when it averaged ~145 km 3 /yr (equivalent to 0.40 mm/yr of sea-level rise, or 13% of the reported rate of 3.1 mm/yr during this recent period).

Geophysical Research Letters↗
USGS2023

Tree-ring derived avalanche frequency and climate associations in a high-latitude, maritime climate

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Alaska↗
USGS1982

Ancient plate boundaries in the Bering Sea region

Plate tectonic models of the Bering Sea suggest that the abyssal Bering Sea Basin is underlain by oceanic crust, a supposition supported by refraction and magnetic data. The oceanic crust is thought to be a remnant of the Kula(?) plate that was isolated within what is now the Bering Sea when the proto-Aleutian arc began to form between the Alaska Peninsula and Kamchatka in late Mesozoic or earliest Tertiary times. Prior to the formation of the Aleutian arc, the Kula(?) plate moved NW, directly underthrusting eastern Siberia; the plate’s eastern edge either obliquely underthrust or slid past the Bering Sea margin along a transform boundary. The Koryak Range in eastern Siberia is composed in part of mélange units that include Palaeozoic and Mesozoic allochthonous blocks juxtaposed within a matrix of Cretaceous sedimentary rocks. Structural trends suggest that these blocks were accreted into the Koryak area from the south along an ancient subduction zone formed by underthrusting of the Kula(?) plate. The base of the Bering Sea continental margin that extends from eastern Siberia to the Alaska Peninsula—the so-called Beringian margin—is underlain by a thick (7–10 km) sedimentary section along the base of the slope. Rocks dredged from the basement exposed farther up the slope (1500–2000 m deep) include shallow-water Upper Jurassic sandstone that is unconformably overlain by shallow-water Eocene to Miocene diatomaceous mudstone. Fauna in the dredge samples indicate that the shelf edge has subsided several kilometres since late Palaeogene time, perhaps in response to the cessation of motion relative to the adjacent oceanic plate and subsequent sediment loading of the oceanic plate. Uplift of the former plate boundary exposed in the Koryak Range occurred principally in late Cenozoic time, and collapse of the adjacent plate boundary, the Beringian margin, began in earliest Tertiary time and has continued to the present. Both tectonic events occurred after the site of active plate collision shifted south to near the present Aleutian Trench. We are uncertain as to why these two ancient, yet adjacent former plate boundaries should behave so differently, i.e. why one area was folded and uplifted while the other was extensionally deformed and subsided, both apparently in response to the cessation of convergent or strike-slip plate motion.

Geological Society, London, Special Publications↗
USGS2004

Preliminary report on the 22 December 2003, M 6.5 San Simeon, California earthquake

The M w 6.5 San Simeon earthquake struck the central California coast on 22 December 2003 at 19:15:56 UTC (11:15:56 am local time.) The epicenter was located 11 km northeast of the town of San Simeon, and 39 km west-northwest of Paso Robles ( Figure 1 ), as reported by the California Integrated Seismic Network (CISN, the California region of the Advanced National Seismic System [ANSS]). The mainshock nucleated at 35.702°N, 121.108°W and a depth of 7.1 km, and the rupture propagated unilaterally to the southeast. The strong directivity of the rupture resulted in a concentration of damage and aftershock activity to the southeast of the hypocenter. The worst earthquake damage occurred in Paso Robles, where two people died in the collapse of an unreinforced masonry building. The accurate and rapid earthquake information provided in near real-time by CISN/ANSS to the Governor's Office of Emergency Services made it possible to focus emergency response in the source area, although the earthquake was felt from San Francisco to Los Angeles. The San Simeon earthquake occurred on a reverse fault striking northwest and most likely dipping to the northeast. Although motion along the Pacific-North America plate boundary in California is dominantly strike-slip, there is a small compressional component through central California. Repeated thrust earthquakes such as the San Simeon event accommodate this compression and build the Coast Ranges. Other recent thrust earthquakes in central California include the 1983 Coalinga (M 6.4) and the 1985 Kettleman Hills (M 6.0) earthquakes. Prior earthquakes in the vicinity of the San Simeon event include a M 5-6 earthquake in 1853, a M 5.7 earthquake in 1906, and the M L 6.2 Bryson earthquake of 1952 ( Figure 1 ) ( McLaren and Savage, 2001 .) The San Simeon earthquake occurred on a previously unknown blind thrust fault. No surface rupture associated with the earthquake has been identified. A number of roads, including State Highway 46, buckled due to the earthquake, but this deformation appears mainly to be failure of road fill due to ground shaking and not the result of tectonic surface rupture. Extrapolation of the fault plane to the surface would roughly align with the surface trace of the Oceanic Fault, but this is thought to be a vertical strike-slip fault. Two models for the kinematics of the region have previously been proposed. The first is a fault-propagation fold model developed by Namson and Davis ( 1990 ) for the Santa Lucia mountains ∼30 km to the southeast of the San Simeon sequence. The mainshock geometry is similar to, although more steeply dipping than, the main blind thrust of this model, implying that this model may be applicable to the San Simeon region as well. The second is the model of McLaren and Savage ( 2001 ), in which the region is dominated by strike-slip faulting with shortening on high-angle reverse faults. This model also may be applicable, although the dip of the San Simeon mainshock is shallower than predicted. The San Simeon earthquake was followed by a vigorous aftershock sequence, with 165 events above M 3 reported by CISN within the first week of the mainshock. Although the event triggered many aftershocks, it did not significantly impact the seismicity rates of other nearby faults such as the San Andreas Fault and the San Simeon-Hosgri fault zone. The only triggered seismicity seems to be a few small events within the mainshock coda at the Geysers geothermal area, north of San Francisco. The San Simeon earthquake did, however, trigger shallow creep on the San Andreas Fault at Parkfield and hydrologic changes in hot springs in Paso Robles.

California↗