Optimal reduction of anthropogenic emissions for air pollution control and the retrieval of emission source from observed pollutants Ⅰ. Application of incomplete adjoint operator

The ultimate solution to anthropogenic air pollution depends on an adjustment and upgrade of industrial and energy structures. Before this process can be completed, reducing the anthropogenic pollutant emissions is an effective measure. This is a problem belonging to “Natural Cybernetics”, i.e., the problem of air pollution control should be solved together with the weather prediction; however, this is very complicated. Considering that heavy air pollution usually occurs in stable weather conditions and that the feedbacks between air pollutants and meteorological changes are insufficient, we propose a simplified natural cybernetics method. Here, an off-line air pollution evolution equation is first solved with data from a given anthropogenic emission inventory under the predicted weather conditions, and then, a related “incomplete adjoint problem” is solved to obtain the optimal reduction of anthropogenic emissions. Usually, such solution is sufficient for satisfying the air quality and economical/ social requirements. However, a better solution can be obtained by iteration after updating the emission inventory with the reduced anthropogenic emissions. Then, this paper discusses the retrieval of the pollutant emission source with a known spatio-temporal distribution of the pollutant concentrations, and a feasible mathematical method to achieve this is proposed. The retrieval of emission source would also help control air pollution.     

A SVD-based ensemble projection algorithm for calculating the conditional nonlinear optimal perturbation

Conditional nonlinear optimal perturbation(CNOP) is an extension of the linear singular vector technique in the nonlinear regime.It represents the initial perturbation that is subjected to a given physical constraint,and results in the largest nonlinear evolution at the prediction time.CNOP-type errors play an important role in the predictability of weather and climate.Generally,when calculating CNOP in a complicated numerical model,we need the gradient of the objective function with respect to the initial perturbations to provide the descent direction for searching the phase space.The adjoint technique is widely used to calculate the gradient of the objective function.However,it is difficult and cumbersome to construct the adjoint model of a complicated numerical model,which imposes a limitation on the application of CNOP.Based on previous research,this study proposes a new ensemble projection algorithm based on singular vector decomposition(SVD).The new algorithm avoids the localization procedure of previous ensemble projection algorithms,and overcomes the uncertainty caused by choosing the localization radius empirically.The new algorithm is applied to calculate the CNOP in an intermediate forecasting model.The results show that the CNOP obtained by the new ensemble-based algorithm can effectively approximate that calculated by the adjoint algorithm,and retains the general spatial characteristics of the latter.Hence,the new SVD-based ensemble projection algorithm proposed in this study is an effective method of approximating the CNOP.     

A variational method for correcting non-systematic errors in numerical weather prediction

A variational method based on previous numerical forecasts is developed to estimate and correct non-systematic component of numerical weather forecast error. In the method, it is assumed that the error is linearly dependent on some combination of the forecast fields, and three types of forecast combination are applied to identifying the forecasting error: 1) the forecasts at the ending time, 2) the combination of initial fields and the forecasts at the ending time, and 3) the combination of the forecasts at...     

Study of statistically correcting model CMAQ-MOS for forecasting regional air quality

Based on analysis of the air pollution observational data at 8 observation sites in Beijing including outer suburbs during the period from September 2004 to March 2005, this paper reveals synchronal and in-phase characteristics in the spatial and temporal variation of air pollutants on a city-proper scale at deferent sites; describes seasonal differences of the pollutant emission influence between the heating and non-heating periods, also significantly local differences of the pollutant emission influence between the urban district and outer suburbs, i.e. the spatial and temporal distribution of air pollutant is closely related with that of the pollutant emission intensity. This study shows that due to complexity of the spatial and temporal distribution of pollution emission sources, the new generation Community Multi-scale Air Quality (CMAQ) model developed by the EPA of USA produced forecasts, as other models did, with a systematic error of significantly lower than observations, albeit the model has better capability than previous models had in predicting the spatial distribution and variation tendency of multi-sort pollutants. The reason might be that the CMAQ adopts average amount of pollutant emission inventory, so that the model is difficult to objectively and finely describe the distribution and variation of pollution emission sources intensity on different spatial and temporal scales in the areas, in which the pollution is to be forecast. In order to correct the systematic prediction error resulting from the average pollutant emission inventory in CMAQ, this study proposes a new way of combining dynamics and statistics and establishes a statistically correcting model CMAQ-MOS for forecasts of regional air quality by utilizing the relationship of CMAQ outputs with corresponding observations, and tests the forecast capability. The investigation of experiments presents that CMAQ-MOS reduces the systematic errors of CMAQ because of the uncertainty of pollution emission inventory and improves the forecast level of air quality. Also this work employed a way of combining point and area forecasting, i.e. taking the products of CMAQ for a center site to forecast air pollution for other sites in vicinity with the scheme of model products "reanalysis" and average over the "area".     

