Horizontal Mixing Research Articles

Modeling the Pathways and Accumulation Patterns of Micro- and Macro-Plastics in the Mediterranean

The Mediterranean is considered a hot-spot for plastic pollution, due to its semi-enclosed nature and heavily populated coastal areas. In the present study, a basin-scale coupled hydrodynamic/particle drift model was used to track the pathways and fate of plastics from major land-based sources (coastal cities and rivers), taking into account of the most important processes (advection, stokes drift, vertical and horizontal mixing, sinking, wind drag, and beaching). A hybrid ensemble Kalman filter algorithm was implemented to correct the near- surface circulation, assimilating satellite data (sea surface height, temperature) in the hydrodynamic model. Different size classes and/or types of both micro- and macroplastics were considered in the model. Biofouling induced sinking was explicitly described, as a possible mechanism of microplastics removal from the surface. A simplified parameterization of size-dependent biofilm growth has been adopted, as a function of bacterial biomass (obtained from a biogeochemical model simulation), being considered a proxy for the biofouling community. The simulated distributions for micro- and macroplastics were validated against available observations, showing reasonable agreement, both in terms of magnitude and horizontal variability. An 8-year simulation was used to identify micro- and macroplastics accumulation patterns in the surface layer, water column, seafloor and beaches. The impact of different processes (vertical mixing, biofouling, and wind/wave drift) was identified through a series of sensitivity experiments. For both micro- and macroplastics, distributions at sea surface were closely related to the adopted sources. The microplastics concentration was drastically reduced away from source areas, due to biofouling induced sinking, with their size distribution dominated by larger (>1 mm) size classes in open sea areas, in agreement with observations. High concentration patches of floating plastics were simulated in convergence areas, characterized by anticyclonic circulation. The distribution of macroplastics on beaches followed the predominant southeastward wind/wave direction. In the water column, a sub-surface maximum in microplastics abundance was simulated, with increasing contribution of smaller particles in deeper layers. Accumulation of microplastics on the seafloor was limited in relatively shallow areas (<500 m), with bottom depth below their relaxation depth due to defouling. The simulated total amount of floating plastics (∼3,760 tonnes) is comparable with estimates from observations.

Tropical Cyclone Intensification Simulated in the Ooyama-type Three-layer Model with a Multilevel Boundary Layer

AbstractThe first successful simulation of tropical cyclone (TC) intensification was achieved with a three-layer model, often named the Ooyama-type three-layer model, which consists of a slab boundary layer and two shallow water layers above. Later studies showed that the use of a slab boundary layer would produce unrealistic boundary layer wind structure and too strong eyewall updraft at the top of TC boundary layer and thus simulate unrealistically rapid intensification compared to the use of a height-parameterized boundary layer. To fully consider the highly height-dependent boundary layer dynamics in the Ooyama-type three-layer model, this study replaced the slab boundary layer with a multilevel boundary layer in the Ooyama-type model and used it to conduct simulations of TC intensification and also compared the simulation with that from the model version with a slab boundary layer. Results show that compared with the simulation with a slab boundary layer, the use of a multilevel boundary layer can greatly improve simulations of the boundary-layer wind structure and the strength and radial location of eyewall updraft, and thus more realistic intensification rate due to better treatments of the surface layer processes and the nonlinear advection terms in the boundary layer. Sensitivity of the simulated TCs to the model configuration and to both horizontal and vertical mixing lengths, sea surface temperature, the Coriolis parameter, and the initial TC vortex structure are also examined. The results demonstrate that this new model can reproduce various sensitivities comparable to those found in previous studies using fully physics models.

Sensitivity of Typhoon Lingling (2019) simulations to horizontal mixing length and planetary boundary layer parameterizations

Sensitivity of Typhoon Lingling (2019) simulations to horizontal mixing length and planetary boundary layer parameterizations

Meridional Variations of C2H2 in Jupiter's Stratosphere From Juno UVS Observations

