SEARCH · USGS Science

Results for “Ocean Modeling”

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 559 records · Page 31Linked to original sources

USGS1998

Evidence for a subsurface ocean on Europa

Ground-based spectroscopy of Jupiter's moon Europa, combined with gravity data, suggests that the satellite has an icy crust roughly 150 km thick and a rocky interior. In addition, images obtained by the Voyager spacecraft revealed that Europa's surface is crossed by numerous intersecting ridges and dark bands (called lineae) and is sparsely cratered, indicating that the terrain is probably significantly younger than that of Ganymede and Callisto. It has been suggested that Europa's thin outer ice shell might be separated from the moon's silicate interior by a liquid water layer, delayed or prevented from freezing by tidal heating; in this model, the lineae could be explained by repetitive tidal deformation of the outer ice shell. However, observational confirmation of a subsurface ocean was largely frustrated by the low resolution (>2 km per pixel) of the Voyager images. Here we present high-resolution (54 m per pixel) Galileo spacecraft images of Europa, in which we find evidence for mobile 'icebergs'. The detailed morphology of the terrain strongly supports the presence of liquid water at shallow depths below the surface, either today or at some time in the past. Moreover, lower- resolution observations of much larger regions suggest that the phenomena reported here are widespread.

Nature↗

USGS1993

Modelling the seasonality of subsurface light and primary production in the Arabian Sea

Seasonal changes in mixed-layer depth and phytoplankton biomass in the Arabian Sea are assessed with climatologies of ship-based hydrographic measurements and ocean-color observations from satellite. At the close of the intermonsoons in November and especially May, the open Arabian Sea resembles the stereotypic, unperturbed tropical ocean, with a thin oligotrophic mixed layer and a pronounced subsurface chlorophyll maximum. Both the northeast and southwest monsoons disrupt this typical tropical hydrography through mixed-layer deepening and eutrophication in the central and northern Arabian Sea. Computations using a spectral model of light penetration suggest that seasonal changes in mixed-layer thickness and phytoplankton concentration result in pronounced fluctuations through the annual cycle in the radiant flux reaching the base of the mixed layer. At the close of the fall and spring intermonsoons the base of the model euphotic zone is in the thermocline across all of the open Arabian Sea. The euphotic zone appears to rise into the mixed layer of the northern Arabian Sea during both the winter and summer monsoons. Strong seasonality in total primary production and its partitioning between the mixed layer and thermocline is predicted byb a photo-synthesis-irradiance model for a site in the western Arabian Sea (14.36 ° N, 57.38° E). Modeled mixed-layer primary production depicts an intense peak for the southwest monsoon and a secondary northeast monsoon peak separated by intermonsoon period of low production. During the fall and spring intermonsoons, in the presence of a subsurface clorophyll maximum, the model estimate of primary production in the thermocline exceeds that in the mixed layer. Our model calculations suggest that the subsurface clorophyll maximum present in the Arabian Sea during the spring intermonsoon is a precursor of the regional, summer, phytoplankton bloom.

Marine Ecology Progress Series↗

USGS2020

Research, monitoring, and evaluation of emerging issues and measures to recover the Snake River Fall Chinook Salmon ESU: January 2019 - December 2019

