Halton iterative partitioning: spatially balanced sampling via partitioning

A new spatially balanced sampling design for environmental surveys is introduced, called Halton iterative partitioning (HIP). The design draws sample locations that are well spread over the study area. Spatially balanced designs are known to be efficient when surveying natural resources because nearby locations tend to be similar. The HIP design uses structural properties of the Halton sequence to partition a resource into nested boxes. Sample locations are then drawn from specific boxes in the partition to ensure spatial diversity. The method is conceptually simple and computationally efficient, draws spatially balanced samples in two or more dimensions and uses standard design-based estimators. Furthermore, HIP samples have an implicit ordering that can be used to define spatially balanced over-samples. This feature is particularly useful when sampling natural resources because we can dynamically add spatially balanced units from the over-sample to the sample as non-target or inaccessible units are discovered. We use several populations to show that HIP sampling draws spatially balanced samples and gives precise estimates of population totals.     

Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

A two-dimensional inviscid model of the gravity-current head by Lagrangian block simulation

A two-dimensional inviscid model of the gravity-current head produced by the release of a relatively small volume of dense fluid from behind a tall lock gate is constructed by Lagrangian block simulation. Three numerical experiments are conducted for the lock’s height-to-length aspect ratios H/L _o  = 8, 4 and 2. The front speeds obtained by the simulations agree with the laboratory observation for a similar range of aspect ratios. The floor velocity in the wake behind these heads is found to be greater than their front speed. The high floor velocity is caused by the impingement of the coherent wake vortex on the floor. It is a condition that permits these gravity-current heads to maintain their structural integrity so that the fine sediments can travel with the head over long distances on the ocean floor. The structural coherence of the current head depends on the lock aspect ratio. The gravity-current head produced by the release from the lock with the highest aspect ratio of H/L _o  = 8 is most coherent and relatively has the greatest floor velocity and the least trailing current behind the head.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Shoaling internal waves may reduce gravity current transport

Gravity currents descending along slopes have typically been studied in quiescent environments, despite the fact that in many geophysical settings there is significant externally driven motion. Here we investigate how the head of a gravity current is influenced by interfacial internal waves at the pycnocline of a two-layer ambient water column. Our experimental measurements show that larger amplitude internal waves, interacting with the gravity current, reduce both the mass transport by the gravity current and its thickness. These results suggest that the ambient internal wave field should be considered when estimating transport by gravity currents in geophysical settings with strong internal waves, such as lakes and the coastal ocean.     

Mobility of free surface in different liquids and its influence on water striders locomotion

Dynamics of the surface layer in different liquids is examined by means of infrared thermography of the surface and simultaneous velocity fields measurements using surface and infrared Particle Image Velocimetry. This technique allows measurements and comparison of two velocity fields—at the surface and at small depth about 50–200 μm. In distilled water the velocity fields at the surface and at small depth exhibit significant dissimilarity. The flow field below the surface is essentially 3D, whereas the surface flow is characterized by vanishing 2D divergence of velocity, indicating predominantly planar motion. In contrast, in ethanol–butanol mixture two velocity fields are well correlated, both corresponding to 3D flow with continuous surface renewal. Thermal patterns, observed at the surface, and the flow field structure in different liquids are associated with different boundary conditions for velocity at the surface. Water surface is seldom renewed, which inhibits heat and mass exchange between the liquid and atmosphere. However, absence of vertical advection also enables organisms to live within the surface layer, to stand and walk on the free surface. This is illustrated by the difficulties a water strider faces on the surface of ultrapure water, which exhibits Marangoni convection.     

Interaction of North Brazil Current rings with the Lesser Antilles Arc and Barbados Island: laboratory experiments and observations

