CAS-ESM2.0 Dataset for the Carbon Dioxide Removal Model Intercomparison Project(CDRMIP)

Understanding the response of the Earth system to varying concentrations of carbon dioxide(CO_(2))is critical for projecting possible future climate change and for providing insight into mitigation and adaptation strategies in the near future.In this study,we generate a dataset by conducting an experiment involving carbon dioxide removal(CDR)—a potential way to suppress global warming—using the Chinese Academy of Sciences Earth System Model version 2.0(CASESM2.0).A preliminary evaluation is provided.The model is integrated from 200–340 years as a 1%yr^(−1) CO_(2) concentration increase experiment,and then to~478 years as a carbon dioxide removal experiment until CO_(2) returns to its original value.Finally,another 80 years is integrated in which CO_(2) is kept constant.Changes in the 2-m temperature,precipitation,sea surface temperature,ocean temperature,Atlantic meridional overturning circulation(AMOC),and sea surface height are all analyzed.In the ramp-up period,the global mean 2-m temperature and precipitation both increase while the AMOC weakens.Values of all the above variables change in the opposite direction in the ramp-down period,with a delayed peak relative to the CO_(2) peak.After CO_(2) returns to its original value,the global mean 2-m temperature is still~1 K higher than in the original state,and precipitation is~0.07 mm d^(–1) higher.At the end of the simulation,there is a~0.5°C increase in ocean temperature and a 1 Sv weakening of the AMOC.Our model simulation produces similar results to those of comparable experiments previously reported in the literature.

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL)is still not well characterized,especially over complex underlying surfaces.Herein,gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL.As for heights above the urban canopy layer(UCL),the relationship between turbulence velocity scale(VTKE)and wind speed(V)was consistent with the“HOckey-Stick”(HOST)theory proposed for a relatively flat area.Four regimes have been identified according to urban nocturnal stable boundary layer.Regime 1 occurs where local shear plays a leading role for weak turbulence under the constraint that the wind speed V<VT(threshold wind speed).Regime 2 is determined by the existence of strong turbulence that occurs when V>VT and is mainly driven by bulk shear.Regime 3 is identified by the existence of moderate turbulence when upside-down turbulence sporadic bursts occur in the presence of otherwise weak turbulence.Regime 4 is identified as buoyancy turbulence,when V>VT,and the turbulence regime is affected by a combination of local wind shear,bulk shear and buoyancy turbulence.The turbulence activities demonstrated a weak thermal stratification dependency in regime 1,for which within the UCL,the turbulence intensity was strongly affected by local wind shear when V<VT.This study further showed typical examples of different stable boundary layers and the variations between turbulence regimes by analyzing the evolution of wind vectors.Partly because of the influence of large-scale motions,the power spectral density of vertical velocity for upsidedown structure showed an increase at low frequencies.The upside-down structures were also characterized by the highest frequency of the stable stratifications in the higher layer.

CAS-ESM2.0 Model Datasets for the CMIP6 Ocean Model Intercomparison Project Phase 1 (OMIP1)

As a member of the Chinese modeling groups,the coupled ocean-ice component of the Chinese Academy of Sciences’Earth System Model,version 2.0(CAS-ESM2.0),is taking part in the Ocean Model Intercomparison Project Phase 1(OMIP1)experiment of phase 6 of the Coupled Model Intercomparison Project(CMIP6).The simulation was conducted,and monthly outputs have been published on the ESGF(Earth System Grid Federation)data server.In this paper,the experimental dataset is introduced,and the preliminary performances of the ocean model in simulating the global ocean temperature,salinity,sea surface temperature,sea surface salinity,sea surface height,sea ice,and Atlantic Meridional Overturning Circulation(AMOC)are evaluated.The results show that the model is at quasi-equilibrium during the integration of 372 years,and performances of the model are reasonable compared with observations.This dataset is ready to be downloaded and used by the community in related research,e.g.,multi-ocean-sea-ice model performance evaluation and interannual variation in oceans driven by prescribed atmospheric forcing.

CAS-ESM2.0 Model Datasets for the CMIP6 Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP)

