Search USGS⌕ Search

SEARCH · Search USGS

Results for “Acoustics”

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 1,009 records · Page 56Linked to original sources
USGS2002

Spatial and temporal distribution of contaminated, effluent-affected sediment on the Palos Verdes margin, southern California

A sedimentary deposit on the continental margin near the Palos Verdes Peninsula, California is comprised of sewage effluent and geologic materials and is contaminated with metals, pesticides (including DDT and associated compounds), and PCBs. The deposit was mapped with subbottom acoustic profilers, and sediment cores were analyzed for geochemical and physical properties to determine the volume of the deposit and the distribution and mass of contaminants. Mapping showed that the deposit ranges up to 60-cm thick, has a total volume exceeding 9 million m3, and covers over 40 km2. Virtually the entire effluent-affected deposit is contaminated with DDT and PCBs. Nearly half of the area of the deposit lies on the continental slope, but 70-75% of the volume of the deposit and total mass of DDT reside on the continental shelf. Analysis of data collected biennially since 1981 by the Sanitation Districts of Los Angeles County show that the mass of DDT has apparently decreased at some stations but has remained essentially constant at others. Temporal changes m mass per unit area of DDT are not statistically significant (at the 90% confidence level) at the most contaminated locations over a 16-yr period. The results of this mapping effort were used as a basis for modeling efforts described elsewhere in this issue.

California↗
USGS2003

Long-term, high-frequency current and temperature measurements along central California: Insights into upwelling/relaxation and internal waves on the inner shelf

Thermistor chains and acoustic Doppler current profilers were deployed at the northern and southern ends of Monterey Bay to examine the thermal and hydrodynamic structure of the inner (h ??? 20 m) shelf of central California. These instruments sampled temperature and current velocity at 2-min intervals over a 13-month period from June 2000 to July 2001. Time series of these data, in conjunction with SST imagery and CODAR sea surface current maps, helped to establish the basic hydrography for Monterey Bay. Analysis of time series data revealed that depth integrated flow at both sites was shore parallel (northwest-southeast) with net flows out of the Bay (northwest). The current and temperature records were dominated by semi-diurnal and diurnal tidal signals that lagged the surface tides by 3 h on average. Over the course of an internal tidal cycle these flows were asymmetric, with the flow during the flooding internal tide to the southeast typically lasting only one-third as long as the flow to the northwest during the ebbing internal tide. The transitions from ebb to flood were rapid and bore-like in nature; they were also marked by rapid increases in temperature and high shear. During the spring and summer, when thermal stratification was high, we observed almost 2000 high-frequency (Tp ??? 4-20 min) internal waves in packets of 8-10 following the heads of these bore-like features. Previous studies along the West Coast of the US have concluded that warm water bores and high-frequency internal waves may play a significant role in the onshore transport of larvae.

Continental Shelf Research↗
USGS2000

Instrumental lahar monitoring at Merapi Volcano, Central Java, Indonesia

More than 50 volcanic debris flows or lahars were generated around Mt Merapi during the first rainy season following the nuees ardentes of 22 November 1994. The rainfalls that triggered the lahars were analyzed, using such instruments as weather radar and telemetered rain gauges. Lahar dynamics were also monitored, using new non-contact detection instrumentation installed on the slopes of the volcano. These devices include real-time seismic amplitude measurement (RSAM), seismic spectral amplitude measurement (SSAM) and acoustic flow monitoring (AFM) systems. Calibration of the various systems was accomplished by field measurements of flow velocities and discharge, contemporaneously with instrumental monitoring. The 1994–1995 lahars were relatively short events, their duration in the Boyong river commonly ranging between 30 min and 1 h 30 min. The great majority (90%) of the lahars was recognized at Kaliurang village between 13:00 and 17:30 h, due to the predominance of afternoon rainfalls. The observed mean velocity of lahar fronts ranged between 1.1 and 3.4 m/s, whereas the peak velocity of the flows varied from 11 to 15 m/s, under the Gardu Pandang viewpoint location at Kaliurang, to 8–10 m/s at a section 500 m downstream from this site. River slopes vary from 28 to 22 m/km at the two sites. Peak discharges recorded in various events ranged from 33 to 360 m 3 /s, with the maximum value of peak discharge 360 m 3 /s, on 20 May 1995. To improve the lahar warning system along Boyong river, some instrumental thresholds were proposed: large and potentially hazardous lahars may be detected by RSAM units exceeding 400, SSAM units exceeding 80 on the highest frequency band, or AFM values greater than 1500 mV on the low-gain, broad-band setting.

