Estimation bias under model selection for distance sampling detection functions

Many simulation studies have examined the properties of distance sampling estimators of wildlife population size. When assumptions hold, if distances are generated from a detection model and fitted using the same model, they are known to perform well. However, in practice, the true model is unknown. Therefore, standard practice includes model selection, typically using model comparison tools like Akaike Information Criterion. Here we examine the performance of standard distance sampling estimators under model selection. We compare line and point transect estimators with distances simulated from two detection functions, hazard-rate and exponential power series (EPS), over a range of sample sizes. To mimic the real-world context where the true model may not be part of the candidate set, EPS models were not included as candidates, except for the half-normal parameterization. We found median bias depended on sample size (being asymptotically unbiased) and on the form of the true detection function: negative bias (up to 15% for line transects and 30% for point transects) when the shoulder of maximum detectability was narrow, and positive bias (up to 10% for line transects and 15% for point transects) when it was wide. Generating unbiased simulations requires careful choice of detection function or very large datasets. Practitioners should collect data that result in detection functions with a shoulder similar to a half-normal and use the monotonicity constraint. Narrow-shouldered detection functions can be avoided through good field procedures and those with wide shoulder are unlikely to occur, due to heterogeneity in detectability.     

Experimental assessment of characteristic turbulent scales in two-phase flow of hydraulic jump: from bottom to free surface

A hydraulic jump is a turbulent shear flow with a free-surface roller. The turbulent flow pattern is characterised by the development of instantaneous three-dimensional turbulent structures throughout the air–water column up to the free surface. The length and time scales of the turbulent structures are key information to describe the turbulent processes, which is of significant importance for the improvement of numerical models and physical measurement techniques. However, few physical data are available so far due to the complexity of the measurement. This paper presents an investigation of a series of characteristic turbulent scales for hydraulic jumps, covering the length and time scales of turbulent flow structures in bubbly flow, on free surface and at the impingement point. The bubbly-flow turbulent scales are obtained for Fr = 7.5 with 3.4 × 10⁴ < Re < 1.4 × 10⁵ in both longitudinal and transverse directions, and are compared with the free-surface scales. The results highlight three-dimensional flow patterns with anisotropic turbulence field. The turbulent structures are observed with different length and time scales respectively in the shear flow region and free-surface recirculation region. The bubbly structures next to the roller surface and the free-surface fluctuation structures show comparable length and time scales, both larger than the scales of vortical structures in the shear flow and smaller than the scales of impingement perimeter at the jump toe. A decomposition of physical signals indicates that the large turbulent scales are related to the unsteady motion of the flow in the upper part of the roller, while the high-frequency velocity turbulence dominates in the lower part of the roller. Scale effects cannot be ignored for Reynolds number smaller than 4 × 10⁴, mainly linked to the formation of large eddies in the shear layer. The present study provides a comprehensive assessment of turbulent scales in hydraulic jump, including the analyses of first data set of longitudinal bubbly-flow integral scales and transverse jump toe perimeter integral scales.     

Modeling and data assessment of longitudinal salinity in a low-gradient estuarine river

Continuous data of vertical-profile salinity were analyzed for four stations located successively upriver in a macrotidal estuary, the lower St. Johns River (Northeast Florida, USA). The data analysis confirmed well-mixed salinity conditions in the river with at most 1.3 ppt of vertical variability at Dames Point (river km 20), where the main variations of salinity are along the longitudinal axis of the river. Given the well-mixed salinity conditions and dominant horizontal structure of salinity variations in the river, we present and apply a barotropic, two-dimensional modeling approach for hydrodynamic-salinity transport simulation in the lower St. Johns River. When properly forced by offshore surge, high-resolution wind fields and freshwater river inflows, the model replicated the salinity measurements remarkably well, including the separation into tidal and sub-tidal components. The data and model results show that, at times, offshore winds and surge can be more influential on longitudinal salinity variations than local winds over the river. We demonstrate the importance of using proper boundary conditions to force the model relative to the minimal sensitivity of the model to parameter adjustment of horizontal mixing and uncertainty-based perturbation of wind and inflow forcings.     

Direct measurements of the drag force over aligned arrays of cubes exposed to boundary-layer flows

