Scalar mixing in an urban canyon

The scalar dynamics within a unit-aspect-ratio street canyon are studied using large-eddy simulation. The key processes of ventilation and mixing are analysed with the canyon-averaged concentration, mean tracer age and variance. The results are sensitive to the source location and can be classified according to the streamline geometry. The canyon-averaged concentrations for the corner vortices, vortex sea and central vortex do not converge to the same value at large times, though the mean decay rates do. The variance measured with respect to the canyon average shows two distinct decay regimes: the early regime reflects large-scale straining and enhanced diffusion across streamlines, while the late regime is associated with escape from the canyon, i.e., ventilation. Analytical predictions for the variance-decay or mixing time scales are verified for the early regime. It is argued that the presence of an open boundary at the roof level suppresses rapid mixing of the scalar field and is responsible for differences with respect to scalar dynamics within closed domains.

     

Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence

Channel confluences at which two channels merge have an important effect on momentum exchange and contaminant diffusion in both natural rivers and artificial canals. In this study, a three-dimensional numerical model, which is based on the Reynolds Averaged Navier–Stokes equations and Reynolds Stress Turbulence model, is applied to simulate and compare flow patterns and contaminant transport processes for different bed morphologies. The results clearly show that the distribution of contaminant concentrations is mainly controlled by the shear layer and two counter-rotating helical cells, which in turn are affected by the discharge ratio and the bed morphology. As the discharge ratio increases, the shear flow moves to the outer bank and the counter-clockwise tributary helical cell caused by flow deflection is enlarged, leading the mixing happens near the outer bank and the mixing layer distorted. The bed morphology can induce shrinkage of the separation zone and increase of the clockwise main channel helical cell, which is initiated by the interaction between the tributary helical cell and the main channel flow and strengthened by the deep scour hole. The bed morphology can also affect the distortion direction of the mixing layer. Both a large discharge ratio and the bed morphology could lead to an increase in mixing intensity.     

Study on the concentration distribution in a trapezoidal open-channel flow with a side discharge

The number and distribution of pollutant concentration in a trapezoidal open channel flow with a side discharge is calculated and effects of the bank gradient are investigated in this paper. A sigma-coordinate water quality numerical model is used to simulate the process of both water and pollutant transportation in the trapezoidal channel open flow. The diffusion coefficient used in the prediction is determined by two methods including constant coefficient and the depth-averaged k-epsilon turbulence closure model. The change of the concentration with the bank gradient is acquired based on the simulation of cases with different bank gradients. An analytical formula is derived by using the mirror image method and related diffusion theories, ignoring the discharge momentum and the influence of the opposite bank. The formula can predict the number and distribution of pollutant concentration with some acceptable errors. The results demonstrate that the bank gradient has great influence on the concentration distribution which will decrease with the increase of the bank gradient approximately following a hyperbolic law.     

Mass exchange in a shallow channel flow with a series of groynes: LES study and comparison with laboratory and field experiments

The exchange of dissolved matter between a straight open channel and a series of shallow embayments present at one of its sides is investigated using large eddy simulation (LES). The direct link between the mechanism of mass exchange and the dynamics of coherent structures is demonstrated. It is shown that for the geometrical configuration considered in the present study, the mass exchange process is very non-uniform over the depth in the vicinity of the channel–embayment interface. Most of the contaminant is ejected from the embayments close to the free surface. The amount of contaminant re-entrained into the embayments situated downstream of the one in which contaminant was introduced is quantified. The mass exchange coefficient predicted by LES does not vary significantly with the embayment rank and is in very good agreement with the one predicted by the model proposed by Weitbrecht et al. (J Hydraul Eng 134(2):173–183, 2008) based on the value of a dimensionless morphometric groyne-field parameter. Field experiments were purposely performed in a natural stream with embayments whose length over width ratios were close to the ratio in the LES study. The concentration fields predicted by LES are compared with video-records of colored dye used to visualize the mass exchange in the field experiment. It is shown that, for both LES and the field experiment, the dominant passage frequency of the eddies inside the interfacial mixing layer is well predicted by the analytical model of Sukhodolov and Sukhodolova (in: Cowen et al (eds) Hydraulic measurements & experimental methods. Proceedings of international conference, Lake Placid, USA, pp 172–177, 2007). The model is then used to scale the time in the LES animations and field video-records showing the temporal evolution of the concentration field. The results of the comparison indicate several similarities in the mixing process, despite the differences in the bathymetry and the large difference in the Reynolds number between LES and the field experiment. This proves the usefulness of performing detailed LES and laboratory studies in well-controlled environments to understand mass-exchange processes around river groyne fields.     