AbstractThe Ultraviolet Spectrograph (UVS) instrument on the Juno mission records far‐ultraviolet reflected sunlight from Jupiter. These spectra are sensitive to the abundances of chemical species in the upper atmosphere and to the distribution of the stratospheric haze layer. We combine observations from the first 30 perijoves of the mission in order to study the meridional distribution of acetylene (C2H2) in Jupiter's stratosphere. We find that the abundance of C2H2 decreases toward the poles by a factor of 2–4, in agreement with previous analyses of mid‐infrared spectra. This result is expected from insolation rates: near the equator, the UV solar flux is higher, allowing more C2H2 to be generated from the UV photolysis of CH4. The decrease in abundance toward the poles suggests that horizontal mixing rates are not rapid enough to homogenize the latitudinal distribution.

Wind- and Wave-driven Ocean Surface Boundary Layer in a Frontal Zone: Roles of Submesoscale Eddies and Ekman-Stokes Transport

AbstractLarge-eddy simulations are used to investigate the influence of a horizontal frontal zone, represented by a stationary uniform background horizontal temperature gradient, on the wind- and wave-driven ocean surface boundary layers. In a frontal zone, the temperature structure, the ageostrophic mean horizontal current, and the turbulence in the ocean surface boundary layer all change with the relative angle among the wind and the front. The net heating and cooling of the boundary layer could be explained by the depth-integrated horizontal advective buoyancy flux, called the Ekman Buoyancy Flux (or the Ekman-Stokes Buoyancy Flux if wave effects are included). However, the detailed temperature profiles are also modulated by the depth-dependent advective buoyancy flux and submesoscale eddies. The surface current is deflected less (more) to the right of the wind and wave when the depth-integrated advective buoyancy flux cools (warms) the ocean surface boundary layer. Horizontal mixing is greatly enhanced by submesoscale eddies. The eddy-induced horizontal mixing is anisotropic and is stronger to the right of the wind direction. Vertical turbulent mixing depends on the superposition of the geostrophic and ageostrophic current, the depth-dependent advective buoyancy flux, and submesoscale eddies.

Composited structure of non‐precipitating shallow cumulus clouds

AbstractThe normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large‐eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water () decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (), and potential temperature () perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power‐law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud‐top downdrafts cannot be captured well.

The Changing Face of Winter: Lessons and Questions From the Laurentian Great Lakes

AbstractAmong its many impacts, climate warming is leading to increasing winter air temperatures, decreasing ice cover extent, and changing winter precipitation patterns over the Laurentian Great Lakes and their watershed. Understanding and predicting the consequences of these changes is impeded by a shortage of winter‐period studies on most aspects of Great Lake limnology. In this review, we summarize what is known about the Great Lakes during their 3–6 months of winter and identify key open questions about the physics, chemistry, and biology of the Laurentian Great Lakes and other large, seasonally frozen lakes. Existing studies show that winter conditions have important effects on physical, biogeochemical, and biological processes, not only during winter but in subsequent seasons as well. Ice cover, the extent of which fluctuates dramatically among years and the five lakes, emerges as a key variable that controls many aspects of the functioning of the Great Lakes ecosystem. Studies on the properties and formation of Great Lakes ice, its effect on vertical and horizontal mixing, light conditions, and biota, along with winter measurements of fundamental state and rate parameters in the lakes and their watersheds are needed to close the winter knowledge gap. Overcoming the formidable logistical challenges of winter research on these large and dynamic ecosystems may require investment in new, specialized research infrastructure. Perhaps more importantly, it will demand broader recognition of the value of such work and collaboration between physicists, geochemists, and biologists working on the world's seasonally freezing lakes and seas.

Eddy Induced Trapping and Homogenization of Freshwater in the Bay of Bengal