The portion of the Snake River fall Chinook salmon Oncorhynchus tshawytscha evolutionary significant unit (ESU) that spawns upstream of Lower Granite Dam transitioned from low to high abundance during 19922019 in association with U.S. Endangered Species Act recovery efforts and other federally mandated actions. This annual report focuses on (1) numeric and habitat use responses by natural- and hatchery-origin spawners, (2) phenotypic and numeric responses by natural-origin juveniles, and (3) use of a small unmanned aerial system (sUAS) to search for fall Chinook salmon redds and carcasses. Spawners have located and used most of the available spawning habitat and that habitat is gradually approaching redd capacity. Timing of spawning and fry emergence have been relatively stable, but effects of density dependence are evident in juvenile life stages. Apparent abundance of juvenile fall Chinook salmon has increased and we noted the following responses: parr dispersal from riverine rearing habitat into Lower Granite Reservoir has become earlier; growth rate (g/d) and dispersal size of parr declined; and passage timing of smolts from the two Snake River reaches has become earlier and downstream movement rate faster. These findings coupled with stock-recruitment analyses presented in this report provide evidence for density-dependence in the Snake River reaches and in Lower Granite Reservoir that was influenced by the expansion of the recovery program. The long-term goal is to use this information in a comprehensive modeling effort to conduct action-effectiveness and uncertainty research and to inform Fish Population, Hydrosystem, Harvest, Hatchery, and Predation and Invasive Species Management Research, Monitoring, and Evaluation (RM&E) progams. In 2019, the U.S. Geological Survey (USGS) shifted survey efforts in the Snake River toward deepwater redd searches and fish collection for parentage-based tagging (PBT) analyses because all unmanned aerial system (UAS) activities were suspended by the Department of the Interior two weeks into the spawning season. We counted 81 deepwater redds at 17 of the 29 sites surveyed. Redd depths averaged 3.6 m. We collected a total of 123 live fall Chinook salmon from 16 unique geographic locations that spanned 55 river kilometers. Forty-six fish were recovered at Eureka Bar (rkm 307.1) and Kirby Creek (rkm 352.0), which accounted for 37% of all collected fish in 2019. Most (73 fish) post-spawned salmon were collected from early to mid-November just after peak spawning. A summary of 2019 PBT results can be found in Appendix A.1. In 2019, we PIT tagged subyearling fall Chinook salmon in both the Snake and Clearwater rivers. In the Snake River, we tagged 410 fish with 8-mm tags, 666 fish with 9-mm tags, and 1,082 fish 12-mm tags. During seining, our recapture rate of previously tagged fish was slightly higher in the lower reach at 11.4% than in the upper reach at 10.6%. In an effort to represent more of the population through tagging, we tagged fish as small as 45 mm with 8-mm tags at one site in the upper reach. This allowed us to increase the number of fish tagged in that reach by 17.7%. An additional 10.6% of collected fish could have been tagged in the lower reach had we used 8-mm tags in that reach. In the Clearwater River, we tagged 2,451 subyearlings and recaptured 260 (10.6%) fish in the river and 66 fish (48 tagged by USGS, 18 tagged by the Nez Perce Tribe) at Lower Granite Dam during October to provide information for growth estimation. Within riverine habitats, growth in both length and mass were higher for fish tagged with 8-mm tags than with 9- and 12-mm tags. Estimated growth in length and mass of subyearlings was higher in Lower Granite Reservoir than in riverine habitats. We adapted existing statistical models used to estimate abundance of steelhead and spring/summer Chinook salmon for fall Chinook salmon passing Lower Granite Dam. Run reconstruction efforts to date at Lower Granite Dam for Snake River fall Chinook salmon have provided estimates of the number of returning adults but with no measure of uncertainty about the estimates. The objective of this study was to estimate the abundance, with uncertainty, of marked (coded-wire tagged CWT or adipose clipped) and unmarked fall Chinook salmon past Lower Granite Dam for return years 20032018. Estimating uncertainty is important for informing the state-space life cycle model (Chapter 5), which incorporates both observation and process uncertainty into parameter estimates. The coefficient of variation (CV) for log- abundance was 1.0% or less in all years, whereas the CV for abundance averaged 4.2% and ranged from 1.4% to 10.4% among years. Over the past five years, we have been developing a two-stage state-space life-cycle model for naturally produced fall Chinook salmon in the Snake River basin. Initial efforts focused on generating juvenile and adult abundance estimates, with estimates of uncertainty, for informing the life-cycle model. In this report we 1) describe the statistical life-cycle model and improvements made since our last report to the Independent Scientific Advisory Board (ISAB), 2) estimate the effects of covariates on key demographic parameters, and 3) use the fitted life- cycle model to simulate population trajectories under hydrosystem actions proposed for the NOAA 2020 Biological Opinion (hereafter, the Proposed Action). Major recent advancements to the model include revised juvenile abundance estimates, the ability to estimate smolt-to-adult return rates (SAR) separately for subyearling and yearling juvenile fall Chinook salmon, and improvements in the observation model for estimating age, sex, and outmigration structure in adult returns. We examined the effect of numerous environmental, hydrosystem, and ocean covariates on key demographic parameters but only a few covariates were significant. For the adult-to-juvenile transition, we found the maximum weekly river flows during the winter egg- incubation period had a significant negative effect on the juvenile outmigrant abundance from that brood year. For subyearling outmigrants, percent spill during the summer had a significant positive effect on SAR and the mean winter PDO (Pacific Decadal Oscillation) had a significant negative effect on SAR. For yearling outmigrants, NPGO (North Pacific Gyre Oscillation) had a significant positive effect on SAR. We used the fitted model to simulate population trajectories under the Proposed Action, and the median 10-year geometric mean abundance was 8,222 female spawners (interquartile range: 2,592 26,714). Overall, the probability of quasi- extinction (probability of falling below 50 female spawners for 4 consecutive years) was low, with only 1.6% of all simulations having a quasi-extinction probability >0.95. Although quasi-extinction probability was low, we did not assess the additional effect of climate change, which would be expected to increase quasi-extinction probability.

Idaho, Oregon, Washington↗

USGS2010