Structure of turbulent plumes from a momentumless source in a smooth bed

The spatial development of a passive scalar plume is studied within the inhomogeneous turbulence of a boundary layer flow in a recirculating laboratory flume with a smooth bed. The source of the scalar is located flush with the bed, and the low-momentum source design is intended to simulate a diffusive-type scalar release. A weakly diffusive fluorescent dye is used as the scalar. Planar laser-induced fluorescence (PLIF) techniques were used to record the structure of the plume at a spatial resolution of 150 μm. The measured structure of the mean concentration field is compared to an analytical solution for shear-free, homogeneous turbulence. The laboratory plume exhibits spatial development in the mean concentration field that deviates from the self-similar behavior predicted by the analytical solution; this deviation is due to the mean shear and inhomogeneity of the turbulence. In particular, the influence of the viscous sublayer on the plume development is seen to be significant. Nonetheless, the analytical solution replicates some of the features seen in the laboratory plume, and the solution suggests methods of reducing the laboratory data even for cases where the results deviate from the analysis. We also examine the spatial development of the root-mean-square (rms) fluctuating concentration field, and use scalar probability density functions to examine the relationship between the mean and fluctuating concentrations.     

Mixing in the three-dimensional wake of an experimental modelled vehicle

Experimental results on tracer gas diffusion within the near wake of a simplified model car (Ahmed model with a rear slant angle of 25°) are presented. Pollutant emission is simulated using heated air injected through a small pipe at one side of the model base. Fine cold wire thermometry is used to measure instantaneous temperature excess and variance of temperature gradient in the near wake. Measurements of the three mean velocity components were made using a laser Doppler anemometers system. Characteristics of the mean and fluctuating temperature fields, time-averaged flow streamlines and scalar dissipation measurements are presented and discussed. The local mixing time scale is determined from the measured mean dissipation rate of temperature variance. Its value shows that micromixing is not a limiting phenomenon for chemical reactions in the near wake.     

Vertical dense jet in flowing current

The discharge of brackish water, as a dense jet in a natural water body, by the osmotic power plants, undergoes complex mixing processes and has significant environmental impacts. This paper focuses on the mixing processes that develop when a dense round jet outfall perpendicularly enters a shallow flowing current. Extensive experimental measurements of both the salinity and the velocity flow fields were conducted to investigate the hydrodynamic jet behavior within the ambient current. Experiments were carried out in a closed circuit flume at the Coastal Engineering Laboratory (LIC) of the Technical University of Bari (Italy). The salinity concentration and velocity fields were analyzed, providing a more thorough knowledge about the main features of the jet behavior within the ambient flow, such as the jet penetration, spreading, dilution, terminal rise height and its impact point with the flume lower boundary. In this study, special attention is given to understand and confirm the conjecture, not yet experimentally demonstrated, of the development and orientation of the jet vortex structures. Results show that the dense jet is almost characterized by two distinct phases: a rapid ascent phase and a gradually descent phase. The measured flow velocity fields definitely confirm the formation of the counter-rotating vortices pair, within the jet cross-section, during both the ascent and descent phases. Nevertheless, the experimental results show that the counter-rotating vortices pair of both phases (ascent and descent) are of opposite rotational direction.     

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

A stochastic Lagrangian micromixing model for the dispersion of reactive scalars in turbulent flows: role of concentration fluctuations and improvements to the conserved scalar theory under non-homogeneous conditions

Several reaction schemes, based on the conserved scalar theory, are implemented within a stochastic Lagrangian micromixing model to simulate the dispersion of reactive scalars in turbulent flows. In particular, the formulation of the reaction-dominated limit (RDL) reaction scheme is here extended to improve the model performance under non-homogeneous conditions (NHRDL scheme). The validation of the stochastic model is obtained by comparison with the available measurements of reactive pollutant concentrations in a grid-generated turbulent flow. This test case describes the dispersion of two atmospheric reactant species (NO and O₃) and their reaction product (NO₂) in an unbounded turbulent flow. Model inter-comparisons are also assessed, by considering the results of state-of-the-art models for pollutant dispersion. The present validation shows that RDL reaction scheme provides a systematic overestimation (relative error of ca. 85% around the centreline) in computing the local reactant consumption/production rate, whereas the NHRDL scheme drastically reduces this gap (relative error lower than 5% around the centreline). In terms of NO₂ production (or reactant consumption), neglecting concentration fluctuations determines overestimations of the product mean of around 100% and a NO₂ local production of one order of magnitude higher than the reference simulation. In terms of standard deviations, the concentration fluctuations of both the passive and reactive scalars are generally of the same order of magnitude or up to 1 or 2 orders of magnitudes higher than the corresponding ensemble mean values, except for the background reactant close to the plume edges. The study highlights the importance of modelling pollutant reactions depending on the instantaneous instead of the mean concentrations of the reactants, thus quantifying the role of the turbulent fluctuations of concentration, in terms of scalar statistics (mean, standard deviation, intensity of fluctuations, skewness and kurtosis of concentration, segregation coefficient, simulated reaction rate). This stochastic particle method represents an efficient numerical technique to solve the convection–diffusion equation for reactive scalars and involves several application fields: micro-scale air quality (urban and street-canyon scales), accidental releases, impact of odours, water quality and fluid flow industrial processes (e.g. combustion).     