AbstractFreshwater from rivers influences Indian summer monsoon rainfall and tropical cyclones by stratifying the upper layer and warming the subsurface ocean in the Bay of Bengal (BoB). Here, we use in situ and satellite data with reanalysis to showcase how river water experiences a significant increase in salinity on sub‐seasonal timescales. This involves the trapping and homogenization of freshwater by a cyclonic eddy in the Bay. Using a specific example from 2015, river water is shown to enter an eddy along its attracting manifolds within a period of 2 weeks. This leads to the formation of a highly stratified subsurface layer within the eddy. When freshest, the eddy has the largest sea‐level anomaly, spins fastest, and supports strong lateral gradients in salinity. Subsequently, observations reveal a progressive increase in salinity inside the eddy within a month. In particular, salty water spirals in, and freshwater is pulled out across the eddy boundary. Lagrangian experiments elucidate this process, whereby horizontal chaotic mixing provides a mechanism for the rapid increase in surface salinity. A salinity budget also suggests that horizontal advection explains much of the change in mixed layer salinity. Further, the adjustment of this freshwater eddy triggers submesoscale dynamics which appear to be an integral part of the process of salinity homogenization. This pathway is distinct from vertical diffusive mixing and is likely to be important for the evolution of salinity in the BoB.

The Mixed Layer Salinity Budget in the Central Equatorial Indian Ocean

AbstractThe oceanic surface mixed layer salinity (MLS) budget of the central and eastern equatorial Indian Ocean during boreal fall and winter is studied using in situ and remote sensing measurements. Budgets on roughly 100 km scale were constructed using data from two DYNAmics of the Madden–Julian Oscillation and two Research Moored Array for African‐Asian‐Australian Monsoon Analysis and Prediction moorings near 79°E during September 2011 to January 2012. The horizontal advective salinity flux plays a significant role in the seasonal variation of equatorial MLS. In boreal fall, the equatorial and 1.5°S MLS increases due to horizontal advection and turbulent mixing, despite the freshening surface flux associated with MJOs. In boreal winter, with larger sub‐monthly variation and uncertainties, the decreasing of equatorial MLS is accounted by freshening zonal advection and surface flux, abated by salty meridional advection; the 1.5°S MLS is explained by the combination of freshening meridional advection and surface flux, and salty zonal advection. Budgets between 2011 and 2015 are investigated using data products from Tropical Rainfall Measuring Mission, Aquarius, Ocean Surface Current Analyses Real‐time, Objectively Analyzed air‐sea Fluxes, and Argo mixed layers over a wider region. The eastward development of the equatorial salinity tongue in the central to eastern Indian Ocean in boreal fall and the westward retreat in boreal winter are largely determined by the equatorial zonal current. The meridional migration of ITCZ rainfall plays a secondary role. In order to improve model prediction skills of MLS changes in the equatorial Indian Ocean, both zonal and meridional salinity advective fluxes, at a spatial scale of 1° longitude and latitude and a timescale less than days, need to be properly simulated.

Grid of pseudo-2D chemistry models for tidally locked exoplanets – I. The role of vertical and horizontal mixing

ABSTRACT The atmospheres of synchronously rotating exoplanets are intrinsically 3D, and fast vertical and horizontal winds are expected to mix the atmosphere, driving the chemical composition out of equilibrium. Due to the longer computation times associated with multidimensional forward models, horizontal mixing has only been investigated for a few case studies. In this paper, we aim to generalize the impact of horizontal and vertical mixing on the chemistry of exoplanet atmospheres over a large parameter space. We do this by applying a sequence of post-processed forward models for a large grid of synchronously rotating gaseous exoplanets, where we vary the effective temperature (between 400 and 2600 K), surface gravity, and rotation rate. We find that there is a dichotomy in the horizontal homogeneity of the chemical abundances. Planets with effective temperatures below 1400 K tend to have horizontally homogeneous, vertically quenched chemical compositions, while planets hotter than 1400 K exhibit large compositional day-night differences for molecules such as CH4. Furthermore, we find that the planet’s rotation rate impacts the planetary climate, and thus also the molecular abundances and transmission spectrum. By employing a hierarchical modelling approach, we assess the relative importance of disequilibrium chemistry on the exoplanet transmission spectrum, and conclude that the temperature has the most profound impact. Temperature differences are also the main cause of limb asymmetries, which we estimate could be observable with the James Webb Space Telescope. This work highlights the value of applying a consistent modelling setup to a broad parameter space in exploratory theoretical research.

