This study describes a methodology combining a recently proposed renormalization inversion technique with a building‐resolving computational fluid dynamics (CFD) approach for source retrieval in the geometrically complex urban regions. It presents the first application of the renormalization inversion approach to estimate an unknown continuous point release in real situations at an urban scale. The renormalization inversion approach is based on an adjoint source‐receptor relationship and is purely deterministic in nature. The source parameters (i.e., source location and release rate) are reconstructed from a finite set of point measurements of concentration acquired from some sensors and the adjoint functions computed from a CFD model fluidyn‐PANACHE that is able to represent the geometric and flow complexity inherent in the urban regions. The inversion procedure is evaluated for a point source reconstruction using measurements from the Mock Urban Setting Test (MUST) field experiment. Source reconstructions are performed for 20 trials of the MUST experiment of a continuous point release in an idealized urban geometry consisting of a regular array of shipping containers. The steady state flow fields are computed by solving the three‐dimensional Reynolds‐averaged Navier‐Stokes equations by using a finite volume scheme. Then, in each MUST trial adjoint functions are obtained and used for the source identification. Inversion results are presented with both synthetic and real measurements in various atmospheric stabilities varying from neutral to stable and very stable conditions. With real concentration measurements, the point source is retrieved within an average Euclidean distance of 14.6 m from the actual source location. The estimated source intensity is overpredicted by an average factor of 1.37 of the true release rate. In a posterior uncertainty analysis with 10% random noise in measurements, it is demonstrated that standard deviation in the location error and release strength, respectively, varies by 5.22 m and ∼21% from their mean value for all 20 trials. A sensitivity analysis shows that the use of nonzero measurements helps in reducing the uncertainties involved in the source reconstruction. The source reconstruction results in various stability conditions exhibit the reliability of the renormalization inversion methodology coupled with the CFD approach in an urban area. The present methodology can be used by emergency regulators as a tool to detect the unknown accidental or deliberated releases in the complex urban environments.

To obtain, over medium term periods, wind speed time series on a site, located in the southern part of the Paris region (France), where long recording are not available, but where nearby meteorological stations provide large series of data, use was made of ANN based models. The performance of these models have been evaluated by using several commonly used statistics such as average absolute error, root mean square error, normalized mean square error and correlation coefficient. Such global criteria are good indicators of the "robustness" of the models but are unable to provide useful information about their "effectiveness" in accurately generating wind speed fluctuations over a wide range of scales. Therefore a complementary wavelet cross coherence analysis has been performed. Wavelet cross coherence, wavelet cross correlation and spectral wavelet cross correlation coefficients, have been calculated and displayed as functions of the equivalent Fourier period. These coefficients provide quantitative measures of the scale-dependence of the model performance. In particular the spectral wavelet cross coherence coefficient can be used to have a rapid and efficient identification of the validity range of the models. The results show that the ANN models employed in this study are only effective in computing large scale fluctuations of large amplitude. To obtain a more representative time series, with much higher resolution, small scale fluctuations have to be simulated by a superimposed statistical model. By combining ANN and statistical models, both the high and the low frequency segments of the wind velocity spectra can be simulated, over a range of several hours, at the target site.

Reconstruction of unknown atmospheric releases using measured concentrations is an ill-posed inverse problem. Due to insufficient measurements and dispersion model uncertainties, reliable interpretation of a retrieved source is limited by lack of resolution, nonuniqueness, and instability in the inverse solution. The study presents an optimality analysis, in terms of resolution, stability, and reliability, of an inverse solution given by a recently proposed inversion technique, called as renormalization. The inversion technique is based on an adjoint source-receptor framework and construction of a weight function which interprets a priori information about the unknown release apparent to the monitoring network. The properties of weight function provide a perfect data resolution, maximum model resolution, and minimum variance (or stability) for the retrieved source. The reliability of the retrieved source is interpreted in view of the information derived from the geometry of the monitoring network. The inversion technique and resolution features are evaluated for a point source reconstruction using measurements from a recent dispersion experiment (Fusion Field Trials 2007) conducted at Dugway Proving Ground, Utah. With the real measurements, the point release is reconstructed within an average distance of 23 m from the true release where the average distance of the nearest receptor from the true source was 32 m. In all the trials, the point release is retrieved within 3-60 m Euclidean distance from their true location. The source strength is retrieved within a factor of 1.5 to the true release mass. The posterior uncertainty in the release parameters is observed to be within 20% of their mean value. The source localization features are resolved to its maximum extent feasible with the design of the monitoring network. The sensitivity studies are conducted to highlight the importance of receptors reporting zero concentration measurements and variations in the resolution features of the source retrieval with respect to the various arrangements of the receptors.