Journal of Volcanology and Geothermal Research↗
USGS2000

Period doubling and other nonlinear phenomena in volcanic earthquakes and tremor

Evidence of subharmonic period-doubling cascades has recently been recognized in seismograms of volcanic tremor from several volcanoes. This phenomenon occurs only in nonlinear systems, and is the commonest route by which such systems change from periodic to chaotic behavior. It is predicted to occur in a model of volcanic tremor excitation by flow-induced vibration, and it might well also occur in other volcano-seismic source process. If the possibility of period doubling is not taken into account in interpreting spectra of tremor and long-period earthquakes, then low-frequency "sub-harmonic" oscillations may be mis-identified as normal modes of a linear acoustic resonator, leading to errors of an order of magnitude or more in inferred magma-body dimensions. This example illustrates the importance of nonlinear phenomena in attempts to understand volcano-seismic phenomena physically. Linear systems are fundamentally incapable of causing earthquakes or exciting tremor, so nonlinearity is essential to any theory of volcano-seismic phenomena. Nonlinear processes are in many respects qualitatively different from linear ones. A few of their characteristics that might be relevant in volcanoes include the possibility: (1) that damping might increase, rather than decrease, oscillation frequencies; and (2) that these frequencies might be functions of the amplitude of oscillation, so that temporal variations in spectral peak frequencies might not be manifestations of changes of conditions within the magmatic system.

Journal of Volcanology and Geothermal Research↗
USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2004

Pressure changes associated with the ascent and bursting of gas slugs in liquid-filled vertical and inclined conduits

At basaltic volcanoes, the sources of long-period and very-long-period seismicity and acoustic signals are frequently described in terms of fluid dynamic processes, in particular the formation and ascent of gas slugs within the magma column and their bursting at the surface. To investigate pressure changes associated with these processes, two-phase flow experiments have been carried out in vertical and inclined pipes with both single gas slugs and a continuously supplied gas phase. The ascent of individual gas slugs is accompanied by strong dynamic pressure variations resulting from the flow of liquid around the slug. These dynamic transients generate sub-static pressures below the ascending slug in viscosity-controlled systems, and produce super-static pressures when the slug reaches the surface and motion ceases in inertia-dominated systems. Conduit inclination promotes a change of regime from bubbly to slug flow and favours an increase in size and velocity of the slugs at the expense of their frequency of occurrence during continuously supplied two-phase flow. The experimental pressure data support previous theoretical analyses of oscillatory sources in ascending slugs as the slugs approach the surface and burst. Pressure oscillations are also observed during the release of gas slugs and in their wake region.

Journal of Volcanology and Geothermal Research↗
USGS2002

Measurement error associated with surveys of fish abundance in Lake Michigan