Planning for Air Quality Concerns of the Future

— In recognition of growing needs for forecasts of air quality and atmospheric deposition to accompany classical weather forecasts, a new generation of atmospheric prediction models is slowly evolving. These share the common feature that atmospheric chemistry will be directly incorporated into advanced forecast schemes. It is argued that in most practical applications, the over-riding need is not for accurate prediction of some quantifiable air quality component, but rather a forecast of the probability of harmful consequences to exposure. In this event, it is not concentrations that need to be forecast, but the probability that concentrations will exceed some predetermined level at which consequences could be harmful. This argument extends from emergency response applications to ecosystem decline.     

Verification of SPE probability forecasts at the Space Environment Prediction Center (SEPC)

In space weather forecasting, forecast verification is necessary so that the forecast quality can be assessed. This paper provides an example of how to choose and devise verification methods and techniques according to different space weather forecast products. Solar proton events (SPEs) are hazardous space weather events, and forecasting them is one of the major tasks of the Space Environment Prediction Center (SEPC) at the National Space Science Center of the Chinese Academy of Sciences. Through analyzing SPE occurrence characteristics, SPE forecast properties, and verification requirements at SEPC, verification methods for SPE probability forecasts are identified, and verification results obtained. Overall, SPE probability forecasts at SEPC exhibit good accuracy, reliability, and discrimination. Compared with climatology and persistence forecasts, the SPE forecasts are more accurate. However, the forecasts for SPE onset days are substantially underestimated and need to be considerably improved.     

Short-term optimal operation of water systems using ensemble forecasts

Short-term water system operation can be realized using Model Predictive Control (MPC). MPC is a method for operational management of complex dynamic systems. Applied to open water systems, MPC provides integrated, optimal, and proactive management, when forecasts are available. Notwithstanding these properties, if forecast uncertainty is not properly taken into account, the system performance can critically deteriorate.Ensemble forecast is a way to represent short-term forecast uncertainty. An ensemble forecast is a set of possible future trajectories of a meteorological or hydrological system. The growing ensemble forecasts’ availability and accuracy raises the question on how to use them for operational management.The theoretical innovation presented here is the use of ensemble forecasts for optimal operation. Specifically, we introduce a tree based approach. We called the new method Tree-Based Model Predictive Control (TB-MPC). In TB-MPC, a tree is used to set up a Multistage Stochastic Programming, which finds a different optimal strategy for each branch and enhances the adaptivity to forecast uncertainty. Adaptivity reduces the sensitivity to wrong forecasts and improves the operational performance.TB-MPC is applied to the operational management of Salto Grande reservoir, located at the border between Argentina and Uruguay, and compared to other methods.     

Updating Real-Time Flood Forecasting Using a Fuzzy Rule-Based Model/Mise à Jour de Prévision de Crue en Temps Réel Grâce à un Modèle à Base de Règles Floues

Abstract

An updating technique is a tool to update the forecasts of mathematical flood forecasting model based on data observed in real time, and is an important element in a flood forecasting model. An error prediction model based on a fuzzy rule-based method was proposed as the updating technique in this work to improve one- to four-hour-ahead flood forecasts by a model that is composed of the grey rainfall model, the grey rainfall—runoff model and the modified Muskingum flow routing model. The coefficient of efficiency with respect to a benchmark is applied to test the applicability of the proposed fuzzy rule-based method. The analysis reveals that the fuzzy rule-based method can improve flood forecasts one to four hours ahead. The proposed updating technique can mitigate the problem of the phase lag in forecast hydrographs, and especially in forecast hydrographs with longer lead times.     