Multi-stage origin of the Coast Range ophiolite, California: Implications for the life cycle of supra-subduction zone ophiolites

The Coast Range ophiolite of California is one of the most extensive ophiolite terranes in North America, extending over 700 km from the northernmost Sacramento Valley to the southern Transverse Ranges in central California. This ophiolite, and other ophiolite remnants with similar mid-Jurassic ages, represent a major but short-lived episode of oceanic crust formation that affected much of western North America. The history of this ophiolite is important for models of the tectonic evolution of western North America during the Mesozoic, and a range of conflicting interpretations have arisen. Current petrologic, geochemical, stratigraphic, and radiometric age data all favor the interpretation that the Coast Range ophiolite formed to a large extent by rapid extension in the forearc region of a nascent subduction zone. Closer inspection of these data, however, along with detailed studies of field relationships at several locales, show that formation of the ophiolite was more complex, and requires several stages of formation. Our work shows that exposures of the Coast Range ophiolite preserve evidence for four stages of magmatic development. The first three stages represent formation of the ophiolite above a nascent subduction zone. Rocks associated with the first stage include ophiolite layered gabbros, a sheeted complex, and volcanic rocks vith arc tholeiitic or (roore rarely) low-K calc-alkaline affinities. The second stage is characterized by intrusive wehrlite-clinopyroxenite complexes, intrusive gabbros, Cr-rich diorites, and volcanic rocks with high-Ca boninitic or tholeiitic ankaramite affinities. The third stage includes diorite and quartz diorite plutons, felsic dike and sill complexes, and calc-alkaline volcanic rocks. The first three stages of ophiolite formation were terminated by the intrusion of mid-ocean ridge basalt dikes, and the eruption of mid-ocean ridge basalt or ocean-island basalt volcanic suites. We interpret this final magmatic event (MORB dikes) to represent the collision of an active spreading ridge. Subsequent reorganization of relative plate motions led to sinistral transpression, along with renewed subduction and accretion of the Franciscan Complex. The latter event resulted in uplift and exhumation of the ophiolite by the process of accretionary uplift.

California↗

USGS2000

Parameterizing century to model cultivated and noncultivated sites in the Loess region of western Iowa

One of the main questions remaining for global science involves the cycle of carbon among the atmosphere, oceans, and land. Scientists are trying to better determine the amount of carbon stored in and transferred between these three locations. This task has become more complex because in recent decades the amount of carbon released into the atmosphere has increased due to the burning of fossil fuels and land-use changes. The amount of this increase is greater than the amount of carbon accumulating in the atmosphere and oceans. Many scientists are studying different terrestrial ecosystems to find this 'missing" carbon. One such project is the Mississippi Basin Carbon Project (MBCP) of the U.S. Geological Survey (USGS). MBCP is studying the soils and sediments of the Mississippi River Basin, with an emphasis on understanding human influences on erosion and thus the movement of carbon within a landscape. One goal of the MBCP is to understand, at the field scale, the key processes of erosion and sedimentation, and thus the movement of carbon, in upland areas. Both field measurements and modeling efforts are being used for this purpose. On the modeling front, the Century Model is being used to describe the historical carbon dynamics for two field sites, an agricultural field and uncultivated prairie, located in the loess region of western Iowa. The objective of these modeling efforts is to recreate the carbon dynamics of the upper slope in each of these watersheds. The upper slope represents the area of a hillslope with the greatest potential erosion. This report describes how Century was parameterized to represent these two sites.

Iowa↗

USGS2011

Onshore and offshore geologic map of the Coal Oil Point area, southern California

