Vertical Mixing Rates Research Articles

The Rapid Response of Southern Ocean Biological Productivity to Changes in Background Small Scale Turbulence

AbstractBackground subsurface vertical mixing rates in the Southern Ocean (SO) are known to vary by an order of magnitude temporally and spatially, due to variability in their generating mechanisms, which include winds and shear instabilities at the surface, and the interaction of tides and lee waves with rough bottom topography. There is great uncertainty in the parameterization of this mixing in coarse resolution Earth System Models (ESM), and in the impact that this has on SO biological productivity on sub decadal timescales. Using a data assimilating biogeochemical ocean model we show that SO phytoplankton productivity is highly sensitive to differences in background diapycnal mixing over short timescales. Changes in the background vertical mixing rates alter key biogeochemical and physical conditions. The greatest changes to the distribution of physical and biogeochemical tracers occur in regions with very strong tracer vertical gradients. A combination of reduced nutrient limitation and reduced light limitation causes a strong increase in SO phytoplankton productivity with higher background mixing. This leads to increased summer carbon export but reduced wintertime export over the mixed layer depth, which could alter the strength of the SO biological carbon pump and atmospheric concentrations on centennial to millennial timescales. This study demonstrates the importance of accurately representing diapycnal mixing in ESM to predict SO biogeochemical dynamics and their broader climatic implications.

Seasonal deoxygenation in a shallow coastal embayment: The role of stratification and implications for water-quality monitoring.

High frequency measurements at the seabed in a shallow semi-enclosed embayment of the Indian Ocean reveal deoxygenation of the bottom water at a rate of 40 mmol O2 m−3 d−1 when seasonal density stratification occurs. The vertical profile of dissolved oxygen below the pycnocline shows concentration decreasing at a constant rate toward the seabed suggesting a benthic sink and an associated downward flux of ∼30 mmol O2 m−2 d−1. This is closely coupled to an increase in CO2 indicative of microbial respiration. The results and nutrient history of the area suggest that the accumulation of organic material in the marine sediments have led to high sediment oxygen demand making the area extremely sensitive to reduced ventilation from the surface caused by low rates of vertical mixing. The data also shows that the existing weekly water-quality monitoring often misses the most severe events. Global warming is expected to worsen the problem with longer and more frequent stratification periods, reduced oxygen solubility and enhanced rates of microbial respiration.

Vertical mixing and horizontal transport unravel phytoplankton blooms in a large riverine reservoir

A precise understanding of the mechanisms causing phytoplankton blooms in reservoirs is still lacking, especially in large riverine reservoirs. To better understand these blooms, the role of the complex hydrodynamics caused by dam operation must be quantified. Here we examine how synergistic hydrodynamic processes, rather than individual metrics, trigger blooms in Xiangxi Bay, a typical tributary bay of the Three Gorges Reservoir, China. We used a 3D ecological-hydrodynamic model, which integrated hydrodynamics with the abiotic factors that limit phytoplankton growth to simulate one whole year (2010). By implementing a scaling criterion, we quantified the contribution of local phytoplankton growth and hydrodynamic processes, including advection transport and vertical mixing, on bloom dynamics. Results indicated vertical mixing was the main process inhibiting blooms in colder months (from October to February) but horizontal advection, which flushed and diluted blooms, was dominant in warmer months (from May to July) when stratification was intense and nutrients were replete. Accordingly, blooms occurred when both vertical mixing and horizontal advection were low. We suggested a potential dam operation strategy to mitigate blooms during stratification, which involves withdrawing the warm surface water from upstream reservoirs to increase horizontal flows in the surface layer. Extending the application of critical turbulence model, our study shows how vertical mixing and horizontal advection rate interact with phytoplankton growth rate to drive blooms in highly dynamic riverine systems.

A Possible Hysteresis in the Arctic Ocean due to Release of Subsurface Heat during Sea Ice Retreat