Modulation of Cross-Isothermal Velocities with ENSO in the Tropical Pacific Cold Tongue

AbstractThe tropical Pacific Ocean cold tongue (CT) plays a major role in the global climate system. The strength of the CT sets the zonal temperature gradient in the Pacific that couples with the atmospheric Walker circulation. This coupling is an essential component of the El Niño–Southern Oscillation (ENSO). The CT is supplied with cold water by the Equatorial Undercurrent that follows the thermocline as it shoals toward the east, adiabatically transporting cold water toward the surface. As the thermocline shoals, its water is transformed through diabatic processes, producing water mass transformation (WMT) that allows water to cross mean isotherms. Here, we examine WMT in the cold-tongue region from a global high-resolution ocean simulation with saved budget terms that close its heat budget exactly. Using the terms of the heat budget, we quantify each individual component of WMT (vertical mixing, horizontal mixing, eddy fluxes, and solar penetration) and find that vertical mixing is the single most important contribution in the thermocline and solar heating dominates close to the surface. Horizontal diffusion is much smaller. During El Niño events, vertical mixing, and hence cross-isothermal flow as a whole, are much reduced, whereas, during La Niña periods, strong vertical mixing leads to strong WMT, thereby cooling the surface. This analysis demonstrates the enhancement of diabatic processes during cold events, which in turn enhances cooling of the CT from below the surface.

Numerical modelling of alongshore variability in waves and wave-driven currents during the morphodynamic change of a laboratory beach

A non-hydrostatic phase-resolving numerical model (SWASH) is used to investigate the wave-driven nearshore hydrodynamics corresponding to new observations from a three-dimensional physical model in a wave basin. In the physical model, waves from a passing storm are simulated by varying the energy and duration of waves generated by a paddle over time and the beach morphology is surveyed in high spatial-resolution using a 3D laser scanner after each segment of the storm. An oblique incident wave angle introduces spatially variability in wave energy along the beach and results in morphologic differences during the storm, with active transport of beach sand by the wave-driven circulation. SWASH is validated using wave and current observations at each stage of the storm, indicating low root-mean-squared errors between model results and wave measurements, where the high-resolution model is able to predict wave shoaling, breaking and generation of the alongshore current. Model-data agreement is not as high for the velocity, suggesting spatial variability in bottom friction and horizontal mixing, and evolution of the bed in the physical model, could explain the differences. The detailed numerical model results indicate that the alongshore current is generated by wave breaking in the nearshore zone, and it increases in speed across a wider surf zone for increasing wave energy during the storm. This is associated with a bottom shear stress that exceeds the critical bed shear stress, with the potential to initiate motion and transport sediment in agreement with the experimental observations.

Deep Eastern Boundary Currents: Idealized Models and Dynamics

AbstractConcentrated poleward flows along eastern boundaries between 2- and 4-km depth in the southeast Pacific, Atlantic, and Indian Oceans have been observed, and appear in data assimilation products and regional model simulations at sufficiently high horizontal resolution, but their dynamics are still not well understood. We study the local dynamics of these deep eastern boundary currents (DEBCs) using idealized GCM simulations, and we use a conceptual vorticity model for the DEBCs to gain additional insights into the dynamics. Over most of the zonal width of the DEBCs, the vorticity balance is between meridional advection of planetary vorticity and vortex stretching, which is an interior-like vorticity balance. Over a thinner layer very close to the eastern boundary, a balance between vorticity tendencies due to friction and stretching that rapidly decay away from the boundary is found. Over the part of the DEBC that is governed by an interior-like vorticity balance, vertical stretching is driven by both the topography and temperature diffusion, while in the thinner boundary layer, it is driven instead by parameterized horizontal temperature mixing. The topographic driving acts via a cross-isobath flow that leads to stretching and thus to vorticity forcing for the concentrated DEBCs.

