Estimated lead (Pb) exposures for a population of urban community gardeners

Urban community gardens provide affordable, locally grown, healthy foods and many other benefits. However, urban garden soils can contain lead (Pb) that may pose risks to human health. To help evaluate these risks, we measured Pb concentrations in soil, vegetables, and chicken eggs from New York City community gardens, and we asked gardeners about vegetable consumption and time spent in the garden. We then estimated Pb intakes deterministically and probabilistically for adult gardeners, children who spend time in the garden, and adult (non-gardener) household members. Most central tendency Pb intakes were below provisional total tolerable intake (PTTI) levels. High contact intakes generally exceeded PTTIs. Probabilistic estimates showed approximately 40 % of children and 10 % of gardeners exceeding PTTIs. Children’s exposure came primarily from dust ingestion and exposure to higher Pb soil between beds. Gardeners’ Pb intakes were comparable to children’s (in µg/day) but were dominated by vegetable consumption. Adult household members ate less garden-grown produce than gardeners and had the lowest Pb intakes. Our results suggest that healthy gardening practices to reduce Pb exposure in urban community gardens should focus on encouraging cultivation of lower Pb vegetables (i.e., fruits) for adult gardeners and on covering higher Pb non-bed soils accessible to young children. However, the common practice of replacement of root-zone bed soil with clean soil (e.g., in raised beds) has many benefits and should also continue to be encouraged.     

A numerical model of periodically forced circulation near the shore of a lake

Numerical calculations for a model of the near-shore circulation in a lake subject to two diurnal forcing mechanisms are presented. The first mechanism is a heating/cooling term in the heat equation representing the daytime heating and nighttime cooling of the diurnal cycle. The second is a periodic surface stress modelling a sea-breeze/gully wind system typical of some coastal regions. The two forcing mechanisms can either act together or against each other depending on their relative phase. The numerical solutions are compared with previously published analytical solutions and used to explore the extra dynamics associated with non-linear effects (specifically advection). The latter dynamics include the formation of gravity currents and unstable density profiles leading to secondary circulation.     

Development and testing of a micro wind tunnel for on-site wind erosion simulations

Wind erosion processes affect soil surfaces across all land uses worldwide. Understanding the spatial and temporal scales of wind erosion is a challenging undertaking because these processes are diverse and highly variable. Wind tunnels provide a useful tool as they can be used to simulate erosion at small spatial scales. Portable wind tunnels are particularly valued because erosion can be simulated on undisturbed soil surfaces in the field. There has been a long history of use of large portable wind tunnels, with consensus that these wind erosion simulation tools can meet real world aerodynamic criteria. However, one consequence of striving to meet aerodynamic reality is that the size of the tunnels has increased, making them logistically difficult to work with in the field and resulting in a tendency to homogenise naturally complex soil surfaces. This homogenisation is at odds with an increasing awareness of the importance that small scale processes have in wind erosion. To address these logistical and surface homogenisation issues we present here the development and testing of a micro wind tunnel (MWT) designed to simulate wind erosion processes at high spatial resolution. The MWT is a duct-type design—0.05 m tall 0.1 m wide and with a 1.0 m working section. The tunnel uses a centrifugal motor to suck air through a flow‐conditioning section, over the working section and then through a sediment collection trap. Simulated wind velocities range from 5 to 18 m s^?1, with high reproducibility. Wind speeds are laterally uniform and values of u _* at the tunnel bed (calculated by measuring the pressure gradients within the MWT) are comparable with those of larger tunnels in which logarithmic profiles can be developed. Saltation sediment can be added. The tunnel can be deployed by a single person and operated on slopes ranging from 0 to 10°. Evidence is presented here that the MWT provides new and useful understanding of the erodibility of rangelands, claypans and ore stockpiles.     

Modeling the interaction between tides and storm surges for the Taiwan coast

An unstructured grid, two-dimensional hydrodynamic model was established and applied to the coast of Taiwan to investigate the tide-surge interaction. Tidal elevations at the open boundaries coupled with a global ocean tidal model and the meteorological conditions using a cyclone model are used to drive the model. The model was calibrated and verified with the observed tidal levels at six tidal stations for seven typhoon events to ascertain the capability and feasibility of the model. The results show reasonable agreement between the simulated and observed tidal levels. The validated model was then applied to probe the influence of tide-surge interaction on phase, water levels, and storm surge height. We found that the tide-surge interaction influenced both the magnitude and timing of the surge, which depended on the typhoon path. The storm surge heights at different tidal stations were significantly influenced by wind stresses and directions. The water level rise due to the storm surge during high tide was greater at neap tide than at spring tide. Changing tidal ranges altered the prediction of the surge enough to induce the changes in peak water levels.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Modeling and experiments of polydisperse particle clouds

