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 127 records · Page 7Linked to original sources
USGS2023

Processes controlling coastal erosion along Cape Cod Bay, MA

Cape Cod Bay, MA, is a semi-enclosed embayment in the northeastern United States, open on the north to the Gulf of Maine. The coastline experiences impacts typically from strong Nor’easter storms that occur in the late fall or winter months, with some sections of this coastline being affected more severely than others. We investigate the processes that cause spatial variability in storm impacts by using geophysical surveys, shoreline-change analysis, and numerical modeling. We simulated the Gulf of Maine and Cape Cod Bay from Jan–April, 2021, using the COAWST modeling system, including ocean, wave, infragravity wave, and sediment transport models. Results identify bay-scale circulation of alongshore sediment fluxes and convergences at regional shoals. Nearshore modeling also revealed zones of increased wave heights that correlate with regions of increased erosion and coastal angle orientation. Modeled and computed shoreline-change have some correlation but the model does not capture all the variability. This overall approach can potentially be used for other coastal locations to identify regions of storm impacts and to manage coastal resources.

Massachusetts↗
USGS2015

Hierarchical stochastic modeling of large river ecosystems and fish growth across spatio-temporal scales and climate models: the Missouri River endangered pallid sturgeon example

We present a hierarchical series of spatially decreasing and temporally increasing models to evaluate the uncertainty in the atmosphere – ocean global climate model (AOGCM) and the regional climate model (RCM) relative to the uncertainty in the somatic growth of the endangered pallid sturgeon (Scaphirhynchus albus). For effects on fish populations of riverine ecosystems, cli- mate output simulated by coarse-resolution AOGCMs and RCMs must be downscaled to basins to river hydrology to population response. One needs to transfer the information from these climate simulations down to the individual scale in a way that minimizes extrapolation and can account for spatio-temporal variability in the intervening stages. The goal is a framework to determine whether, given uncertainties in the climate models and the biological response, meaningful inference can still be made. The non-linear downscaling of climate information to the river scale requires that one realistically account for spatial and temporal variability across scale. Our down- scaling procedure includes the use of fixed/calibrated hydrological flow and temperature models coupled with a stochastically parameterized sturgeon bioenergetics model. We show that, although there is a large amount of uncertainty associated with both the climate model output and the fish growth process, one can establish significant differences in fish growth distributions between models, and between future and current climates for a given model.

Missouri River basin↗
USGS1984

Upwarp of anomalous asthenosphere beneath the Rio Grande rift

Continental rifts are possible analogues of mid-ocean ridges, although major plate tectonic features are less clearly observed1. Current thermal models of mid-ocean ridges2-4 consist of solid lithospheric plates overlying the hotter, less viscous asthenosphere, with plate thickness increasing away from the ridge axis. The lithospheric lower boundary lies at or near the melting point isotherm, so that at greater depths higher temperatures account for lower viscosity, lower seismic velocities and possibly partial melting. Upwarp of this boundary at the ridge axis concentrates heat there, thus lowering densities by expansion and raising the sea floor to the level of thermal isostatic equilibrium. At slow spreading ridges, a major central graben forms owing to the mechanics of magma injection into the crust5. Topography, heat flow, gravity and seismic studies support these models. On the continents, a low-velocity channel has been observed, although it is poorly developed beneath ancient cratons6-9. Plate tectonic models have been applied to continental basins and margins10-12, but further similarities to the oceanic models remain elusive. Topographic uplift is often ascribed to Airy type isostatic compensation caused by crustal thickening, rather than thermal compensation in the asthenosphere. Here we discuss the Rio Grande rift, in southwestern United States. Teleseismic P-wave residuals show that regional uplift is explained by asthenosphere uplift rather than crustal thickening. ?? 1984 Nature Publishing Group.

Nature↗
USGS2007

Implications of the 26 December 2004 Sumatra-Andaman earthquake on tsunami forecast and assessment models for great subduction-zone earthquakes

