USGS ScienceSearch

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.

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

At least 37 records · Page 2Linked to original sources

Air-sea interactions during strong winter extratropical storms

A high-resolution, regional coupled atmosphere–ocean model is used to investigate strong air–sea interactions during a rapidly developing extratropical cyclone (ETC) off the east coast of the USA. In this two-way coupled system, surface momentum and heat fluxes derived from the Weather Research and Forecasting model and sea surface temperature (SST) from the Regional Ocean Modeling System are exchanged via the Model Coupling Toolkit. Comparisons are made between the modeled and observed wind velocity, sea level pressure, 10 m air temperature, and sea surface temperature time series, as well as a comparison between the model and one glider transect. Vertical profiles of modeled air temperature and winds in the marine atmospheric boundary layer and temperature variations in the upper ocean during a 3-day storm period are examined at various cross-shelf transects along the eastern seaboard. It is found that the air–sea interactions near the Gulf Stream are important for generating and sustaining the ETC. In particular, locally enhanced winds over a warm sea (relative to the land temperature) induce large surface heat fluxes which cool the upper ocean by up to 2 °C, mainly during the cold air outbreak period after the storm passage. Detailed heat budget analyses show the ocean-to-atmosphere heat flux dominates the upper ocean heat content variations. Results clearly show that dynamic air–sea interactions affecting momentum and buoyancy flux exchanges in ETCs need to be resolved accurately in a coupled atmosphere–ocean modeling framework.

Ocean Dynamics

Investigating bedload transport under asymmetrical waves using a coupled ocean-wave model

Transport by asymmetrical wave motions plays a key role in cross-shore movement of sand, which is important for bar migration, exchange through tidal inlets, and beach recovery after storms. We have implemented a modified version of the SANTOSS formulation in the three-dimensional open-source Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling framework. The calculation of bedload transport requires inputs that include: water depth, bulk wave statistics (significant wave height, wave period, wave orbital velocity), current velocity at the edge of the wave boundary layer, and sediment density and grain size. While the coupled ocean-wave model computes water density, depth, and bulk wave statistics, we implement a method to calculate current velocity assuming a log profile and using the Madsen formulations for wave-current bottom boundary layer flows. We investigate the sensitivity of the calculation of near-bottom current velocity to model choices and its influence on cross-shore bedload transport. Results are compared to available numerical experiments using coupled fluid and discrete element model (CFDEM) simulations.

Conference Paper