Abstract The Arctic Ocean is characterized by an ice-covered layer of cold and relatively fresh water above layers of warmer and saltier water. It is estimated that enough heat is stored in these deeper layers to melt all the Arctic sea ice many times over, but they are isolated from the surface by a stable halocline. Current vertical mixing rates across the Arctic Ocean halocline are small, due in part to sea ice reducing wind–ocean momentum transfer and damping internal waves. However, recent observational studies have argued that sea ice retreat results in enhanced mixing. This could create a positive feedback whereby increased vertical mixing due to sea ice retreat causes the previously isolated subsurface heat to melt more sea ice. Here, we use an idealized climate model to investigate the impacts of such a feedback. We find that an abrupt “tipping point” can occur under global warming, with an associated hysteresis window bounded by saddle-node bifurcations. We show that the presence and magnitude of the hysteresis are sensitive to the choice of model parameters, and the hysteresis occurs for only a limited range of parameters. During the critical transition at the bifurcation point, we find that only a small percentage of the heat stored in the deep layer is released, although this is still enough to lead to substantial sea ice melt. Furthermore, no clear relationship is apparent between this change in heat storage and the level of hysteresis when the parameters are varied.

Marine Plastic Drift from the Mekong River to Southeast Asia

Southeast Asia is the world’s most polluted area in terms of marine plastics. The Mekong River is one of the largest rivers in the area, and ranked as somewhere between the 8th- and 11th-biggest contributor to plastics in the world’s oceans. Here, we investigate how microplastics drift from the Mekong river to Southeast Asia, and which coastlines are most exposed. We identify potential factors (wind drift, rivers, vertical mixing and sinking rates) that affect plastic drift in the region using the OpenDrift model with realistic wind and ocean currents for simulations between three months (summer and winter) and 15 months. We find that the seasonal drift is influenced by the monsoon systems and that most of the plastics strand in the Philippines and Indonesia. In addition, the role of wind drift is significant in strong winds. Vertical mixing and sinking rates are unknowns that affect the relative importance of wind drift (near the surface) and ocean currents. Simulations with different terminal velocities show that, unsurprisingly, the higher the terminal velocities are, the closer they deposit to the source. In light of the large uncertainties in sinking rates, we find that the plastic distribution has large uncertainties, but is clearly seasonal and influenced by wind, vertical mixing, river discharge and sinking rates. The Philippines and Indonesia are found to have the coastlines that are most exposed to plastic pollution from the Mekong river. This study shows that simulations of marine plastic drift are very variable, depending on many factors and assumptions. However, it provides more detailed information on marine plastic pollution in Southeast Asia, and hopefully helps authorities take more practical actions.

Aerosol deposition to the boreal forest in the vicinity of the Alberta Oil Sands

Abstract. Measurements of size-resolved aerosol concentration and fluxes were made in a forest in the Athabasca Oil Sands Region (AOSR) of Alberta, Canada, in August 2021 with the aim of investigating (a) particle size distributions from different sources, (b) size-resolved particle deposition velocities, and (c) the rate of vertical mixing in the canopy. Particle size distributions were attributed to different sources determined by wind direction. Air mixed with smokestack plumes from oil sands processing facilities had higher number concentrations with peak number at diameters near 70 nm. Aerosols from the direction of open-pit mine faces showed number concentration peaks near 150 nm and volume distribution peaks near 250 nm (with secondary peaks near 600 nm). Size-resolved deposition fluxes were calculated which show good agreement with previous measurements and a recent parameterization. There is a minimum deposition velocity of vd=0.02 cm s−1 for particles of 80 nm diameter; however, there is a large amount of variation in the measurements, and this value is not significantly different from zero in the 68 % confidence interval. Finally, gradient measurements of aerosol particles (with diameters <1 µm) demonstrated nighttime decoupling of air within and above the forest canopy, with median lag times at night of up to 40 min and lag times between 2 and 5 min during the day. Aerosol mass fluxes (diameters <1 µm) determined using flux–gradient methods (with different diffusion parameterizations) underestimate the flux magnitude relative to eddy covariance flux measurements when averaged over the nearly 1-month measurement period. However, there is significant uncertainty in the averages determined using the flux–gradient method.

L-band spectroscopy of young brown dwarfs

