Search USGS⌕ Search

Geology topics

J.C. Warner

Publications and source records attributed to J.C. Warner.

11 recordsLinked to original sources
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↗
USGS2009

Coastal change along the shore of northeastern South Carolina— The South Carolina Coastal Erosion Study

The U.S. Geological Survey, in cooperation with the South Carolina Sea Grant Consortium, conducted a 7-year, multidisciplinary study of coastal erosion in northeastern South Carolina. Shoreline behavior along the coast of Long Bay is dictated by waves, tidal currents, and sediment supply that act within the overall constraints of the regional geologic setting. Beaches are thin ribbons of sand that sit on top of layered sedimentary rocks, which have been deeply eroded by rivers and coastal processes over millions of years. Offshore of the beaches, these sedimentary rocks are exposed as hardgrounds over large expanses of shallow seafloor and are locally overlain by a discontinuous veneer of sandy sediment generally less than 1 m thick. Rates of shoreline retreat largely depend on the geologic framework of the shoreface that is being excavated by ocean processes. Mainland-attached beaches have remained relatively stable, whereas barrier islands have experienced large shifts in shoreline position. In this sediment-limited region, erosion of the shoreface and inner shelf probably contributes a significant amount of new material to the beach system. Oceanographic studies and numerical modeling show that sediment transport varies along the coast, depending on the type and travel path of storms that impact Long Bay, but the long-term net transport direction is generally from north to south. Changes in storm activity that might accompany climate change, coupled with anticipated increases in sea-level rise, are expected to strongly affect low-lying, heavily developed areas of the coast.

South Carolina↗
USGS2009

High resolution near-bed observations in winter near Cape Hatteras, North Carolina

The U.S. Geological Survey (USGS) Coastal and Marine Science Center in Woods Hole, Massachusetts, is leading an effort to understand the regional sediment dynamics along the coastline of North and South Carolina. As part of the Carolinas Coastal Change Processes Project, a geologic framework study in June of 2008 by the Woods Hole Coastal and Marine Science Center's Sea Floor Mapping Group focused on the seaward limit of Diamond Shoals and provided high resolution bathymetric data, surficial sediment characteristics, and subsurface geologic stratigraphy. These data also provided unprecedented guidance to identify deployment locations for tripods and moorings to investigate the processes that control sediment transport at Diamond Shoals. Equipment was deployed at three sites from early January, 2009 through early May, 2009: north and south of the shoals at 15 m depth, and at the tip at 24 m depth. Many strong storm systems were recorded during that time period. Mounted on the tripods were instruments to measure surface waves, pressure, current velocity, bottom turbulence, suspended-sediment profiles, and sea-floor sand-ripple bedforms. Many instruments were designed and programmed to sample in high resolution in time and space, as fast as 8 Hz hourly bursts and as small as 6 cm bin sizes in near bottom profiles. A second tripod at the north site also held a visual camera system and sonar imaging system which document seafloor bedforms. The region is known for its dynamics, and one of the tripods tipped over towards the end of the experiment. A preliminary look at the data suggests the region is characterized by high energy. Raw data from a burst recorded at the south site on Mar. 26 th show instantaneous flow speed at 150 cm/s at 0.5 m above the seabed. This paper reports preliminary highlights of the observations, based on raw data, and lessons learned from a deployment of large tripod systems in such a dynamic location. ??2009 MTS.

Conference Paper↗
USGS2008

Development of a three-dimensional, regional, coupled wave, current, and sediment-transport model

We are developing a three-dimensional numerical model that implements algorithms for sediment transport and evolution of bottom morphology in the coastal-circulation model Regional Ocean Modeling System (ROMS v3.0), and provides a two-way link between ROMS and the wave model Simulating Waves in the Nearshore (SWAN) via the Model-Coupling Toolkit. The coupled model is applicable for fluvial, estuarine, shelf, and nearshore (surfzone) environments. Three-dimensional radiation-stress terms have been included in the momentum equations, along with effects of a surface wave roller model. The sediment-transport algorithms are implemented for an unlimited number of user-defined non-cohesive sediment classes. Each class has attributes of grain diameter, density, settling velocity, critical stress threshold for erosion, and erodibility constant. Suspended-sediment transport in the water column is computed with the same advection-diffusion algorithm used for all passive tracers and an additional algorithm for vertical settling that is not limited by the CFL criterion. Erosion and deposition are based on flux formulations. A multi-level bed framework tracks the distribution of every size class in each layer and stores bulk properties including layer thickness, porosity, and mass, allowing computation of bed morphology and stratigraphy. Also tracked are bed-surface properties including active-layer thickness, ripple geometry, and bed roughness. Bedload transport is calculated for mobile sediment classes in the top layer. Bottom-boundary layer submodels parameterize wave-current interactions that enhance bottom stresses and thereby facilitate sediment transport and increase bottom drag, creating a feedback to the circulation. The model is demonstrated in a series of simple test cases and a realistic application in Massachusetts Bay.