In fisheries, imprecise measurements in catch data from surveys add uncertainty to the results of fishery stock assessments. The USGS Great Lakes Science Center (GLSC) began to survey the fall fish community of Lake Michigan in 1962 with bottom trawls. The measurement error was evaluated at the level of individual tows for nine fish species collected in this survey by applying a measurementerror regression model to replicated trawl data. It was found that the estimates of measurement-error variance ranged from 0.37 (deepwater sculpin, Myoxocephalus thompsoni ) to 1.23 (alewife, Alosa pseudoharengus ) on a logarithmic scale corresponding to a coefficient of variation = 66 to 156%. The estimates appeared to increase with the range of temperature occupied by the fish species. This association may be a result of the variability in the fall thermal structure of the lake. The estimates may also be influenced by other factors, such as pelagic behavior and schooling. Measurement error might be reduced by surveying the fish community during other seasons and/or by using additional technologies, such as acoustics. Measurement-error estimates should be considered when interpreting results of assessments that use abundance information from USGS-GLSC surveys of Lake Michigan and could be used if the survey design was altered. This study is the first to report estimates of measurement-error variance associated with this survey.

Lake Michigan↗
USGS2002

Flow distribution in selected branches of St. Clair and Detroit rivers

St. Clair and Detroit rivers, which are connecting channels between Lake Huron and Lake Erie in the Great Lakes basin, form part of the boundary between the state of Michigan and the province of Ontario. In 13 reaches, this flow divides locally around islands and dikes to form 31 branches. This study develops a set of simple linear regression equations for computing expected flow proportions in branches, generally as a function of the total flow within the reach. The equations are based on 533 acoustic Doppler current profiler measurements of flow obtained between 1996 and 2000. Root-mean-square errors of these regressions range from 0.00323 to 0.0895. In seven upstream reaches where flow is known because of flow specifications at the boundaries of the waterway and continuity constraints, the uncertainties of the flow proportions can be used to directly infer the uncertainties of the corresponding flows. In six downstream reaches, the uncertainties of flows are determined by both the uncertainties of the flow proportions and the uncertainties of the total flow in the reach. For these reaches, Monte Carlo simulations quantify the ratios of total uncertainty to flow proportion uncertainty, which range from 1.0026 to 13.984. To facilitate routine calculation, polynomial regression equations are developed to approximate these ratios as a function of flow. Results provide a mechanism for computing the magnitudes and uncertainties of steady-state flows within selected branches of the connecting channels by specifying inflows at the headwaters of St. Clair River, seven intervening tributaries, and Lake St. Clair.

Michigan, Ontario↗
USGS2004

Landscape scale measures of steelhead (Oncorhynchus mykiss) bioenergetic growth rate potential in Lake Michigan and comparison with angler catch rates

The relative quality of a habitat can influence fish consumption, growth, mortality, and production. In order to quantify habitat quality, several authors have combined bioenergetic and foraging models to generate spatially explicit estimates of fish growth rate potential (GRP). However, the capacity of GRP to reflect the spatial distributions of fishes over large areas has not been fully evaluated. We generated landscape scale estimates of steelhead ( Oncorhynchus mykiss ) GRP throughout Lake Michigan for 1994-1996, and used these estimates to test the hypotheses that GRP is a good predictor of spatial patterns of steelhead catch rates. We used surface temperatures (measured with AVHRR satellite imagery) and acoustically measured steelhead prey densities (alewife, Alosa pseudoharengus ) as inputs for the GRP model. Our analyses demonstrate that potential steelhead growth rates in Lake Michigan are highly variable in both space and time. Steelhead GRP tended to increase with latitude, and mean GRP was much higher during September 1995, compared to 1994 and 1996. In addition, our study suggests that landscape scale measures of GRP are not good predictors of steelhead catch rates throughout Lake Michigan, but may provide an index of interannual variation in system-wide habitat quality.

Lake Michigan↗
USGS2022

Planning and operations of the Hydrate 01 Stratigraphic Test Well, Prudhoe Bay Unit, Alaska North Slope

