Assimilation of GOES-R Geostationary Lightning Mapper Flash Extent Density Data in GSI 3DVar, EnKF, and Hybrid En3DVar for the Analysis and Short-Term Forecast of a Supercell Storm Case

Capabilities to assimilate Geostationary Operational Environmental Satellite “R-series ”(GOES-R) Geostationary Lightning Mapper(GLM) flash extent density(FED) data within the operational Gridpoint Statistical Interpolation ensemble Kalman filter(GSI-EnKF) framework were previously developed and tested with a mesoscale convective system(MCS) case. In this study, such capabilities are further developed to assimilate GOES GLM FED data within the GSI ensemble-variational(EnVar) hybrid data assimilation(DA) framework. The results of assimilating the GLM FED data using 3DVar, and pure En3DVar(PEn3DVar, using 100% ensemble covariance and no static covariance) are compared with those of EnKF/DfEnKF for a supercell storm case. The focus of this study is to validate the correctness and evaluate the performance of the new implementation rather than comparing the performance of FED DA among different DA schemes. Only the results of 3DVar and pEn3DVar are examined and compared with EnKF/DfEnKF. Assimilation of a single FED observation shows that the magnitude and horizontal extent of the analysis increments from PEn3DVar are generally larger than from EnKF, which is mainly caused by using different localization strategies in EnFK/DfEnKF and PEn3DVar as well as the integration limits of the graupel mass in the observation operator. Overall, the forecast performance of PEn3DVar is comparable to EnKF/DfEnKF, suggesting correct implementation.

The effect of glutathione on glucosinolate biosynthesis through the sulfur assimilation pathway in pakchoi associated with the growth conditions

Glucosinolates(GSLs) are a group of nitrogen-and sulfur-containing secondary metabolites, synthesized primarily in members of the Brassicaceae family, that play an important role in food flavor, plant antimicrobial activity, resistance to insect attack, stress tolerance, and human anti-cancer effects. As a sulfur-containing compound, glutathione has a strong connection with GSLs biosynthesis as a sulfur donor or redox system, and exists in reduced(glutathione;GSH) and oxidized(glutathione disulfide;GSSG) forms. However, the mechanism of GSH regulating GSLs biosynthesis remainds unclear. Hence, the exogenous therapy to pakchoi under normal growth condition and sulfur deficiency condition were conducted in this work to explore the relevant mechanism. The results showed that exogenous application of buthionine sulfoximine, an inhibitor of GSH synthesis, decreased the transcript levels of GSLs synthesis-related genes and transcription factors, as well as sulfur assimilation-related genes under the normal growth condition. Application of exogenous GSH inhibited the expression of GSLs synthesis-and sulfur assimilation-related genes under the normal condition, while the GSLs biosynthesis and the sulfur assimilation pathway were activated by exogenous application of GSH when the content of GSH in vivo of plants decreased owing to sulfur deficiency. Moreover,exogenous application of GSSG increased the transcript levels of GSLs synthesis-and sulfur assimilation-related genes under the normal growth condition and under sulfur deficiency. The present work provides new insights into the molecular mechanisms of GSLs biosynthesis underlying glutathione regulation.

Application of the finite analytic numerical method to a flowdependent variational data assimilation

An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.

Reconstructing urban wind flows for urban air mobility using reduced-order data assimilation