Results from different tsunami forecasting and hazard assessment models are compared with observed tsunami wave heights from the 26 December 2004 Indian Ocean tsunami. Forecast models are based on initial earthquake information and are used to estimate tsunami wave heights during propagation. An empirical forecast relationship based only on seismic moment provides a close estimate to the observed mean regional and maximum local tsunami runup heights for the 2004 Indian Ocean tsunami but underestimates mean regional tsunami heights at azimuths in line with the tsunami beaming pattern (e.g., Sri Lanka, Thailand). Standard forecast models developed from subfault discretization of earthquake rupture, in which deep- ocean sea level observations are used to constrain slip, are also tested. Forecast models of this type use tsunami time-series measurements at points in the deep ocean. As a proxy for the 2004 Indian Ocean tsunami, a transect of deep-ocean tsunami amplitudes recorded by satellite altimetry is used to constrain slip along four subfaults of the M >9 Sumatra–Andaman earthquake. This proxy model performs well in comparison to observed tsunami wave heights, travel times, and inundation patterns at Banda Aceh. Hypothetical tsunami hazard assessments models based on end- member estimates for average slip and rupture length (M w 9.0–9.3) are compared with tsunami observations. Using average slip (low end member) and rupture length (high end member) (M w 9.14) consistent with many seismic, geodetic, and tsunami inversions adequately estimates tsunami runup in most regions, except the extreme runup in the western Aceh province. The high slip that occurred in the southern part of the rupture zone linked to runup in this location is a larger fluctuation than expected from standard stochastic slip models. In addition, excess moment release (∼9%) deduced from geodetic studies in comparison to seismic moment estimates may generate additional tsunami energy, if the exponential time constant of slip is less than approximately 1 hr. Overall, there is significant variation in assessed runup heights caused by quantifiable uncertainty in both first-order source parameters (e.g., rupture length, slip-length scaling) and spatiotemporal complexity of earthquake rupture.

Bulletin of the Seismological Society of America↗
USGS2011

Pliocene Model Intercomparison Project (PlioMIP): experimental design and boundary conditions (Experiment 2)

The Palaeoclimate Modelling Intercomparison Project has expanded to include a model intercomparison for the mid-Pliocene warm period (3.29 to 2.97 million yr ago). This project is referred to as PlioMIP (the Pliocene Model Intercomparison Project). Two experiments have been agreed upon and together compose the initial phase of PlioMIP. The first (Experiment 1) is being performed with atmosphere-only climate models. The second (Experiment 2) utilises fully coupled ocean-atmosphere climate models. Following on from the publication of the experimental design and boundary conditions for Experiment 1 in Geoscientific Model Development, this paper provides the necessary description of differences and/or additions to the experimental design for Experiment 2.

Geoscientific Model Development↗
USGS2018

Thermal erosion of a permafrost coastline: Improving process-based models using time-lapse photography

Coastal erosion rates locally exceeding 30 m y-1 have been documented along Alaska's Beaufort Sea coastline, and a number of studies suggest that these erosion rates have accelerated as a result of climate change. However, a lack of direct observational evidence has limited our progress in quantifying the specific processes that connect climate change to coastal erosion rates in the Arctic. In particular, while longer ice-free periods are likely to lead to both warmer surface waters and longer fetch, the relative roles of thermal and mechanical (wave) erosion in driving coastal retreat have not been comprehensively quantified. We focus on a permafrost coastline in the northern National Petroleum Reserve-Alaska (NPR-A), where coastal erosion rates have averaged 10-15 m y-1 over two years of direct monitoring. We take advantage of these extraordinary rates of coastal erosion to observe and quantify coastal erosion directly via time-lapse photography in combination with meteorological observations. Our observations indicate that the erosion of these bluffs is largely thermally driven, but that surface winds play a crucial role in exposing the frozen bluffs to the radiatively warmed seawater that drives melting of interstitial ice. To first order, erosion in this setting can be modeled using formulations developed to describe iceberg deterioration in the open ocean. These simple models provide a conceptual framework for evaluating how climate-induced changes in thermal and wave energy might influence future erosion rates in this setting.

Arctic, Antarctic, and Alpine Research↗
USGS2012

Sea ice melting in the marginal ice zone