In the event of an accidental or intentional contaminant release in the atmosphere, it is imperative, for managing emergency response, to diagnose the release parameters of the source from measured data. Reconstruction of the source information exploiting measured data is called an inverse problem. To solve such a problem, several techniques are currently being developed. The first part of this paper provides a detailed description of one of them, known as the renormalization method. This technique, proposed by Issartel (2005), has been derived using an approach different from that of standard inversion methods and gives a linear solution to the continuous Source Term Estimation (STE) problem. In the second part of this paper, the discrete counterpart of this method is presented. By using matrix notation, common in data assimilation and suitable for numerical computing, it is shown that the discrete renormalized solution belongs to a family of well-known inverse solutions (minimum weighted norm solutions), which can be computed by using the concept of generalized inverse operator. It is shown that, when the weight matrix satisfies the renormalization condition, this operator satisfies the criteria used in geophysics to define good inverses. Notably, by means of the Model Resolution Matrix (MRM) formalism, we demonstrate that the renormalized solution fulfils optimal properties for the localization of single point sources. Throughout the article, the main concepts are illustrated with data from a wind tunnel experiment conducted at the Environmental Flow Research Centre at the University of Surrey, UK.

The study highlights a theoretical comparison and various interpretations of a recent inversion technique, called renormalization, developed for the reconstruction of unknown tracer emissions from their measured concentrations. The comparative interpretations are presented in relation to the other inversion techniques based on principle of regularization, Bayesian, minimum norm, maximum entropy on mean, and model resolution optimization. It is shown that the renormalization technique can be interpreted in a similar manner to other techniques, with a practical choice of a priori information and error statistics, while eliminating the need of additional constraints. The study shows that the proposed weight matrix and weighted Gram matrix offer a suitable deterministic choice to the background error and measurement covariance matrices, respectively, in the absence of statistical knowledge about background and measurement errors. The technique is advantageous since it (i) utilizes weights representing a priori information apparent to the monitoring network, (ii) avoids dependence on background source estimates, (iii) improves on alternative choices for the error statistics, (iv) overcomes the colocalization problem in a natural manner, and (v) provides an optimally resolved source reconstruction. A comparative illustration of source retrieval is made by using the real measurements from a continuous point release conducted in Fusion Field Trials, Dugway Proving Ground, Utah.

No abstract

In case of a release of a hazardous material (e.g., a chemical or a biological agent) in the atmosphere, estimation of the source from concentration observations (provided by a network of sensors) is a challenging inverse problem known as the atmospheric source term estimation (STE) problem. This study emphasizes a method, known in the literature as the renormalization inversion technique, for addressing this problem. This method provides a solution that has been interpreted as a weighted minimal norm solution and can be computed in terms of a generalized inverse of the sensitivity matrix of the sensors. This inverse is constructed by using an appropriate diagonal weight matrix whose components fulfill the so-called renormalizing conditions. The main contribution of this paper is that it proposes a new compact algorithm (it requires less than 15 lines of MATLAB code) to obtain, in a fast and efficient way, those optimal weights. To show that the algorithm, based on the properties of the resolution matrix, matches the requirements of emergency situations, analysis of the computational complexity and memory requirements is included. Some numerical experiments are also reported to show the efficiency of the algorithm.

This article highlights an optimization‐based approach to estimate the main characteristics (i.e., location and intensity) that describe the source of a release of a hazardous material (e.g., a chemical or biological agent) in the atmosphere. The method, which uses, as input data, concentration measurements provided by a network of sensors (and some additional meteorological data), is first formalized in the general framework of parametric source term estimation inverse problems. Some simple and new geometric interpretations of the method are proposed. Next, it is shown that the approach can also provide a realistic understanding of the uncertainties in the estimation of the characteristics of the source. Then, a strategy for an optimal deployment of mobile sensors in hazardous materials operations is proposed. To conclude, the source estimation approach and the network design strategy are evaluated with data from a wind tunnel experiment conducted at the Environmental Flow Research Center at the University of Surrey, UK.