ABSTRACT We present a L-band (2.98–3.96 $\mu$m) spectroscopic study of eight young L dwarfs with spectral types ranging from L2 to L7. Our spectra (λ/Δλ ≈ 250–600) were collected using the Gemini near-infrared spectrograph. We first examine the young L-band spectral sequence, most notably analysing the evolution of the Q-branch of methane absorption feature at 3.3 $\mu$m. We find the Q-branch feature first appears between L3 and L6, as previously seen in older field dwarfs. Secondly, we analyse how well various atmospheric models reproduce the Lband and published near-IR (0.7–2.5 $\mu$m) spectra of our objects by fitting five different grids of model spectra to the data. Best-fit parameters for the combined near-IR and L-band data are compared to best-fit parameters for just the near-IR data, isolating the impact that the addition of the L band has on the results. This addition notably causes a ∼100 K drop in the best-fit effective temperature. Also, when clouds and a vertical mixing rate (Kzz) are included in the models, thick clouds, and higher Kzz values are preferred. Five of our objects also have previously published effective temperatures and surface gravities derived using evolutionary models, age estimates, and bolometric luminosities. Comparing model spectra matching these parameters to our spectra, we find disequilibrium chemistry and clouds are needed to match these published effective temperatures. Three of these objects are members of AB Dor, allowing us to show the temperature dependence of the Q-branch of methane.

The Origin and Fate of Antarctic Intermediate Water in the Southern Ocean

Abstract Using observationally based hydrographic and eddy diffusivity datasets, a volume budget analysis is performed to identify the main mechanisms governing the spatial and seasonal variability of Antarctic Intermediate Water (AAIW) within the density range γn = (27.25–27.7) kg m−3 in the Southern Ocean. The subduction rates and water mass transformation rates by mesoscale and small-scale turbulent mixing are estimated. First, Ekman pumping upwells the dense variety of AAIW into the mixed layer south of the Polar Front, which can be advected northward by Ekman transport into the subduction regions of lighter-variety AAIW and Subantarctic Mode Water (SAMW). The subduction of light AAIW occurs mainly by lateral advection in the southeast Pacific and Drake Passage as well as eddy-induced flow between the Subantarctic and Polar Fronts. The circumpolar-integrated total subduction yields from −5 to 19 Sv (1 Sv ≡ 106 m3 s−1) of AAIW volume loss. Second, the diapycnal transport from subducted SAMW into the AAIW layer is predominantly by mesoscale mixing (2–13 Sv) near the Subantarctic Front and vertical mixing in the South Pacific, while AAIW is further replenished by transformation from Upper Circumpolar Deep Water by vertical mixing (1–10 Sv). Last, 3–14 Sv of AAIW are exported out of the Southern Ocean. Our results suggest that the distribution of AAIW is set by its formation due to subduction and mixing, and its circulation eastward along the ACC and northward into the subtropical gyres. The volume budget analysis reveals strong seasonal variability in the rate of subduction, vertical mixing, and volume transport driving volume change within the AAIW layer. The nonzero volume budget residual suggests that more observations are needed to better constrain the estimate of geostrophic flow and mesoscale and small-scale mixing diffusivities.

Effect of nutrient supply on cell size evolution of marine phytoplankton.

The variation of nutrient supply not only leads to the differences in the phytoplankton biomass and primary productivity but also induces the long-term phenotypic evolution of phytoplankton. It is widely accepted that marine phytoplankton follows Bergmann's Rule and becomes smaller with climate warming. Compared with the direct effect of increasing temperature, the indirect effect via nutrient supply is considered to be an important and dominant factor in the reduction of phytoplankton cell size. In this paper, a size-dependent nutrient-phytoplankton model is developed to explore the effects of nutrient supply on the evolutionary dynamics of functional traits associated with phytoplankton size. The ecological reproductive index is introduced to investigate the impacts of input nitrogen concentration and vertical mixing rate on the persistence of phytoplankton and the distribution of cell size. In addition, by applying the adaptive dynamics theory, we study the relationship between nutrient input and the evolutionary dynamics of phytoplankton. The results show that input nitrogen concentration and vertical mixing rate have significant effects on the cell size evolution of phytoplankton. Specifically, cell size tends to increase with the input nutrient concentration, as does the diversity of cell sizes. In addition, a single-peaked relationship between vertical mixing rate and cell size is observed. When the vertical mixing rate is too low or too high, only small individuals are dominant in the water column. When the vertical mixing rate is moderate, large individuals can coexist with small individuals, so the diversity of phytoplankton is elevated. We predict that reduced intensity of nutrient input due to climate warming will lead to a trend towards smaller cell size and will reduce the diversity of phytoplankton.