Advancements in uncrewed aircrafts and communications technologies have led to a wave of interest and investment in unmanned aircraft systems(UASs)and urban air mobility(UAM)vehicles over the past decade.To support this emerging aviation application,concepts for UAS/UAM traffic management(UTM)systems have been explored.Accurately characterizing and predicting the microscale weather conditions,winds in particular,will be critical to safe and efficient operations of the small UASs/UAM aircrafts within the UTM.This study implements a reduced order data assimilation approach to reduce discrepancies between the predicted urban wind speed with computational fluid dynamics(CFD)Reynolds-averaged Navier Stokes(RANS)model with real-world,limited and sparse observations.The developed data assimilation system is UrbanDA.These observations are simulated using a large eddy simulation(LES).The data assimilation approach is based on the time-independent variational framework and uses space reduction to reduce the memory cost of the process.This approach leads to error reduction throughout the simulated domain and the reconstructed field is different than the initial guess by ingesting wind speeds at sensor locations and hence taking into account flow unsteadiness in a time when only the mean flow quantities are resolved.Different locations where wind sensors can be installed are discussed in terms of their impact on the resulting wind field.It is shown that near-wall locations,near turbulence generation areas with high wind speeds have the highest impact.Approximating the model error with its principal mode provides a better agreement with the truth and the hazardous areas for UAS navigation increases by more than 10%as wind hazards resulting from buildings wakes are better simulated through this process.

Assimilation of FY-3D MWTS-Ⅱ Radiance with 3D Precipitation Detection and the Impacts on Typhoon Forecasts

Precipitation detection is an essential step in radiance assimilation because the uncertainties in precipitation would affect the radiative transfer calculation and observation errors.The traditional precipitation detection method for microwave only detects clouds and precipitation horizontally,without considering the three-dimensional distribution of clouds.Extending precipitation detection from 2D to 3D is expected to bring more useful information to the data assimilation without using the all-sky approach.In this study,the 3D precipitation detection method is adopted to assimilate Microwave Temperature Sounder-2(MWTS-Ⅱ)onboard the Fengyun-3D,which can dynamically detect the channels above precipitating clouds by considering the near-real-time cloud parameters.Cycling data assimilation and forecasting experiments for Typhoons Lekima(2019)and Mitag(2019)are carried out.Compared with the control experiment,the quantity of assimilated data with the 3D precipitation detection increases by approximately 23%.The quality of the additional MWTS-Ⅱradiance data is close to the clear-sky data.The case studies show that the average root-mean-square errors(RMSE)of prognostic variables are reduced by 1.7%in the upper troposphere,leading to an average reduction of4.53%in typhoon track forecasts.The detailed diagnoses of Typhoon Lekima(2019)further show that the additional MWTS-Ⅱradiances brought by the 3D precipitation detection facilitate portraying a more reasonable circulation situation,thus providing more precise structures.This paper preliminarily proves that 3D precipitation detection has potential added value for increasing satellite data utilization and improving typhoon forecasts.

Assimilation of the FY-4A AGRI Clear-Sky Radiance Data in a Regional Numerical Model and Its Impact on the Forecast of the“21·7”Henan Extremely Persistent Heavy Rainfall

Assimilation of the Advanced Geostationary Radiance Imager(AGRI)clear-sky radiance in a regional model is performed.The forecasting effectiveness of the assimilation of two water vapor(WV)channels with conventional observations for the“21·7”Henan extremely heavy rainfall is analyzed and compared with a baseline test that assimilates only conventional observations in this study.The results show that the 24-h cumulative precipitation forecast by the assimilation experiment with the addition of the AGRI exceeds 500 mm,compared to a maximum value of 532.6 mm measured by the national meteorological stations,and that the location of the maximum precipitation is consistent with the observations.The results for the short periods of intense precipitation processes are that the simulation of the location and intensity of the 3-h cumulative precipitation is also relatively accurate.The analysis increment shows that the main difference between the two sets of assimilation experiments is over the ocean due to the additional ocean observations provided by FY-4A,which compensates for the lack of ocean observations.The assimilation of satellite data adjusts the vertical and horizontal wind fields over the ocean by adjusting the atmospheric temperature and humidity,which ultimately results in a narrower and stronger WV transport path to the center of heavy precipitation in Zhengzhou in the lower troposphere.Conversely,the WV convergence and upward motion in the control experiment are more dispersed;therefore,the precipitation centers are also correspondingly more dispersed.