The National Energy Technology Laboratory, the Japan Oil, Gas and Metals National Corporation, and the U.S. Geological Survey are leading an effort to conduct an extended gas hydrate production test in northern Alaska. The proposed production test required the drilling of an initial stratigraphic test well (STW) to confirm the geologic conditions of the proposed test site. This well was completed in January 2019 in cooperation with the Prudhoe Bay Unit Working Interest Owners. The Prudhoe Bay Unit Hydrate-01 STW was spudded on 10-December-2018. Downhole data acquisition was completed on 25-December-2018, and the rig was released on 01-January-2019. The Hydrate-01 STW was drilled in two sections, including the surface hole that was drilled to a depth of 2248 ft measured depth (MD) (685 m MD) and cased, and the production hole section that was drilled to a depth of 3558 ft MD (1084 m MD) and also cased. A thermally chilled mineral-oil-based mud was used in the main (production) hole section of the well to maintain wellbore stability and quality of the wellbore acquired data. The primary wellbore data were acquired using logging-while-drilling tools. A sidewall pressure core system was also deployed to gather grain size and other data needed for the design of the future production test wells. In addition to confirming the geologic conditions at the test site, the Hydrate-01 STW was designed to serve as a monitoring well during future field operations. Therefore, two sets of fiber-optic cables, each including a bundled distributed acoustic sensor (DAS) and a distributed temperature sensor (DTS), were clamped to the outside of the production casing and cemented in place. In March 2019, the project team acquired three-dimensional (3D) DAS vertical seismic profiling data in the Hydrate-01 STW. Temperature surveys were also acquired with the DTS as deployed in the Hydrate-01 STW during the completion of the well and nearly continuously since March-2019.

Alaska↗
USGS2023

Wide-area debris field and seabed characterization of a deep ocean dump site surveyed by autonomous underwater vehicles

Disposal of industrial and hazardous waste in the deep ocean was a pervasive global practice near coastlines in the 20th century. Uncertainty in the quantity, location, and contents of dumped materials at historical disposal sites underscores ongoing risks to marine ecosystems and human health. This study presents analysis of a 150-km2 wide-area sidescan sonar survey conducted in March 2021 with two autonomous underwater vehicles (AUVs) at an offshore dumpsite in San Pedro Basin, California. Previous camera surveys located 60 barrels and other debris; sediment analysis in the region showed varying concentrations of the insecticidal chemical Dichlorodiphenyl-trichloroethane (DDT), for which an estimated 350-700 metric tons was discarded in the San Pedro Basin between 1947 and 1961. A lack of primary historical documents from Montrose Chemical or California Salvage Company that specify the DDT acid waste disposal method has contributed to the ambiguity surrounding whether the dumping method was via bulk discharge or containerized units. Barrels and debris visually detected during previous surveys by Remotely Operated Vehicles (ROVs) Jason and SuBastian and AUV Sentry served as the basis for ground truth testing and object classification using size and acoustic intensity characteristics. Target detection and classification algorithms were developed using image and signal processing techniques that resulted in the identification of over 74,000 debris targets within the survey region. Statistical, spectral, and machine learning techniques are used to characterize seabed variability and for bottom-type classification. These analytical techniques combined with the survey capabilities of AUVs provide a framework for efficient mapping and characterization of uncharted deep-water disposal sites.

California↗
USGS1997

Review of electric and magnetic fields accompanying seismic and volcanic activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

Surveys in Geophysics↗
USGS2001

Fish species and community distributions as proxies for sea-floor habitat distributions: the Stellwagen Bank National Marine Sanctuary example (northwest Atlantic, Gulf Of Maine)

Defining the habitats of fishes and associated fauna on outer continental shelves is problematic given the paucity of data on the actual types and distributions of seafloor habitats. However many regions have good data on the distributions of fishes from resource surveys or catch statistics because of the economic importance of the fisheries. Fish distribution data (species or communities) have been used as a proxy for the distribution of habitats to develop precautionary conservation strategies for habitat protection (e.g., marine protected areas, fishing gear restrictions). In this study we assessed the relationships between the distributions of fish communities and species derived from trawl survey data with the spatial distribution of sediment types determined by sampling and acoustic reflectance derived from multibeam sonar surveys in Stellwagen Bank National Marine Sanctuary. Fish communities were correlated with reflectance values but all communities did not occur in unique sediment types. This suggests that use of community distributions as proxies for habitats should include the caveat that a greater number of communities within an area could indicate a greater range of habitat types. Single species distributions showed relationships between abundance and reflectance values. Trawl catches with low abundances had wide variations in reflectance values while those with high abundances had narrower ranges indicating habitat affinities. Significant non-random frequency-dependent relationships were observed for 17 of 20 species although only 12 of 20 species had significant relationships based on rank correlation. These results suggest that species distributions based on trawl survey data can be used as proxies for the distribution of seafloor habitats. Species with known habitat associations can be used to infer habitat requirements of co-occurring species and can be used to identify a range of habitat types.

