Gravity currents with double stratification: a numerical and analytical investigation

We consider high-Reynolds-number Boussinesq gravity current and intrusion systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a rectangular lock; the height ratio of the fluids $H$ , the stratification parameter of the ambient $S$ , and the internal stratification parameter of the current, $\sigma $ , are quite general. We perform two-dimensional Navier–Stokes simulation and compare the results with those of a previously-published one-layer shallow-water model. The results provide insights into the behavior of the system and enhance the confidence in the approximate model while also revealing its limitations. The qualitative predictions of the model are confirmed, in particular: (1) there is an initial “slumping” stage of propagation with constant speed $u_N$ , after which $u_N$ decays with time; (2) for fixed $H$ and $S$ , the increase of $\sigma $ causes a slower propagation of the current; (3) for some combinations of the parameters $H,S, \sigma $ the fluid released from the lock lacks initially (or runs out quickly of) buoyancy “driving power” in the horizontal direction, and does not propagate like a gravity current. There is also a fair quantitative agreement between the predictions of the model and the simulations concerning the spread of the current.     

Gravity currents and intrusions of stratified fluids into a stratified ambient

We consider high-Reynolds-number Boussinesq gravity currents and intrusions systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a rectangular lock; the height ratio of the fluids H, and the stratification parameter of the ambient S, are quite general. We develop a one-layer shallow-water (SW) model which is an extension of previously used and tested formulations for currents and intrusions of constant density. The internal stratification enters as a new dimensionless parameter, s ? [0,1]{\sigma \in [0,1]} . Analytical results are obtained for the initial “slumping” stage during which the speed of propagation is constant, and finite-difference solutions are presented for the more general time-dependent motion. Overall, this is a versatile and robust self-contained prediction tool, which reduces smoothly to the classical case when σ = 0. We show that, in general, the speed of propagation decreases when the internal stratification becomes more pronounced (σ increases). An interesting non-expected behavior was detected: when the stratification of the ambient is weak and moderate then the height of the current decreases with σ, but the opposite occurs when the stratification of the ambient is strong (S ≈ 1, including the case of an intrusion). Moreover, when the stratification of the ambient is strong a current with internal stratification may “run out” of driving power. We also consider the Benjamin-type steady state current with internal linear stratification in a non-stratified ambient, and show that an analytical solution exists, and that the maximal thickness decreases to below half-channel depth when σ increases.     

Gravity currents with tailwaters in Boussinesq and non-Boussinesq systems: two-layer shallow-water dam-break solutions and Navier–Stokes simulations

We consider the dam-break initial stage of propagation of a gravity current of density $\rho _{c}$ released from a lock (reservoir) of height $h_0$ in a channel of height $H$ . The channel contains two-layer stratified fluid. One layer, called the “tailwater,” is of the same density as the current and is of thickness $h_T (< h_0)$ , and the other layer, called the “ambient,” is of different density $\rho _{a}$ . Both Boussinesq ( $\rho _{c}/\rho _{a}\approx 1$ ) and non-Boussinesq systems are investigated. By assuming a large Reynolds number, we can model the flow with the two-layer shallow-water approximation. Due to the presence of the tailwater, the “jump conditions” at the front of the current are different from the classical Benjamin formula, and in some circumstances (clarified in the paper) the interface of the current matches smoothly with the horizontal interface of the tailwater. Using the method of characteristics, analytical solutions are derived for various combinations of the governing parameters. To corroborate the results, two-dimensional direct numerical Navier–Stokes simulations are used, and comparisons for about 80 combinations of parameters in the Boussinesq and non-Boussinesq domains are performed. The agreement of speed and height of the current is very close. We conclude that the model yields self-contained and fairly accurate analytical solutions for the dam-break problem under consideration. The results provide reliable insights into the influence of the tailwater on the propagation of the gravity current, for both heavy-into-light and light-into-heavy motions. This is a significant extension of the classical gravity-current theory from the particular $h_T=0$ point to the $h_T > 0$ domain.     