Vegetation cover influence on the morphology and migration patterns of river mouths

The present study focuses on the effects of vegetation cover changes on the dynamic morphology of seven southeastern Mediterranean river mouths. The methodology used comprised monitoring and mapping by GIS techniques, with data derived from historic aerial photographs, which were applied in the investigation of the morphological spatial and temporal migration patterns of the mouths, and subsequent analysis of the vegetation cover changes influencing them. Vegetation cover adjacent to river mouths influences river mouth morphology through five primary mechanisms: a) bank vegetation; b) dune advancement toward the shoreline; c) changes in the beach??s micro-topography; as well as d) long-term continuous channel migration through permanent vegetation patches; and e) channel switching through permanent vegetation patches. The five mechanisms are part of a system of interactions between channel water flow and fluvial processes; coastal sediment transport and coastal processes; and the evolution of plant communities. In the interplay between these factors they all affect and are being affected by one another. In many river mouths artificial channel diversion is often needed due to uncontrolled channel migration. It is demonstrated that vegetation cover can serve as a mean of ??soft?? channel regulation. Therefore, a better understanding of the five influencing mechanisms may aid in controlling and managing river mouth migration patterns. The study contributes to the knowledge about bank vegetation as a tool of ??soft?? channel regulation and thus can contribute to the improvement of coastal zone management.     

Trajectory of a jet in crossflow in a channel bend

Mixing in rivers is an important issue with many applications in water quality and water resource management. Mixing of effluents with ambient river water is especially important, particularly in river bends, where secondary circulation complicates the mixing process. By comparing measured trajectories from dye tests to velocimetry data measured with an acoustic Doppler velocimeter, this paper models the trajectory of a jet in an open channel bend using a modified formula for a jet trajectory in a straight crossflow. The original formula is shown to be insufficient for modeling the trajectory in the bend. Modifications are proposed using the position of the centre of the main secondary circulation cell to account for the bend effects. In the absence of secondary circulation, the modified formula reduces to the original formula. Once the secondary circulation has developed, the proposed formula is shown to have better residuals, lower root mean squared error, and higher \(R^2\) than the original formula.     

The Information Content of a Scalar Plume – A Plume Tracing Perspective

The ability of many animals and insects to track a plume to its source is a particularly impressive feat when the fluid dynamics is considered. Inspired by this observation this research seeks to identify the information in a passive scalar plume suitable for developing robust and efficient plume tracing algorithms. The subject of this study is a scalar plume emanating from a point source in a turbulent boundary layer which has been modeled in a laboratory facility built specifically for this purpose. A coupled PIV-LIF technique is used to measure the velocity and scalar field in a time resolved fashion. This data set is analyzed and the convergence rates of five single-point statistics, suitable as kernels of plume tracing algorithms, are investigated. The experimental data shows that the scalar fluctuations over long downstream distances from the source are characterized by filamentary structures that lead to relatively slow convergence rates for any statistic that is based on mean concentrations. The scalar intermittency, however, converges rapidly toward its true value, in fact converging to a testable hypothesis for source location direction faster than the time scale of the larger scale plume meander.     

Modelling of concentrations along a moving observer in an inhomogeneous plume. Biological application: model of odour-mediated insect flights