A parallel Eulerian–Lagrangian localized adjoint method

We describe an approach to parallelize Eulerian–Lagrangian localized adjoint methods such that no errors are introduced compared to the sequential case. This parallelization approach fully captures the hyperbolic features of the underlying problem. It uses an overlapping domain decomposition technique, and does not involve the introduction of artificial boundary conditions between subdomains. Implementation details on different parallel architectures are discussed.     

Nowcasting with INCA During SNOW-V10

A Central-European nowcasting system which has been developed for use in mountainous terrain is tested in the Whistler/Vancouver area as part of the SNOW-V10 experiment. The integrated nowcasting through comprehensive analysis system provides hourly updated gridded forecasts of temperature, humidity, and wind, as well as precipitation forecasts which are updated every 15 min. It is based on numerical weather prediction (NWP) output and real-time surface weather station and radar data. Verification of temperature, relative humidity, and wind against surface stations shows that forecast errors are significantly reduced in the nowcasting range compared to those of the driving NWP model. The main contribution to the improvement comes from the implicit bias correction due to use of the latest observations. Relative humidity shows the longest lasting effect, with >50 % reduction of mean absolute error up to +4 h. For temperature and wind speed this percentage is reached after +2 and +3 h, respectively. Two cases of precipitation nowcasting are discussed and verified qualitatively.     

An explicit four-dimensional variational data assimilation method based on the proper orthogonal decomposition: Theoretics and evaluation

The proper orthogonal decomposition (POD) method is used to construct a set of basis functions for spanning the ensemble of data in a certain least squares optimal sense. Compared with the singular value decomposition (SVD), the POD basis functions can capture more energy in the forecast ensemble space and can represent its spatial structure and temporal evolution more effectively. After the analysis variables are expressed by a truncated expansion of the POD basis vectors in the ensemble space, the control variables appear explicitly in the cost function, so that the adjoint model, which is used to derive the gradient of the cost function with respect to the control variables, is no longer needed. The application of this new technique significantly simplifies the data assimilation process. Several assimilation experiments show that this POD-based explicit four-dimensional variational data assimilation method performs much better than the usual ensemble Kalman filter method on both enhancing the assimilation precision and reducing the computation cost. It is also better than the SVD-based explicit four-dimensional assimilation method, especially when the forecast model is not perfect and the forecast error comes from both the noise of the initial filed and the uncertainty of the forecast model. Supported by the National Natural Science Foundation of China (Grant No. 40705035), National High Technology Research and Development Program of China (Grant No. 2007AA12Z144), Knowledge Innovation Project of Chinese Academy of Sciences (Grant Nos. KZCX2-YW-217 and KZCX2-YW-126-2), and National Basic Research Program of China (Grant No. 2005CB321704)     

An explicit four-dimensional variational data assimilation method

A new data assimilation method called the explicit four-dimensional variational (4DVAR) method is proposed. In this method, the singular value decomposition (SVD) is used to construct the orthogonal basis vectors from a forecast ensemble in a 4D space. The basis vectors represent not only the spatial structure of the analysis variables but also the temporal evolution. After the analysis variables are ex-pressed by a truncated expansion of the basis vectors in the 4D space, the control variables in the cost function appear explicitly, so that the adjoint model, which is used to derive the gradient of cost func-tion with respect to the control variables, is no longer needed. The new technique significantly simpli-fies the data assimilation process. The advantage of the proposed method is demonstrated by several experiments using a shallow water numerical model and the results are compared with those of the conventional 4DVAR. It is shown that when the observation points are very dense, the conventional 4DVAR is better than the proposed method. However, when the observation points are sparse, the proposed method performs better. The sensitivity of the proposed method with respect to errors in the observations and the numerical model is lower than that of the conventional method.     