Vertical Mixing and Heat Fluxes Conditioned by a Seismically Imaged Oceanic Front

The southwest Atlantic gyre connects several distinct water masses, which means that this oceanic region is characterized by a complex frontal system and enhanced water mass modification. Despite its significance, the distribution and variability of vertical mixing rates have yet to be determined for this system. Specifically, potential conditioning of mixing rates by frontal structures, in this location and elsewhere, is poorly understood. Here, we analyze vertical seismic (i.e., acoustic) sections from a three-dimensional survey that straddles a major front along the northern portion of the Brazil-Falkland Confluence. Hydrographic analyses constrain the structure and properties of water masses. By spectrally analyzing seismic reflectivity, we calculate spatial and temporal distributions of the dissipation rate of turbulent kinetic energy, ε, of diapycnal mixing rate, K, and of vertical diffusive heat flux, FH. We show that estimates of ε, K, and FH are elevated compared to regional and global mean values. Notably, cross-sectional mean estimates vary little over a 6 week period whilst smaller scale thermohaline structures appear to have a spatially localized effect upon ε, K, and FH. In contrast, a mesoscale front modifies ε and K to a depth of 1 km, across a region of O(100) km. This front clearly enhances mixing rates, both adjacent to its surface outcrop and beneath the mixed layer, whilst also locally suppressing ε and K to a depth of 1 km. As a result, estimates of FH increase by a factor of two in the vicinity of the surface outcrop of the front. Our results yield estimates of ε, K and FH that can be attributed to identifiable thermohaline structures and they show that fronts can play a significant role in water mass modification to depths of 1 km.

Winter Dynamics in an Epishelf Lake: Quantitative Mixing Estimates and Ice Shelf Basal Channel Considerations

AbstractMilne Ice Shelf is located at the mouth of Milne Fiord (82.6°N, 81.0°W), on Ellesmere Island, Nunavut. This floating ice feature is attached to both sides of the fjord. During the melt season, the ice shelf acts as a floating barrier preventing surface runoff in Milne Fiord from flowing freely to the ocean. This results in a permanent layer of freshwater that floats on top of the seawater of the fjord, commonly known as an epishelf lake. Few studies address the physical characteristics of epishelf lakes. Here, winter data from a mooring installed in Milne Fiord epishelf lake during 2011–2019 is analyzed in the framework of a one‐dimensional model in order to study mixing in the upper water column and infer the characteristics of a basal channel in the ice shelf. The results show that vertical mixing rates are higher in the epishelf lake than in the seawater below. Estimation of the Richardson number using a geostrophic balance approach reveals that enhanced mixing in the epishelf lake is associated with horizontal temperature gradients. In contrast to most studies on basal channels in Greenland and Antarctica, the results presented here show the basal channel of Milne Ice Shelf is apparently not evolving (melting) rapidly.

Radium-228-derived ocean mixing and trace element inputs in the South Atlantic

