USGS Science⌕ Search

SEARCH · USGS Science

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 37 records · Page 2Linked to original sources
USGS2008

Storm-driven sediment transport in Massachusetts Bay

Massachusetts Bay is a semi-enclosed embayment in the western Gulf of Maine about 50 km wide and 100 km long. Bottom sediment resuspension is controlled predominately by storm-induced surface waves and transport by the tidal- and wind-driven circulation. Because the Bay is open to the northeast, winds from the northeast ('Northeasters') generate the largest surface waves and are thus the most effective in resuspending sediments. The three-dimensional oceanographic circulation model Regional Ocean Modeling System (ROMS) is used to explore the resuspension, transport, and deposition of sediment caused by Northeasters. The model transports multiple sediment classes and tracks the evolution of a multilevel sediment bed. The surficial sediment characteristics of the bed are coupled to one of several bottom-boundary layer modules that calculate enhanced bottom roughness due to wave-current interaction. The wave field is calculated from the model Simulating WAves Nearshore (SWAN). Two idealized simulations were carried out to explore the effects of Northeasters on the transport and fate of sediments. In one simulation, an initially spatially uniform bed of mixed sediments exposed to a series of Northeasters evolved to a pattern similar to the existing surficial sediment distribution. A second set of simulations explored sediment-transport pathways caused by storms with winds from the northeast quadrant by simulating release of sediment at selected locations. Storms with winds from the north cause transport southward along the western shore of Massachusetts Bay, while storms with winds from the east and southeast drive northerly nearshore flow. The simulations show that Northeasters can effectively transport sediments from Boston Harbor and the area offshore of the harbor to the southeast into Cape Cod Bay and offshore into Stellwagen Basin. This transport pattern is consistent with Boston Harbor as the source of silver found in the surficial sediments of Cape Cod Bay and Stellwagen Basin.

Massachusetts↗
USGS2014

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↗
USGS2019

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↗
USGS2011

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↗
USGS2007

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↗
USGS1996

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↗
USGS2007

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↗
USGS2012

The structural and geochemical evolution of the continental crust: Support for the oceanic plateau model of continental growth

The problem of the origin of the continental crust can be resolved into two fundamental questions: (1) the location and mechanisms of initial mantle extraction of the primitive crust and (2) the processes by which this primitive crust is converted into the continental crust that presently exists. We know that Archean continental crust is compositionally distinct from younger continental crust. Archean magmatism was dominantly bimodal, mafic thoeleiitic plus dacitic, heavy rare earth element depleted, in contrast to the dominantly unimodal, roughly andesitic calc‐alkaline magmatism on younger crust [ Taylor and McLennan , 1985; Condie , 1989]. The problem is whether these compositional differences are primarily due to different mechanisms of crustal extraction from the mantle or to different mechanisms of differentiation and alteration of newly formed continental crust.

Reviews of Geophysics↗
USGS2021

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↗
USGS2012

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↗
USGS2010

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↗
USGS2014

Modeling future scenarios of light attenuation and potential seagrass success in a eutrophic estuary

Estuarine eutrophication has led to numerous ecological changes, including loss of seagrass beds. One potential cause of these losses is a reduction in light availability due to increased attenuation by phytoplankton. Future sea level rise will also tend to reduce light penetration and modify seagrass habitat. In the present study, we integrate a spectral irradiance model into a biogeochemical model coupled to the Regional Ocean Model System (ROMS). It is linked to a bio-optical seagrass model to assess potential seagrass habitat in a eutrophic estuary under future nitrate loading and sea-level rise scenarios. The model was applied to West Falmouth Harbor, a shallow estuary located on Cape Cod (Massachusetts) where nitrate from groundwater has led to eutrophication and seagrass loss in landward portions of the estuary. Measurements of chlorophyll, turbidity, light attenuation, and seagrass coverage were used to assess the model accuracy. Mean chlorophyll based on uncalibrated in-situ fluorometry varied from 28 μg L −1 at the landward-most site to 6.5 μg L −1 at the seaward site, while light attenuation ranged from 0.86 to 0.45 m -1 . The model reproduced the spatial variability in chlorophyll and light attenuation with RMS errors of 3.72 μg L −1 and 0.07 m -1 respectively. Scenarios of future nitrate reduction and sea-level rise suggest an improvement in light climate in the landward basin with a 75% reduction in nitrate loading. This coupled model may be useful to assess habitat availability changes due to eutrophication and sediment resuspension and fully considers spatial variability on the tidal timescale.