Wind tunnel measurements of the total drag force for aligned arrays of cubes exposed to two different boundary-layer flows at three flow velocities are discussed. The drag force for eight different building packing densities λ _p (from 0.028 to 1) is measured with a standard load cell generating a novel dataset. Different λ _p are reproduced by increasing the number of buildings on the same lot area; this represents a real situation that an urban planner is faced with when a lot area of a given (fixed) size is allocated to the development of new built areas. It is assumed that the surrounding terrain is uniform and there is a transition from a given roughness (smooth) to a new roughness (rough). The approaching flow will adjust itself over the new surface within a distance that in general may be larger than the horizontal length covered by the array. We investigate the region where the flow adjustment occurs. The wide range of packing densities allowed us to analyse in detail the evolution of the drag force. The drag force increases with increasing packing densities until it reaches a maximum at an intermediate packing density (λ _p  = 0.25 in our case) followed by a slight decrease at larger packing densities. The value of the drag force depends on the flow adjustment along the array which is evaluated by introducing the parameter “drag area” to retrieve information about the drag distribution at different λ _p . Results clearly suggest a change of the distribution of the drag force, which is found to be relatively uniform at low packing densities, while most of the force acts on first rows of the arrays at large packing densities. The drag area constitutes the basis for the formulation of a new adjustment length scale defined as the ratio between the volume of the air within the array and the drag area. The proposed adjustment length scale automatically takes into account the change in drag distribution along the array for a better parameterization of urban effects in dispersion models.     

A mathematical model for type II profile of concentration distribution in turbulent flows

This paper presents a mathematical model to investigate type II profile of suspension concentration distribution (i.e., the concentration profile where the maximum concentration appears at some distance above the bed surface) in a steady, uniform turbulent flow through open-channels. Starting from the mass and momentum conservation equations of two-phase flow, a theoretical model has been derived. The distribution equation is derived considering the effects of fluid lift force, drag force, particle inertia, particle–particle interactions, particle velocity fluctuations and drift diffusion. The equation is solved numerically and is compared with available experimental data as well as with other models existing in the literature. Good agreement between the observed value and computed result, and minimum error in comparison to other models indicate that the present model can be applied in predicting particle concentration distribution for type II profile for a wide range of flow conditions. The proposed model is also able to show the transition from type I profile to type II profile.     

Comparison of analytical descriptions of nocturnal low-level jets within the Ekman model framework

This study focuses on the inertial oscillation aspect of the nocturnal low-level jet (NLLJ). In the context of the Ekman model solutions, conceptual NLLJ inertial oscillation analytical frameworks proposed by Blackadar in 1957 and Shapiro and Fedorovich and van de Wiel et al. in 2010 are compared. Considering a NLLJ produced via direct numerical simulation over flat terrain with no baroclinic influence as a reference case, the deficiencies of each framework in representing a realistic NLLJ are assessed. The Blackadar theory results in unrealistic wind profiles near the surface. While extensions of Blackadar’s framework by Shapiro and Fedorovich and van de Wiel et al. produce more realistic NLLJs, the simpler approach taken by van de Wiel et al. does not describe the NLLJ wind hodograph at later times sufficiently in qualitative terms.     

Influence of the secondary motions on pollutant mixing in a meandering open channel flow

This paper presents large eddy simulation of turbulent flow in a meandering open channel with smooth wall and rectangular cross-section. The Reynolds number based on the channel height is 40,000 and the aspect ratio of the cross-section is 4.48. The depth-averaged mean stream-wise velocity agree well to experimental measurements. In this specific case, two interacting cells are formed that swap from one bend to the other. Transport and mixing of a pollutant is analysed using three different positions of release, e.g. on the inner bank, on the outer bank and on the centre of the cross section. The obtained depth-average mean concentration profiles are reasonably consistent with available experimental data. The role of the secondary motions in the mixing processes is the main focus of the discussion. It is found that the mixing when the scalar is released on the centre of the cross-section is stronger and faster than the mixing of the scalar released on the sides. When the position of release is close to a bank side, the mixing is weaker and a clear concentration of scalar close to the corresponding side-wall can be observed in both cases.     

Arsenic and other trace elements in groundwater and human urine in Ha Nam province,the Northern Vietnam: contamination characteristics and risk assessment