Geologic maps that span the shoreline and include both onshore and offshore areas are potentially valuable tools that can lead to a more in depth understanding of coastal environments. Such maps can contribute to the understanding of shoreline change, geologic hazards, both offshore and along-shore sediment and pollutant transport. They are also useful in assessing geologic and biologic resources. Several intermediate-scale (1:100,000) geologic maps that include both onshore and offshore areas (herein called onshore-offshore geologic maps) have been produced of areas along the California coast (see Saucedo and others, 2003; Kennedy and others, 2007; Kennedy and Tan, 2008), but few large-scale (1:24,000) maps have been produced that can address local coastal issues. A cooperative project between Federal and State agencies and universities has produced an onshore-offshore geologic map at 1:24,000 scale of the Coal Oil Point area and part of the Santa Barbara Channel, southern California (fig. 1). As part of the project, the U.S. Geological Survey (USGS) and the California Geological Survey (CGS) hosted a workshop (May 2nd and 3rd, 2007) for producers and users of coastal map products (see list of participants) to develop a consensus on the content and format of onshore-offshore geologic maps (and accompanying GIS files) so that they have relevance for coastal-zone management. The USGS and CGS are working to develop coastal maps that combine geospatial information from offshore and onshore and serve as an important tool for addressing a broad range of coastal-zone management issues. The workshop was divided into sessions for presentations and discussion of bathymetry and topography, geology, and habitat products and needs of end users. During the workshop, participants reviewed existing maps and discussed their merits and shortcomings. This report addresses a number of items discussed in the workshop and details the onshore and offshore geologic map of the Coal Oil Point area. Results from this report directly address issues raised in the California Ocean Protection Act (COPA) Five Year Strategic Plan. For example, one of the guiding principles of the COPA five-year strategic plan is to 'Recognize the interconnectedness of the land and the sea, supporting sustainable uses of the coast and ensuring the health of ecosystems.' Results from this USGS report directly connect the land and sea with the creation of both a seamless onshore and offshore digital terrain model (DTM) and geologic map. One of the priority goals (and objectives) of the COPA plan is to 'monitor and map the ocean environment to provide data about conditions and trends.' Maps within this report provide land and sea geologic information for mapping and monitoring nearshore sediment processes, pollution transport, and sea-level rise and fall.

California↗

USGS1984

Models for the deposition of Mesozoic-Cenozoic fine-grained organic-carbon-rich sediment in the deep sea

The widespread occurrence of organic-carbon-rich strata (‘black shales’) in certain portions of Jurassic, Cretaceous and Cenozoic sequences has been well-documented from Deep Sea Drilling Project sites in the Atlantic and Pacific Oceans and from sequences, now exposed on land, originally deposited in the Tethyan ocean. These ancient black shales usually have been explained by analogy with examples of modern deep-sea sediments in which organic matter locally is preserved by (1) increasing the supply of organic matter, (2) increasing the rate of sedimentation, and/or (3) decreasing the oxygen content of the bottom water. However, detailed examination of many black shales reveals characteristics that cannot be explained by simple local models, including: their approximate coincidence in time globally; their occurrence in a variety of different environments, including open oxygenated oceans, restricted basins, deep and shallow water; their interbedding with organic-carbonpoor strata which often dominate a so-called black shale sequence; their deposition by pelagic, hemipelagic, turbiditic and other processes; and the variations in type and amount of organic matter that occur even within the same sequence. A more complex model for the origin of black shales therefore appears most appropriate, in which the cyclic preservation of organic matter depends on the interplay of the three main variables, namely supply of organic matter, sedimentation rate, and deep-water oxygenation, each of which varies independently to some extent. The variation and relative importance of these parameters in individual basins and widespread black shale deposition in general are linked globally and temporally by changes in global sea-level, climate and related changes in oceanic circulation. An important and often overlooked factor for the supply of organic matter to deep-basin sediments is the frequency and magnitude of redepositional processes. The interplay of these variables is discussed in relation to the middle Cretaceous and Cenozoic organic-carbon-rich strata, in particular, which show marked differences in the relative importance of the different variables.

Geological Society of London Special Publications↗

USGS2010

Marine magnetic anomalies

Marine magnetic data have been available for many years from all of the world's oceans, and their contribution to marine geophysics and geology is profound. These data, for example, have allowed charting the age of the ocean floor, reconstruction of the geologic history of the major ocean basins, development of a Cenozoic and Mesozoic timescale of geomagnetic reversals, and speculation on the processes of sea‐floor spreading. Research on these and similar problems actively continued during this quadrennial, but here we discuss only a few topics in which we believe the most significant advances have been made during the last four years: the source of marine magnetic anomalies, the geomagnetic time‐scale, high‐amplitude anomalies, and studies of back‐arc basins. Source of Anomalies Studies of magnetic surveys over continental areas can often be constrained by the magnetic properties of samples collected from below the survey. Until recently, model studies of marine magnetic anomalies have enjoyed the luxury of an inaccessible source. Except for samples scraped from the ocean floor by dredging [e.g., Irving , 1970] and inferences drawn from magnetic studies of ophiolite sequences [e.g., Vine and Moores , 1972], very little information was available for constraining studies of marine magnetic anomalies.

Reviews of Geophysics↗

USGS2017

Retrospective analysis of seasonal ocean growth rates of two sea winter Atlantic Salmon in eastern Maine using historic scales