The second version of the Chinese Academy of Sciences Earth System Model(CAS-ESM2.0)is participating in the Flux-Anomaly-Forced Model Intercomparison Project(FAFMIP)experiments in phase 6 of the Coupled Model Intercomparison Project(CMIP6).The purpose of FAFMIP is to understand and reduce the uncertainty of ocean climate changes in response to increased CO2 forcing in atmosphere-ocean general circulation models(AOGCMs),including the simulations of ocean heat content(OHC)change,ocean circulation change,and sea level rise due to thermal expansion.FAFMIP experiments(including faf-heat,faf-stress,faf-water,faf-all,faf-passiveheat,faf-heat-NA50pct and faf-heat-NA0pct)have been conducted.All of the experiments were integrated over a 70-year period and the corresponding data have been uploaded to the Earth System Grid Federation data server for CMIP6 users to download.This paper describes the experimental design and model datasets and evaluates the preliminary results of CAS-ESM2.0 simulations of ocean climate changes in the FAFMIP experiments.The simulations of the changes in global ocean temperature,Atlantic Meridional Overturning Circulation(AMOC),OHC,and dynamic sea level(DSL),are all reasonably reproduced.

Ocean Response to a Climate Change Heat-Flux Perturbation in an Ocean Model and Its Corresponding Coupled Model

State-of-the-art coupled general circulation models(CGCMs)are used to predict ocean heat uptake(OHU)and sealevel change under global warming.However,the projections of different models vary,resulting in high uncertainty.Much of the inter-model spread is driven by responses to surface heat perturbations.This study mainly focuses on the response of the ocean to a surface heat flux perturbation F,as prescribed by the Flux-Anomaly-Forced Model Intercomparison Project(FAFMIP).The results of ocean model were compared with those of a CGCM with the same ocean component.On the global scale,the changes in global mean temperature,ocean heat content(OHC),and steric sea level(SSL)simulated in the OGCM are generally consistent with CGCM simulations.Differences in changes in ocean temperature,OHC,and SSL between the two models primarily occur in the Arctic and Atlantic Oceans(AA)and the Southern Ocean(SO)basins.In addition to the differences in surface heat flux anomalies between the two models,differences in heat exchange between basins also play an important role in the inconsistencies in ocean climate changes in the AA and SO basins.These discrepancies are largely due to both the larger initial value and the greater weakening change of the Atlantic meridional overturning circulation(AMOC)in CGCM.The greater weakening of the AMOC in the CGCM is associated with the atmosphere–ocean feedback and the lack of a restoring salinity boundary condition.Furthermore,differences in surface salinity boundary conditions between the two models contribute to discrepancies in SSL changes.

A high-top version of IAP-AGCM:Preliminary assessment and sensitivity IAP-AGCM

Extending the atmospheric model top to high altitude is important for simulation of upper atmospheric phenomena,such as the stratospheric quasi-biennial oscillation.The high-top version of the Institute of Atmospheric Physics Atmospheric General Circulation Model with 91 vertical layers(IAP-AGCML91)extends to the mesopause at about 0.01 hPa(~80 km).The high-top model with a fully resolved stratosphere is found to simulate a warmer stratosphere than the low-top version,except near the South Pole,thus reducing its overall cold bias in the stratosphere,and significantly in the upper stratosphere.This sensitivity is shown to be consistent with two separate mechanisms:larger shortwave heating and larger poleward stratospheric meridional eddy heat flux in the hightop model than in the low-top model.Results indicate a significant influence of vertical resolution and model top on climate simulations in IAP-AGCM.

Simulation of the QBO in IAP-AGCM:Analysis of momentum budget

The quasi-biennial oscillation(QBO),a dominant mode of the equatorial stratospheric(~100–1 hPa)variability,is known to impact tropospheric circulation in the middle and high latitudes.Yet,its realistic simulation in general circulation models remains a challenge.The authors examine the simulated QBO in the 69-layer version of the Institute of Atmospheric Physics Atmospheric General Circulation Model(IAP-AGCML69)and analyze its momentum budget.The authors find that the QBO is primarily caused by parameterized gravity-wave forcing due to tropospheric convection,but the downward propagation of the momentum source is significantly offset by the upward advection of zonal wind by the equatorial upwelling in the stratosphere.Resolved-scale waves act as a positive contribution to the total zonal wind tendency of the QBO over the equator with comparable magnitude to the gravity-wave forcing in the upper stratosphere.Results provide insights into the mechanism of the QBO and possible causes of differences in models.

China’s EarthLab-Forefront of Earth System Simulation Research

This article introduces“EarthLab”,a major new Earth system numerical simulation facility developed in China.EarthLab is a numerical simulation system for a physical climate system,an environmental system,an ecological system,a solid earth system,and a space weather system as a whole with a high-performance scientific computing platform.EarthLab consists of five key elements-namely:a global earth numerical simulation system,a regional high-precision simulation system,a supercomputing support and management system,a database,data assimilation and visualization system,and a high-performance computing system for earth sciences.EarthLab helps to study the atmosphere,hydrosphere,cryosphere,lithosphere,and biosphere,as well as their interactions,to improve the accuracy of predictions by integrating simulations and observations,and to provide a scientific foundation for major issues such as national disaster prevention and mitigation.The construction and operation of EarthLab will involve close cooperation with joint contributions and shared benefits.

