This study examines the performance of coupling the deterministic four-dimensional variational assimilation system （4DVAR） with an ensemble Kalman filter （EnKF） to produce a superior hybrid approach for data assim...This study examines the performance of coupling the deterministic four-dimensional variational assimilation system （4DVAR） with an ensemble Kalman filter （EnKF） to produce a superior hybrid approach for data assimilation. The coupled assimilation scheme （E4DVAR） benefits from using the state-dependent uncertainty provided by EnKF while taking advantage of 4DVAR in preventing filter divergence： the 4DVAR analysis produces posterior maximum likelihood solutions through minimization of a cost function about which the ensemble perturbations are transformed, and the resulting ensemble analysis can be propagated forward both for the next assimilation cycle and as a basis for ensemble forecasting. The feasibility and effectiveness of this coupled approach are demonstrated in an idealized model with simulated observations. It is found that the E4DVAR is capable of outperforming both 4DVAR and the EnKF under both perfect- and imperfect-model scenarios. The performance of the coupled scheme is also less sensitive to either the ensemble size or the assimilation window length than those for standard EnKF or 4DVAR implementations.展开更多
Environmental systems including our atmosphere oceans, biological… etc. can be modeled by mathematical equations to estimate their states. These equations can be solved with numerical methods. Initial and boundary co...Environmental systems including our atmosphere oceans, biological… etc. can be modeled by mathematical equations to estimate their states. These equations can be solved with numerical methods. Initial and boundary conditions are needed for such of these numerical methods. Predication and simulations for different case studies are major sources for the great importance of these models. Satellite data from different wide ranges of sensors provide observations that indicate system state. So both numerical models and satellite data provide estimation of system states, and between the different estimations it is required the best estimate for system state. Assimilation of observations in numerical weather models with data assimilation techniques provide an improved estimate of system states. In this work, highlights on the mathematical perspective for data assimilation methods are introduced. Least square estimation techniques are introduced because it is considered the basic mathematical building block for data assimilation methods. Stochastic version of least square is included to handle the error in both model and observation. Then the three and four dimensional variational assimilation 3dvar and 4dvar respectively will be handled. Kalman filters and its derivatives Extended, (KF, EKF, ENKF) and hybrid filters are introduced.展开更多
Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation,and is called groundwater contaminant source identification.This is a kind of typical gro...Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation,and is called groundwater contaminant source identification.This is a kind of typical groundwater inverse problem,and the solution is usually ill-posed.Especially considering the spatial variability of hydraulic conductivity field,the identification process is more challenging.In this paper,the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model(MT3DMS)and a data assimilation method(Iterative local update ensemble smoother,ILUES).In addition,Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction.In practical problems,the geostatistical method is usually used to characterize the hydraulic conductivity field,and only the contaminant source information is inversely calculated in the identification process.In this study,the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field.The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site,and ILUES has good performance in solving high-dimensional parameter inversion problems.展开更多
A method has been presented to improve ensemble forecast by utilizing these initial members generated by four-dimensional variational data assimilation (4-D VDA), to conquer limitation of those initial members generat...A method has been presented to improve ensemble forecast by utilizing these initial members generated by four-dimensional variational data assimilation (4-D VDA), to conquer limitation of those initial members generated by Monte Carlo forecast (MCF) or lagged average forecast (LAF). This method possesses significant statistical characteristic of MCF, and by virtue of LAF that contains multi-time information and its initial members are harmonic with展开更多
Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational...Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.展开更多
An observation localization scheme is introduced into an ensemble-based three-dimensional variational （3DVar） assimilation method based on the singular value decomposition technique （SVD-En3DVar） to im- prove assi...An observation localization scheme is introduced into an ensemble-based three-dimensional variational （3DVar） assimilation method based on the singular value decomposition technique （SVD-En3DVar） to im- prove assimilation skill. A point-by-point analysis technique is adopted in which the weight of each obser- vation decreases with increasing distance between the analysis point and the observation point. A set of numerical experiments, in which simulated Doppler radar data are assimilated into the Weather Research and Forecasting （WRF） model, is designed to test the scheme. The results are compared with those ob- tained using the original global and local patch schemes in SVD-En3DVar, neither of which includes this type of observation localization. The observation localization scheme not only eliminates spurious analysis increments in areas of missing data, but also avoids the discontinuous analysis fields that arise from the local patch scheme. The new scheme provides better analysis fields and a more reasonable short-range rainfall forecast than the original schemes. Additional forecast experiments that assimilate real data from i0 radars indicate that the short-term precipitation forecast skill can be improved by assimilating radar data and the observation localization scheme provides a better forecast than the other two schemes.展开更多
Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere,ocean,and continental water sources.The terrestrial water storage changes,detected as gravity chan...Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere,ocean,and continental water sources.The terrestrial water storage changes,detected as gravity changes by the Gravity Recovery and Climate Experiment(GRACE) satellites,are mainly caused by precipitation,evapotranspiration,river transportation and downward infiltration processes.In this study,a land data assimilation system LDAS-G was developed to assimilate the GRACE terrestrial water storage(TWS) data into the Community Land Model(CLM3.5) using the POD-based ensemble four-dimensional variational assimilation method PODEn4 DVar,disaggregating the GRACE large-scale terrestrial water storage changes vertically and in time,and placing constraints on the simulation of vertical hydrological variables to improve land surface hydrological simulations.The ideal experiments conducted at a single point and assimilation experiments carried out over China by the LDAS-G data assimilation system showed that the system developed in this study improved the simulation of land surface hydrological variables,indicating the potential of GRACE data assimilation in large-scale land surface hydrological research and applications.展开更多
Background error covariance(BEC)plays an essential role in variational data assimilation.Most variational data assimilation systems still use static BEC.Actually,the characteristics of BEC vary with season,day,and eve...Background error covariance(BEC)plays an essential role in variational data assimilation.Most variational data assimilation systems still use static BEC.Actually,the characteristics of BEC vary with season,day,and even hour of the background.National Meteorological Center-based diurnally varying BECs had been proposed,but the diurnal variation characteristics were gained by climatic samples.Ensemble methods can obtain the background error characteristics that suit the samples in the current moment.Therefore,to gain more reasonable diurnally varying BECs,in this study,ensemble-based diurnally varying BECs are generated and the diurnal variation characteristics are discussed.Their impacts are then evaluated by cycling data assimilation and forecasting experiments for a week based on the operational China Meteorological Administration-Beijing system.Clear diurnal variation in the standard deviation of ensemble forecasts and ensemble-based BECs can be identified,consistent with the diurnal variation characteristics of the atmosphere.The results of one-week cycling data assimilation and forecasting show that the application of diurnally varying BECs reduces the RMSEs in the analysis and 6-h forecast.Detailed analysis of a convective rainfall case shows that the distribution of the accumulated precipitation forecast using the diurnally varying BECs is closer to the observation than using the static BEC.Besides,the cycle-averaged precipitation scores in all magnitudes are improved,especially for the heavy precipitation,indicating the potential of using diurnally varying BEC in operational applications.展开更多
In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moi...In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moisture profile by assimilating surface soil moisture observations. The Arou observation station, located in the upper reaches of the Heihe River in northwestern China, was selected to test the proposed method. Three numeric experiments were designed and performed to analyze the influence of uncertainties in model parameters, atmospheric forcing, and the model's physical mechanics on soil moisture estimates. Several assimilation schemes based on the ensemble Kalman filter(En KF), ensemble Kalman smoother(En KS), and dual En KF(DEn KF) were also compared in this study. The results demonstrate that soil moisture and soil properties can be simultaneously estimated by state-parameter estimation methods, which can provide more accurate estimation of soil moisture than traditional filter methods such as En KF and En KS. The estimation accuracy of the model parameters decreased with increasing error sources. DEn KS outperformed DEn KF in estimating soil moisture in most cases, especially where few observations were available. This study demonstrates that the DEn KS approach is a useful and practical way to improve soil moisture estimation.展开更多
针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance...针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance tomography,ERT)法采集的ERT观测数据,实现对污染源源强和渗透系数场的联合反演。以此为基础设计3组数值算例,比较不同类型观测数据对反演精度的影响。研究结果表明:融合ERT数据的ES-MDA算法对模型参数的反演精度更高,并且将ERT数据和传统的质量浓度与水头观测数据相结合,能进一步优化反演结果。展开更多
基金sponsored by the U.S. National Science Foundation (Grant No.ATM0205599)the U.S. Offce of Navy Research under Grant N000140410471Dr. James A. Hansen was partially supported by US Offce of Naval Research (Grant No. N00014-06-1-0500)
文摘This study examines the performance of coupling the deterministic four-dimensional variational assimilation system （4DVAR） with an ensemble Kalman filter （EnKF） to produce a superior hybrid approach for data assimilation. The coupled assimilation scheme （E4DVAR） benefits from using the state-dependent uncertainty provided by EnKF while taking advantage of 4DVAR in preventing filter divergence： the 4DVAR analysis produces posterior maximum likelihood solutions through minimization of a cost function about which the ensemble perturbations are transformed, and the resulting ensemble analysis can be propagated forward both for the next assimilation cycle and as a basis for ensemble forecasting. The feasibility and effectiveness of this coupled approach are demonstrated in an idealized model with simulated observations. It is found that the E4DVAR is capable of outperforming both 4DVAR and the EnKF under both perfect- and imperfect-model scenarios. The performance of the coupled scheme is also less sensitive to either the ensemble size or the assimilation window length than those for standard EnKF or 4DVAR implementations.
文摘Environmental systems including our atmosphere oceans, biological… etc. can be modeled by mathematical equations to estimate their states. These equations can be solved with numerical methods. Initial and boundary conditions are needed for such of these numerical methods. Predication and simulations for different case studies are major sources for the great importance of these models. Satellite data from different wide ranges of sensors provide observations that indicate system state. So both numerical models and satellite data provide estimation of system states, and between the different estimations it is required the best estimate for system state. Assimilation of observations in numerical weather models with data assimilation techniques provide an improved estimate of system states. In this work, highlights on the mathematical perspective for data assimilation methods are introduced. Least square estimation techniques are introduced because it is considered the basic mathematical building block for data assimilation methods. Stochastic version of least square is included to handle the error in both model and observation. Then the three and four dimensional variational assimilation 3dvar and 4dvar respectively will be handled. Kalman filters and its derivatives Extended, (KF, EKF, ENKF) and hybrid filters are introduced.
基金supported by the Fundamental Research Funds for the Central Universities(No.22120190013)National Natural Science Foundation of China(No.41807187)
文摘Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation,and is called groundwater contaminant source identification.This is a kind of typical groundwater inverse problem,and the solution is usually ill-posed.Especially considering the spatial variability of hydraulic conductivity field,the identification process is more challenging.In this paper,the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model(MT3DMS)and a data assimilation method(Iterative local update ensemble smoother,ILUES).In addition,Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction.In practical problems,the geostatistical method is usually used to characterize the hydraulic conductivity field,and only the contaminant source information is inversely calculated in the identification process.In this study,the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field.The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site,and ILUES has good performance in solving high-dimensional parameter inversion problems.
文摘A method has been presented to improve ensemble forecast by utilizing these initial members generated by four-dimensional variational data assimilation (4-D VDA), to conquer limitation of those initial members generated by Monte Carlo forecast (MCF) or lagged average forecast (LAF). This method possesses significant statistical characteristic of MCF, and by virtue of LAF that contains multi-time information and its initial members are harmonic with
基金This work was jointly supported by the National Natural Science Foundation of China(Grant Nos.41975137,42175012,and 41475097)the National Key Research and Development Program(Grant No.2018YFF0300103).
文摘Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.
基金Supported by the Open Project Fund of the State Key Laboratory of Severe Weather of Chinese Academy of Meteorological Sciences, National Natural Science Foundation of China (40875063 and 41275102)Fundamental Research Fund for Central Universities of China (lzujbky-2010-9)
文摘An observation localization scheme is introduced into an ensemble-based three-dimensional variational （3DVar） assimilation method based on the singular value decomposition technique （SVD-En3DVar） to im- prove assimilation skill. A point-by-point analysis technique is adopted in which the weight of each obser- vation decreases with increasing distance between the analysis point and the observation point. A set of numerical experiments, in which simulated Doppler radar data are assimilated into the Weather Research and Forecasting （WRF） model, is designed to test the scheme. The results are compared with those ob- tained using the original global and local patch schemes in SVD-En3DVar, neither of which includes this type of observation localization. The observation localization scheme not only eliminates spurious analysis increments in areas of missing data, but also avoids the discontinuous analysis fields that arise from the local patch scheme. The new scheme provides better analysis fields and a more reasonable short-range rainfall forecast than the original schemes. Additional forecast experiments that assimilate real data from i0 radars indicate that the short-term precipitation forecast skill can be improved by assimilating radar data and the observation localization scheme provides a better forecast than the other two schemes.
基金supported by the National Natural Science Foundation of China(Grant Nos.41075062,91125016)the National Basic Research Program of China(Grants Nos.2010CB951001,2010CB428403)
文摘Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere,ocean,and continental water sources.The terrestrial water storage changes,detected as gravity changes by the Gravity Recovery and Climate Experiment(GRACE) satellites,are mainly caused by precipitation,evapotranspiration,river transportation and downward infiltration processes.In this study,a land data assimilation system LDAS-G was developed to assimilate the GRACE terrestrial water storage(TWS) data into the Community Land Model(CLM3.5) using the POD-based ensemble four-dimensional variational assimilation method PODEn4 DVar,disaggregating the GRACE large-scale terrestrial water storage changes vertically and in time,and placing constraints on the simulation of vertical hydrological variables to improve land surface hydrological simulations.The ideal experiments conducted at a single point and assimilation experiments carried out over China by the LDAS-G data assimilation system showed that the system developed in this study improved the simulation of land surface hydrological variables,indicating the potential of GRACE data assimilation in large-scale land surface hydrological research and applications.
基金This work was jointly sponsored by the National Natural Science Foundation of China[grant number 42075148]the Outreach Projects of the State Key Laboratory of Severe Weather[grant number 2021LASWA08]+1 种基金the Outreach Projects of the Key Laboratory of Meteorological Disaster[grant number KLME202209]supported by the High-Performance Computing Center of Nanjing University of Information Science and Technology(NUIST).
文摘Background error covariance(BEC)plays an essential role in variational data assimilation.Most variational data assimilation systems still use static BEC.Actually,the characteristics of BEC vary with season,day,and even hour of the background.National Meteorological Center-based diurnally varying BECs had been proposed,but the diurnal variation characteristics were gained by climatic samples.Ensemble methods can obtain the background error characteristics that suit the samples in the current moment.Therefore,to gain more reasonable diurnally varying BECs,in this study,ensemble-based diurnally varying BECs are generated and the diurnal variation characteristics are discussed.Their impacts are then evaluated by cycling data assimilation and forecasting experiments for a week based on the operational China Meteorological Administration-Beijing system.Clear diurnal variation in the standard deviation of ensemble forecasts and ensemble-based BECs can be identified,consistent with the diurnal variation characteristics of the atmosphere.The results of one-week cycling data assimilation and forecasting show that the application of diurnally varying BECs reduces the RMSEs in the analysis and 6-h forecast.Detailed analysis of a convective rainfall case shows that the distribution of the accumulated precipitation forecast using the diurnally varying BECs is closer to the observation than using the static BEC.Besides,the cycle-averaged precipitation scores in all magnitudes are improved,especially for the heavy precipitation,indicating the potential of using diurnally varying BEC in operational applications.
基金supported by the Natural National Science Foundation of China(Grant Nos.91325106&41271358)the Hundred Talent Program of the Chinese Academy of Sciences(Grant No.29Y127D01)+1 种基金the Cross-disciplinary Collaborative Teams Program for ScienceTechnology and Innovation of the Chinese Academy of Sciences
文摘In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moisture profile by assimilating surface soil moisture observations. The Arou observation station, located in the upper reaches of the Heihe River in northwestern China, was selected to test the proposed method. Three numeric experiments were designed and performed to analyze the influence of uncertainties in model parameters, atmospheric forcing, and the model's physical mechanics on soil moisture estimates. Several assimilation schemes based on the ensemble Kalman filter(En KF), ensemble Kalman smoother(En KS), and dual En KF(DEn KF) were also compared in this study. The results demonstrate that soil moisture and soil properties can be simultaneously estimated by state-parameter estimation methods, which can provide more accurate estimation of soil moisture than traditional filter methods such as En KF and En KS. The estimation accuracy of the model parameters decreased with increasing error sources. DEn KS outperformed DEn KF in estimating soil moisture in most cases, especially where few observations were available. This study demonstrates that the DEn KS approach is a useful and practical way to improve soil moisture estimation.
文摘针对未知的污染场地,为了准确估计污染物运移模型的参数,提出一种基于多重数据同化集合平滑器(ensemble smoother with multiple data assimilation,ES-MDA)算法的地下水模型参数反演方法,通过融合由高密度电阻率(electrical resistance tomography,ERT)法采集的ERT观测数据,实现对污染源源强和渗透系数场的联合反演。以此为基础设计3组数值算例,比较不同类型观测数据对反演精度的影响。研究结果表明:融合ERT数据的ES-MDA算法对模型参数的反演精度更高,并且将ERT数据和传统的质量浓度与水头观测数据相结合,能进一步优化反演结果。