A Study on the Performance of the NCMRWF Analysis and Forecasting System During Asian Summer Monsoon: Thermodynamic Aspects

—The thermodynamic characteristics of the Asian summer monsoon are examined with a global analysis-forecast system. In this study, we investigated the large-scale balances of heat and moisture by making use of operational analyses as well as forecast fields for June, July and August (JJA), 1994. Apart from elucidating systematic errors in the temperature and moisture fields, the study expounds the influence of these errors on the large-scale budgets of heat and moisture over the monsoon region. The temperature forecasts of the model delineate predominant cooling in the middle and lower tropospheres over the monsoon region. Similarly, the moisture forecasts evince a drying tendency in the lower troposphere. However, certain sectors of moderate moistening exist over the peninsular India and adjoining oceanic sectors of the Arabian Sea and Bay of Bengal.¶The broad features of the large-scale heat and moisture budgets represented by the analysis/forecast fields indicate good agreement with the observed aspects of the summer monsoon circulation. The model forecasts fail to retain the analyzed atmospheric variability in terms of the mean circulation, which is indicated by underestimation of various terms of heat and moisture budgets with an increase in the forecast period. Further, the forecasts depict an anomalous diabatic cooling layer in the lower middle troposphere of the monsoon region which inhibits vertical transfer of heat and moisture from the mixed layer of the atmospheric boundary layer to the middle troposphere. In effect, the monsoon circulation is considerably weakened with an increase in the forecast period. The treatment of shallow convection and the use of interactive clouds in the model can reduce the cooling bias considerably.     

An Application of the Multi-Physics Ensemble Kalman Filter to Typhoon Forecast

This study examines the roles of the multi-physics approach in accounting for model errors for typhoon forecasts with the local ensemble transform Kalman filter (LETKF). Experiments with forecasts of Typhoon Conson (2010) using the weather research and forecasting (WRF) model show that use of the WRF’s multiple physical parameterization schemes to represent the model uncertainties can help the LETKF provide better forecasts of Typhoon Conson in terms of the forecast errors, the ensemble spread, the root mean square errors, the cross-correlation between mass and wind field as well as the coherent structure of the ensemble spread along the storm center. Sensitivity experiments with the WRF model show that the optimum number of the multi-physics ensemble is roughly equal to the number of combinations of different physics schemes assigned in the multi-physics ensemble. Additional idealized experiments with the Lorenz 40-variable model to isolate the dual roles of the multi-physics ensemble in correcting model errors and expanding the local ensemble space show that the multi-physics approach appears to be more essential in augmenting the local rank representation of the LETKF algorithm rather than directly accounting for model errors during the early cycles. The results in this study suggest that the multi-physics approach is a good option for short-range forecast applications with full physics models in which the spinup of the ensemble Kalman filter may take too long for the ensemble spread to capture efficiently model errors and cross-correlations among model variables.     

Identification of hydraulic conductivity distributions in density dependent flow fields of submarine groundwater discharge modeling using adjoint-state sensitivities

A mathematical optimal control method is developed to identify a hydraulic conductivity distribution in a density dependent flow field. Using a variational method, the adjoint partial differential equations are obtained for the density- dependent state equations used for the saline aquifer water flow. The adjoint equations are numerically solved in through a finite difference method. The developed method is applied to identify the hydraulic conductivity distribution through the numerical solution of an optimal control problem. To demonstrate the effectiveness of the optimal control method, three numerical experiments are conducted with artificial observation data. The results indicate that the developed method has the potential to accurately identify the hydraulic conductivity distribution in a saline water aquifer flow system.     

Adjoint method for the optimum planning of industrial pollutant sources

The deterioration of air quality is threatening the life and health of people. Scientists in China and other countries have done a great deal of research work on the details of air pollution and the methods of preven-tion and control during the past decades. Up to now, most of the achievements are concentrated on the techniques of controlling pollutant sources and the programs of reduction, which focus on the improve-ment of air quality and the restoration of environment. The techniques of con…     