We develop a stochastic model for the time-evolution of scalar concentrations and temporal gradients in concentration experienced by observers moving within inhomogeneous plumes that are dispersing within turbulent flows. In this model, scalar concentrations and their gradients evolve jointly as a Markovian process. Underlying the model formulation is a natural generalisation of Thomson’s well mixed condition [Thomson DJ (1987) J Fluid Mech 180:529–556]. As a consequence model outputs are necessarily compatible with statistical properties of scalars observed in experiment that are used here as model input. We then use the model to examine how insects aloft within the atmospheric boundary-layer can locate odour sources by modulating their flight patterns in response to odour cues. Mechanisms underlying odour-mediated flights have been studied extensively at laboratory-scale but an understanding of these flights over landscape scales is still lacking. Insect flights are simulated by combining the stochastic model with a simple model of insect olfactory response. These simulations show the strong influence of wind speed on the distributions of the times taken by insects to locate the source. In accordance with experimental observations [Baker TC, Vickers NJ (1997) In: Insect pheromone research: new directions, pp 248–264; Mafra-Neto A, Cardé RT (1994) Nature 369:142–144], flight patterns are predicted to become straighter and shorter, and source location is predicted to become more likely as the mean wind speed increases. The most probable arrival time to the source decreases with the mean wind speed. It is shown that scale-free movement patterns arising from olfactory-driven foraging stem directly from the power-law distribution of concentration excursion times above/below a threshold level and are robust with respect to variations in Reynolds number. Flight lengths are well represented by a power law distribution in agreement with the observed patterns of foraging bumblebees [Heinrich B (1979) Oecologia 40(3):235–245].     

Instantaneous transport of a passive scalar in a turbulent separated flow

The results of large-eddy simulations of flow and transient solute transport over a backward facing step and through a 180° bend are presented. The simulations are validated successfully in terms of hydrodynamics and tracer transport with experimental velocity data and measured residence time distribution curves confirming the accuracy of the method. The hydrodynamics are characterised by flow separation and subsequent recirculation in vertical and horizontal directions and the solute dispersion process is a direct response to the significant unsteadiness and turbulence in the flow. The turbulence in the system is analysed and quantified in terms of power density spectra and covariance of velocity fluctuations. The injection of an instantaneous passive tracer and its dispersion through the system is simulated. Large-eddy simulations enable the resolution of the instantaneous flow field and it is demonstrated that the instabilities of intermittent large-scale structures play a distinguished role in the solute transport. The advection and diffusion of the scalar is governed by the severe unsteadiness of the flow and this is visualised and quantified. The analysis of the scalar mass transport budget quantifies the mechanisms controlling the turbulent mixing and reveals that the mass flux is dominated by advection.     

Influence of planform geometry and momentum ratio on thermal mixing at a stream confluence with a concordant bed

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.     

Prediction of the upper tail of concentration distributions of a continuous point source release in urban environments

The peak values observed in a measured concentration time series of a dispersing gaseous pollutant released continuously from a point source in urban environments, and the hazard level associated with them, demonstrate the necessity of predicting the upper tail of concentration distributions. For the prediction of concentration distributions statistical models are preferably employed which provide information about the probability of occurrence. In this paper a concentration database pertaining to a field experiment is used for the selection of the statistical distribution. The inverses of the gamma cumulative distribution function (cdf) for 75th–99th percentiles of concentration are found to be more consistent with the experimental data than those of the log-normal distribution. The experimental values have been derived from measured high frequency time series by sorting first the concentrations and then finding the concentration which corresponds to each probability. Then the concentration mean and variance that are predicted with Computational Fluid Dynamics-Reynolds Averaged Navier–Stokes (RANS) methodology are used to construct the gamma distribution. The proposed model (“RANS-gamma”) is included in the framework of a computational code (ADREA-HF) suitable for simulating the dispersion of airborne pollutants over complex geometries. The methodology is validated by comparing the inverses of the model cdfs with the observed ones from two wind tunnel experiments. The evaluation is performed in the form of validation metrics such as the fractional bias, the normalized mean square error and the factor-of-two percentage. From the above comparisons it is concluded that the overall model performance for the present cases is satisfactory.     

A σ-coordinate transport model coupled with rotational boussinesq-type equations

This paper describes a σ-coordinate scalar transport model coupled with a Boussinesq-type hydrodynamic model. The Boussinesq model has the ability to calculate both three-dimensional velocity distributions and the water surface motion. To capture ‘dispersion’ processes in open channel flow, horizontal vorticity effects induced by a bottom shear stress are included in the Boussinesq model. Thus, a reasonable representation of vertical flow structure can be captured in shallow and wavy flow fields. To solve the coupled Boussinesq and scalar transport system, a finite-volume method, based on a Godunov-type scheme with the HLL Riemann solver, is employed. Basic advection and advection–diffusion numerical tests in a non-rectangular domain were carried out and the computed results show good agreement with analytic solutions. With quantitative comparisons of dispersion experiments in an open channel, it is verified that the proposed coupled model is appropriate for both near and far field scalar transport predictions. From numerical simulations in the surf zone, physically reasonable results showing expected vertical variation are obtained.     

