An official website of the United States government
Abstract In this study, we report on turbulent mixing observed during the annual stratification cycle in the hypolimnetic waters of Lake Michigan (USA), highlighting stratified, convective, and transitional mixing periods. Measurements were collected using a combination of moored instruments and microstructure profiles. Observations during the stratified summer showed a shallow, wind‐driven surface mixed layer (SML) with locally elevated dissipation rates in the thermocline () potentially associated with internal wave shear. Below the thermocline, turbulence was weak () and buoyancy‐suppressed (< 8.5), with low hypolimnetic mixing rates () limiting benthic particle delivery. During the convective winter period, a diurnal cycle of radiative convection was observed over each day of measurement, where temperature overturns were directly correlated with elevated turbulence levels throughout the water column (;). A transitional mixing period was observed for spring conditions when surface temperatures were near the temperature of maximum density (TMD3.98) and the water column began to stably stratify. While small temperature gradients allowed strong mixing over the transitional period (), hypolimnetic velocity shear was overwhelmed by weakly stable stratification (;), limiting the development of the SML. These results highlight the importance of radiative convection for breaking down weak hypolimnetic stratification and driving energetic, full water column mixing during a substantial portion of the year (>100 days at our sample site). Ongoing surface water warming in the Laurentian Great Lakes is significantly reducing the annual impact of convective mixing, with important consequences for nutrient cycling, primary production, and benthic‐pelagic coupling.
more »
« less
Global Contrasts Between Oceanic Cycling of Cadmium and Phosphate
Abstract Cadmium (Cd) is a trace metal whose distribution in the ocean bears a remarkable resemblance to the nutrient phosphate (PO43−). This resemblance has led to the use of Cd as a proxy for ocean nutrient cycling in paleoceanographic applications, but the processes governing the cycling of Cd in the modern ocean remain unclear. In this study, we use previously published Cd observations and an Artificial Neural Network to produce a dissolved Cd climatology that reproduces the observed subtle deviations between the Cd anddistributions. We use the Cd andclimatologies, along with an ocean circulation inverse model, to diagnose the biogeochemical sources and sinks of dissolved Cd and. Our calculations reveal that dissolved Cd, like, is removed in the surface ocean and has a source in the subsurface, consistent with the simultaneous incorporation of Cd andinto sinking organic particles. However, there are also contrasts between the cycling of dissolved Cd andIn particular, thesurface export ratio varies 8‐fold across latitudes, reaching highest values in the iron‐limited sub‐Antarctic Southern Ocean. This depletes Cd relative toin the low‐latitude thermocline while adding excess Cd to deep waters by the regeneration of Cd‐enriched particles. Also, Cd tends to regenerate slightly deeper thanin the subsurface ocean, and theregeneration ratio reaches a maximum at 700–1,500 m. These contrasts are responsible for a slight concavity in therelationship and should be considered when interpreting paleoceanographic Cd records.
more »
« less
- Award ID(s):
- 1658392
- PAR ID:
- 10364887
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Global Biogeochemical Cycles
- Volume:
- 35
- Issue:
- 6
- ISSN:
- 0886-6236
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The UV finiteness found in calculations of the 4‐point amplitude insupergravity at loop orderhas not been explained, which motivates our study of the relevant superspace invariants and on‐shell superamplitudes for bothand. The local 4‐point superinvariants forare expected to have nonlinear completions whose 6‐point amplitudes have non‐vanishing SSL's (soft scalar limits), violating the behavior required of Goldstone bosons. For, we find atthat local 6‐point superinvariant and superamplitudes, which might cancel these SSL's, do not exist. This rules out the candidate 4‐point counterterm and thus gives a plausible explanation of the observedfiniteness. However, atwe construct a local 6‐point superinvariant with non‐vanishing SSL's, so the SSL argument does not explain the observedUV finiteness. Forsupergravity there are no 6‐point invariants at eitheror 4, so the SSL argument predicts UV finiteness.more » « less
-
Abstract Estimates of turbulence kinetic energy (TKE) dissipation rate (ε) are key in understanding how heat, gas, and other climate‐relevant properties are transferred across the air‐sea interface and mixed within the ocean. A relatively new method involving moored pulse‐coherent acoustic Doppler current profilers (ADCPs) allows for estimates ofεwith concurrent surface flux and wave measurements across an extensive length of time and range of conditions. Here, we present 9 months of moored estimates ofεat a fixed depth of 8.4 m at the Stratus mooring site (20°S, 85°W). We find that turbulence regimes are quantified similarly using the Obukhov length scaleand the newer Langmuir stability length scale, suggesting that ocean‐side friction velocityimplicitly captures the influence of Langmuir turbulence at this site. This is illustrated by a strong correlation between surface Stokes driftandthat is likely facilitated by the steady Southeast trade winds regime. In certain regimes,, whereis the von Kármán constant andis instrument depth, and surface buoyancy flux capture our estimates ofwell, collapsing data points near unity. We find that a newer Langmuir turbulence scaling, based onand, scalesεwell at times but is overall less consistent than. Monin‐Obukhov similarity theory (MOST) relationships from prior studies in a variety of aquatic and atmospheric settings largely agree with our data in conditions where convection and wind‐driven current shear are both significant sources of TKE, but diverge in other regimes.more » « less
-
Abstract In field observations from a sinuous estuary, the drag coefficientbased on the momentum balance was in the range of, much greater than expected from bottom friction alone.also varied at tidal and seasonal timescales.was greater during flood tides than ebbs, most notably during spring tides. The ebb tidewas negatively correlated with river discharge, while the flood tideshowed no dependence on discharge. The large values ofare explained by form drag from flow separation at sharp channel bends. Greater water depths during flood tides corresponded with increased values of, consistent with the expected depth dependence for flow separation, as flow separation becomes stronger in deeper water. Additionally, the strength of the adverse pressure gradient downstream of the bend apex, which is indicative of flow separation, correlated withduring flood tides. Whilegenerally increased with water depth,decreased for the highest water levels that corresponded with overbank flow. The decrease inmay be due to the inhibition of flow separation with flow over the vegetated marsh. The dependence ofduring ebbs on discharge corresponds with the inhibition of flow separation by a favoring baroclinic pressure gradient that is locally generated at the bend apex due to curvature‐induced secondary circulation. This effect increases with stratification, which increases with discharge. Additional factors may contribute to the high drag, including secondary circulation, multiple scales of bedforms, and shallow shoals, but the observations suggest that flow separation is the primary source.more » « less
-
Abstract We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating theslender body inverse problem, where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius, we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff. For both the truncated and‐regularized approximations, we obtain rigorous‐based convergence to the PDE solution as: A fiber velocity withregularity givesconvergence, while a fiber velocity with at leastregularity yieldsconvergence. Moreover, we determine the dependence of the‐regularized error estimate on the regularization parameter.more » « less