Machine Learning With Data Assimilation and Uncertainty Quantification for Dynamical Systems:A Review

Data assimilation(DA)and uncertainty quantification(UQ)are extensively used in analysing and reducing error propagation in high-dimensional spatial-temporal dynamics.Typical applications span from computational fluid dynamics(CFD)to geoscience and climate systems.Recently,much effort has been given in combining DA,UQ and machine learning(ML)techniques.These research efforts seek to address some critical challenges in high-dimensional dynamical systems,including but not limited to dynamical system identification,reduced order surrogate modelling,error covariance specification and model error correction.A large number of developed techniques and methodologies exhibit a broad applicability across numerous domains,resulting in the necessity for a comprehensive guide.This paper provides the first overview of state-of-the-art researches in this interdisciplinary field,covering a wide range of applications.This review is aimed at ML scientists who attempt to apply DA and UQ techniques to improve the accuracy and the interpretability of their models,but also at DA and UQ experts who intend to integrate cutting-edge ML approaches to their systems.Therefore,this article has a special focus on how ML methods can overcome the existing limits of DA and UQ,and vice versa.Some exciting perspectives of this rapidly developing research field are also discussed.Index Terms-Data assimilation(DA),deep learning,machine learning(ML),reduced-order-modelling,uncertainty quantification(UQ).

Increasing nitrogen absorption and assimilation ability under mixed NO_(3)^(–)and NH_(4)^(+) supply is a driver to promote growth of maize seedlings

Compared with sole nitrate (NO_(3)^(-)) or sole ammonium (NH_(4)^(+)) supply,mixed nitrogen (N) supply may promote growth of maize seedlings.Previous study suggested that mixed N supply not only increased photosynthesis rate,but also enhanced leaf growth by increasing auxin synthesis to build a large sink for C and N utilization.However,whether this process depends on N absorption is unknown.Here,maize seedlings were grown hydroponically with three N forms (NO_(3)^(-)only,75/25 NO_(3)^(-)/NH_(4)^(+) and NH_(4)^(+) only).The study results suggested that maize growth rate and N content of shoots under mixed N supply was little different to that under sole NO_(3)^(-)supply at 0–3 d,but was higher than under sole NO_(3)^(-)supply at 6–9 d.^(15)N influx rate under mixed N supply was greater than under sole NO_(3)^(-) or NH_(4)^(+) supply at 6–9 d,although NO_(3)^(-) and NH_(4)^(+) influx under mixed N supply were reduced compared to sole NO_(3)^(-) and NH_(4)^(+) supply,respectively.qRT-PCR determination suggested that the increased N absorption under mixed N supply may be related to the higher expression of NO_(3)^(-) transporters in roots,such as ZmNRT1.1A,ZmNRT1.1B,ZmNRT1.1C,ZmNRT1.2 and ZmNRT1.3,or NH_(4)^(+) absorption transporters,such as Zm AMT1.1A,especially the latter.Furthermore,plants had higher nitrate reductase (NR)glutamine synthase (GS) activity and amino acid content under mixed N supply than when under sole NO_(3)^(-) supply.The experiments with inhibitors of NR reductase and GS synthase further confirmed that N assimilation ability under mixed N supply was necessary to promote maize growth,especially for the reduction of NO_(3)^(-) by NR reductase.This research suggested that the increased processes of NO_(3)^(-)and NH_(4)^(+) assimilation by improving N-absorption ability of roots under mixed N supply may be the main driving force to increase maize growth.

Prediction of landslide block movement based on Kalman filtering data assimilation method