Mixing layer in open-channel junction flows

When two open-channel flows merge in a three-branch subcritical junction, a mixing layer appears at the interface between the two inflows. If the width of the downstream channel is equal to the width of each inlet channel, this mixing layer is accelerated and is curved due to the junction geometry. The present work is dedicated to simplified geometries, considering a flat bed and a \(90^{\circ }\) angle where two configurations with different momentum ratios are tested. Due to the complex flow pattern in the junction, the so-called Serret–Frenet frame-axis based on the local direction of the velocity must be employed to characterize the flow pattern and the mixing layer as Cartesian and cylindrical frame-axes are not adapted. The analysis reveals that the centerline of the mixing layer, defined as the location of maximum Reynolds stress and velocity gradient, fairly fits the streamline separating at the upstream corner, even though a slight shift of the mixing layer towards the center of curvature is observed. The shape of the mixing layer appears to be strongly affected by the streamwise acceleration and the complex lateral confinement due to the side walls and the corners of the junction, leading to a streamwise increase of the mean velocity along the centerline and a decrease of the velocity difference. This results in a specific streamwise evolution of the mixing layer width, which reaches a plateau in the downstream region of the junction. Finally, the evaluation of the terms in the Reynolds-Averaged-Navier–Stokes equations reveals that the streamwise and normal acceleration and the pressure gradient remain dominant, which is typical of accelerated and rotational flows.     

Ants can sort their brood without a gaseous template

A fundamental issue of collective intelligence is whether the collective pattern or process is based on environmental information that explicitly codes for it or arises through self-organization of the individuals. Sometimes, these alternatives occur together. Adaptive systems may also be capable of utilizing different types of mechanism under different conditions. Sendova-Franks et al. (Anim Behav 68:1095–1106, 2004) demonstrated evidence for a self-organization mechanism of brood sorting in the ant Temnothorax albipennis, where the brood are sorted in a series of bands or concentric annuli that increase in size with distance from the colony centre. The work by Cox and Blanchard (J Theor Biol 204:223-238, 2000) suggests an alternative or complementary mechanism whereby the brood pattern is specified by the template of a CO₂ gradient. Here, we test for a gaseous template as a necessary condition for brood sorting. Under the experimental condition, we pumped the air out of the nest continuously to prevent the accumulation of any gaseous substances. We compared the brood pattern between the experimental and control conditions according to four characteristics: mean distance from centre, mean nearest-neighbour distance, shape and area. Under the experimental condition, the order of brood types according to the first two characteristics was the same as in the control. The area of the brood pattern was smaller, and its shape elongated under the experimental condition. As expected on the basis of these differences, mean distance from centre was greater and mean nearest-neighbour distance was smaller under the experimental condition (although not statistically significantly) and by the expected amount. We found evidence that ants avoid placing brood in the strongest airflow stream. This could explain the reduced area and elongated shape of the brood pattern under the experimental condition. We conclude that a gaseous template is not a necessary condition for brood sorting. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Estimation method to identify scalar point source in turbulent flow based on Taylor’s diffusion theory

We introduce a new approach to diffusion-source estimation for quick identification of the unknown source, based on Taylor’s diffusion theory for turbulent transport of passive scalar from a fixed point source. In order to evaluate the method, we used planar laser-induced fluorescence to measure the concentration field of fluorescent dye in water flowing in a channel. We considered two kinds of datasets: basis data and observed data. The former is used to determine the basis functions characterizing the streamwise dependence of variances for three statistics: the mean concentration, root-mean-square (RMS) of fluctuations in the concentration, and RMS of the temporal gradient of the fluctuating concentration. Consistent with Taylor’s theory, we found that the lateral distribution of each statistic was basically Gaussian, and their standard deviations increased as a function of the square root of the distance from the emitted point. Based on these facts, a basis function can be formulated and expected to be valid for estimation of unknown sources. Source estimation was performed with the observed data, which corresponded to limited available information about the concentration from an unknown point source. We confirmed a good prediction accuracy of the proposed method with an averaged bias as small as the turbulent integral scale. Better precision was achieved by employing several statistics simultaneously. In this case, the standard deviation of the estimated source position was assessed at 14 % of the mean distance between the source and measurement points, after 100 source-estimate trials with different datasets. The methodology tested in this paper is expected to be applicable more general and complex environmental diffusion issues involving anisotropic turbulent dispersion, and space–time variable mainstream systems; but its versatility in such systems is currently under investigation.