Abstract. Trace elements (TEs) play important roles as micronutrients in modulating marine productivity in the global ocean. The South Atlantic around 40∘ S is a prominent region of high productivity and a transition zone between the nitrate-depleted subtropical gyre and the iron-limited Southern Ocean. However, the sources and fluxes of trace elements to this region remain unclear. In this study, the distribution of the naturally occurring radioisotope 228Ra in the water column of the South Atlantic (Cape Basin and Argentine Basin) has been investigated along a 40∘ S zonal transect to estimate ocean mixing and trace element supply to the surface ocean. Ra-228 profiles have been used to determine the horizontal and vertical mixing rates in the near-surface open ocean. In the Argentine Basin, horizontal mixing from the continental shelf to the open ocean shows an eddy diffusion of Kx=1.8±1.4 (106 cm2 s−1) and an integrated advection velocity w=0.6±0.3 cm s−1. In the Cape Basin, horizontal mixing is Kx=2.7±0.8 (107 cm2 s−1) and vertical mixing Kz = 1.0–1.7 cm2 s−1 in the upper 600 m layer. Three different approaches (228Ra diffusion, 228Ra advection, and 228Ra/TE ratio) have been applied to estimate the dissolved trace element fluxes from the shelf to the open ocean. These approaches bracket the possible range of off-shelf fluxes from the Argentine Basin margin to be 4–21 (×103) nmol Co m−2 d−1, 8–19 (×104) nmol Fe m−2 d−1 and 2.7–6.3 (×104) nmol Zn m−2 d−1. Off-shelf fluxes from the Cape Basin margin are 4.3–6.2 (×103) nmol Co m−2 d−1, 1.2–3.1 (×104) nmol Fe m−2 d−1, and 0.9–1.2 (×104) nmol Zn m−2 d−1. On average, at 40∘ S in the Atlantic, vertical mixing supplies 0.1–1.2 nmol Co m−2 d−1, 6–9 nmol Fe m−2 d−1, and 5–7 nmol Zn m−2 d−1 to the euphotic zone. Compared with atmospheric dust and continental shelf inputs, vertical mixing is a more important source for supplying dissolved trace elements to the surface 40∘ S Atlantic transect. It is insufficient, however, to provide the trace elements removed by biological uptake, particularly for Fe. Other inputs (e.g. particulate or from winter deep mixing) are required to balance the trace element budgets in this region.

Oxygen Seasonality, Utilization Rate, and Impacts of Vertical Mixing in the Eighteen Degree Water Region of the Sargasso Sea as Observed by Profiling Biogeochemical Floats

AbstractSeasonal oxygen structure and utilization in the Sargasso Sea are characterized using nine profiling floats with oxygen 2021 sensors (years 2005–2008), deployed in an Eighteen Degree Water (EDW) experiment (CLIMODE). During autumn‐winter when the mixed layer is deepening, oxygen increases from the surface to the base of the EDW at 400 m. During spring‐summer, oxygen decreases except between the seasonal pycnocline and compensation depth, creating the seasonal shallow oxygen maximum layer (SOMax) with oxygen production of 0.04 μmol kg−1·day−1. In the underlying seasonal oxygen minimum (SOMin), the oxygen utilization rate (OUR) is 0.10 μmol kg−1·day−1, decreasing with depth to 0.04 μmol kg−1·day−1 in the EDW. Remineralization in May to August is double that of August to November. The Sargasso Sea is a net carbon producer; estimated annual export production from the top 100–250 m is 2.9 mol C m−2 and from the top 400 m is 4.2 mol C m−2. Below the EDW, oxygen decreases seasonally at the same time as in the EDW, indicating remineralization down to 700 m. However, on isopycnals in this deeper layer, oxygen increases during May to September, likely due to lateral advection from nonlocal surface outcrops. Summer shoaling of these isopycnals creates this paradox. The complex vertical oxygen structure in the upper 200 m enables important vertical diffusive flux that modifies the OUR calculated from oxygen change. Ignoring mixing underestimates maximum remineralization by 19% and underestimates maximum net production by 88%. However, vertical mixing is negligible in the deeper layers, so the associated total integrated remineralization error is 5%–9%.

Radium isotope ratios as a tool to characterise nutrient dynamics in a variably stratified temperate fjord

The effect of freshwater discharge on solute vertical mixing and horizontal transfer rates is evaluated in a variably stratified temperate fjord (Killary Harbour, Ireland) using the activity ratios of 224Ra and 223Ra in the water column. With high river discharge (>1.7 × 106 m3 d−1), surface 224Ra/223Ra activity ratios (ARs) decreased from 18.0 ± 7.1 to 8.6 ± 0.9 from the inner to the outer fjord and the system becomes stratified, resulting in clear differences of 224Ra/223Ra ARs between surface and bottom layers. When river discharge was low (<0.6 × 106 m3 d−1), 224Ra/223Ra ARs dipped from 19.6 ± 5.1 to 8.0 ± 1.8 in the inner system before increasing again towards the outlet to the sea up to 11.3 ± 3.6, and the water column was well mixed: the difference in 224Ra/223Ra ARs and salinity between surface and bottom layers was minimal. Longitudinal seaward decreases of surface water 224Ra/223Ra ARs within the inner fjord independently of river discharge suggested that the contribution of benthic-pelagic coupling to surface solute inventories is small in the inner section of the system. Conversely, the increase of 224Ra/223Ra ARs during low discharge conditions indicated that the slower flushing of the system, coupled with a well-mixed water column, facilitated solute transfers from deeper waters to the surface in the outer fjord. The effect of discharge on the residence time of water in the inner fjord was also determined using radium ages: longitudinal water flushing timescales are quickly reduced in the inner part of the bay with increasing discharge. We show that the degree of vertical mixing is an important driver of nutrient ratios in the system. Retention of soluble reactive phosphorus (SRP) and Ra in deep water and a higher molar N:P ratio in surface water were more likely under stratified conditions. Stratification may therefore amplify the effect of variable watershed inputs on the estuarine N:P ratio in Killary Harbour.