Compared with the study of single point motion of landslides,studying landslide block movement based on data from multiple monitoring points is of great significance for improving the accurate identification of landslide deformation.Based on the study of landslide block,this paper regarded the landslide block as a rigid body in particle swarm optimization algorithm.The monitoring data were organized to achieve the optimal state of landslide block,and the 6-degree of freedom pose of the landslide block was calculated after the regularization.Based on the characteristics of data from multiple monitoring points of landslide blocks,a prediction equation for the motion state of landslide blocks was established.By using Kalman filtering data assimilation method,the parameters of prediction equation for landslide block motion state were adjusted to achieve the optimal prediction.This paper took the Baishuihe landslide in the Three Gorges reservoir area as the research object.Based on the block segmentation of the landslide,the monitoring data of the Baishuihe landslide block were organized,6-degree of freedom pose of block B was calculated,and the Kalman filtering data assimilation method was used to predict the landslide block movement.The research results showed that the proposed prediction method of the landslide movement state has good prediction accuracy and meets the expected goal.This paper provides a new research method and thinking angle to study the motion state of landslide block.

A data assimilation-based forecast model of outer radiation belt electron fluxes

Because radiation belt electrons can pose a potential threat to the safety of satellites orbiting in space,it is of great importance to develop a reliable model that can predict the highly dynamic variations in outer radiation belt electron fluxes.In the present study,we develop a forecast model of radiation belt electron fluxes based on the data assimilation method,in terms of Van Allen Probe measurements combined with three-dimensional radiation belt numerical simulations.Our forecast model can cover the entire outer radiation belt with a high temporal resolution(1 hour)and a spatial resolution of 0.25 L over a wide range of both electron energy(0.1-5.0 MeV)and pitch angle(5°-90°).On the basis of this model,we forecast hourly electron fluxes for the next 1,2,and 3 days during an intense geomagnetic storm and evaluate the corresponding prediction performance.Our model can reasonably predict the stormtime evolution of radiation belt electrons with high prediction efficiency(up to~0.8-1).The best prediction performance is found for~0.3-3 MeV electrons at L=~3.25-4.5,which extends to higher L and lower energies with increasing pitch angle.Our results demonstrate that the forecast model developed can be a powerful tool to predict the spatiotemporal changes in outer radiation belt electron fluxes,and the model has both scientific significance and practical implications.

The effect of model errors in ensemble sequential assimilation of geomagnetic field

Earth’s magnetic field is generated in the fluid outer core through the dynamo process.Over the last decade,data assimilation has been used to retrieve the core dynamics and predict the evolution of the geomagnetic field.The presence of model errors in the geomagnetic data assimilation is inevitable because current numerical geodynamo models are still far from realistic core dynamics.In this paper,we investigate the effect of model errors in geomagnetic data assimilation based on ensemble Kalman filter(EnKF).We construct two dynamo models with different control parameters but exhibiting similar force balance and magnetic morphology at the core surface.We then use one dynamo model to generate synthetic observations and the other as the forward model in EnKF.Our test experiments show that the EnKF approach with the pre-setting model errors can nevertheless recover large-scale core surface flow and make a rough short-term(5-year)prediction.However,the data assimilation in the presence of model errors cannot keep improving the core state even though new observations are available.Motivated by the planned Macao Science Satellite-1,which is expected to provide improved internal geomagnetic field model,we also perform a test experiment using synthetic observations up to spherical harmonic degree l=18.Our results indicate that high-resolution observations are crucial in reconstructing small scale flow.

Satellite All-sky Infrared Radiance Assimilation:Recent Progress and Future Perspectives

Satellite infrared(IR)sounder and imager measurements have become one of the main sources of data used by data assimilation systems to generate initial conditions for numerical weather prediction(NWP)models and atmospheric analysis/reanalysis.This paper reviews the development of satellite IR data assimilation in NWP in recent years,especially the assimilation of all-sky satellite IR observations.The major challenges and future directions are outlined and discussed.

Incremental 4D-VAR assimilation scheme based on Lorenz model

