USGS Science⌕ Search

Geology topics

J.C. Warner

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

16 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

Sediment dispersal in the northwestern Adriatic Sea

Sediment dispersal in the Adriatic Sea was evaluated using coupled three-dimensional circulation and sediment transport models, representing conditions from autumn 2002 through spring 2003. The calculations accounted for fluvial sources, resuspension by waves and currents, and suspended transport. Sediment fluxes peaked during southwestward Bora wind conditions that produced energetic waves and strengthened the Western Adriatic Coastal Current. Transport along the western Adriatic continental shelf was nearly always to the south, except during brief periods when northward Sirocco winds reduced the coastal current. Much of the modeled fluvial sediment deposition was near river mouths, such as the Po subaqueous delta. Nearly all Po sediment remained in the northern Adriatic. Material from rivers that drain the Apennine Mountains traveled farther before deposition than Po sediment, because it was modeled with a lower settling velocity. Fluvial sediment delivered to areas with high average bed shear stress was more highly dispersed than material delivered to more quiescent areas. Modeled depositional patterns were similar to observed patterns that have developed over longer timescales. Specifically, modeled Po sediment accumulation was thickest near the river mouth with a very thin deposit extending to the northeast, consistent with patterns of modern sediment texture in the northern Adriatic. Sediment resuspended from the bed and delivered by Apennine Rivers was preferentially deposited on the northern side of the Gargano Peninsula, in the location of thick Holocene accumulation. Deposition here was highest during Bora winds when convergences in current velocities and off-shelf flux enhanced delivery of material to the midshelf. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
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

Using the Model Coupling Toolkit to couple earth system models

Continued advances in computational resources are providing the opportunity to operate more sophisticated numerical models. Additionally, there is an increasing demand for multidisciplinary studies that include interactions between different physical processes. Therefore there is a strong desire to develop coupled modeling systems that utilize existing models and allow efficient data exchange and model control. The basic system would entail model "1" running on "M" processors and model "2" running on "N" processors, with efficient exchange of model fields at predetermined synchronization intervals. Here we demonstrate two coupled systems: the coupling of the ocean circulation model Regional Ocean Modeling System (ROMS) to the surface wave model Simulating WAves Nearshore (SWAN), and the coupling of ROMS to the atmospheric model Coupled Ocean Atmosphere Prediction System (COAMPS). Both coupled systems use the Model Coupling Toolkit (MCT) as a mechanism for operation control and inter-model distributed memory transfer of model variables. In this paper we describe requirements and other options for model coupling, explain the MCT library, ROMS, SWAN and COAMPS models, methods for grid decomposition and sparse matrix interpolation, and provide an example from each coupled system. Methods presented in this paper are clearly applicable for coupling of other types of models. ?? 2008 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software↗
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↗
USGS2004

Tidal oscillation of sediment between a river and a bay: A conceptual model

A conceptual model of fine sediment transport between a river and a bay is proposed, based on observations at two rivers feeding the same bay. The conceptual model consists of river, transitional, and bay regimes. Within the transitional regime, resuspension, advection, and deposition create a mass of sediment that oscillates landward and seaward. While suspended, this sediment mass forms an estuarine turbidity maximum. At slack tides this sediment mass temporarily deposits on the bed, creating landward and seaward deposits. Tidal excursion and slack tide deposition limit the range of the sediment mass. To verify this conceptual model, data from two small tributary rivers of San Pablo Bay are presented. Tidal variability of suspended-sediment concentration markedly differs between the landward and seaward deposits, allowing interpretation of the intratidal movement of the oscillating sediment mass. Application of this model in suitable estuaries will assist in numerical model calibration as well as in data interpretation. A similar model has been applied to some larger-scale European estuaries, which bear a geometric resemblance to the systems analyzed in this study. ?? 2004 Elsevier Ltd. All rights reserved.

California↗
USGS2004

Floodtide pulses after low tides in shallow subembayments adjacent to deep channels

In shallow waters surface gravity waves (tides) propagate with a speed proportional to the square root of water depth (c=g(h+η)). As the ratio of free surface displacement to mean depth (η/h) approaches unity the wave will travel noticeably faster at high tide than at low tide, creating asymmetries in the tidal form. This physical process is explained analytically by the increased significance of friction and the nonlinear terms in the continuity and momentum equations. In a tidal system comprising a shallow bay adjacent to a deeper channel, tidal asymmetries will be more prevalent in the shallow bay. Thus strong barotropic gradients can be generated between the two, producing rapid accelerations of currents into the bay (relative to other bay tidal processes) and create a maximum peak in the flood tide that we describe as a floodtide pulse. These floodtide pulses can promote a landward flux of suspended-sediment into the bay. In Grizzly Bay (part of northern San Francisco Bay, USA), field observations verify the occurrence of floodtide pulses during the lowest low tides of the year. No pulses were observed in neighboring Honker Bay, which has an average depth ~30 cm greater than Grizzly Bay. Numerical simulations of northern San Francisco Bay using realistic bathymetry demonstrated that floodtide pulses occurred in Grizzly Bay but not in Honker Bay, consistent with the observations. Both observations and numerical simulations show that floodtide pulses promote a landward flux of sediment into Grizzly Bay. Numerical simulations of an idealized bay-channel system quantify the importance of mean depth and friction in creating these floodtide pulses.

California↗
USGS2003

Tidal truncation and barotropic convergence in a channel network tidally driven from opposing entrances

Residual circulation patterns in a channel network that is tidally driven from entrances on opposite sides are controlled by the temporal phasing and spatial asymmetry of the two forcing tides. The Napa/Sonoma Marsh Complex in San Francisco Bay, CA, is such a system. A sill on the west entrance to the system prevents a complete tidal range at spring tides that results in tidal truncation of water levels. Tidal truncation does not occur on the east side but asymmetries develop due to friction and off-channel wetland storage. The east and west asymmetric tides meet in the middle to produce a barotropic convergence zone that controls the transport of water and sediment. During spring tides, tidally averaged water-surface elevations are higher on the truncated west side. This creates tidally averaged fluxes of water and sediment to the east. During neap tides, the water levels are not truncated and the propagation speed of the tides controls residual circulation, creating a tidally averaged flux in the opposite direction. ?? 2003 Elsevier Science B.V. All rights reserved.

California↗