The heat and salt flux boundary conditions together with the freezing curve relationship are a necessary component of any ice-sea water thermodynamic model. A neutral two-layer oceanic planetary boundary layer model that incorporates these boundary conditions gives the following results: The interfacial salinity is within 10% of the far-field salinity for conditions commonly encountered in the MIZ and depends only on the turbulent Lewis number and the far-field temperature and salinity. The predicted melt rates agree with the limited field observations, of the order of 1 m day −1 . The Obukov lengths, determined from the predicted interface conditions and melt rates, are generally much greater than the Ekman layer thicknesses; hence, the surface buoyancy flux has little effect on the turbulence in the planetary boundary layer.

Journal of Geophysical Research - Oceans↗
USGS2009

Projected climate-induced faunal change in the Western Hemisphere

Climate change is predicted to be one of the greatest drivers of ecological change in the coming century. Increases in temperature over the last century have clearly been linked to shifts in species distributions. Given the magnitude of projected future climatic changes, we can expect even larger range shifts in the coming century. These changes will, in turn, alter ecological communities and the functioning of ecosystems. Despite the seriousness of predicted climate change, the uncertainty in climate-change projections makes it difficult for conservation managers and planners to proactively respond to climate stresses. To address one aspect of this uncertainty, we identified predictions of faunal change for which a high level of consensus was exhibited by different climate models. Specifically, we assessed the potential effects of 30 coupled atmosphere-ocean general circulation model (AOGCM) future-climate simulations on the geographic ranges of 2954 species of birds, mammals, and amphibians in the Western Hemisphere. Eighty percent of the climate projections based on a relatively low greenhouse-gas emissions scenario result in the local loss of at least 10% of the vertebrate fauna over much of North and South America. The largest changes in fauna are predicted for the tundra, Central America, and the Andes Mountains where, assuming no dispersal constraints, specific areas are likely to experience over 90% turnover, so that faunal distributions in the future will bear little resemblance to those of today. ?? 2009 by the Ecological Society of America.

Ecology↗
USGS1988

Descriptive and grade-tonnage models of volcanogenic manganese deposits in oceanic environments; a modification

Four types of volcanogenic manganese deposits, distinguished on the basis of geologic, geochemical, and geophysical characteristics, appear to result from a combination of volcanic and hydrothermal processes related to hot-spring activity in oceanic environments. We compare these four desposit types, here called the Franciscan, Cuban, Olympic Peninsula, and Cyprus, with respect to host rocks, associated rocks, minerals, deposit shape, dimensions, volume, tonnage, grade, and mineral-deposit density (number of deposits per unit area). Franciscan-type deposits occur in obducted oceanic ridge and backarc marginal-basin environments, are associated with chert, shale, and graywacke aroun the margins of mafic volcanic centers, and have a median tonnage of 450 t and median grades of 36 weight percent Mn and less than 5.1 weight percent Fe. Cuban-type deposits occur in island-arc environments, are associated with tuff and limestone around domal structures or intrusions inferred to be volcanic centers, and have a median tonnage of 6,400 t and median grades of 39 weight percent Mn and less than 4.4 weight percent Fe. Olympic Peninsula-type deposits occur in obducted oceanic midplate settings, are associated with argillaceous limestone, argillite, and graywacke around mafic volcanic centers (seamounts or islands), and have a median tonnage of 340 t and median grades of 35 weight percent Mn and less than 6.5 weight percent Fe. Cyprus-type deposits occur in the same tectonic environments as Franciscan type but are associated with basalt, marl, chalk, silt, and chert off the ridge-axis position and have a median tonnage of 41,000 t and median grades of 33 weight percent Fe and 8 weight percent Mn. All these deposits are thin ellipsoids, concordant to the host rocks, but Cyprus-and Cuban-type deposits are larger than Franciscan- and Olympic Peninsula-type deposits. Except for Cyprus-type deposits, which are manganiferous iron (umber) deposits composed of hydrated iron and manganese oxides, all volcanogenic manganese deposits contain manganese oxides, silicates, and carbonates. Mineral-deposit densities, along with grade and tonnage information, are useful for estimating the number, size, andgrades of these deposits in resource assessments.

Bulletin↗
USGS2009

Investigating the impact of surface wave breaking on modeling the trajectories of drifters in the northern Adriatic Sea during a wind-storm event