Four-dimensional variational(4D-VAR) data assimilation method is a perfect data assimilation solution in theory, but the computational issue is quite difficult in operational implementation.The incremental 4D-VAR assimilation scheme is set up in order to reduce the computational cost. It is shown that the accuracy of the observations, the length of the assimilation window and the choice of the first guess have an important influence on the assimilation outcome through the contrast experiment. Compared with the standard 4D-VAR assimilation scheme, the incremental 4D-VAR assimilation scheme shows its advantage in the computation speed through an assimilation experiment.

An experimenal analysis for the impact of 2D variation assimilation of satellite data on typhoon track simulation

A series of test simulations are performed to evaluate the impact of satellite-derived meteorological data on numerical typhoon track prediction. Geostationary meteorological satellite (GMS-5) and NOAA’s TIROS operational vertical sounder (TOVS) observations are used in the experiments. A two-dimensional variation assimilation scheme is developed to assimilate the satellite data directly into the Penn State-NCAR nonhydrostatic meteorological model (MM5). Three-dimensional objective analyses fields based on T213 results and routine observations are employed as the background fields of the initialization. The comparisons of the simulated typhoon tracks are also carried out, which correspond respectively to the initialization scheme with two-dimensional variation(2D-Var), three-dimensional observational nudging and direct assimilation of satellite data. It is found that, comparing with the experiments without satellite data assimilation, the first two assimilation schemes lead to significant improvements on typhoon track prediction. Track errors reduce by 18% at 12 h for 2D-Var and from about 16% at 24 h to about 35% at 48 h for observational nudging. The simulated results based on assimilating different kinds of satellite data are also compared .

Wind Speed and Altitude Dependent AMDAR Observational Error and Its Impacts on Data Assimilation and Forecasting

Aircraft Meteorological Data Relay(AMDAR)observations have been widely used in numerical weather prediction(NWP)because of its high spatiotemporal resolution.The observational error of AMDAR is influenced by aircraft flight altitude and atmospheric condition.In this study,the wind speed and altitude dependent observational error of AMDAR is estimated.The statistical results show that the temperature and the observational error in wind speeds slightly decrease as altitude increases,and the observational error in wind speed increases as wind speed increases.Pseudo single AMDAR observation assimilation tests demonstrate that the wind speed and altitude dependent observational error can provide more reasonable analysis increment.Furthermore,to assess the performance of wind speed and altitude dependent observational error on data assimilation and forecasting,two-month 3-hourly cycling data assimilation and forecast experiments based on the Weather Research and Forecasting Model(WRF)and its Data Assimilation system(WRFDA)are performed for the period during 1 September-31 October,2017.The results of the two-month 3-hourly cycling experiments indicate that new observational error improves analysis and forecast of wind field and geo-potential height,and has slight improvements on temperature.The Fractions Skill Score(FSS)of the 6-h accumulated precipitation shows that new wind speed and altitude dependent observational error leads to better precipitation forecast skill than the default observational error in the WRFDA does.

The bZIP transcription factor MdHY5 regulates anthocyanin accumulation and nitrate assimilation in apple

The basic leucine zipper(bZIP)transcription factor HY5 plays a multifaceted role in plant growth and development.Here the apple MdHY5 gene was cloned based on its homology with Arabidopsis HY5.Expression analysis demonstrated that MdHY5 transcription was induced by light and abscisic acid treatments.Electrophoretic mobility shift assays and transient expression assays subsequently showed that MdHY5 positively regulated both its own transcription and that of MdMYB10 by binding to E-box and G-box motifs,respectively.Furthermore,we obtained transgenic apple calli that overexpressed the MdHY5 gene,and apple calli coloration assays showed that MdHY5 promoted anthocyanin accumulation by regulating expression of the MdMYB10 gene and downstream anthocyanin biosynthesis genes.In addition,the transcript levels of a series of nitrate reductase genes and nitrate uptake genes in both wild-type and transgenic apple calli were detected.In association with increased nitrate reductase activities and nitrate contents,the results indicated that MdHY5 might be an important regulator in nutrient assimilation.Taken together,these results indicate that MdHY5 plays a vital role in anthocyanin accumulation and nitrate assimilation in apple.