A model for polydisperse particle clouds has been developed in this study. We extended the monodisperse particle cloud model of Lai et al. (Environ Fluid Mech 13(5):435–463, 2013) to the case of polydisperse particles. The particle cloud is first considered to be a thermal or buoyant vortex ring, with the thermal induced velocity field modeled by an expanding spherical Hill’s vortex. The buoyancy of the composite thermal is assumed to be the sum of buoyancy contributed by the all particles inside the thermal. Individual particles (of different particle properties) in the cloud are then tracked by the particle tracking equation using the computed induced velocity field. The turbulent dispersion effect is also accounted for by using a random walk model. Experiments of polydisperse particle clouds were carried out to validate the model. The agreement between model predictions and experiments was reasonable. We further validate our model by comparing it with the LES study of Wang et al. (J Hydraul Eng ASCE 141(7):06015006, 2014). The limitations of our model are then discussed with reference to the comparison. Overall, although some flow details are not captured by our model, the simplicity and generality of the model makes it useful in engineering applications.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

Diffusional mass transfer coefficient at the water–sediment interface for wind-induced flow in very shallow lagoons

Very shallow lagoons that are a few centimeters deep are common in the arid Andes of Northern Chile, Argentina, Bolivia and Perú. The dynamics of these lagoons are dominated by the water–sediment interface (WSI) and strong afternoon winds. Although many studies have examined the diffusional mass transfer coefficients (k _t ) of open channel flows, estimates for wind-induced flows are still unknown. The aim of this article is to propose and validate an analytical expression for computing k _t at the WSI for wind-induced flow. The laboratory measurements were conducted in a wind tunnel with a water tank of variable depth located at its downwind end. Natural muddy sediments were placed in the middle of the tank so that the dissolved oxygen (DO) was consumed in the sediments. The diffusional mass transfer coefficient that characterizes the DO uptake in the sediment was obtained from DO micro-profiles measured with an OX-25 Unisense microelectrode. Water velocity profiles were measured with a 2D side-view Sontek acoustic doppler velocimetry (ADV), and the wind shear velocity was computed based on wind velocity profiles that were measured with an Extech hot-wire anemometer. A total of 16 experiments were conducted with different water depths and wind shear stresses. The constants required by the model were determined from these experiments, and the analytical expression was successfully validated by the laboratory observations. The analytical expression obtained for computing k _t was also validated with field observations that were conducted in October, 2012, in Salar del Huasco, Northern Chile (20.274° S, 68.883° W, 3800 m above sea level). The comparison between the observed and predicted values of k _t provides a determination coefficient of r ² = 0.48 and a p value < 0.01. The results show that the value of k _t for wind-induced flow is proportional to the wind shear velocity and the inverse of the Reynolds number of the wind-induced current.     

A numerical study of hypothetical storm surge and coastal inundation for AILA cyclone in the Bay of Bengal

The head Bay region bordering the Bay of Bengal is highly vulnerable to tropical cyclones. Catastrophic risks from storm surge and associated inundation are quite high due to high population density in coastal areas, socio-economic conditions, and shallow bathymetry. It features the world’s largest deltaic system comprising of ‘Sunderbans’ bordered by West Bengal and Bangladesh. In a geomorphologic sense, the head Bay region is a low-lying belt comprising several barrier islands and river drainage systems, numerous tidal creeks, and mud flats having a high risk for widespread inundation. In addition, the high tidal range together with low-lying topography leads to high risk and vulnerability from storm surge inundation. During May 2009, a severe cyclonic storm Aila struck West Bengal causing enormous destruction to life and property along coastal belts of West Bengal and Bangladesh. It was the strongest pre-monsoon cyclone in the past two decades that had landfall in West Bengal. This work reports on a numerical study for hypothetical storm surge and associated inundation from Aila using the ADCIRC model. The study covers a comprehensive qualitative analysis on water level elevation and onshore inundation for West Bengal and Bangladesh regions. The estimated peak storm surge was about 4 m in the Sunderban region that propagated into all major riverine systems, inundating the river banks as well the inland areas. Numerical simulations indicate an average inland penetration distance of 350 m with a maximum of 600 m at various coastal locations in West Bengal and Bangladesh. The study emphasizes the need and importance of inundation modeling system required for emergency preparedness and disaster management.     

Smoothing of land use maps for trend and change detection in urbanization

Urban sprawl and its evolution over relatively short periods of time demands that we develop statistical tools to make best use of the routinely produced land use data from satellites. An efficient smoothing framework to estimate spatial patterns in binary raster maps derived from land use datasets is developed and presented in this paper. The framework is motivated by the need to model urbanization, specifically urban sprawl, and also its temporal evolution. We frame the problem as estimation of a probability of urbanization surface and use Bayesian P-splines as the tool of choice. Once such a probability map is produced, with associated uncertainty, we develop exploratory tools to identify regions of significant change across space and time. The proposal is used to study urbanisation and its development around the city of Bologna, Emilia Romagna, Italy, using land use data from the Cartography Archive of Emilia Romagna Region for the period 1976–2008.