Massachusetts↗
USGS1999

Analysis of stress-induced oval fractures in a borehole at Deep Sea Drilling Project Site 504, eastern equatorial Pacific

Deep Sea Drilling Project (DSDP) Hole 504B is located in the eastern equatorial Pacific Ocean and extends to a total depth of 2111 m beneath the seafloor (mbsf). Several acoustic televiewer logs have been obtained in this well during successive stages of drilling, and the resulting digital images have revealed numerous oval-shaped fractures seemingly etched into the borehole wall. A theoretical examination of these stress-induced features identifies a unique and ephemeral set of stress distributions and magnitudes that are necessary for their production. Consequently, the ovals provide a basis for quantifying the magnitudes and orientations of the maximum and minimum horizontal principal stresses, S H and S h at this site. Vertical, truncated breakouts and horizontal tensile fractures define the spatial boundaries of the ovals. Explicit criteria for their occurrence are combined with estimates for various physical properties of the rock to yield a range of possible values for the horizontal principal stresses. The conspicuous oval geometry is completed by a curved fracture that joins the vertical and horizontal components. Its degree of curvature is delineated by the modified Griffith failure criterion and is directly related to the principal stress difference ( S H − S h ). Matching a series of type curves corresponding to specific values for ( S H − S h ) with the actual undistorted well bore images allows the magnitude of the stress difference to be further constrained. With a value for ( S H − S h ) of 45 ± 5 MPa the individual magnitudes of S H and S h are determined more precisely. Final estimates for the horizontal principal stresses in DSDP Hole 504B at a depth of 1200 mbsf are 141 MPa ≤ S H ≤ 149 MPa and 91 MPa ≤ S h ≤ 109 MPa. Stress magnitudes derived from this approach rely heavily upon the values of a variety of physical properties, and complementary laboratory measurements performed on relevant rock samples provide critical information. Uncertainties in estimating these properties translate into less precise predictions of principal stresses.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Measurements of Reynolds stress profiles in unstratified tidal flow

In this paper we present a method for measuring profiles of turbulence quantities using a broadband acoustic doppler current profiler (ADCP). The method follows previous work on the continental shelf and extends the analysis to develop estimates of the errors associated with the estimation methods. ADCP data was collected in an unstratified channel and the results of the analysis are compared to theory. This comparison shows that the method provides an estimate of the Reynolds stresses, which is unbiased by Doppler noise, and an estimate of the turbulent kinetic energy (TKE) which is biased by an amount proportional to the Doppler noise. The noise in each of these quantities as well as the bias in the TKE match well with the theoretical values produced by the error analysis. The quantification of profiles of Reynolds stresses simultaneous with the measurement of mean velocity profiles allows for extensive analysis of the turbulence of the flow. In this paper, we examine the relation between the turbulence and the mean flow through the calculation of u * , the friction velocity, and C d the coefficient of drag. Finally, we calculate quantities of particular interest in turbulence modeling and analysis, the characteristic lengthscales, including a lengthscale which represents the stream-wise scale of the eddies which dominate the Reynolds stresses.

Journal of Geophysical Research C: Oceans↗
USGS1999

Estimates of bottom roughness length and bottom shear stress in South San Francisco Bay, California