Impact of the Assimilation Frequency of Radar Data with the ARPS 3DVar and Cloud Analysis System on Forecasts of a Squall Line in Southern China

Assimilation configurations have significant impacts on analysis results and subsequent forecasts. A squall line system that occurred on 23 April 2007 over southern China was used to investigate the impacts of the data assimilation frequency of radar data on analyses and forecasts. A three-dimensional variational system was used to assimilate radial velocity data,and a cloud analysis system was used for reflectivity assimilation with a 2-h assimilation window covering the initial stage of the squall line. Two operators of radar reflectivity for cloud analyses corresponding to single-and double-moment schemes were used. In this study, we examined the sensitivity of assimilation frequency using 10-, 20-, 30-, and 60-min assimilation intervals. The results showed that analysis fields were not consistent with model dynamics and microphysics in general;thus, model states, including dynamic and microphysical variables, required approximately 20 min to reach a new balance after data assimilation in all experiments. Moreover, a 20-min data assimilation interval generally produced better forecasts for both single-and double-moment schemes in terms of equitable threat and bias scores. We conclude that a higher data assimilation frequency can produce a more intense cold pool and rear inflow jets but does not necessarily lead to a better forecast.

An Examination of the Predictability of Tropical Cyclone Genesis in High-Resolution Coupled Models with Dynamically Downscaled Coupled Data Assimilation Initialization

Predicting tropical cyclone(TC)genesis is of great societal importance but scientifically challenging.It requires fineresolution coupled models that properly represent air−sea interactions in the atmospheric responses to local warm sea surface temperatures and feedbacks,with aid from coherent coupled initialization.This study uses three sets of highresolution regional coupled models(RCMs)covering the Asia−Pacific(AP)region initialized with local observations and dynamically downscaled coupled data assimilation to evaluate the predictability of TC genesis in the West Pacific.The APRCMs consist of three sets of high-resolution configurations of the Weather Research and Forecasting−Regional Ocean Model System(WRF-ROMS):27-km WRF with 9-km ROMS,and 9-km WRF with 3-km ROMS.In this study,a 9-km WRF with 9-km ROMS coupled model system is also used in a case test for the predictability of TC genesis.Since the local sea surface temperatures and wind shear conditions that favor TC formation are better resolved,the enhanced-resolution coupled model tends to improve the predictability of TC genesis,which could be further improved by improving planetary boundary layer physics,thus resolving better air−sea and air−land interactions.

Analysis of influence of observation operator on sequential data assimilation through soil temperature simulation with common land model

An observation operator is a bridge linking the system state vector and observations in a data assimilation system. Despite its importance, the degree to which an observation operator influences the performance of data assimilation methods is still poorly understood. This study aimed to analyze the influences of linear and nonlinear observation operators on the sequential data assimilation through soil temperature simulation using the unscented particle filter(UPF) and the common land model. The linear observation operator between unprocessed simulations and observations was first established. To improve the correlation between simulations and observations, both were processed based on a series of equations. This processing essentially resulted in a nonlinear observation operator. The linear and nonlinear observation operators were then used along with the UPF in three assimilation experiments: an hourly in situ soil surface temperature assimilation, a daily in situ soil surface temperature assimilation, and a moderate resolution imaging spectroradiometer(MODIS) land surface temperature(LST) assimilation. The results show that the filter improved the soil temperature simulation significantly with the linear and nonlinear observation operators. The nonlinear observation operator improved the UPF's performance more significantly for the hourly and daily in situ observation assimilations than the linear observation operator did, while the situation was opposite for the MODIS LST assimilation. Because of the high assimilation frequency and data quality, the simulation accuracy was significantly improved in all soil layers for hourly in situ soil surface temperature assimilation, while the significant improvements of the simulation accuracy were limited to the lower soil layers for the assimilation experiments with low assimilation frequency or low data quality.