Substantial declines of anadromous Atlantic Salmon Salmo salar have occurred throughout its range, with many populations at the southern extent of the distribution currently extirpated or endangered. While both one sea winter (1SW) and two sea winter (2SW) spawner numbers for the North American stocks have declined since the 1950s, the decline has been most severe in 2SW spawners. The first months at sea are considered a period of high mortality. However, early ocean mortality alone cannot explain the more pronounced decline of 2SW spawners, suggesting that the second year at sea may be more critical than previously thought. Atlantic Salmon scales collected by anglers and the state agency from 1946 to 2013 from five rivers in eastern Maine were used to estimate smolt age and ocean age of returning adults. Additionally, seasonal growth rates of maiden 2SW spawners were estimated using intercirculi measurements and linear back-calculation methods. Generalized linear mixed models (Gaussian family, log link function) were used to investigate the influence of average sea surface temperature, accumulated thermal units, the Atlantic Multidecadal Oscillation (AMO) and North Atlantic Oscillation indices, smolt age, smolt length, postsmolt growth, and river of origin on growth rate during the oceanic migration of North American Atlantic Salmon. Results suggest that different factors influence salmon growth throughout their oceanic migration, and previous growth can be a strong predictor of future size. Growth was negatively impacted by the phase of the AMO, which has been linked to salmon abundance trends, in early spring following the postsmolt period. This is likely when the 1SW and 2SW stock components separate, and our results suggest that this period may be of interest in future work examining the disproportionate decline in 2SW spawners.

Maine↗

USGS2022

Phenological variation in spring migration timing of adult alewife (Alosa pseudoharengus) in coastal Massachusetts

The timing of biological events in plants and animals, such as migration and reproduction, is shifting due to climate change. Anadromous fishes are particularly susceptible to these shifts as they are subject to strong seasonal cycles when transitioning between marine and freshwater habitats to spawn. We used linear models to determine the extent of phenological shifts in adult Alewife Alosa pseudoharengus as they migrated from ocean to freshwater environments during spring to spawn at 12 sites along the northeastern USA. We also evaluated broadscale oceanic and atmospheric drivers that trigger their movements from offshore to inland habitats, including sea surface temperature, North Atlantic Oscillation index, and Gulf Stream index. Run timing metrics of initiation, median (an indicator of peak run timing), end, and duration were found to vary among sites. Although most sites showed negligible shifts towards earlier timing, statistically significant changes were detected in three systems. Overall, winter sea surface temperature, spring and fall transition dates, and annual run size were the strongest predictors of run initiation and median dates, while a combination of within-season and seasonal-lag effects influenced run end and duration timing. Disparate results observed across the 12 spawning runs suggest that regional environmental processes were not consistent drivers of phenology and local environmental and ecological conditions may be more important. Additional years of data to extend time series and monitoring of Alewife timing and movements in nearshore habitats may provide important information about staging behaviors just before adults transition between ocean and freshwater habitats.

Massachusetts↗

USGS2015

Heat flux from magmatic hydrothermal systems related to availability of fluid recharge

Magmatic hydrothermal systems are of increasing interest as a renewable energy source. Surface heat flux indicates system resource potential, and can be inferred from soil CO 2 flux measurements and fumarole gas chemistry. Here we compile and reanalyze results from previous CO 2 flux surveys worldwide to compare heat flux from a variety of magma-hydrothermal areas. We infer that availability of water to recharge magmatic hydrothermal systems is correlated with heat flux. Recharge availability is in turn governed by permeability, structure, lithology, rainfall, topography, and perhaps unsurprisingly, proximity to a large supply of water such as the ocean. The relationship between recharge and heat flux interpreted by this study is consistent with recent numerical modeling that relates hydrothermal system heat output to rainfall catchment area. This result highlights the importance of recharge as a consideration when evaluating hydrothermal systems for electricity generation, and the utility of CO 2 flux as a resource evaluation tool.

Journal of Volcanology and Geothermal Research↗

USGS2025

Post-wildfire sediment fluxes and turbidity plumes in a coastal-draining watershed

Coastal watersheds impacted by wildfires experience higher erosion resulting in increased sediment delivery to the ocean that alters limiting factors (i.e., light) for marine organisms. With increasing wildfire magnitude and severity, it is critical to explore changes in riverine discharges to the ocean to assess cascading hazards associated with wildfires. In situ data, remotely sensed turbidity data, and hydrological model (Soil and Water Assessment Tool “SWAT”) simulations have been adapted to capture and investigate fire-related land use change impacts on Malibu Creek, California, USA. Modifying SWAT land cover inputs using burn severity data had minimal impact on simulations, requiring additional parameterization for acceptable model performance. Remotely sensed turbidity, in situ discharge, rating curve sediment loads, and SWAT simulated discharge and sediment loads increased following the Woolsey Fire. When compared to in situ and rating curve data in similar non-fire water years, the 2019 Woolsey Fire water year in situ discharge was 1.8 times higher, SWAT simulated discharges were 1.4–1.7 times higher, and rating curve sediment load was 1.3 times higher. However, the SWAT simulated sediment loads were slightly lower (0.8–0.9 times) than rating curve sediment loads in similar non-fire water years. Mean coastal turbidity increased to 18.2 Formazin Nephelometric Unit (FNU) during the first storm post-fire (mean background value of 4.3 FNU). Synergies between methods demonstrated rapid coastal sediment exports (remote sensing) and ongoing erosion post-fire (SWAT). These data are essential to understanding fire-related marine ecological changes and implementing effective management and conservation initiatives.