An accurate numerical prediction of the oceanic upper layer velocity is a demanding requirement for many applications at sea and is a function of several near-surface processes that need to be incorporated in a numerical model. Among them, we assess the effects of vertical resolution, different vertical mixing parameterization (the so-called Generic Length Scale -GLS- set of k-??, k-??, gen, and the Mellor-Yamada), and surface roughness values on turbulent kinetic energy (k) injection from breaking waves. First, we modified the GLS turbulence closure formulation in the Regional Ocean Modeling System (ROMS) to incorporate the surface flux of turbulent kinetic energy due to wave breaking. Then, we applied the model to idealized test cases, exploring the sensitivity to the above mentioned factors. Last, the model was applied to a realistic situation in the Adriatic Sea driven by numerical meteorological forcings and river discharges. In this case, numerical drifters were released during an intense episode of Bora winds that occurred in mid-February 2003, and their trajectories compared to the displacement of satellite-tracked drifters deployed during the ADRIA02-03 sea-truth campaign. Results indicted that the inclusion of the wave breaking process helps improve the accuracy of the numerical simulations, subject to an increase in the typical value of the surface roughness z 0 . Specifically, the best performance was obtained using ?? CH = 56,000 in the Charnok formula, the wave breaking parameterization activated, k-?? as the turbulence closure model. With these options, the relative error with respect to the average distance of the drifter was about 25% (5.5 km/day). The most sensitive factors in the model were found to be the value of ?? CH enhanced with respect to a standard value, followed by the adoption of wave breaking parameterization and the particular turbulence closure model selected. ?? 2009 Elsevier Ltd.

Ocean Modelling↗
USGS2021

Digital Twin Earth - Coasts: Developing a fast and physics-informed surrogate model for coastal floods via neural operators

Developing fast and accurate surrogates for physics-based coastal and ocean mod- els is an urgent need due to the coastal flood risk under accelerating sea level rise, and the computational expense of deterministic numerical models. For this purpose, we develop the first digital twin of Earth coastlines with new physics-informed machine learning techniques extending the state-of-art Neural Operator. As a proof-of-concept study, we built Fourier Neural Operator (FNO) surrogates on the simulations of an industry-standard coastal and ocean model – Nucleus for Euro- pean Modelling of the Ocean (NEMO). The resulting FNO surrogate accurately predicts the sea surface height in most regions while achieving upwards of 45x acceleration of NEMO. We delivered an open-source CoastalTwin platform in an end-to-end and modular way, to enable easy extensions to other simulations and ML-based surrogate methods. Our results and deliverable provide a promising approach to massively accelerate coastal dynamics simulators, which can enable scientists to efficiently execute many simulations for decision-making, uncertainty quantification, and other research activities.

Conference Paper↗
USGS2011

Sea surface temperatures of the mid-Piacenzian Warm Period: A comparison of PRISM3 and HadCM3

It is essential to document how well the current generation of climate models performs in simulating past climates to have confidence in their ability to project future conditions. We present the first global, in-depth comparison of Pliocene sea surface temperature (SST) estimates from a coupled ocean–atmosphere climate model experiment and a SST reconstruction based on proxy data. This enables the identification of areas in which both the climate model and the proxy dataset require improvement. In general, the fit between model-produced SST anomalies and those formed from the available data is very good. We focus our discussion on three regions where the data–model anomaly exceeds 2 °C. 1) In the high latitude North Pacific, a systematic model error may result in anomalies that are too cold. Also, the deeper Pliocene thermocline may cause disagreement along the California margin; either the upwelling in the model is too strong or the modeled thermocline is not deep enough. 2) In the North Atlantic, the model predicts cooling in the center of a data-based warming trend that steadily increases with latitude from + 1.5 °C to >+ 6 °C. The discrepancy may arise because the modeled North Atlantic Current is too zonal compared to reality, which is reinforced by the lowering of the altitude of the Pliocene Western Cordillera Mountains. In addition, the model's use of modern bathymetry in the higher latitudes may have led the model to underestimate the northward penetration of warmer surface water into the Arctic. 3) Finally, though the data and model show good general agreement across most of the Southern Ocean, a few locations show offsets due to the modern land–sea mask used in the model. Additional considerations could account for many of the modest data–model anomalies, such as differences between calibration climatologies, the oversimplification of the seasonal cycle, and differences between SST proxies (i.e. seasonality and water depth). New SST estimates from data-sparse and regionally important areas will greatly enhance our ability to judge model performance.