Numerical simulations of intrusive gravity currents interacting with a bottom-mounted obstacle in a continuously stratified ambient

In this study, the flow dynamics of intrusive gravity currents past a bottom-mounted obstacle were investigated using highly resolved numerical simulations. The propagation dynamics of a classic intrusive gravity current was first simulated in order to validate the numerical model with previous laboratory experiments. A bottom-mounted obstacle with a varying non-dimensional height of \(\tilde{D}=D/H\), where D is the obstacle height and H is the total flow depth, was then added to the problem in order to study the downstream flow pattern of the intrusive gravity current. For short obstacles, the intrusion re-established itself downstream without much distortion. However, for tall obstacles, the downstream flow was found to be a joint effect of horizontal advection, overshoot-springback phenomenon, and associated Kelvin-Helmholtz instabilities. Analysis of the numerical results show that the relationship between the downstream propagation speed and the obstacle height can be subdivided into three regimes: (1) a retarding regime (\(\tilde{D}\) \(\approx \) 0–0.3) where a 30 % increase in obstacle height leads to a 20 % reduction in propagation speed, simply due to the obstacle’s retarding effect; (2) an impounding regime (\(\tilde{D}\) \(\approx \) 0.3–0.6) where the additional 30 % increase in obstacle height only leads to a further (negligible) 5 % reduction in propagation speed, due to the accelerating effect of upstream impoundment and downstream enhanced mixing; and (3) a choking regime (\(\tilde{D}\) \(\approx \) 0.6–1.0) where the propagation speed is dramatically reduced due to the dominance of the obstacle’s blocking effect. The obstacle thickness was found to be irrelevant in determining the downstream propagation speed at least for the parameter range explored in this study. The present work highlights the significance of topographic effects in stratified flows with horizontal pressure forcing.     

Thin-layer models for gravity currents in channels of general cross-section area,a review

We present a brief review of the recent investigations on gravity currents in horizontal channels with non-rectangular cross-section area (such as triangle, \(\bigvee \)-valley, circle/semi-circle, trapezoid) which occur in nature (e.g., rivers) and constructed environment (tunnels, reservoirs, canals). To be specific, we discuss the propagation of a gravity current (GC) in a horizontal channel along the horizontal coordinate x, with gravity g acting in the \(-z\) direction, and y the horizontal–lateral coordinate. The bottom and top of the channel are at \(z=0,H\). The “standard” problem is concerned with 2D flow in a channel with rectangular (or laterally unbounded) cross-section area (CSA). Recent investigations have successfully extended the standard knowledge to the channels of CSA given by the quite general \(-f_1(z)\le y \le f_2(z)\) for \(0 \le z \le H\). This includes the practical \(\bigvee \)-valley, triangle, circle/semi-circle and trapezoid; these geometries may be in “up” or “down” setting with respect to gravity, e.g., \(\bigtriangleup \) and \(\bigtriangledown \). The major objective of the extended theory is to predict the height of the interface \(z=h(x,t)\) and the velocity (averaged over the CSA) u(x, t), where t is time; the prediction includes the speed and position of the nose \(u_N(t), x_N(t)\). We show that the motion is governed by a set of simplified equations, called “model,” that provides versatile and insightful solutions and trends. The emphasis in on a high-Reynolds-number current whose motion is dominated by buoyancy–inertia balance; in particular a GC released from a lock, which also contains general effects such as front and internal jumps (shocks), and reflected bore. We discuss two-layer, one-layer, and box models; Boussinesq and non-Boussinesq systems; compositional and particle-driven cases; and the effect of stratification of the ambient fluid. The models are self-contained, and admit realistic initial and boundary conditions. The governing equations are amenable to analytical solutions in some special circumstances. Some salient features of the buoyancy-viscous regime, and the estimate for the length at which transition to this regime takes place, are also presented. Some experimental support to the theory, and open questions for further investigations, are also mentioned. The major conclusions are (1) The CSA geometry has significant influence on the motion of the GC; and (2) The new theory is a useful, very significant, extension of the standard two-dimensional GC problem. The standard current is just a particular case, \(f_{1,2} =\) constants, among many other covered by the new theory.     