California↗

USGS2010

Two-dimensional time dependent hurricane overwash and erosion modeling at Santa Rosa Island

A 2DH numerical, model which is capable of computing nearshore circulation and morphodynamics, including dune erosion, breaching and overwash, is used to simulate overwash caused by Hurricane Ivan (2004) on a barrier island. The model is forced using parametric wave and surge time series based on field data and large-scale numerical model results. The model predicted beach face and dune erosion reasonably well as well as the development of washover fans. Furthermore, the model demonstrated considerable quantitative skill (upwards of 66% of variance explained, maximum bias - 0.21 m) in hindcasting the post-storm shape and elevation of the subaerial barrier island when a sheet flow sediment transport limiter was applied. The prediction skill ranged between 0.66 and 0.77 in a series of sensitivity tests in which several hydraulic forcing parameters were varied. The sensitivity studies showed that the variations in the incident wave height and wave period affected the entire simulated island morphology while variations in the surge level gradient between the ocean and back barrier bay affected the amount of deposition on the back barrier and in the back barrier bay. The model sensitivity to the sheet flow sediment transport limiter, which served as a proxy for unknown factors controlling the resistance to erosion, was significantly greater than the sensitivity to the hydraulic forcing parameters. If no limiter was applied the simulated morphological response of the barrier island was an order of magnitude greater than the measured morphological response.

Coastal Engineering↗

USGS2007

Two lithospheric profiles across southern California derived from gravity and seismic data

We present two detailed 2-D density transects for the crust and uppermost mantle across southern California using a linear gravity inversion technique. This technique parameterizes the crust and upper mantle as a set of blocks that are based on published geologic and seismic models. Each block can have a range of densities that are constrained where possible by borehole measurements, seismic velocities, and petrologic data. To further constrain the models, it is assumed that the lithosphere is close to isostatic equilibrium at both ends of the profiles, in the deep ocean and east of the Mojave Desert. We calculate the lithostatic pressure variations field for the whole cross section to rule out the geophysically insignificant solutions. In the linear equation, ρ = a + bV ( V , seismic P-wave velocity; ρ , density), which approximates the mantle density–velocity ( ρ – V ) relationship, different coefficients for b were evaluated. Lower coefficients ( b < 0.2) correspond to an almost purely thermally perturbed mantle, while higher coefficients ( b > 0.3) imply that other effects, such as composition and/or metamorphic changes, play an important role in the mantle. Density models were constructed with the coefficient b ranging from 0 to 0.6. The results indicate that a high b value in the mantle ρ – V relationship is associated with less dense crust in the Mojave block and more dense crust in the Catalina schist block. In the less dense Mojave block, the average density of the whole crust is ∼2.75 g/cm 3 , while that of the lower crust is ∼2.72 g/cm 3 . These densities imply a high silica content in the crust, and a minor fraction of basic rock in the lower crust, or perhaps the absence of a basaltic layer altogether. By comparison, the average density of a typical continental stable platform is ∼2.85 g/cm 3 . Models with higher b coefficients (0.5–0.6) are characterized by a large isostatic imbalance. On the other hand, lower b values (0–0.2) require a consolidated whole crust density in the Mojave Desert of ∼2.78 g/cm 3 , and a lower crust density of ∼2.89 g/cm 3 with mostly basaltic composition. This contradicts the observed, lower V p / V s -ratio in the Mojave Desert associated with mostly felsic and low-density crust. Models with lower b coefficients (0.1–0.2) are characterized by an absence of local Airy compensation beneath the San Gabriel Mountains at the LARSE-1 profile. These, and other non-gravity arguments, suggest optimal solutions to the mantle ρ – V relation of b ∼ 0.2–0.4. This, in turn, means that both thermal and petrological effects occur inside the downwelling of the uppermost mantle high velocity body located beneath the Transverse Ranges. During the development of this mantle downwelling, the basaltic layer of the Mojave block was likely eroded and pulled down into the high velocity body. Those basaltic fragments may have been transformed into eclogites, and this metamorphic change implies a higher b -coefficient density–velocity relationship than would be expected for a purely thermal process.

California↗

USGS2005

Effects of pressure on aqueous chemical equilibria at subzero temperatures with applications to Europa