Implementation and modification of a three-dimensional radiation stress formulation for surf zone and rip-current applications

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and appears to be more appropriate for sigma coordinates in very shallow waters where wave breaking conditions dominate. The improvements of the modified model are shown through simulations of several cases that include: (a) obliquely incident spectral waves on a planar beach; (b) obliquely incident spectral waves on a natural barred beach (DUCK'94 experiment); (c) alongshore variable offshore wave forcing on a planar beach; (d) alongshore varying bathymetry with constant offshore wave forcing; and (e) nearshore barred morphology with rip-channels. Quantitative and qualitative comparisons to previous analytical, numerical, laboratory studies and field measurements show that the modified model replicates surf zone recirculation patterns (onshore drift at the surface and undertow at the bottom) more accurately than previous formulations based on radiation stress (Haas and Warner, 2009). The results of the model and test cases are further explored for identifying the forces operating in rip current development and the potential implication for sediment transport and rip channel development. Also, model analysis showed that rip current strength is higher when waves approach at angles of 5° to 10° in comparison to normally incident waves.

Coastal Engineering

Coupling alongshore variations in wave energy to beach morphologic change using the SWAN wave model at Ocean Beach, San Francisco, CA

Coastal managers have faced increasing pressure to manage their resources wisely over the last century as a result of heightened development and changing environmental forcing. It is crucial to understand seasonal changes in beach volume and shape in order to identify areas vulnerable to accelerated erosion. Shepard (1950) was among the first to quantify seasonal beach cycles. Sonu and Van Beek (1971) and Wright et al. (1985) described commonly occurring beach states. Most studies utilize widest spaced 2-D cross shore profiles or shorelines extracted from aerial photographs (e.g. Winant et al. 1975; Aubrey, 1979, Aubrey and Ross, 1985; Larson and Kraus, 1994; Jimenez et al., 1977; Lacey and Peck, 1998; Guillen et al., 1999; Norcorss et al., 2002) to analyzed systematic changes in beach evolution. But with the exception of established field stations, such as Duck, NC (Birkemeier and Mason, 1984), ans Hazaki Oceanographical Research Station (HORS) in Japan (Katoh, 1997), there are very few beach change data sets with high temporal and spatial resolutions (e.g. Dail et al., 2000; Ruggiero et al., 2005; Yates et al., in press). Comprehensive sets of nearshore morphological data and local in situ measurements outside of these field stations are very rare and virtually non-existent high-energy coasts. Studied that have attempted to relate wave statistics to beach morphology change require some knowledge of the nearshore wave climate, and have had limited success using offshore measurement (Sonu and Van Beek, 1971; Dail et al., 2000). The primary objective of this study is to qualitatively compare spatially variable nearshore wave predictions to beach change measurements in order to understand the processes responsible for a persistent erosion 'hotspot' at Ocean Beach, San Francisco, CA. Local wave measurements are used to calibrate and validate a wave model that provides nearshore wave prediction along the beach. The model is run for thousands of binned offshore wave conditions to help isolate the effects of offshore wave direction and period on nearshore wave predictions. Alongshore varying average beach change statistics are computed at specific profile locations from topographic beach surveys and lidar data. The study area is located in the San Francisco Bight in central California. Ocean Beach is a seven kilometer long north-south trending sandy coastline located just south of the entrance to the San Francisco Bay Estuary (Figure 1). It contains an erosion hotspot in the southern part of the beach which has resulted in damage to local infrastructure and is the cause of continued concern. A wide range of field data collection and numerical modeling efforts have been focused here as part of the United States Geological Survey's (USGS) San Francisco Bight Coastal Processes Study, which began in October 2003 and represents the first comprehensive study of coastal processes at the mouth of San Francisco Bay. Ocean Beach is exposed to very strong tidal flows, with measured currents often in excess of 1 m/s at the north end of the beach. Current profiler measurements indicate that current magnitudes are greater in the northern portion of the beach, while wave energy is greater in the southern portion where erosion problems are greatest (Barnard et al., 2007). The sub-aerial beach volume fluctuates seasonally over a maximum envelope of 400,000 m 3 for the seven kilometer stretch (Barnard et al, 2007). The wave climate in the region is dominated by an abundance of low frequency energy (greater than 20 s period) and prevailing northwest incident wave angles. The application of a wave model to the region is further complicated by the presence of the Farallon Islands 40 kilometers west, and a massive ebb tidal delta at the mouth of San Francisco Bay (~150 km 2 ), which creates complicated refraction patterns as wave energy moves from offshore Ocean Beach; however the cost and threat of the energetic nearshore environment have limited the temporal and spatial resolution of these measurements. Applying numerical models to predict wave and current patterns along the beach can help supplement the filed data that exists and provide opportunities to make prediction about the impacts of changing environmental forcing.

California

Potential role of vegetation feedback in the climate sensitivity of high-latitude regions: A case study at 6000 years B.P.

Previous climate model simulations have shown that the configuration of the Earth's orbit during the early to mid-Holocene (approximately 10–5 kyr) can account for the generally warmer-than-present conditions experienced by the high latitudes of the northern hemisphere. New simulations for 6 kyr with two atmospheric/mixed-layer ocean models (Community Climate Model, version 1, CCMl, and Global ENvironmental and Ecological Simulation of Interactive Systems, version 2, GENESIS 2) are presented here and compared with results from two previous simulations with GENESIS 1 that were obtained with and without the albedo feedback due to climate-induced poleward expansion of the boreal forest. The climate model results are summarized in the form of potential vegetation maps obtained with the global BIOME model, which facilitates visual comparisons both among models and with pollen and plant macrofossil data recording shifts of the forest-tundra boundary. A preliminary synthesis shows that the forest limit was shifted 100–200 km north in most sectors. Both CCMl and GENESIS 2 produced a shift of this magnitude. GENESIS 1 however produced too small a shift, except when the boreal forest albedo feedback was included. The feedback in this case was estimated to have amplified forest expansion by approximately 50%. The forest limit changes also show meridional patterns (greatest expansion in central Siberia and little or none in Alaska and Labrador) which have yet to be reproduced by models. Further progress in understanding of the processes involved in the response of climate and vegetation to orbital forcing will require both the deployment of coupled atmosphere-biosphere-ocean models and the development of more comprehensive observational data sets.

Global Biogeochemical Cycles

Bora event variability and the role of air-sea feedback

A two-way interacting high resolution numerical simulation of the Adriatic Sea using the Navy Coastal Ocean Model (NCOM) and Coupled Ocean/ Atmosphere Mesoscale Prediction System (COAMPS??) was conducted to improve forecast momentum and heat flux fields, and to evaluate surface flux field differences for two consecutive bora events during February 2003. (COAMPS?? is a registered trademark of the Naval Research Laboratory.) The strength, mean positions and extensions of the bora jets, and the atmospheric conditions driving them varied considerably between the two events. Bora 1 had 62% stronger heat flux and 51% larger momentum flux than bora 2. The latter displayed much greater diurnal variability characterized by inertial oscillations and the early morning strengthening of a west Adriatic barrier jet, beneath which a stronger west Adriatic ocean current developed. Elsewhere, surface ocean current differences between the two events were directly related to differences in wind stress curl generated by the position and strength of the individual bora jets. The mean heat flux bias was reduced by 72%, and heat flux RMSE reduced by 30% on average at four instrumented over-water sites in the two-way coupled simulation relative to the uncoupled control. Largest reductions in wind stress were found in the bora jets, while the biggest reductions in heat flux were found along the north and west coasts of the Adriatic. In bora 2, SST gradients impacted the wind stress curl along the north and west coasts, and in bora 1 wind stress curl was sensitive to the Istrian front position and strength. The two-way coupled simulation produced diminished surface current speeds of ???12% over the northern Adriatic during both bora compared with a one-way coupled simulation. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans

Climate impacts on the Gulf of Maine ecosystem: A review of observed and expected changes in 2050 from rising temperatures

The Gulf of Maine has recently experienced its warmest 5-year period (2015–2020) in the instrumental record. This warming was associated with a decline in the signature subarctic zooplankton species, Calanus finmarchicus . The temperature changes have also led to impacts on commercial species such as Atlantic cod ( Gadus morhua ) and American lobster ( Homarus americanus ) and protected species including Atlantic puffins ( Fratercula arctica ) and northern right whales ( Eubalaena glacialis ). The recent period also saw a decline in Atlantic herring ( Clupea harengus ) recruitment and an increase in novel harmful algal species, although these have not been attributed to the recent warming. Here, we use an ensemble of numerical ocean models to characterize expected ocean conditions in the middle of this century. Under the high CO 2 emissions scenario (RCP8.5), the average temperature in the Gulf of Maine is expected to increase 1.1°C to 2.4°C relative to the 1976–2005 average. Surface salinity is expected to decrease, leading to enhanced water column stratification. These physical changes are likely to lead to additional declines in subarctic species including C. finmarchicus , American lobster, and Atlantic cod and an increase in temperate species. The ecosystem changes have already impacted human communities through altered delivery of ecosystem services derived from the marine environment. Continued warming is expected to lead to a loss of heritage, changes in culture, and the necessity for adaptation.

Maine, Massachusetts, New Brunswick, New Hampshire

GFDL's ESM2 global coupled climate-carbon Earth System Models. Part I: physical formulation and baseline simulation characteristics

We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory's previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.

Journal of Climate

Effect of roughness formulation on the performance of a coupled wave, hydrodynamic, and sediment transport model

A variety of algorithms are available for parameterizing the hydrodynamic bottom roughness associated with grain size, saltation, bedforms, and wave–current interaction in coastal ocean models. These parameterizations give rise to spatially and temporally variable bottom-drag coefficients that ostensibly provide better representations of physical processes than uniform and constant coefficients. However, few studies have been performed to determine whether improved representation of these variable bottom roughness components translates into measurable improvements in model skill. We test the hypothesis that improved representation of variable bottom roughness improves performance with respect to near-bed circulation, bottom stresses, or turbulence dissipation. The inner shelf south of Martha’s Vineyard, Massachusetts, is the site of sorted grain-size features which exhibit sharp alongshore variations in grain size and ripple geometry over gentle bathymetric relief; this area provides a suitable testing ground for roughness parameterizations. We first establish the skill of a nested regional model for currents, waves, stresses, and turbulent quantities using a uniform and constant roughness; we then gauge model skill with various parameterization of roughness, which account for the influence of the wave-boundary layer, grain size, saltation, and rippled bedforms. We find that commonly used representations of ripple-induced roughness, when combined with a wave–current interaction routine, do not significantly improve skill for circulation, and significantly decrease skill with respect to stresses and turbulence dissipation. Ripple orientation with respect to dominant currents and ripple shape may be responsible for complicating a straightforward estimate of the roughness contribution from ripples. In addition, sediment-induced stratification may be responsible for lower stresses than predicted by the wave–current interaction model.

Ocean Modelling

Predicting larval alewife transport in Lake Michigan using hydrodynamic and Lagrangian particle dispersion models

Several species of fish in large lakes and marine environments have a pelagic larval stage, and are subject to variable transport that can ultimately regulate survival and recruitment success. Alewife, Alosa pseudoharengus , are subject to transport by complex coastal currents during their pelagic larval stage (~ 30 d). We assessed backward-trajectory simulations, consisting of a Lagrangian particle dispersion model linked to the Finite Volume Community Ocean Model, to estimate likely hatch locations of aged larval alewife collected from locations on both the eastern and western sides of Lake Michigan during July 2015. We used four deployments of three satellite-tracked drifter buoys in coastal waters to assess model skill in estimating the origin of a drifter from its final location. We found that the trajectories of drifters varied greatly, depending on wind events and associated coastal transport processes, including upwelling/downwelling and coastal jet currents. In 2 of 12 cases, the backward trajectory simulations failed to predict the drifter origin, associated with transport of 170 km in a narrow coastal jet current. In the remaining 10 cases, the known drifter origin was within 3.5 km of the spatial patch of predicted possible origins for a scenario of horizontal diffusivity (188 m 2 s −1 ) consistent with the offshore model grid resolution. Modeled backward trajectories estimated that alewife originated from the same side of the lake where they were collected, within ~ 100 km of the collection site. Our paper demonstrates the utility of hydrodynamic models to estimate a region of origin for aged larval fish.

Lake Michigan

Towards a comprehensive water quality modeling of Barnegat Bay: Development of ROMS to WASP Coupler

The Regional Ocean Modeling System (ROMS) has been coupled with the Water Quality Analysis Simulation Program (WASP) to be used in a comprehensive analysis of water quality in Barnegat Bay, New Jersey. The coupler can spatially aggregate hydrodynamic information in ROMS cells into larger WASP segments. It can also be used to resample ROMS output at a finer temporal scale to meet WASP time-stepping requirements. The coupler aggregates flow components, temperature, and salinity in ROMS output for input to WASP via a hydrodynamic linkage file. The coupler was tested initially with idealized cases designed to verify the water mass balance and conservation of constituent mass using one-to-one and one-to-many connectivity options between segments. A realistic example from the Toms River embayment, a subdomain of Barnegat Bay, was used to demonstrate the functionality of the coupling. A WASP eutrophication model accounting for dissolved oxygen (DO), nitrogen, and constant phytoplankton concentrations was applied to explore the distribution and trends in DO and nitrogen in the embayment for the period of July–August 2012. Results of DO modeling indicate satisfactory agreement with measurements collected at in-bay stations and also indicate that this coupled approach, despite substantial differences in spatiotemporal discretization between the models, provides adequate predictive capabilities.

New Jersey

West Florida shelf circulation and temperature budget for the 1999 spring transition

Mid-latitude continental shelves undergo a spring transition as the net surface heat flux changes from cooling to warming. Using in situ data and a numerical circulation model we investigate the circulation and temperature budget on the West Florida Continental Shelf (WFS) for the spring transition of 1999. The model is a regional adaptation of the primitive equation, Princeton Ocean Model forced by NCEP reanalysis wind and heat flux fields and by river inflows. Based on agreements between the modeled and observed fields we use the model to draw inferences on how the surface momentum and heat fluxes affect the seasonal and synoptic scale variability. We account for a strong southeastward current at mid-shelf by the baroclinic response to combined wind and buoyancy forcing, and we show how this local forcing leads to annually occurring cold and low salinity tongues. Through term-by-term analyses of the temperature budget we describe the WFS temperature evolution in spring. Heat flux largely controls the seasonal transition, whereas ocean circulation largely controls the synoptic scale variability. These two processes, however, are closely linked. Bottom topography and coastline geometry are important in generating regions of convergence and divergence. Rivers contribute to the local hydrography and are important ecologically. Along with upwelling, river inflows facilitate frontal aggregation of nutrients and the spring formation of a high concentration chlorophyll plume near the shelf break (the so-called ‘Green River’) coinciding with the cold, low salinity tongues. These features originate by local, shelf-wide forcing; the Loop Current is not an essential ingredient.

Gulf of Mexico

The GFDL global atmosphere and land model AM4.0/LM4.0: 2. Model description, sensitivity studies, and tuning strategies

In Part 2 of this two‐part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part 1. Part 2 provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.

Journal of Advances in Modeling Earth Systems

Circulation and effluent dilution modeling in Massachusetts Bay : model implementation, verification and results

A three-dimensional hydrodynamic model was developed as part of a cooperative U.S. Geological Survey/Massachusetts Water Resources Authority program to study contaminated sediment accumulation and transport in Massachusetts Bay. This report details the development of the model and assesses how well the model represents observed currents and water properties in the bay. It also summarizes circulation and comparative effluent dilution simulations from existing and future Boston sewage outfalls over a three-year period from October 1, 1989 to December 31, 1992. The ECOM-si model, a semi-implicit version of the Blumberg and Mellor (1987) Estuarine, Coastal and Ocean Model, is shown to reproduce many of the important hydrodynamical features of Massachusetts Bay: the seasonal evolution of the pycnocline, the mean flow pattern, and the strength of sub-tidal current fluctuations. Throughout the simulation period, during both vertically well-mixed and stratified conditions, the seasonal statistics of observed currents are well-represented by the model. The model is therefore appropriate for studying the average dilution of sewage effluent and other continuously discharged substances over seasonal time scales. The ability of the model to reproduce individual flow events varies with season and location within the bay. Flow events during unstratified conditions in western Massachusetts Bay are particularly well-represented, indicating that the model is appropriate for studying processes such as the transport of suspended material from the future outfall site due to winter storms. Individual flow events during stratified conditions and in the offshore Stellwagen Bank region, however, are less well-represented due to small length scales (caused by upwelling and river discharge events) coupled with insufficient data to specify open boundary forcing from the Gulf of Maine. Thus while the model might be used to answer issues such as the frequency with which Gulf of Maine river plumes visit the new outfall site, attempting to predict whether a particular plume would visit the outfall site could be problematic. Comparative simulations of effluent discharged from the existing and future Boston outfalls show that the region of relatively high effluent concentrations (1 part effluent to 200 parts sea water) is significantly smaller with the future outfall and is limited to Western Massachusetts Bay during both unstratified and stratified seasons. The region of even higher concentration (1 part effluent to 50 parts sea water) that covers much of Boston Harbor with the existing outfall is non-existent in the future outfall simulation. Additional simulations of chlorination plant failure predict that the offshore location of the future outfall will lead to dramatically lower levels of pathogens at area beaches.

Open-File Report

Modeling extreme water levels in the Salish Sea: The importance of including remote sea level anomalies for application in hydrodynamic simulations

Extreme water-level recurrence estimates for a complex estuary using a high-resolution 2D model and a new method for estimating remotely generated sea level anomalies (SLAs) at the model boundary have been developed. The hydrodynamic model accurately resolves the dominant physical processes contributing to extreme water levels across the Washington State waters of the Salish Sea, including the relative contribution of remote SLA and other non-tidal residual processes that drive extreme water levels above the predicted tide. The model’s predictions have errors of less than 15 cm (<5% of 3–4 m tidal range) at eight tide gauge locations across the model domain. The influence of remote SLAs at the seaward boundary of the model was implemented using a multivariate regression of readily available and locally relevant wind, sea surface temperature, and pressure anomaly data, combined with El Niño Index data (R 2 = 0.76). The hydrodynamic model simulations using the remote SLA predictor compared well with simulations using the widely used data-assimilative global ocean model HYCOM SLA data (root mean square difference of 5.5 cm). Extreme water-level recurrence estimates with and without remote SLA show that remote forcing accounts for 50–60% of the total water level anomaly observed along Salish Sea shorelines. The resulting model simulations across decadal timescales provide estimates of extreme water level recurrence across the Salish Sea, capturing climate variability important to long-term coastal hazard planning. This approach has widespread applications for other complex estuarine systems.

Water

The role of model dynamics in ensemble Kalman filter performance for chaotic systems

The ensemble Kalman filter (EnKF) is susceptible to losing track of observations, or ‘diverging’, when applied to large chaotic systems such as atmospheric and ocean models. Past studies have demonstrated the adverse impact of sampling error during the filter’s update step. We examine how system dynamics affect EnKF performance, and whether the absence of certain dynamic features in the ensemble may lead to divergence. The EnKF is applied to a simple chaotic model, and ensembles are checked against singular vectors of the tangent linear model, corresponding to short-term growth and Lyapunov vectors, corresponding to long-term growth. Results show that the ensemble strongly aligns itself with the subspace spanned by unstable Lyapunov vectors. Furthermore, the filter avoids divergence only if the full linearized long-term unstable subspace is spanned. However, short-term dynamics also become important as nonlinearity in the system increases. Non-linear movement prevents errors in the long-term stable subspace from decaying indefinitely. If these errors then undergo linear intermittent growth, a small ensemble may fail to properly represent all important modes, causing filter divergence. A combination of long and short-term growth dynamics are thus critical to EnKF performance. These findings can help in developing practical robust filters based on model dynamics.

Tellus, Series A: Dynamic Meteorology and Oceanogr

Sedex hydrothermal systems triggered upheavals to marine chemistry and mass extinctions: Applications for ore genesis research and mineral exploration

New USGS research reveals that the discharge of metal-rich brine that formed sedex deposits on ancient seafloors had profound effects on global ocean chemistry and biologic evolution. For example, brine expulsion caused short-duration positive excursions (“spikes”) in the global marine Sr isotope record. While these spikes are unexplained by conventional oceanic models, our chronostratigraphic correlations of major sedex events, combined with mass balance evidence and oceanographic modeling, confirm that the flux of radiogenic Sr from sedex brines during ore formation was sufficient to cause these once enigmatic 87Sr/86Sr spikes. Recognition that the timing of peak 87 Sr/ 86 Sr spikes correlates exactly with global δ 13 C (and δ 18 O) spikes, climate change, deposition of metal-rich black shales and ironstones, metal-induced malformation (teratology) of marine organisms, and mass extinctions, establishes a causal relationship between sedex deposits and these dramatic events in earth history. The relationships among these features are not fully understood. However, our new model demonstrates that the flux of key biolimiting nutrients and metals contained in sedex brines exceeds that of the total modern riverine flux to the ocean. Undoubtably, these immense nutrient fluxes spurred ocean eutrophication, which, ultimately, through a series of positive feedback mechanisms, may be a previously unrecognized trigger of global ocean anoxic events (OAEs) that produced these chemical and biological perturbations. A derivative result from this integrative research is the recognition that OAEs resulted in the formation of “bathtub rim” deposits at redox boundaries along continental margins that concentrated various redox sensitive critical minerals. For example, we have identified midcontinent phosphorite deposits that contains heavy REE grades and tonnages that rival any REE deposit in the world. The recognition that sedex-forming fluid expulsion events are recorded in the global marine isotopic, geologic, and biological records, defines a new approach to the study of and exploration for sedex deposits. Traditional ore genesis research, coupled with chronostratigraphic correlation and high-resolution 87 Sr/ 86 Sr isotope chemostratigraphy can be used to answer long-standing questions about geologic processes responsible for formation of these extraordinary deposits. This approach allows us to constrain, for the first time, the age, duration, and fluxes of fluids and metals vented into the ocean by these giant hydrothermal systems. Accordingly, the fact that large mineralizing events are recorded in the marine sedimentary record opens the tantalizing prospect that we have the ability to conduct effective resource assessments and define prospective basins anywhere in the world. This innovative approach allows for identification of favorable stratigraphic ages and basins and remote evaluation of the size (and, thus, the mineral potential) of undiscovered mineral deposits. This methodology could be applied on regional basin-wide assessments, to evaluate sedimentary basin prospectivity, resource favorability of specific horizons therein, and to the evaluation of the potential of early-stage prospects.

Geological Society of Nevada Newsletter

Development of a submerged aquatic vegetation growth model in the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST v3.4) model