The interaction of North Brazil Current (NBC) rings with the Lesser Antilles Arc (LAA) and the Barbados Island (BI) is addressed by experimental modeling and observations. Our results compare well with previous experimental results and numerical simulations. Several sizes, intensities and two different vorticity profiles (non-isolated and initially isolated vortices) were tested. Three regimes were found namely: (1) the vortex surrounds the BI and its translational motion (TM) stops North of BI; (2) the vortex passes through the corridor between the LAA and the BI by reducing its size; and (3) The vortex stopped at the entrance to the corridor South of the BI. Isolated vortices were prone to stopped North of the BI. Apparently the intensity in the outer vorticity ring has an influence on the fate of the NBC ring. Non-isolated vortices can also stop its TM North of the BI because when in \(\beta\) plane they develop an outer ring similar to the isolated vortices. From these results we conclude that intense and big NBC rings are likely to stop its TM North of the BI. Medium and moderate vortices stops its TM South of the BI and they reduce their size until they are able to pass through the corridor between the LAA and the BI. Mild vortices of all sizes stop South of the corridor, close to the BI and the LAA. Drifter trajectories and Sea Surface Height altimetry confirm the results.     

Human health impacts of exposure to metals through extreme consumption of fish from the Colombian Caribbean Sea

The health risks of metal exposure due to the high consumption of fish were assessed for a riverine population living on the Caribbean coast of Colombia. The concentrations of metals (Cd, Cr, Cu, Hg, Ni, Pb and Zn) in the edible tissues of fish were determined and used for risk assessment. The daily fish consumption of residents (n = 95) was as high as 283, 366 and 469 g/day in children (CH), women of childbearing age (WCHA) and the remaining population groups (RP), respectively. The estimation of the potential risk (HQ) indicated that there was no health risk from most of the metals, because they did not exceed their related reference doses, with values of HQ < 1. Although the concentrations of Pb and Hg were not particularly high in fish (<0.2 µg/g), their possible health effects for vulnerable groups are of great concern due to the extremely high fish intake. The Pb intake for all groups was higher than the lower confidence limit of the benchmark dose for nephrotoxicity and neurodevelopmental effects in children. The weekly intake of methylmercury was also elevated, with values approximately 3, 2 and 1.5 times the provisional tolerable weekly intake for CH, WCHA and RP, respectively. Moreover, higher Hg levels were found in top predators, whereas maximum levels for other metals were found in bottom-feeding fish. This study highlights that an accurate data of daily intake, a continuous monitoring of metals in fish and their related fish consumption advisories to protect subsistence fishing communities are recommended in a local and worldwide context.     

Hydrogeochemistry study and groundwater quality assessment in the north of Isfahan,Iran

This study presents the groundwater quality assessment in the north of Isfahan, Iran. In the study area, assessment and measurement of groundwater hydrochemical parameters such as pH, total dissolved solids (TDS), electrical conductivity (EC), sodium absorption ratio (SAR), total hardness, major cations (K⁺, Na⁺, Ca²⁺ and Mg²⁺) and major anions (Cl^?, \({\text{HCO}}_{ 3}^{ - } ,{\text{CO}}_{3}^{2 - }\) and \({\text{SO}}_{4}^{2 - }\)) concentrations were performed. Accordingly, the 66 water samples from different locations were collected during April and May 2015. Water samples collected in the field were analyzed in the laboratory for cations and anions using the standard methods. In this research, the analytical results of physiochemical parameters of groundwater were compared with the standard guideline values as recommended by the world health organization (WHO) for drinking and public health purposes. The pH values of groundwater samples varied from 7.05 to 8.95 with a mean of 7.78, indicating a neutral to slightly alkaline water. TDS values showed that 14% of the samples exceeds the desirable limit given by WHO. EC values varied from 213 to 4320 µS/cm, while 23% of the samples were more than the standard limit. Gibbs diagram had shown that 90% of the samples in the study area fall in the rock weathering zone, and this means that chemical weathering of rock-forming minerals is the main factor controlling the water chemistry in the study area. Irrigation suitability and risk assessment of groundwater are evaluated by measuring EC, %Na, SAR and RSC. According to the dominant cations and anions, five types of water were identified in the water samples: Ca-HCO₃, Ca-SO₄, Na-Cl, Na-HCO₃ and Na-SO₄. The results show that the majority of samples (30 samples, 45%) belongs to the mixed Na-SO₄ water type. Correlation analysis and principal component analysis was used to identify the relationship between ions and physicochemical parameters. Results indicated that 18 stations of the study area had the best quality and can be used for irrigation and drinking purposes in the future.