Pressure plays a critical role in controlling aqueous geochemical processes in deep oceans and deep ice. The putative ocean of Europa could have pressures of 1200 bars or higher on the seafloor, a pressure not dissimilar to the deepest ocean basin on Earth (the Mariana Trench at 1100 bars of pressure). At such high pressures, chemical thermodynamic relations need to explicitly consider pressure. A number of papers have addressed the role of pressure on equilibrium constants, activity coefficients, and the activity of water. None of these models deal, however, with processes at subzero temperatures, which may be important in cold environments on Earth and other planetary bodies. The objectives of this work were to (1) incorporate a pressure dependence into an existing geochemical model parameterized for subzero temperatures (FREZCHEM), (2) validate the model, and (3) simulate pressure-dependent processes on Europa. As part of objective 1, we examined two models for quantifying the volumetric properties of liquid water at subzero temperatures: one model is based on the measured properties of supercooled water, and the other model is based on the properties of liquid water in equilibrium with ice. The relative effect of pressure on solution properties falls in the order: equilibrium constants(K) > activity coefficients (??) > activity of water (aw). The errors (%) in our model associated with these properties, however, fall in the order: ?? > K > aw. The transposition between K and ?? is due to a more accurate model for estimating K than for estimating ??. Only activity coefficients are likely to be significantly in error. However, even in this case, the errors are likely to be only in the range of 2 to 5% up to 1000 bars of pressure. Evidence based on the pressure/temperature melting of ice and salt solution densities argue in favor of the equilibrium water model, which depends on extrapolations, for characterizing the properties of liquid water in electrolyte solutions at subzero temperatures, rather than the supercooled water model. Model-derived estimates of mixed salt solution densities and chemical equilibria as a function of pressure are in reasonably good agreement with experimental measurements. To demonstrate the usefulness of this low-temperature, high-pressure model, we examined two hypothetical cases for Europa. Case 1 dealt with the ice cover of Europa, where we asked the question: How far above the putative ocean in the ice layer could we expect to find thermodynamically stable brine pockets that could serve as habitats for life? For a hypothetical nonconvecting 20 km icy shell, this potential life zone only extends 2.8 km into the icy shell before the eutectic is reached. For the case of a nonconvecting icy shell, the cold surface of Europa precludes stable aqueous phases (habitats for life) anywhere near the surface. Case 2 compared chemical equilibria at 1 bar (based on previous work) with a more realistic 1460 bars of pressure at the base of a 100 km Europan ocean. A pressure of 1460 bars, compared to 1 bar, caused a 12 K decrease in the temperature at which ice first formed and a 11 K increase in the temperature at which MgSO4. 12H2O first formed. Remarkably, there was only a 1.2 K decrease in the eutectic temperatures between 1 and 1460 bars of pressure. Chemical systems and their response to pressure depend, ultimately, on the volumetric properties of individual constituents, which makes every system response highly individualistic. Copyright ?? 2005 Elsevier Ltd.

Geochimica et Cosmochimica Acta↗

USGS2009

Far field tsunami simulations of the 1755 Lisbon earthquake: Implications for tsunami hazard to the U.S. East Coast and the Caribbean

The great Lisbon earthquake of November 1st, 1755 with an estimated moment magnitude of 8.5-9.0 was the most destructive earthquake in European history. The associated tsunami run-up was reported to have reached 5-15??m along the Portuguese and Moroccan coasts and the run-up was significant at the Azores and Madeira Island. Run-up reports from a trans-oceanic tsunami were documented in the Caribbean, Brazil and Newfoundland (Canada). No reports were documented along the U.S. East Coast. Many attempts have been made to characterize the 1755 Lisbon earthquake source using geophysical surveys and modeling the near-field earthquake intensity and tsunami effects. Studying far field effects, as presented in this paper, is advantageous in establishing constraints on source location and strike orientation because trans-oceanic tsunamis are less influenced by near source bathymetry and are unaffected by triggered submarine landslides at the source. Source location, fault orientation and bathymetry are the main elements governing transatlantic tsunami propagation to sites along the U.S. East Coast, much more than distance from the source and continental shelf width. Results of our far and near-field tsunami simulations based on relative amplitude comparison limit the earthquake source area to a region located south of the Gorringe Bank in the center of the Horseshoe Plain. This is in contrast with previously suggested sources such as Marqu??s de Pombal Fault, and Gulf of C??diz Fault, which are farther east of the Horseshoe Plain. The earthquake was likely to be a thrust event on a fault striking ~ 345?? and dipping to the ENE as opposed to the suggested earthquake source of the Gorringe Bank Fault, which trends NE-SW. Gorringe Bank, the Madeira-Tore Rise (MTR), and the Azores appear to have acted as topographic scatterers for tsunami energy, shielding most of the U.S. East Coast from the 1755 Lisbon tsunami. Additional simulations to assess tsunami hazard to the U.S. East Coast from possible future earthquakes along the Azores-Iberia plate boundary indicate that sources west of the MTR and in the Gulf of Cadiz may affect the southeastern coast of the U.S. The Azores-Iberia plate boundary west of the MTR is characterized by strike-slip faults, not thrusts, but the Gulf of Cadiz may have thrust faults. Southern Florida seems to be at risk from sources located east of MTR and South of the Gorringe Bank, but it is mostly shielded by the Bahamas. Higher resolution near-shore bathymetry along the U.S. East Coast and the Caribbean as well as a detailed study of potential tsunami sources in the central west part of the Horseshoe Plain are necessary to verify our simulation results. ?? 2008 Elsevier B.V.