Computers & Geosciences↗
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↗
USGS2008

Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems.

Journal of Computational Physics↗
USGS2007

Sensitivity of estuarine turbidity maximum to settling velocity, tidal mixing, and sediment supply

Estuarine turbidity maximum, numerical modeling, settling velocity, stratification The spatial and temporal distribution of suspended material in an Estuarine Turbidity Maxima (ETM) is primarily controlled by particle settling velocity, tidal mixing, shear-stress thresholds for resuspension, and sediment supply. We vary these parameters in numerical experiments of an idealized two-dimensional (x-z) estuary to demonstrate their affects on the development and retention of particles in an ETM. Parameters varied are the settling velocity (0.01, 0.1, and 0.5 mm/s), tidal amplitude (0.4 m 12 hour tide and 0.3 to 0.6 m 14 day spring neap cycle), and sediment availability (spatial supply limited or unlimited; and temporal supply as a riverine pulse during spring vs. neap tide). Results identify that particles with a low settling velocity are advected out of the estuary and particles with a high settling velocity provide little material transport to an ETM. Particles with an intermediate settling velocity develop an ETM with the greatest amount of material retained. For an unlimited supply of sediment the ETM and limit of salt intrusion co-vary during the spring neap cycle. The ETM migrates landward of the salt intrusion during spring tides and seaward during neap tides. For limited sediment supply the ETM does not respond as an erodible pool of sediment that advects landward and seaward with the salt front. The ETM is maintained seaward of the salt intrusion and controlled by the locus of sediment convergence in the bed. For temporal variability of sediment supplied from a riverine pulse, the ETM traps more sediment if the pulse encounters the salt intrusion at neap tides than during spring tides. ?? 2007 Elsevier B.V. All rights reserved.

Proceedings in Marine Science↗
USGS2005

Numerical modeling of an estuary: A comprehensive skill assessment

Numerical simulations of the Hudson River estuary using a terrain-following, three-dimensional model (Regional Ocean Modeling System (ROMS)) are compared with an extensive set of time series and spatially resolved measurements over a 43 day period with large variations in tidal forcing and river discharge. The model is particularly effective at reproducing the observed temporal variations in both the salinity and current structure, including tidal, spring neap, and river discharge-induced variability. Large observed variations in stratification between neap and spring tides are captured qualitatively and quantitatively by the model. The observed structure and variations of the longitudinal salinity gradient are also well reproduced. The most notable discrepancy between the model and the data is in the vertical salinity structure. While the surface-to-bottom salinity difference is well reproduced, the stratification in the model tends to extend all the way to the water surface, whereas the observations indicate a distinct pycnocline and a surface mixed layer. Because the southern boundary coindition is located near the mouth the estuary, the salinity within the domain is particularly sensitive to the specification of salinity at the boundary. A boundary condition for the horizontal salinity gradient, based on the local value of salinity, is developed to incorporate physical processes beyond the open boundary not resolved by the model. Model results are sensitive to the specification of the bottom roughness length and vertical stability functions, insofar as they influence the intensity of vertical mixing. The results only varied slightly between different turbulence closure methods of k-??, k-??, and k-kl. Copyright 2005 by the American Geophysical Union.

New York↗
USGS2005

Performance of four turbulence closure models implemented using a generic length scale method

A two-equation turbulence model (one equation for turbulence kinetic energy and a second for a generic turbulence length-scale quantity) proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235] is implemented in a three-dimensional oceanographic model (Regional Oceanographic Modeling System; ROMS v2.0). These two equations, along with several stability functions, can represent many popular turbulence closures, including the k-kl (Mellor-Yamada Level 2.5), k-??, and k-?? schemes. The implementation adds flexibility to the model by providing an unprecedented range of turbulence closure selections in a single 3D oceanographic model and allows comparison and evaluation of turbulence models in an otherwise identical numerical environment. This also allows evaluation of the effect of turbulence models on other processes such as suspended-sediment distribution or ecological processes. Performance of the turbulence models and sediment-transport schemes is investigated with three test cases for (1) steady barotropic flow in a rectangular channel, (2) wind-induced surface mixed-layer deepening in a stratified fluid, and (3) oscillatory stratified pressure-gradient driven flow (estuarine circulation) in a rectangular channel. Results from k-??, k-??, and gen (a new closure proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235]) are very similar for these cases, but the k-kl closure results depend on a wall-proximity function that must be chosen to suit the flow. Greater variations appear in simulations of suspended-sediment concentrations than in salinity simulations because the transport of suspended-sediment amplifies minor variations in the methods. The amplification is caused by the added physics of a vertical settling rate, bottom stress dependent resuspension, and diffusive transport of sediment in regions of well mixed salt and temperature. Despite the amplified sensitivity of sediment to turbulence models in the estuary test case, the four closures investigated here all generated estuarine turbidity maxima that were similar in their shape, location, and concentrations.

Ocean Modelling↗