A New Horizontal Mixing-Length Formulation for Numerical Simulations of Tropical Cyclones

AbstractA new physically based horizontal mixing-length formulation is introduced and evaluated in the Hurricane Weather and Research Forecasting (HWRF) Model. Recent studies have shown that the structure and intensity of tropical cyclones (TCs) simulated by numerical models are sensitive to horizontal mixing length in the parameterization of horizontal diffusion. Currently, many numerical models including the operational HWRF Model formulate the horizontal mixing length as a fixed fraction of grid spacing or a constant value, which is not realistic. To improve the representation of the horizontal diffusion process, the new formulation relates the horizontal mixing length to local wind and its horizontal gradients. The resulting horizontal mixing length and diffusivity are much closer to those derived from field measurements. To understand the impact of different mixing-length formulations, we analyze the evolutions of an idealized TC simulated by the HWRF Model with the new formulation and with the current formulation (i.e., constant values) of horizontal mixing length. In two real-case tests, the HWRF Model with the new formulation produces the intensity and track forecasts of Hurricanes Harvey (2017) and Lane (2018) that are much closer to observations. Retrospective runs of hundreds of forecast cycles of multiple hurricanes show that the mean errors in intensity and track simulated by HWRF with the new formulation can be reduced approximately by 10%.

Direct Estimates of Turbulent Mixing in the Indonesian Archipelago and Its Role in the Transformation of the Indonesian Throughflow Waters

AbstractTurbulent mixing enhanced in the Indonesian Archipelago (IA) is thought to influence the coupled ocean‐climate system through sea surface temperature cooling and water‐mass transformation. The present study describes the first extensive microstructure measurements conducted in the IA in 2019. A total of 120 profiles of turbulent energy dissipation rates showed that strong turbulent mixing was extremely localized in a few narrow tidal straits, where large density overturns and high vertical diffusivity associated with downslope jets were found. Far away from mixing hotspots, in contrast, observed vertical diffusivities were similar to typical open ocean values. Microstructure measurements together with a high‐resolution numerical simulation suggest that significant water‐mass transformation occurred within very short distances, O (10 km), in these narrow tidal straits through not only vertical mixing but also horizontal mixing enhanced by submesoscale topographic eddies.

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.

Short-Term Response of Chlorophyll a Concentration Due to Intense Wind and Freshwater Peak Episodes in Estuaries: The Case of Fangar Bay (Ebro Delta)

Estuaries and coastal bays are areas of large spatio-temporal variability in physical and biological variables due to environmental factors such as local wind, light availability, freshwater inputs or tides. This study focuses on the effect of strong wind events and freshwater peaks on short-term chlorophyll a (Chl a) concentration distribution in the small-scale and microtidal, Fangar Bay (Ebro Delta, northwestern Mediterranean). The hydrodynamics of this bay are primarily driven by local wind episodes modulated by stratification in the water column. Results based on field-campaign observations and Sentinel-2 images revealed that intense wind episodes from both NW (offshore) and NE-E (onshore) caused an increase in the concentration of surface Chl a. The mechanisms responsible were horizontal mixing and the bottom resuspension (also linked to the breakage of the stratification) that presumably resuspended Chl a containing biomass (i.e., micropyhtobentos) and/or incorporated nutrients into the water column. On the other hand, sea-breeze was not capable of breaking up the stratification, so the chlorophyll a concentration did not change significantly during these episodes. It was concluded that the mixing produced by the strong winds favoured an accumulation of Chl a concentration, while the stratification that causes a positive estuarine circulation reduced this accumulation. However, the spatial-temporal variability of the Chl a concentration in small-scale estuaries and coastal bays is quite complex due to the many factors involved and deserve further intensive field campaigns and additional numerical modelling efforts.