Analysis of Pollutant Entrainment from Localized Sources in a Street Network

The propagation of a pollutant emitted from localized sources both within and above a regular street network is studied by analyzing data from direct numerical simulations of passive scalar dispersion. Two wind directions are considered, corresponding to aligned and oblique flow with respect to the street axes. Particular attention is paid to the role of entrainment of the scalar into the urban canopy from an elevated source and re-entrainment of material originally released further upstream from a ground source. The variation of concentration differences and vertical fluxes between the streets and the air above as a function of distance reveals important differences between the rate of lateral and vertical mixing for the two sources. Detrainment and entrainment need a longer fetch to equilibrate for the elevated source than for the ground source. There are large differences between the advection and detrainment velocities for the aligned and oblique cases, so that a change in wind direction could affect ventilation efficiency considerably. Time scales associated with different dispersion processes are computed and the time of first appearance of the scalar from the onset of release in different streets is mapped. It is shown that re-entrainment can provide a shortcut dispersion pathway for reaching certain parts of the network. This is particularly striking in the case of oblique flow, when material can be transferred by entrainment up to twice as rapidly as by advection. Taken together, these results highlight the overall message that vertical exchange is a two-way process and that entrainment needs to be considered in the context of emergency response as well as urban ventilation.

Role of internal tide mixing in keeping the deep Andaman Sea warmer than the Bay of Bengal

Vertical profiles of temperature obtained from various hydrographic datasets show that deep waters (below 1,200 m) in the Andaman Sea are warmer (about 2 °C) than that of the Bay of Bengal. As a result, the biochemical properties in the deep waters also exhibit significant differences between these two basins. Higher temperature in the deep waters of Andaman Sea compared to the BoB had been widely attributed to the enclosed nature of the Andaman Sea. In this study, we show that strong tidal energy dissipation in the Andaman Sea also plays an important role in maintaining the higher temperatures in the deep waters. Dissipation rates inferred from the hydrographic data and internal tide energy budget suggests that the rate of vertical mixing in the Andaman Sea is about two-orders of magnitude larger than that in the Bay of Bengal. This elevated internal tide induced vertical mixing results in the efficient transfer of heat into the deeper layers, which keeps the deep Andaman Sea warm. Numerical experiments conducted using a high-resolution setup of Regional Ocean Modelling System (ROMS) further confirm the effect of tidal mixing in the Andaman Sea.

Circulation and Thermohaline Variability of the Hanna Shoal Region on the Northeastern Chukchi Sea Shelf

AbstractWe analyzed velocity and hydrographic data from 23 moorings in the northeast Chukchi Sea from 2011 to 2014. In most years the eastern side of Hanna Shoal was strongly stratified year‐round, while weakly stratified regions prevailed on the shelf south and west of the Shoal. Stratification differences cause differential vertical mixing rates, which in conjunction with advection of different bottom water properties resulted in seasonally varying along‐isobath density gradients. In agreement with numerical models, we find that bottom waters flow anticyclonically around the Shoal. Whereas most of the shelf responded barotropically to wind‐forcing, there was a strong baroclinic component to the flow field northeast of Hanna Shoal, resulting in no net vertically integrated transport on average. In contrast there is a net eastward transport from west of the Shoal, which implies convergence north of the Shoal. Convergence and along‐isobath density gradients may foster cross‐shelf exchange north of Hanna Shoal. Modal analyses indicate that the shelf south of the Shoal and Barrow Canyon responded coherently to local and remote winds, whereas the wind‐current response around Hanna Shoal was less coherent. Barotropic topographic waves, of ~3‐day period, were generated episodically northeast of the Shoal and propagate clockwise around Hanna Shoal, but are blocked from entering Barrow Canyon and are possibly scattered by the horizontally sheared flow and converging isobaths on the western side of the Shoal. Analysis of water properties on the western side of Hanna Shoal suggests that these include contributions from the western and southern portions of the Chukchi Sea.