Marine Geology↗

USGS2008

Complex faulting associated with the 22 December 2003 M w 6.5 San Simeon California, earthquake, aftershocks and postseismic surface deformation

We use data from two seismic networks and satellite interferometric synthetic aperture radar (InSAR) imagery to characterize the 22 December 2003 M w 6.5 San Simeon earthquake sequence. Absolute locations for the mainshock and nearly 10,000 aftershocks were determined using a new three-dimensional (3D) seismic velocity model; relative locations were obtained using double difference. The mainshock location found using the 3D velocity model is 35.704° N, 121.096° W at a depth of 9.7±0.7 km. The aftershocks concentrate at the northwest and southeast parts of the aftershock zone, between the mapped traces of the Oceanic and Nacimiento fault zones. The northwest end of the mainshock rupture, as defined by the aftershocks, projects from the mainshock hypocenter to the surface a few kilometers west of the mapped trace of the Oceanic fault, near the Santa Lucia Range front and the >5 mm postseismic InSAR imagery contour. The Oceanic fault in this area, as mapped by Hall (1991) , is therefore probably a second-order synthetic thrust or reverse fault that splays upward from the main seismogenic fault at depth. The southeast end of the rupture projects closer to the mapped Oceanic fault trace, suggesting much of the slip was along this fault, or at a minimum is accommodating much of the postseismic deformation. InSAR imagery shows ∼72 mm of postseismic uplift in the vicinity of maximum coseismic slip in the central section of the rupture, and ∼48 and ∼45 mm at the northwest and southeast end of the aftershock zone, respectively. From these observations, we model a ∼30-km-long northwest-trending northeast-dipping mainshock rupture surface—called the mainthrust—which is likely the Oceanic fault at depth, a ∼10-km-long southwest-dipping backthrust parallel to the mainthrust near the hypocenter, several smaller southwest-dipping structures in the southeast, and perhaps additional northeast-dipping or subvertical structures southeast of the mainshock plane. Discontinuous backthrust features opposite the mainthrust in the southeast part of the aftershock zone may offset the relic Nacimiento fault zone at depth. The InSAR data image surface deformation associated with both aseismic slip and aftershock production on the mainthrust and the backthrusts at the northwest and southeast ends of the aftershock zone. The well-defined mainthrust at the latitude of the epicenter and antithetic backthrust illuminated by the aftershock zone indicate uplift of the Santa Lucia Range as a popup block; aftershocks in the southeast part of the zone also indicate a popup block, but it is less well defined. The absence of backthrust features in the central part of the zone suggests range-front uplift by fault-propagation folding, or backthrusts in the central part were not activated during the mainshock.

California↗

USGS2015

Robust global ocean cooling trend for the pre-industrial Common Era

The oceans mediate the response of global climate to natural and anthropogenic forcings. Yet for the past 2,000 years — a key interval for understanding the present and future climate response to these forcings — global sea surface temperature changes and the underlying driving mechanisms are poorly constrained. Here we present a global synthesis of sea surface temperatures for the Common Era ( CE ) derived from 57 individual marine reconstructions that meet strict quality control criteria. We observe a cooling trend from 1 to 1800 CE that is robust against explicit tests for potential biases in the reconstructions. Between 801 and 1800 CE , the surface cooling trend is qualitatively consistent with an independent synthesis of terrestrial temperature reconstructions, and with a sea surface temperature composite derived from an ensemble of climate model simulations using best estimates of past external radiative forcings. Climate simulations using single and cumulative forcings suggest that the ocean surface cooling trend from 801 to 1800 CE is not primarily a response to orbital forcing but arises from a high frequency of explosive volcanism. Our results show that repeated clusters of volcanic eruptions can induce a net negative radiative forcing that results in a centennial and global scale cooling trend via a decline in mixed-layer oceanic heat content.

Nature Geoscience↗