3D1D hydro-nucleosynthesis simulations – I. Advective–reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs

ABSTRACT We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective–reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with [Fe/H] = −2.6, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the 12C(p, γ)13N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the XH profiles between the upstream and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the upstream and downstream. For example, 89Kr and 90Kr with half-lives of $3.18\,\,\mathrm{\mathrm{min}}$ and $32.3\,\,\mathrm{\mathrm{s}}$ differ by a factor 2–10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050.

The Influence of Meltwater on Phytoplankton Blooms Near the Sea‐Ice Edge

AbstractPhytoplankton blooms occur annually at the sea‐ice edge throughout the Arctic during the spring melt period. Our study considers how these spring blooms may depend on sea‐ice meltwater, focusing on the role of horizontal mixing and advection. We extend the classic Fisher reaction‐diffusion equation to consider a time‐ and space‐varying death rate that represents the role of meltwater in the system. Our results indicate that blooms peak at a characteristic distance from the ice edge where (i) meltwater is concentrated enough to stratify the upper ocean such that the phytoplankton are confined near the surface and (ii) phytoplankton have been exposed to sufficient sunlight to allow for optimized growth. The results reproduce key characteristics of a large bloom observed in Fram Strait in May 2019. Our findings support the idea that sea‐ice meltwater is of central importance in setting the spatial patterns of Arctic phytoplankton blooms.

Comparing single‐phase, non‐Newtonian approaches with experimental results: Validating flume‐scale mud and debris flow in HEC‐RAS

AbstractPost‐fire debris flows and tailing impoundment failures destroy lives and property. These geologic hazards – and other similar processes – fall on a continuum between classic Newtonian flood analyses and geotechnical stability analyses. The US Army Corps of Engineers (USACE) is developing a non‐Newtonian library (DebrisLib) that includes a suite of rheological and clastic approaches to hyper‐concentrated, mudflow and debris flow dynamics. The Hydrologic Engineering Center (HEC) has implemented these non‐Newtonian methods into the widely used, public‐domain open‐channel hydraulics and morphodynamic software, HEC‐RAS (river analysis system). This work presents part of the verification and validation of these non‐Newtonian approaches, applying several rheological equations to published laboratory results high‐concentration flume experiments.This study tested the linear Bingham model as well as the turbulent and Bagnold quadratic terms of the O'Brien equation. HEC‐RAS also includes the non‐linear Herschel–Bulkley (HB) approach, which quantifies shear‐thickening or shear‐thinning processes. The study used these non‐Newtonian models in HEC‐RAS to simulate 10 of Parsons et al.’s (2001) flume experiments, which measured the snout and plug velocity of fluids with high solid concentrations (Cv = 68–74%) and a broad range of material gradations (d50 = 0.05–1 mm, d15 = 0.006–0.1 mm). The experiments also measured and back‐calculated Bingham and HB parameters of the materials, finding HB powers between 0.45 and 1.25 (i.e. fluids that are dilatant, pseudo‐plastic and visco‐plastic).The rheological models incorporated into DebrisLib and implemented in HEC‐RAS reproduce experimental data well for most experiments. The Bingham model generated a plug velocity root‐mean‐square error (RMSE) of 0.21 m/s using standard flow parameters and Parsons et al.’s calibrated parameters, a substantial improvement over the unmodified shallow water flow equations (RMSE 0.77 m/s). Experiments with strong snout effects tended to generate higher residuals, especially in the snout velocity. The RMSE associated with the O'Brien equation was larger with the Parsons et al. fit parameters, but similar (0.23 m/s) with measured parameters. The turbulent parameter was the largest (often the dominant) parameter in most O'Brien simulations, with the dispersive stress only proving significant for the coarsest material. DebrisLib had to use a modified version of the dispersive term to simulate these concentrations. Both the 2D depth‐averaged shallow water equation (SWE) and diffusion wave equation were used to simulate the experiments. The best results were obtained with the SWE with horizontal mixing.