Tracking the global flows of atmospheric moisture and associated uncertainties

Abstract. Many processes in hydrology and Earth system science relate to continental moisture recycling, the contribution of terrestrial evaporation to precipitation. For example, the effects of land-cover changes on regional rainfall regimes depend on this process. To study moisture recycling, a range of moisture-tracking models are in use that are forced with output from atmospheric models but differ in various ways. They can be Eulerian (grid-based) or Lagrangian (trajectory-based), have two or three spatial dimensions, and rely on a range of other assumptions. Which model is most suitable depends not only on the purpose of the study but also on the quality and resolution of the data with which it is forced. Recently, the high-resolution ERA5 reanalysis data set has become the state of the art, paving the way for a new generation of moisture-tracking models. However, it is unclear how the new data can best be used to obtain accurate estimates of atmospheric moisture flows. Here we develop a set of moisture-tracking models forced with ERA5 data and systematically test their performance regarding continental evaporation recycling ratio, distances of moisture flows, and “footprints” of evaporation from seven point sources across the globe. We report simulation times to assess possible trade-offs between accuracy and speed. Three-dimensional Lagrangian models were most accurate and ran faster than Eulerian versions for tracking water from single grid cells. The rate of vertical mixing of moisture in the atmosphere was the greatest source of uncertainty in moisture tracking. We conclude that the recently improved resolution of atmospheric reanalysis data allows for more accurate moisture tracking results in a Lagrangian setting, but that considerable uncertainty regarding turbulent mixing remains. We present an efficient Lagrangian method to track atmospheric moisture flows from any location globally using ERA5 reanalysis data and make the code for this model, which we call UTrack-atmospheric-moisture, publicly available.

The Influence of Submesoscales and Vertical Mixing on the Export of Sinking Tracers in Large-Eddy Simulations

AbstractWe use idealized large-eddy simulations (LES) and a simple analytical theory to study the influence of submesoscales on the concentration and export of sinking particles from the mixed layer. We find that restratification of the mixed layer following the development of submesoscales reduces the rate of vertical mixing which, in turn, enhances the export rate associated with gravitational settling. For a neutral tracer initially confined to the mixed layer, subinertial (submesoscale) motions enhance the downward tracer flux, consistent with previous studies. However, the sign of the advective flux associated with the concentration of sinking particles reverses, indicating reentrainment into the mixed layer. A new theory is developed to model the gravitational settling flux when the particle concentration is nonuniform. The theory broadly agrees with the LES results and allows us to extend the analysis to a wider range of parameters.

Trace element and Sr isotopic composition of bottom and near-surface oceanic water in the southern region of the Emperor Ridge

During the 86th voyage of the research vessel “Akademik M.A. Lavrentiev” (July–August 2019), near-bottom and near-surface water samples (at a depth of 10 m) were taken in the immediate vicinity and above the Koko, Jingu, Nintoku, and Suiko seamounts in the southern part of the Emperor seamount chain. Trace element concentrations (including REE) in these samples are higher than in the IAPSO Standard Atlantic Seawater. The distribution and mutual correlations of elements involved in biochemical processes indicate an active role of plankton in the accumulation of these elements, as well as an increase in its amount with increasing water temperature. 87Sr/86Sr ratios in water samples vary from 0.70891 to 0.70917. Moreover, near-bottom water, taken farther south of 44°N, is more radiogenic than near-surface. This difference can be explained by the dilution of the surface layer with meteoric water, because at the time of sample collection over the Koko, Jingu, and Nintoku seamounts it was raining daily, while no rainfall was observed over the Suiko seamount during the sample collection and the week before it. The rapid restoration of isotopic equilibrium indicates a high rate of vertical mixing of ocean water, which was indicated by the earlier studies of oceanic water chemical composition.