Towards improved post‐processing of hydrologic forecast ensembles

Forecast ensembles of hydrological and hydrometeorologial variables are prone to various uncertainties arising from climatology, model structure and parameters, and initial conditions at the forecast date. Post‐processing methods are usually applied to adjust the mean and variance of the ensemble without any knowledge about the uncertainty sources. This study initially addresses the drawbacks of a commonly used statistical technique, quantile mapping (QM), in bias correction of hydrologic forecasts. Then, an auxiliary variable, the failure index (γ), is proposed to estimate the ineffectiveness of the post‐processing method based on the agreement of adjusted forecasts with corresponding observations during an analysis period prior to the forecast date. An alternative post‐processor based on copula functions is then introduced such that marginal distributions of observations and model simulations are combined to create a multivariate joint distribution. A set of 2500 hypothetical forecast ensembles with parametric marginal distributions of simulated and observed variables are post‐processed with both QM and the proposed multivariate post‐processor. Deterministic forecast skills show that the proposed copula‐based post‐processing is more effective than the QM method in improving the forecasts. It is found that the performance of QM is highly correlated with the failure index, unlike the multivariate post‐processor. In probabilistic metrics, the proposed multivariate post‐processor generally outperforms QM. Further evaluation of techniques is conducted for river flow forecast of Sprague River basin in southern Oregon. Results show that the multivariate post‐processor performs better than the QM technique; it reduces the ensemble spread and is a more reliable approach for improving the forecast. Copyright © 2012 John Wiley & Sons, Ltd.     

Variational data assimilation in the transport of sediment in river

The variational method of data assimilation is used to solve an inverse problem in the transport of sediment in river, which plays an important role in the change of natural environment. The cost function is defined to measure the error between model predictions and field observations. The adjoint model of IAP river sedimentation model is created to obtain the gradient of the cost function with respect to control variables. The initial conditions are taken as the control variables; their optimal values can be retrieved by minimizing the cost function with limited memory quasi-Newton method (LMQN). The results show that the adjoint method approach can successfully make the model prediction well fit the simulated observations. And it is expected to use this method to solve other inverse problems of river sedimentation. But some numerical problems need to be discussed before applying to real river data. Project partially supported by the State Key Laboratory of Numerical Modelling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences     

大别山库区降水预报性能评估及应用对策

对降水预报进行性能评估及应用对策研究可以更好地发挥降水预报在水库调度中的决策支持作用.基于大别山库区近10 a汛期(2007—2016年5月1日—9月30日)24~168 h共7个预见期降水预报和地面降水观测资料,采用正确率、TS评分、概率统计、ROC曲线以及CTS等方法评估大别山库区降水预报性能,并以响洪甸水库为重点研究区域分析降水预报在水库调度中的应用对策.结果表明:1)大别山库区各量级的降水预报都有正预报技巧;24~72 h预见期降水预报的TS评分较高且空报率、漏报率也较低,具有较高的预报性能;但96 h及以上预见期降水预报性能明显下降,中雨以上量级空报率、漏报率较大,特别是对大暴雨及其以上量级的降水预报性能显著下降.2)大别山库区预报降水量级与实况降水量级基本符合,预报降水量级大于等于实况降水量级的概率超过75%;虽然降水预报量级上呈现出过度预报的现象,但降水过程预报对水库调度仍有较好的应用价值,应用时要考虑到降水预报量级可能存在偏差.3)转折性天气预报96 h及以上预见期CTS评分较低,但72 h以内预见期的性能明显改进,尤其是24 h预见期CTS评分也提高到了38.2%;水库调度可从长预见期的降水预报获取降水过程及其可能发生转折的信息,根据短预见期的降水预报进行调度方案调整.     

Fundamental Control Functions and Error Analysis in Variational Data Assimilation

The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find the initial condition function. The equation for the error of the optimal solution (analysis) is derived through the errors of the input data (background and observation errors). The numerical algorithm is developed to compute the sensitivity coefficients for the analysis error using the fundamental control functions. Application to the variational data assimilation problem for a model of ocean thermodynamics is considered.