Palaeogeography, Palaeoclimatology, Palaeoecology↗
USGS2007

Simulation and observations of annual density banding in skeletons of Montastraea (Cnidaria: Scleractinia) growing under thermal stress associated with ocean warming

We present a model of annual density banding in skeletons of Montastraea coral species growing under thermal stress associated with an ocean-warming scenario. The model predicts that at sea-surface temperatures (SSTs) <29??C, high-density bands (HDBs) are formed during the warmest months of the year. As temperature rises and oscillates around the optimal calcification temperature, an annual doublet in the HDB (dHDB) occurs that consists of two narrow HDBs. The presence of such dHDBs in skeletons of Montastraea species is a clear indication of thermal stress. When all monthly SSTs exceed the optimal calcification temperature, HDBs form during the coldest, not the warmest, months of the year. In addition, a decline in mean-annual calcification rate also occurs during this period of elevated SST. A comparison of our model results with annual density patterns observed in skeletons of M. faveolata and M. franksi, collected from several localities in the Mexican Caribbean, indicates that elevated SSTs are already resulting in the presence of dHDBs as a first sign of thermal stress, which occurs even without coral bleaching. ?? 2007, by the American Society of Limnology and Oceanography, Inc.

Limnology and Oceanography↗
USGS2001

Predictions of hydrothermal alteration within near-ridge oceanic crust from coordinated geochemical and fluid flow models

Coordinated geochemical and hydrological calculations guide our understanding of the composition, fluid flow patterns, and thermal structure of near-ridge oceanic crust. The case study presented here illustrates geochemical and thermal changes taking place as oceanic crust ages from 0.2 to 1.0 Myr. Using a finite element code, we model fluid flow and heat transport through the upper few hundred meters of an abyssal hill created at an intermediate spreading rate. We use a reaction path model with a customized database to calculate equilibrium fluid compositions and mineral assemblages of basalt and seawater at 500 bars and temperatures ranging from 150 to 400??C. In one scenario, reaction path calculations suggest that volume increases on the order of 10% may occur within portions of the basaltic basement. If this change in volume occurred, it would be sufficient to fill all primary porosity in some locations, effectively sealing off portions of the oceanic crust. Thermal profiles resulting from fluid flow simulations indicate that volume changes along this possible reaction path occur primarily within the first 0.4 Myr of crustal aging. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2015

Projected wave conditions in the Eastern North Pacific under the influence of two CMIP5 climate scenarios

Hindcast and 21st century winds, simulated by General Circulation Models (GCMs), were used to drive global- and regional-scale spectral wind-wave generation models in the Pacific Ocean Basin to assess future wave conditions along the margins of the North American west coast and Hawaiian Islands. Three-hourly winds simulated by four separate GCMs were used to generate an ensemble of wave conditions for a recent historical time-period (1976&ndash;2005) and projections for the mid and latter parts of the 21st century under two radiative forcing scenarios (RCP 4.5 and RCP 8.5), as defined by the fifth phase of the Coupled Model Inter-comparison Project (CMIP5) experiments. Comparisons of results from historical simulations with wave buoy and ERA-Interim wave reanalysis data indicate acceptable model performance of wave heights, periods, and directions, giving credence to generating projections. Mean and extreme wave heights are projected to decrease along much of the North American west coast. Extreme wave heights are projected to decrease south of &sim;50&deg;N and increase to the north, whereas extreme wave periods are projected to mostly increase. Incident wave directions associated with extreme wave heights are projected to rotate clockwise at the eastern end of the Aleutian Islands and counterclockwise offshore of Southern California. Local spatial patterns of the changing wave climate are similar under the RCP 4.5 and RCP 8.5 scenarios, but stronger magnitudes of change are projected under RCP 8.5. Findings of this study are similar to previous work using CMIP3 GCMs that indicates decreasing mean and extreme wave conditions in the Eastern North Pacific, but differ from other studies with respect to magnitude and local patterns of change. This study contributes toward a larger ensemble of global and regional climate projections needed to better assess uncertainty of potential future wave climate change, and provides model boundary conditions for assessing the impacts of climate change on coastal systems.

Ocean Modelling↗
USGS1987

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering↗