CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate ...CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.展开更多
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 Interp...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.展开更多
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 ac...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.展开更多
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...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.展开更多
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 ...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).展开更多
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 ob...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.展开更多
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 ...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.展开更多
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 det...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.展开更多
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 photosynthes...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.展开更多
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 th...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.展开更多
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 f...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.展开更多
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 landsl...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.展开更多
This paper presents an attempt at assimilating clear-sky FY-4A Advanced Geosynchronous Radiation Imager(AGRI)radiances from two water vapor channels for the prediction of three landfalling typhoon events over the West...This paper presents an attempt at assimilating clear-sky FY-4A Advanced Geosynchronous Radiation Imager(AGRI)radiances from two water vapor channels for the prediction of three landfalling typhoon events over the West Pacific Ocean using the 3DVar data assimilation(DA)method along with the WRF model.A channel-sensitive cloud detection scheme based on the particle filter(PF)algorithm is developed and examined against a cloud detection scheme using the multivariate and minimum residual(MMR)algorithm and another traditional cloud mask–dependent cloud detection scheme.Results show that both channel-sensitive cloud detection schemes are effective,while the PF scheme is able to reserve more pixels than the MMR scheme for the same channel.In general,the added value of AGRI radiances is confirmed when comparing with the control experiment without AGRI radiances.Moreover,it is found that the analysis fields of the PF experiment are mostly improved in terms of better depicting the typhoon,including the temperature,moisture,and dynamical conditions.The typhoon track forecast skill is improved with AGRI radiance DA,which could be explained by better simulating the upper trough.The impact of assimilating AGRI radiances on typhoon intensity forecasts is small.On the other hand,improved rainfall forecasts from AGRI DA experiments are found along with reduced errors for both the thermodynamic and moisture fields,albeit the improvements are limited.展开更多
It is not reasonable that one can only use the adjoint of model in data assimilation. The simulated numerical experiment shows that for the tidal model, the result of the adjoint of equation is almost the same as that...It is not reasonable that one can only use the adjoint of model in data assimilation. The simulated numerical experiment shows that for the tidal model, the result of the adjoint of equation is almost the same as that of the adjoint of model: the averaged absolute difference of the amplitude between observations and simulation is less than 5.0 cm and that of the phase-lag is less than 5.0°. The results are both in good agreement with the observed M2 tide in the Bohai Sea and the Yellow Sea. For comparison, the traditional methods also have been used to simulate M2 tide in the Bohai Sea and the Yellow Sea. The initial guess values of the boundary conditions are given first, and then are adjusted to acquire the simulated results that are as close as possible to the observations. As the boundary conditions contain 72 values, which should be adjusted and how to adjust them can only be partially solved by adjusting them many times. The satisfied results are hard to acquire even gigantic efforts are done. Here, the automation of the treatment of the open boundary conditions is realized. The method is unique and superior to the traditional methods. It is emphasized that if the adjoint of equation is used, tedious and complicated mathematical deduction can be avoided. Therefore the adjoint of equation should attract much attention.展开更多
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.展开更多
This paper mainly discusses the openness and assimilation of the words system from the viewpoint of loan words between English and Chinese.Both the history of loan words and the ways of borrowing words in English and ...This paper mainly discusses the openness and assimilation of the words system from the viewpoint of loan words between English and Chinese.Both the history of loan words and the ways of borrowing words in English and Chinese have shown us the strong abilities of languages.展开更多
In order to solve the difficult problem of typhoon track prediction due to the sparsity of conventional data over the tropical ocean, in this paper, the No. 0205 typhoon Rammasun of 4-6 July 2002 is studied and an exp...In order to solve the difficult problem of typhoon track prediction due to the sparsity of conventional data over the tropical ocean, in this paper, the No. 0205 typhoon Rammasun of 4-6 July 2002 is studied and an experiment of the typhoon track prediction is made with the direct use of the Advanced TIROS-N Operational Vertical Sounder (ATOVS) microwave radiance data in three-dimensional variational data assimilation. The prediction result shows that the experiment with the ATOVS microwave radiance data can not only successfully predict the observed fact that typhoon Rammasun moves northward and turns right, but can also simulate the action of the fast movement of the typhoon, which cannot be simulated with only conventional radiosonde data. The skill of the typhoon track prediction with the ATOVS microwave radiance data is much better than that without the ATOVS data. The typhoon track prediction of the former scheme is consistent in time and in location with the observation. The direct assimilation of展开更多
The Panzhihua intrusion in southwest China is part of the Emeishan Large Igneous Province and host of a large Fe-Ti-V ore deposit.During emplacement of the main intrusion,multiple generations of mafic dykes invaded ca...The Panzhihua intrusion in southwest China is part of the Emeishan Large Igneous Province and host of a large Fe-Ti-V ore deposit.During emplacement of the main intrusion,multiple generations of mafic dykes invaded carbonate wall rocks,producing a large contact aureole.We measured the oxygen-isotope composition of the intrusions,their constituent minerals,and samples of the country rock.Magnetite and plagioclase from Panzhihua intrusion haveδ18O values that are consistent with magmatic equilibrium, and formed from magmas withδ18O values that were 1-2‰higher than expected in a mantle-derived magma.The unmetamorphosed country rock has highδ18O values,ranging from 13.2‰(sandstone) to 24.6-28.6‰(dolomite).The skarns and marbles from the aureole have lowerδ18O andδ13C values than their protolith suggesting interaction with fluids that were in exchange equilibrium with the adjacent mafic magmas and especially the numerous mafic dykes that intruded the aureole.This would explain the alteration ofδ18O of the dykes which have significantly higher values than expected for a mantle-derived magma.Depending on the exactδ18O values assumed for the magma and contaminant, the amount of assimilation required to produce the elevatedδ18O value of the Panzhihua intrusion was between 8 and 13.7 wt.%,assuming simple mixing.The exact mechanism of contamination is unclear but may involve a combination of assimilation of bulk country rock,mixing with a melt of the country rock and exchange with CO2-rich fluid derived from decarbonation of the marls and dolomites.These mechanisms,particularly the latter,were probably involved in the formation of the Fe-Ti-V ores.展开更多
基金supported by the General Project of Top-Design of Multi-Scale Nature-Social ModelsData Support and Decision Support System for NSFC Carbon Neutrality Major Project(42341202)the Basic Scientific Research Fund of the Chinese Academy of Meteorological Sciences(2021Z014)。
文摘CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.
基金supported by NOAA JTTI award via Grant #NA21OAR4590165, NOAA GOESR Program funding via Grant #NA16OAR4320115provided by NOAA/Office of Oceanic and Atmospheric Research under NOAA-University of Oklahoma Cooperative Agreement #NA11OAR4320072, U.S. Department of Commercesupported by the National Oceanic and Atmospheric Administration (NOAA) of the U.S. Department of Commerce via Grant #NA18NWS4680063。
文摘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.
基金funded by the National Natural Science Foundation of China (Grant Nos.31972394 and 31501748)。
文摘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.
基金The National Key Research and Development Program of China under contract Nos 2022YFC3104804,2021YFC3101501,and 2017YFC1404103the National Programme on Global Change and Air-Sea Interaction of China under contract No.GASI-IPOVAI-04the National Natural Science Foundation of China under contract Nos 41876014,41606039,and 11801402.
文摘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.
基金the support of the Leverhulme Centre for Wildfires,Environment and Society through the Leverhulme Trust(RC-2018-023)Sibo Cheng,César Quilodran-Casas,and Rossella Arcucci acknowledge the support of the PREMIERE project(EP/T000414/1)+5 种基金the support of EPSRC grant:PURIFY(EP/V000756/1)the Fundamental Research Funds for the Central Universitiesthe support of the SASIP project(353)funded by Schmidt Futures–a philanthropic initiative that seeks to improve societal outcomes through the development of emerging science and technologiesDFG for the Heisenberg Programm Award(JA 1077/4-1)the National Natural Science Foundation of China(61976120)the Natural Science Key Foundat ion of Jiangsu Education Department(21KJA510004)。
文摘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).
基金supported by the National Key R&D Program of China(Grant Nos.2017YFC1501803 and 2017YFC1502102)。
文摘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.
基金supported by the National Natural Science Foundation of China (Grant Nos. 42025404, 42188101, and 42241143)the National Key R&D Program of China (Grant Nos. 2022YFF0503700 and 2022YFF0503900)+1 种基金the B-type Strategic Priority Program of the Chinese Academy of Sciences (Grant No. XDB41000000)the Fundamental Research Funds for the Central Universities (Grant No. 2042022kf1012)
文摘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.
基金jointly sponsored by the National Key Research and Development Program of China(Grant Nos.2018YFC1506701 and 2017YFC1502102)the National Natural Science Foundation of China(Grant No.41675102)。
文摘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.
基金supported by the National Basic Research Program of China (2015CB150402)the National Natural Science Foundation of China (31672221 and 31421092)the Science Foundation for Young Scholars of Tobacco Research Institute of Chinese Academy of Agricultural Sciences (2022C03 and 20211302)。
文摘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.
文摘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.
基金supported by the Macao Foundation and preresearch project on Civil Aerospace Technologies of CNSA(D020308,D020303)the Macao Science and Technology Development Fund(0001/2019/A1)the National Natural Science Foundation of China(41904066,42142034)。
文摘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.
基金supported by National Natural Science Foundation of China(Grant Nos.42090054,52027814 and 41772376)the Open Fund of the Technology Innovation Center for Automated Geological Disaster Monitoring,Ministry of Natural Resources(Grant No.2022058014)。
文摘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.
基金primarily supported by the Chinese National Natural Science Foundation of China(Grant No. G42192553)Open Fund of Fujian Key Laboratory ofSevere Weather and Key Laboratory of Straits Severe Weather(Grant No. 2023KFKT03)+6 种基金the Open Project Fund of China Meteorological Administration Basin Heavy Rainfall Key Laboratory(Grant No. 2023BHR-Y20)the Open Fund of the State Key Laboratory of Remote Sensing Science (Grant No. OFSLRSS202321)the Program of Shanghai Academic/Technology Research Leader(Grant No. 21XD1404500)the Shanghai Typhoon Research Foundation (Grant No. TFJJ202107)the Chinese National Natural Science Foundation of China (Grant No. G41805016)the National Meteorological Center Foundation (Grant No. FY-APP-2021.0207)the High Performance Computing Center of Nanjing University of Information Science&Technology for their support of this work
文摘This paper presents an attempt at assimilating clear-sky FY-4A Advanced Geosynchronous Radiation Imager(AGRI)radiances from two water vapor channels for the prediction of three landfalling typhoon events over the West Pacific Ocean using the 3DVar data assimilation(DA)method along with the WRF model.A channel-sensitive cloud detection scheme based on the particle filter(PF)algorithm is developed and examined against a cloud detection scheme using the multivariate and minimum residual(MMR)algorithm and another traditional cloud mask–dependent cloud detection scheme.Results show that both channel-sensitive cloud detection schemes are effective,while the PF scheme is able to reserve more pixels than the MMR scheme for the same channel.In general,the added value of AGRI radiances is confirmed when comparing with the control experiment without AGRI radiances.Moreover,it is found that the analysis fields of the PF experiment are mostly improved in terms of better depicting the typhoon,including the temperature,moisture,and dynamical conditions.The typhoon track forecast skill is improved with AGRI radiance DA,which could be explained by better simulating the upper trough.The impact of assimilating AGRI radiances on typhoon intensity forecasts is small.On the other hand,improved rainfall forecasts from AGRI DA experiments are found along with reduced errors for both the thermodynamic and moisture fields,albeit the improvements are limited.
文摘It is not reasonable that one can only use the adjoint of model in data assimilation. The simulated numerical experiment shows that for the tidal model, the result of the adjoint of equation is almost the same as that of the adjoint of model: the averaged absolute difference of the amplitude between observations and simulation is less than 5.0 cm and that of the phase-lag is less than 5.0°. The results are both in good agreement with the observed M2 tide in the Bohai Sea and the Yellow Sea. For comparison, the traditional methods also have been used to simulate M2 tide in the Bohai Sea and the Yellow Sea. The initial guess values of the boundary conditions are given first, and then are adjusted to acquire the simulated results that are as close as possible to the observations. As the boundary conditions contain 72 values, which should be adjusted and how to adjust them can only be partially solved by adjusting them many times. The satisfied results are hard to acquire even gigantic efforts are done. Here, the automation of the treatment of the open boundary conditions is realized. The method is unique and superior to the traditional methods. It is emphasized that if the adjoint of equation is used, tedious and complicated mathematical deduction can be avoided. Therefore the adjoint of equation should attract much attention.
基金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.
文摘This paper mainly discusses the openness and assimilation of the words system from the viewpoint of loan words between English and Chinese.Both the history of loan words and the ways of borrowing words in English and Chinese have shown us the strong abilities of languages.
文摘In order to solve the difficult problem of typhoon track prediction due to the sparsity of conventional data over the tropical ocean, in this paper, the No. 0205 typhoon Rammasun of 4-6 July 2002 is studied and an experiment of the typhoon track prediction is made with the direct use of the Advanced TIROS-N Operational Vertical Sounder (ATOVS) microwave radiance data in three-dimensional variational data assimilation. The prediction result shows that the experiment with the ATOVS microwave radiance data can not only successfully predict the observed fact that typhoon Rammasun moves northward and turns right, but can also simulate the action of the fast movement of the typhoon, which cannot be simulated with only conventional radiosonde data. The skill of the typhoon track prediction with the ATOVS microwave radiance data is much better than that without the ATOVS data. The typhoon track prediction of the former scheme is consistent in time and in location with the observation. The direct assimilation of
基金The project benefited from a PROCORE Hong Kong-France exchange grant to Arndt and Zhou and a grant from the US National Science Foundation
文摘The Panzhihua intrusion in southwest China is part of the Emeishan Large Igneous Province and host of a large Fe-Ti-V ore deposit.During emplacement of the main intrusion,multiple generations of mafic dykes invaded carbonate wall rocks,producing a large contact aureole.We measured the oxygen-isotope composition of the intrusions,their constituent minerals,and samples of the country rock.Magnetite and plagioclase from Panzhihua intrusion haveδ18O values that are consistent with magmatic equilibrium, and formed from magmas withδ18O values that were 1-2‰higher than expected in a mantle-derived magma.The unmetamorphosed country rock has highδ18O values,ranging from 13.2‰(sandstone) to 24.6-28.6‰(dolomite).The skarns and marbles from the aureole have lowerδ18O andδ13C values than their protolith suggesting interaction with fluids that were in exchange equilibrium with the adjacent mafic magmas and especially the numerous mafic dykes that intruded the aureole.This would explain the alteration ofδ18O of the dykes which have significantly higher values than expected for a mantle-derived magma.Depending on the exactδ18O values assumed for the magma and contaminant, the amount of assimilation required to produce the elevatedδ18O value of the Panzhihua intrusion was between 8 and 13.7 wt.%,assuming simple mixing.The exact mechanism of contamination is unclear but may involve a combination of assimilation of bulk country rock,mixing with a melt of the country rock and exchange with CO2-rich fluid derived from decarbonation of the marls and dolomites.These mechanisms,particularly the latter,were probably involved in the formation of the Fe-Ti-V ores.