The contamination characteristics of arsenic and other trace elements in groundwater and the potential risks of arsenic from the groundwater were investigated. Elevated contamination of arsenic, barium and manganese was observed in tube-well water of two villages (Chuyen Ngoai and Chau Giang) in Ha Nam province in the Northern Vietnam. Concentrations of As in the groundwater ranged from 12.8 to 884 µg/L with mean values in Chuyen Ngoai and Chau Giang were 614.7 and 160.1 µg/L, respectively. About 83 % of these samples contained As concentrations exceeding WHO drinking water guideline of 10 μg/L. The mean values of Mn and Ba in groundwater from Chuyen Ngoai and Chau Giang were 300 and 657 μg/L and 650 and 468 μg/L, respectively. The mean value of Ba concentration in groundwater in both Chuyen Ngoai and Chau Giang was about 22 % of the samples exceeded the WHO guideline (700 µg/L). Arsenic concentrations in human urine of residents from Chuyen Ngoai and Chau Giang were the range from 8.6 to 458 µg/L. The mean values of Mn and Ba in human urine of local people from Chuyen Ngoai were 46.9 and 62.8 μg/L, respectively, while those in people from Chau Giang were 25.9 and 45.9 μg/L, respectively. The average daily dose from ingesting arsenic for consuming both untreated and treated groundwater is from 0.02 to 11.5 and 0.003 to 1.6 μg/kg day, respectively. Approximately, 57 % of the families using treated groundwater and 64 % of the families using untreated groundwater could be affected by elevated arsenic exposure.     

Health risk assessment of potentially harmful elements and dietary minerals from vegetables irrigated with untreated wastewater,Pakistan

In the developing world, vegetables are commonly grown in suburban areas irrigated with untreated wastewater containing potentially harmful elements (PHEs). In Pakistan, there is no published work on the bioaccessibility aspect of PHEs and dietary minerals (DMs) in sewage-irrigated soil or the vegetables grown on such soils in Pakistan. Several industrial districts of Pakistan were selected for assessment of the risk associated with the ingestion of vegetables grown over sewage-irrigated soils. Both the total and bioaccessible fraction of PHEs (Cd, Co, Cr, Ni, and Pb) and DMs (Fe, Cu, Mn, Zn, Ca, Mg, and I) in soils and vegetable samples were measured. The concentrations of these PHEs and DMs in sewage-irrigated and control soils were below published upper threshold limits. However, compared to control soils, sewage irrigation over the years decreased soil pH (7.7 vs 8.1) and enhanced dissolved organic carbon (1.8 vs 0.8 %), which could enhance the phyto-availability of PHEs and DMs to crops. Of the PHEs and DMs, the highest transfer factor (soil to plant) was noted for Cd and Ca, respectively. Concentrations of PHEs in most of the sewage-irrigated vegetables were below the published upper threshold limits, except for Cd in the fruiting portion of eggplant and bell pepper (0.06–0.08 mg/kg Cd, dry weight) at three locations in Gujarat and Kasur districts. The bioaccessible fraction of PHEs can reduce the context of dietary intake measurements compared to total concentrations, but differences between both measurements were not significant for Cd. Since the soils of the sampled districts are not overly contaminated compared to control sites, vegetables grown over sewage-irrigated soils would provide an opportunity to harvest mineral-rich vegetables potentially providing consumers 62, 60, 12, 104, and 63 % higher dietary intake of Cu, Mn, Zn, Ca, and Mg, respectively. Based on Fe and vanadium correlations in vegetables, it is inferred that a significant proportion of total dietary Fe intake could be contributed by soil particles adhered to the consumable portion of vegetables. Faecal sterol ratios were used to identify and distinguish the source of faecal contamination in soils from Gujranwala, Gujarat, and Lahore districts, confirming the presence of human-derived sewage biomarkers at different stages of environmental alteration. A strong correlation of some metals with soil organic matter concentration was observed, but none with sewage biomarkers.     

Efficacy of woody biomass and biochar for alleviating heavy metal bioavailability in serpentine soil

Crops grown in metal-rich serpentine soils are vulnerable to phytotoxicity. In this study, Gliricidia sepium (Jacq.) biomass and woody biochar were examined as amendments on heavy metal immobilization in a serpentine soil. Woody biochar was produced by slow pyrolysis of Gliricidia sepium (Jacq.) biomass at 300 and 500 °C. A pot experiment was conducted for 6 weeks with tomato (Lycopersicon esculentum L.) at biochar application rates of 0, 22, 55 and 110 t ha^?1. The CaCl₂ and sequential extractions were adopted to assess metal bioavailability and fractionation. Six weeks after germination, plants cultivated on the control could not survive, while all the plants were grown normally on the soils amended with biochars. The most effective treatment for metal immobilization was BC500-110 as indicated by the immobilization efficiencies for Ni, Mn and Cr that were 68, 92 and 42 %, respectively, compared to the control. Biochar produced at 500 °C and at high application rates immobilized heavy metals significantly. Improvements in plant growth in biochar-amended soil were related to decreasing in metal toxicity as a consequence of metal immobilization through strong sorption due to high surface area and functional groups.