Massachusetts↗
USGS2018

Chapter 5: Impact of sedimentary-exhalative hydrothermal systems on marine chemistry and mass extinctions: Applications for ore genesis research and mineral exploration

Times of metal-rich brine discharge into ancient ocean basins, associated with the formation of sedimentary-exhalative (sedex) Zn-Pb-Ba ore deposits, coincided with short-duration positive excursions (“spikes”) in the global marine Sr isotope record. While these spikes are unexplained by conventional oceanic models, chronostratigraphic correlations, combined with mass balance evidence and oceanographic modeling, suggest that the flux of radiogenic Sr from sedex brines during ore formation is sufficient to explain these previously enigmatic 87 Sr/ 86 Sr spikes. We review existing 87 Sr/ 86 Sr data and present new data as verification of these global 87 Sr/ 86 Sr spikes and their correlations with the formation of giant sedex ore deposits. Major events include an 1 ×10 −4 (~0.7078–~0.7079) excursion contemporaneous with formation of the Rammelsberg deposit at ~389 Ma; spikes on the order of 1 to 3 × 10 −4 , coeval with formation of the Meggen deposit at ~381 Ma, several ore deposits in the Macmillan Pass district at ~379 to 375 Ma, and the Silvermines deposits at ~352 Ma; and two >6 × 10 −4 spikes coincident with formation of the giant Navan deposit at ~346 Ma and Red Dog deposits at ~337 Ma. Moreover, the timing of peak 87 Sr/ 86 Sr spikes correlates with global δ 13 C and δ 18 O spikes, deposition of metal-rich black shales and ironstones, metal-induced malformation (teratology) of marine organisms, and mass extinctions. The relationships among these features were poorly understood, but our new model explains how the flux of key biolimiting nutrients and metals contained in sedex brines, demonstrably equivalent to or exceeding that of the total modern riverine flux to the ocean, spurred ocean eutrophication, which, ultimately, through a series of positive feedback mechanisms, may have triggered global chemical and biological events. If, as we hypothesize, sedex hydrothermal systems are recorded in the global marine isotopic, geologic, and biological records, our findings define a new approach to the study of and exploration for sedex deposits. We demonstrate that fluid inclusion solute chemistry and isotopic and stratigraphic studies of sedex deposits, 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 provides evidence for the age, duration, and fluxes of fluids and metals vented into the ocean by these giant hydrothermal systems. Accordingly, the marine 87 Sr/ 86 Sr curve constitutes a global exploration tool that could be applied to assess the mineral potential of sedimentary basins. To illustrate the potential of this tool to identify favorable stratigraphic ages and basins with potential for undiscovered giant sedex deposits, we highlight several spikes, on par with those characteristic of the Red Dog and Navan deposits, which have not been correlated with known metal deposits. Given these strong temporal correlations, mass balance estimates, and results of ocean chemistry modeling, our study suggests that further work is warranted to determine the extent to which periodic venting of hydrothermal basinal brines into the ocean has influenced the evolution of marine chemistry. Ultimately, these global signatures can be applied to the study of and exploration for sedex deposits.

Book chapter↗
USGS2022

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↗
USGS2017

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↗
USGS2002

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↗
USGS1984

Numerical noise in ocean and estuarine models

Approximate methods for solving the shallow water equations may lead to solutions exhibiting large fictitious, numerically-induced oscillations. The analysis of the discrete dispersion relation and modal solutions of small wavelengths provides a powerful technique for assessing the sensitivity of alternative numerical schemes to irregular data which may lead to such oscillatory numerical noise. For those schemes where phase speed vanishes at a finite wavenumber or there are multiple roots for wavenumber, oscillation modes can exist which are uncoupled from the dynamics of the problem. The discrete modal analysis approach is used here to identify two classes of spurious oscillation modes associated respectively with the two different asymptotic limits corresponding to estuarine and large scale ocean models. The analysis provides further insight into recent numerical results for models which include large spatial scales and Coriolis acceleration.

Advances in Water Resources↗