On the front conditions for gravity currents in channels of general cross-section

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel with general cross-section whose width is \(f(z), 0 \le z\le H\), and the gravity acceleration g acts in \(-z\) direction. (The classical rectangular cross-section is covered by the particular case \(f(z) =\) const.) We assume a two-layer system of homogeneous fluids of constant densities \(\rho _{c}\) (current, of height \(h < H \)) and smaller \(\rho _{a}\) (ambient, filling the remaining part of the channel). We focus attention on the calculation and assessment of the nose Froude-number condition \(Fr = U/(g' h)^{1/2}\); here U is the speed of propagation of the current and \(g' = (\rho _{c}/\rho _{a}-1) g\) is the reduced gravity. We first revisit the steady-state current, and derive compact insightful expressions of Fr and energy dissipation as a function of \(\varphi \) (\(=\) area fraction occupied by the current in the cross-section). We show that the head loss \(\delta _0\) on the stagnation line is formally a degree of freedom in the determination of \(Fr(\varphi )\), and we clarify the strong connections with the head loss \(\delta \) in the ambient fluid, and with the overall rate of dissipation \(\dot{{\mathcal{D}}}\). We demonstrate that the closure \(\delta _0 = 0\) [suggested by Benjamin (J Fluid Mech 31, 209–248, 1968) for the rectangular cross-section] produces in general the smallest Fr for a given \(\varphi \); the results are valid for a significant range \([0, \varphi _{\max }]\), in which the current is dissipative, except for the point \(\varphi _{\max }\) where \(\delta = \dot{{\mathcal{D}}} = 0\). We show that imposing the closure \(\delta = \dot{{\mathcal{D}}} = 0\), which corresponds to an energy-conserving or non-dissipative current, produces in general unacceptable restrictions of the range of validity, and large values of Fr; in particular, deep currents (\(\varphi < 0.3\) say) must be excluded because they are inherently dissipative. On the other hand, the compromise closure \(\delta (\varphi ) =\delta _0(\varphi )\) produces the simple \(Fr(\varphi ) = \sqrt{2}(1 - \varphi )\) formula whose values and dissipation properties are very close, and the range of validity is identical, to these obtained with Benjamin’s closure (moreover, we show that this corresponds to circulation-conservation solutions). The results are illustrated for practical cross-section geometries (rectangle, \(\Delta \) and \(\nabla \) triangle, circle, and the general power-law \(f(z) = b z ^\alpha \) (\(b>0, \alpha \ge 0, 0< z \le H\)). Next, we investigate the connection of the steady-state results with the time-dependent current, and show that in a lock-released current the rate of dissipation of the system is equal to, or larger than, that obtained for Fr corresponding to the conditions at the nose of the current. The results and insights of this study cover a wide range of cross-section geometry and apply to both Boussinesq and non-Boussinesq systems; they reveal a remarkable robustness of Fr as a function of \(\varphi \).     

A general solution of Benjamin-type gravity current in a channel of non-rectangular cross-section

We consider the steady-state propagation of a high-Reynolds-number gravity current in a horizontal channel along the horizontal coordinate x. The bottom and top of the channel are at z =?0, H, and the cross-section is given by the quite general form ?f ₁(z) ≤?y ≤?f ₂(z) for 0 ≤?z ≤?H, where f _1,2 are piecewise continuous functions and f ₁ +?f ₂ >?0 for ${z \in(0,H)}$ . The interface of the current is horizontal, the (maximum) thickness is h, its density is ρ _c . The reduced gravity g′ =?|ρ _c /ρ _a ? 1|g (where ${- g\hat{z}}$ is the gravity acceleration and ρ _a the density of the ambient) drives the current with speed U into the stationary ambient fluid. We show that the dimensionless Fr =?U/(g′ h)^1/2, the rate of energy dissipation (scaled with the rate of pressure work), and the dimensionless head-loss Δ/h, can be expressed by compact formulas which involve three integrals over the cross-section areas of the current and ambient. By some standard manipulations these integrals are simplified into quite simple line-integrals of the shape-function of the channel, f(z) =?f ₁(z) +?f ₂(z), and of z f(z). This theory applies to Boussinesq and non-Boussinesq currents of “heavy” (bottom) and “light” (top) type. The classical results of Benjamin (J Fluid Mech 31:209–248, 1968) for a rectangular channel are fully recovered. We also recover the Fr results of Marino and Thomas (J Fluid Eng 131(5):051201, 2009) for channels of shape y =?±b z ^α (where b, α are positive constants); in addition, we consider the energy dissipation of these flows. The results provide insights into the effect of the cross-section shape on the behavior of the steady-state current, in quite general cases, for both heavy-into-light and light-into-heavy fluid systems, Boussinesq and non-Boussinesq. In particular, we show that a very deep current displays ${Fr = \sqrt{2}}$ , and is dissipative; the value of Fr and rate of dissipation (absolute value) decrease when the thickness of the current increases. However, in general, energy considerations restrict the thickness of the current by a clear-cut condition of the form h/H ≤?a _max ?< 1.     

Waves and turbulence in katabatic winds

The measurements taken during the Vertical Transport and Mixing Experiment (VTMX, October, 2000) on a northeastern slope of Salt Lake Valley, Utah, were used to calculate the statistics of velocity fluctuations in a katabatic gravity current in the absence of synoptic forcing. The data from ultrasonic anemometer-thermometers placed at elevations 4.5 and 13.9 m were used. The contributions of small-scale turbulence and waves were isolated by applying a high-pass digital (Elliptical) filter, whereupon the filtered quantities were identified as small-scale turbulence and the rest as internal gravity waves. Internal waves were found to play a role not only at canonical large gradient Richardson numbers $(\overline{\hbox {Ri}_\mathrm{g} } >1)$ , but sometimes at smaller values $(0.1 < \overline{\hbox {Ri}_\mathrm{g}}<1)$ , in contrast to typical observations in flat-terrain stable boundary layers. This may be attributed, at least partly, to (critical) internal waves on the slope, identified by Princevac et al. [1], which degenerate into turbulence and help maintain an active internal wave field. The applicability of both Monin-Obukhov (MO) similarity theory and local scaling to filtered and unfiltered data was tested by analyzing rms velocity fluctuations as a function of the stability parameter z/L, where L is the Obukhov length and z the height above the ground. For weaker stabilities, $\hbox {z/L}<1$ , the MO similarity and local scaling were valid for both filtered and unfiltered data. Conversely, when $\hbox {z/L}>1$ , the use of both scaling types is questionable, although filtered data showed a tendency to follow local scaling. A relationship between z/L and $\overline{\hbox {Ri}_\mathrm{g} }$ was identified. Eddy diffusivities of momentum $\hbox {K}_\mathrm{M}$ and heat $\hbox {K}_\mathrm{H}$ were dependent on wave activities, notably when $\overline{\hbox {Ri}_\mathrm{g} } > 1$ . The ratio $\hbox {K}_{\mathrm{H}}/\hbox {K}_{\mathrm{M}}$ dropped well below unity at high $\overline{\hbox {Ri}_\mathrm{g} }$ , in consonance with previous laboratory stratified shear layer measurements as well as other field observations.     

Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

Direct numerical simulations of boundary condition effects on the propagation of density current in wall-bounded and open channels

The propagation of density current under different boundary conditions is investigated using high resolution direct numerical simulations (DNS). A revised Kleiser and Schumann influence-matrix method is used to treat the general Robin type velocity boundary conditions and the related “tau” error corrections in the numerical simulations. Comparison of the simulation results reveals that the boundary conditions change the turbulent flow field and therefore the propagation of the front. This paper mainly focuses on the effects of boundary conditions and initial depth of the dense fluid. The differences in energy dissipation and overall front development in wall-bounded and open channels are examined. Through DNS simulations, it is evident that with the decrease of initial release depth ratio ( $D/H$ ), the effect of the top boundary becomes less important. In wall-bounded channels, there are three distinctive layers in the vertical distribution of energy dissipation corresponding to the contributions from bottom wall, interface, and top wall, respectively. In open channels, there are only two layers with the top one missing due to the shear free nature of the boundary. It is found that the energy dissipation distribution in the bottom layer is similar for cases with the same $D/H$ ratio regardless the top boundary condition. The simulation results also reveal that for low Reynolds number cases, the energy change due to concentration diffusion cannot be neglected in the energy budget. To reflect the real dynamics of density current, the dimensionless Froude number and Reynolds number should be defined using the release depth $D$ as the length scale.     

Initial mixing of inclined dense jet in perpendicular crossflow

A comprehensive experimental investigation for an inclined ( $60^{\circ }$ to vertical) dense jet in perpendicular crossflow—with a three-dimensional trajectory—is reported. The detailed tracer concentration field in the vertical cross-section of the bent-over jet is measured by the laser-induced fluorescence technique for a wide range of jet densimetric Froude number $Fr$ and ambient to jet velocity ratios $U_r$ . The jet trajectory and dilution determined from a large number of cross-sectional scalar fields are interpreted by the Lagrangian model over the entire range of jet-dominated to crossflow-dominated regimes. The mixing during the ascent phase of the dense jet resembles that of an advected jet or line puff and changes to a negatively buoyant thermal on descent. It is found that the mixing behavior is governed by a crossflow Froude number $\mathbf{F} = U_r Fr$ . For $\mathbf{F} < 0.8$ , the mixing is jet-dominated and governed by shear entrainment; significant detrainment occurs and the maximum height of rise $Z_{max}$ is under-predicted as in the case of a dense jet in stagnant fluid. While the jet trajectory in the horizontal momentum plane is well-predicted, the measurements indicate a greater rise and slower descent. For $\mathbf{F} \ge 0.8$ the dense jet becomes significantly bent-over during its ascent phase; the jet mixing is dominated by vortex entrainment. For $\mathbf{F} \ge 2$ , the detrainment ceases to have any effect on the jet behavior. The jet trajectory in both the horizontal momentum and buoyancy planes are well predicted by the model. Despite the under-prediction of terminal rise, the jet dilution at a large number of cross-sections covering the ascent and descent of the dense jet are well-predicted. Both the terminal rise and the initial dilution for the inclined jet in perpendicular crossflow are smaller than those of a corresponding vertical jet. Both the maximum terminal rise $Z_{max}$ and horizontal lateral penetration $Y_{max}$ follow a $\mathbf{F}^{-1/2}$ dependence in the crossflow-dominated regime. The initial dilution at terminal rise follows a $S \sim \mathbf{F}^{1/3}$ dependence.     

Data depth for the uniform distribution

Given a set $X$ of $k$ points and a point $z$ in the $n$ -dimensional euclidean space, the Tukey depth of $z$ with respect to $X$ , is defined as $m/k$ , where $m$ is the minimum integer such that $z$ is not in the convex hull of some set of $k-m$ points of $X$ . If $z$ belongs to the closed region $B$ delimited by an ellipsoid, define the continuous depth of $z$ with respect to $B$ as the quotient $V(z)/\text{ Vol }(B)$ , where $V(z)$ is the minimum volume of the intersection of $B$ with the halfspaces defined by any hyperplane passing through $z$ , and $\text{ Vol }(B)$ is the volume of $B$ . We consider $z$ a random variable and prove that, if $z$ is uniformly distributed in $B$ , the continuous depth of $z$ with respect to $B$ has expected value $1/2^{n+1}$ . This result implies that if $z$ and $X$ are uniformly distributed in $B$ , the expected value of Tukey depth of $z$ with respect to $X$ converges to $1/2^{n+1}$ as the number of points $k$ goes to infinity. These findings have applications in ecology, namely within the niche theory, where it is useful to explore and characterize the distribution of points inside species niche.     

Fitting a parabola to growth data of fishes and some applications to fisheries

Denoting a fish length or weight at age t by X _t , a reference age by t _m , and the corresponding fish length or weight by X _m , the relation between age and length or weight may be described by a parabola as follows: $$\left| {X_t } \right. - X_m \left| = \right.a + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ or $$X_t = A + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ where a, b and c are constants. Each of the above Eqs. describes one curve at ages older than t _m and another one at younger ages, which is made possible by means of the transformation of t to (|t-t _m |). In 2 cases out of 10, the parabola takes the form of a cubic equation. Evidence is given that, as the growth data become fewer, the better fit of the parabola or cubic equation will probably be less in comparison to the von Bertalanffy equation (1938, 1949) as developed by Beverton and Holt (1957) and the power-growth equation (Rafail, 1971), and vice versa. This growth equation is used to derive models for estimating the optimum age and yield for fish populations.     

A simple and precise method for fitting a von Bertalanffy growth curve

The parameter K of the von Bertalanffy equation, as developed by Beverton and Holt (1957), is first estimated by the relation $$\log _e \left( {dL_t /dt} \right) = A - Kt$$ where dLt/dt denotes growth increments per a unit of age, t denotes age, and A is a constant. The K estimate is used to evaluate L∞; $$L_\infty = \left( {e^K \sum\limits_2^n {L_t - \sum\limits_1^{n - 1} {L_t } } } \right)/\left( {n - 1} \right)\left( {e^K - 1} \right)$$ The L∞ estimate is used to estimate t _o, and to obtain a better estimate for K; $$\log _e \left( {1 - L_t /L_\infty } \right) = - Kt + Kt_0 $$ The K estimate may be used to obtain another estimate for L∞. Solved examples show that a single iteration is sufficient to obtain fitted equations which are, on the average, as precise as equations fitted by the least squares method shown by Tomlinson and Abramson (1961). This new method can be used, with a slight modification, for the second equation given above, if growth data have unequal age intervals. The variance of K, t _o and log _e L∞ can be estimated by applying the simple methods used in the case of straight-line relationships.     

A flexible spatio-temporal model for air pollution with spatial and spatio-temporal covariates

The development of models that provide accurate spatio-temporal predictions of ambient air pollution at small spatial scales is of great importance for the assessment of potential health effects of air pollution. Here we present a spatio-temporal framework that predicts ambient air pollution by combining data from several different monitoring networks and deterministic air pollution model(s) with geographic information system covariates. The model presented in this paper has been implemented in an R package, SpatioTemporal, available on CRAN. The model is used by the EPA funded Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air) to produce estimates of ambient air pollution; MESA Air uses the estimates to investigate the relationship between chronic exposure to air pollution and cardiovascular disease. In this paper we use the model to predict long-term average concentrations of \(\text {NO}_{x}\) in the Los Angeles area during a 10 year period. Predictions are based on measurements from the EPA Air Quality System, MESA Air specific monitoring, and output from a source dispersion model for traffic related air pollution (Caline3QHCR). Accuracy in predicting long-term average concentrations is evaluated using an elaborate cross-validation setup that accounts for a sparse spatio-temporal sampling pattern in the data, and adjusts for temporal effects. The predictive ability of the model is good with cross-validated \(R^2\) of approximately \(0.7\) at subject sites. Replacing four geographic covariate indicators of traffic density with the Caline3QHCR dispersion model output resulted in very similar prediction accuracy from a more parsimonious and more interpretable model. Adding traffic-related geographic covariates to the model that included Caline3QHCR did not further improve the prediction accuracy.     

Partitioning of \alpha and \beta diversity using hierarchical Bayesian modeling of species distribution and abundance

Diversity partitioning is becoming widely used to decompose the total number of species recorded in an area or region \((\gamma )\) into the average number of species within samples \((\alpha )\) and the average difference in species composition \((\beta )\) among samples. Single-value metrics of \(\alpha \) and \(\beta \) diversity are popular because they may be applied at multiple scales and because of their ease in computation and interpretation. Studies thus far, however, have emphasized observed diversity components or comparisons to randomized, null distributions. In addition, prediction of \(\alpha \) and \(\beta \) components using environmental or spatial variables has been limited to more extensive data sets because multiple samples are required to estimate single \(\alpha \) and \(\beta \) components. Lastly, observed diversity components do not incorporate variation in detection probabilities among species or samples. In this study, we used hierarchical Bayesian models of species abundances to provide predictions of \(\alpha \) and \(\beta \) components in species richness and composition using environmental and spatial variables. We illustrate our approach using butterfly data collected from 26 grassland remnants to predict spatially nested patterns of \(\alpha \) and \(\beta \) based on the predicted counts of butterflies. Diversity partitioning using a Bayesian hierarchical model incorporated variation in detection probabilities by butterfly species and habitat patches, and provided prediction intervals for \(\alpha \) and \(\beta \) components using environmental and spatial variables.     

Physical and biological controls of algal blooms in the Río de la Plata

Coupled three-dimensional hydrodynamic and ecological numerical simulations were used to investigate the role of transport, stagnation zones and dispersion on inter-annual blooms of the diatom Aulacoseira sp. in the vicinity of the drinking water intakes of the Buenos Aires city (Argentina) in the upper Río de la Plata. Three different summer events were analyzed. First, a mild biomass bloom year (2006–2007), second, a high biomass bloom year (2007–2008) and third, a “normal” no bloom year (2009–2010). Simulated water height, water temperature, suspended solids and chlorophyll \(a\) concentrations patterns compared well with field data. Results revealed that the advection of phytoplankton cells via inflows to the Río de la Plata triggered Aulacoseira sp. blooms in the domain. In addition, excessive growth observed near the drinking water intakes, along the Argentinean margin, were associated with long retention times (stagnant region) and weak horizontal dispersion. Increased concentrations of suspended solids in the water column, in response to re-suspension events, did not prevent the blooms, however, were found to also play a key role in controlling the rate of phytoplankton growth. Finally, a non-dimensional parameter, R, that considers phytoplankton patch size, e-folding growth and dispersion time scales is shown to determine the potential bloom occurrences, as well as bloom intensity; R values higher than 5.7 suggest intense phytoplankton growth. For the mild biomass bloom year, \(R = 7.5\) , for the high biomass bloom year, \(R = 11\) and for the “normal” no bloom year \(R= 0.4\) .     

Quantifying a significance of sediment particle size to hyporheic sedimentary oxygen demand with a permeable stream bed

A mechanistic model of sedimentary oxygen demand (SOD) for hyporheic flow is presented. The permeable sediment bed, e.g. sand or fine gravel, is considered with hydraulic conductivity in the range $0.1 < K < 20$  cm/s. Hyporheic pore water flow is induced by pressure fluctuations at the sediment/water interface due to near-bed turbulent coherent motions. A 2-D advection–diffusion equation is linked to the pore water flow model to simulate the effect of advection–dispersion driven by interstitial flow on oxygen transfer through the permeable sediment. Microbial oxygen uptake in the sediment is expressed as a function of the microbial growth rate, and is related to the sediment properties, i.e. the grain diameter $(d_{s})$ and porosity $(\phi )$ . The model describes the significance of sediment particle size to oxygen transfer through the sediment and microbial oxygen uptake: With increasing grain diameter $(d_{s})$ , the hydraulic conductivity $(K)$ increases so does the oxygen transfer rate, while particle surface area per volume (the available surface area for colonization by biofilms) decreases reducing the microbial oxygen uptake rate. Simulation results show that SOD increases as the hydraulic conductivity $(K)$ increases before a threshold has been reached. After that, SOD diminishes with the increment of the hydraulic conductivity $(K)$ .