A field investigation of the hydrodynamics and the resuspension and transport of particulate matter in a bottom boundary layer was carried out in South San Francisco Bay (South Bay), California, during March-April 1995. Using broadband acoustic Doppler current profilers, detailed measurements of turbulent mean velocity distribution within 1.5 m above bed have been obtained. A global method of data analysis was used for estimating bottom roughness length z o and bottom shear stress (or friction velocities u * ). Field data have been examined by dividing the time series of velocity profiles into 24-hour periods and independently analyzing the velocity profile time series by flooding and ebbing periods. The global method of solution gives consistent properties of bottom roughness length z o and bottom shear stress values (or friction velocities u * ) in South Bay. Estimated mean values of z o and u * for flooding and ebbing cycles are different. The differences in mean z o and u * are shown to be caused by tidal current flood-ebb inequality, rather than the flooding or ebbing of tidal currents. The bed shear stress correlates well with a reference velocity; the slope of the correlation defines a drag coefficient. Forty-three days of field data in South Bay show two regimes of z o (and drag coefficient) as a function of a reference velocity. When the mean velocity is >25–30 cm s −1 , the ln z o (and thus the drag coefficient) is inversely proportional to the reference velocity. The cause for the reduction of roughness length is hypothesized as sediment erosion due to intensifying tidal currents thereby reducing bed roughness. When the mean velocity is <25–30 cm s −1 , the correlation between z o and the reference velocity is less clear. A plausible explanation of scattered values of z o under this condition may be sediment deposition. Measured sediment data were inadequate to support this hypothesis, but the proposed hypothesis warrants further field investigation.

Journal of Geophysical Research C: Oceans↗
USGS2000

Measuring stream discharge by non-contact methods: A proof-of-concept experiment

This report describes an experiment to make a completely non-contact open-channel discharge measurement. A van-mounted, pulsed doppler (10GHz) radar collected surface-velocity data across the 183-m wide Skagit River, Washington at a USGS streamgaging station using Bragg scattering from short waves produced by turbulent boils on the surface of the river. Surface velocities were converted to mean velocities for 25 sub-sections by assuming a normal open-channel velocity profile (surface velocity times 0.85). Channel cross-sectional area was measured using a 100 MHz ground-penetrating radar antenna suspended from a cableway car over the river. Seven acoustic doppler current profiler discharge measurements and a conventional current-meter discharge measurement were also made. Three non-contact discharge measurements completed in about a 1-hour period were within 1 % of the gaging station rating curve discharge values. With further refinements, it is thought that open-channel flow can be measured reliably by non-contact methods.

Geophysical Research Letters↗
USGS2001

Creation of residual flows in a partially stratified estuary

The creation of residual flows in estuaries is examined using acoustic Doppler current profiler data sets from northern San Francisco Bay. The data sets are analyzed using principal component analysis to examine the temporal variability of the flows which create the residual circulation. It is seen that in this periodically and partially stratified estuary the residual flows are created through a series of pulses with strong variability at the 24‐hour timescale, through the interaction of shear, stratification and mixing. This interaction is captured through the use of a dimensionless number, the horizontal Richardson number ( Ri x ), which is developed to examine the local balance between the stratifying and destratifying forces at the tidal timescale. It is seen that Ri x is a valuable parameter in predicting the onset of the residual‐creating events, with a threshold value of ≈3 on ebb tides. This critical value is argued to be a threshold, above which the stratification and shear flow create a feedback effect, each further intensifying the other. This feedback results in a highly variable exchange flow which creates the estuarine residual in intermittent pulses rather than as a steady flow. Although typically attributed to baroclinic forcing, an argument is made that these pulses of residual‐creating exchange flow could be created by barotropic forcing in the presence of variable stratification which is asymmetric between flood and ebb tides. This result poses a great challenge for turbulence modeling, as the timing and magnitude of stratification and shear must be correctly simulated on the tidal timescale in order to reproduce the effects seen in the data sets presented.

Journal of Geophysical Research C: Oceans↗