The coupled biophysical interactions between submerged aquatic vegetation (SAV), hydrodynamics (currents and waves), sediment dynamics, and nutrient cycling have long been of interest in estuarine environments. Recent observational studies have addressed feedbacks between SAV meadows and their role in modifying current velocity, sedimentation, and nutrient cycling. To represent these dynamic processes in a numerical model, the presence of SAV and its effect on hydrodynamics (currents and waves) and sediment dynamics was incorporated into the open-source Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) model. In this study, we extend the COAWST modeling framework to account for dynamic changes of SAV and associated epiphyte biomass. Modeled SAV biomass is represented as a function of temperature, light, and nutrient availability. The modeled SAV community exchanges nutrients, detritus, dissolved inorganic carbon, and dissolved oxygen with the water-column biogeochemistry model. The dynamic simulation of SAV biomass allows the plants to both respond to and cause changes in the water column and sediment bed properties, hydrodynamics, and sediment transport (i.e., a two-way feedback). We demonstrate the behavior of these modeled processes through application to an idealized domain and then apply the model to a eutrophic harbor where SAV dieback is a result of anthropogenic nitrate loading and eutrophication. These cases demonstrate an advance in the deterministic modeling of coupled biophysical processes and will further our understanding of future ecosystem change.

Geoscientific Model Development