The source of Double ITCZ induced by the SST bias over the tropical western Pacific as reflected in CAS-ESM2 Model

“Double ITCZ”is a common precipitation bias over the tropical Pacific in current climate models and Earth system models,but the reasons for its formation are still worth exploring and discussing.In this study,we adopted the second version of Chinese Academy of Sciences Earth System Model(CAS-ESM2),by comparing a set of sea surface temperature(SST)bias correction experiments over the tropical Pacific,to explore the possible mechanism of SST bias on inducing the“Double ITCZ”from the perspective of the climatic annual mean bias in the coupled model.We revealed that the simulated climatic annual mean SST bias over the tropical Pacific can affect the bias of latent heat flux through the saturation specific humidity,and the bias of latent heat flux can further affect that of vertical velocity of humid air by the condensation release mechanism,and finally modulate the simulated bias in precipitation.Furthermore,through the inter-comparison between different experiments,it is found that the source of Double ITCZ might mainly come from the annual mean SST bias over the tropical western Pacific through the proposed air-sea coupled process of“SST-saturated specific humidity-latent heat flux-vertical velocity-precipitation”,indicating a possible way on reducing the couple biases in models over the tropical Pacific to improve the accuracy of CAS-ESM2 for climate simulation.

An Overview of the Dynamic Framework in Earth-System Model and Its Well-Posedness

The well-posedness of the dynamic framework in earth-system model(ESM for short)is a common issue in earth sciences and mathematics.In this paper,the authors first introduce the research history and fundamental roles of the well-posedness of the dynamic framework in the ESM,emphasizing the three core components of ESM,i.e.,the atmospheric general circulation model(AGCM for short),land-surface model(LSM for short)and oceanic general circulation model(OGCM for short)and their couplings.Then,some research advances made by their own research group are outlined.Finally,future research prospects are discussed.

Stochastic method to determine the scale and anomalous diffusion of gusts in a windy atmospheric boundary layer

In the atmospheric boundary layer, especially during strong wind period, the coherent structures are obvious and related to the direct interaction of the air masses with the ground. In this paper, we used the observation data during dust weather in Northwest Gansu to study the coherent structure and their ‘‘anomalous diffusion''. The structures in the atmospheric boundary layer included turbulent fluctuations and gusty wind disturbances, and could be denoted as ‘‘critical events' '. Their fractal dimensions were expressed by the complex index l of waiting times. Although the complex index can indicate the ability of the system to generate coherent structures, it has a strong dependence on the threshold marking the‘‘critical events' '. Hence, the continuous time random walk method was used to analyze the coherent structures. The scaling law of anomalous diffusion of coherent structures was obtained, and the diffusion scaling exponent H that indicated the ability of diffusion of different structures was analyzed. The exponents changed with structure scales which were affected by velocities and heights. At small scales, it was almost isotropic, and at large scales, the coherent structures were obvious and the diffusion was anomalous.

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

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

Optimal reduction of anthropogenic emissions for air pollution control and the retrieval of emission source from observed pollutants II: Iterative optimization using a positive-negative discriminant

In part Ⅰ of this paper series, the application of an incomplete adjoint operator to calculate the optimal reduction of the total emissions S was suggested. This paper, part Ⅱ of the series, focuses on calculating the reduction of the anthropogenic emission source S. As aSc(the source term due to the chemical reaction in the atmosphere) is a complex function of the pollutant concentration, we propose an iterative optimization method using a positive-negative discriminant to obtain Sfrom he aS. Tconvergence of the iterations is also proven in this paper.

Optimal reduction of anthropogenic emissions for air pollution control and the retrieval of emission source from observed pollutantsⅢ:Emission source inversion using a double correction iterative method

Using the incomplete adjoint operator method in part I of this series of papers,the total emission source S can be retrieved from the pollutant concentrationsρob obtained from the air pollution monitoring network.This paper studies the problem of retrieving anthropogenic emission sources from S.Assuming that the natural source Sn is known,and as the internal source Sc due to chemical reactions is a function of pollutant concentrations,if the chemical reaction equations are complete and the parameters are accurate,Sc can be calculated directly fromρob,and then Sa can be obtained from S.However,if the chemical reaction parameters(denoted asγ)are insufficiently accurate,bothγand Sc should be corrected.This article proposes a"double correction iterative method"to retrieve Sc and correctγand proves that this iterative method converges.