A Neural-network-based Alternative Scheme to Include Nonhydrostatic Processes in an Atmospheric Dynamical Core

Here,a nonhydrostatic alternative scheme(NAS)is proposed for the grey zone where the nonhydrostatic impact on the atmosphere is evident but not large enough to justify the necessity to include an implicit nonhydrostatic solver in an atmospheric dynamical core.The NAS is designed to replace this solver,which can be incorporated into any hydrostatic models so that existing well-developed hydrostatic models can effectively serve for a longer time.Recent advances in machine learning(ML)provide a potential tool for capturing the main complicated nonlinear-nonhydrostatic relationship.In this study,an ML approach called a neural network(NN)was adopted to select leading input features and develop the NAS.The NNs were trained and evaluated with 12-day simulation results of dry baroclinic-wave tests by the Weather Research and Forecasting(WRF)model.The forward time difference of the nonhydrostatic tendency was used as the target variable,and the five selected features were the nonhydrostatic tendency at the last time step,and four hydrostatic variables at the current step including geopotential height,pressure in two different forms,and potential temperature,respectively.Finally,a practical NAS was developed with these features and trained layer by layer at a 20-km horizontal resolution,which can accurately reproduce the temporal variation and vertical distribution of the nonhydrostatic tendency.Corrected by the NN-based NAS,the improved hydrostatic solver at different horizontal resolutions can run stably for at least one month and effectively reduce most of the nonhydrostatic errors in terms of system bias,anomaly root-mean-square error,and the error of the wave spatial pattern,which proves the feasibility and superiority of this scheme.

The Asymmetric Effects of El Ni?o and La Ni?a on the East Asian Winter Monsoon and Their Simulation by CMIP5 Atmospheric Models

E1 Nifio-Southem Oscillation （ENSO） events significantly affect the year-by-year variations of the East Asian winter monsoon （EAWM）. However, the effect of La Nifia events on the EAWM is not a mirror image of that of E1 Nifio events. Although the EAWM becomes generally weaker during El Nifio events and stronger during La Nifia winters, the enhanced precipitation over the southeastern China and warmer surface air temperature along the East Asian coastline during E1 Nifio years are more significant. These asymmetric effects are caused by the asymmetric longitudinal positions of the western North Pacific （WNP） anticyclone during El Nifio events and the WNP cyclone during La Nifia events; specifically, the center of the WNP cyclone during La Nifia events is westward-shifted relat- ive to its El Nifio counterpart. This central-position shift results from the longitudinal shift of remote E1 Nifio and La Nifia anomalous heating, and asymmetry in the amplitude of local sea surface temperature anomalies over the WNP. However, such asymmetric effects of ENSO on the EAWM are barely reproduced by the atmospheric models of Phase 5 of the Coupled Model Intercomparison Project （CMIP5）, although the spatial patterns of anomalous circula- tions are reasonably reproduced. The major limitation of the CMIP5 models is an overestimation of the anomalous WNP anticyclone/cyclone, which leads to stronger EAWM rainfall responses. The overestimated latent heat flux an- omalies near the South China Sea and the northern WNP might be a key factor behind the overestimated anomalous circulations.

Boreal Winter Rainfall Anomaly over the Tropical Indo-Pacific and Its Effect on Northern Hemisphere Atmospheric Circulation in CMIP5 Models

Experimental outputs of 11 Atmospheric Model Intercomparison Project （AMIP） models from phase 5 of the Coupled Model Intercomparison Project （CMIP5） are analyzed to assess the atmospheric circulation anomaly over Northern Hemisphere induced by the anomalous rainfall over tropical Pacific and Indian Ocean during boreal winter.The analysis shows that the main features of the interannual variation of tropical rainfall anomalies,especially over the Central Pacific （CP） （5°S-5°N,175°E-135°W） and Indo-western Pacific （IWP） （20°S-20°N,110°-150°E） are well captured in all the CMIP5/AMIP models.For the IWP and western Indian Ocean （WIO） （10°S-10°N,45°-75°E）,the anomalous rainfall is weaker in the 11 CMIP5/AMIP models than in the observation.During El Ni（n）o/La Ni（n）a mature phases in boreal winter,consistent with observations,there are geopotential height anomalies known as the Pacific North American （PNA） pattern and Indo-western Pacific and East Asia （IWPEA） pattern in the upper troposphere,and the northwestern Pacific anticyclone （cyclone） （NWPA） in the lower troposphere in the models.Comparison between the models and observations shows that the ability to simulate the PNA and NWPA pattern depends on the ability to simulate the anomalous rainfall over the CP,while the ability to simulate the IWPEA pattern is related to the ability to simulate the rainfall anomaly in the IWP and WIO,as the SST anomaly is same in AMIP experiments.It is found that the tropical rainfall anomaly is important in modeling the impact of the tropical Indo-Pacific Ocean on the extratropical atmospheric circulation anomaly.

New Developments on Existence and Uniqueness of Solutions to Some Models in Atmospheric Dynamics

This survey is concerned with the new developments on existence and uniqueness of solutions of some basic models in atmospheric dynamics, such as two-and three-dimensional quasi-geostrophic models and three-dimensional balanced model. The main aim of this paper is to introduce some results about the global and local (with respect to time) existence of solutions given by the authors in recent years, but others' important contributions and the literature on this subject are also quoted. We discuss briefly the relationships among the existence and uniqueness, physical instability and computational instability. In the appendixes, some key mathematical techniques in obtaining our results are presented, which are of vital importance to other problems in geophysical fluid dynamics as well.

ON THE VERTICAL DISCRETIZATION OF ATMOSPHERIC MODELS—PART I:EXISTING PROBLEMS,DATA PROCESSING AND CHOICE OF VERTICAL COORDINATE AND BASIS FUNCTIONS

Discussions are carried out on the vertical discretization of current atmospheric models.It is pointed out that there exist problems in the integration of the hydrostatic equation and the computation of vertical advection,vertical diffusion and so on.Then some possible ways for solving or alleviating them are suggested.Finally,the choice of vertical coordinate and basis functions is discussed.

Improved Diurnal Cycle of Precipitation on Land in a Global Non-Hydrostatic Model Using a Revised NSAS Deep Convective Scheme

In relatively coarse-resolution atmospheric models,cumulus parameterization helps account for the effect of subgridscale convection,which produces supplemental rainfall to the grid-scale precipitation and impacts the diurnal cycle of precipitation.In this study,the diurnal cycle of precipitation was studied using the new simplified Arakawa-Schubert scheme in a global non-hydrostatic atmospheric model,i.e.,the Yin-Yang-grid Unified Model for the Atmosphere.Two new diagnostic closures and a convective trigger function were suggested to emphasize the job of the cloud work function corresponding to the free tropospheric large-scale forcing.Numerical results of the 0.25-degree model in 3-month batched real-case simulations revealed an improvement in the diurnal precipitation variation by using a revised trigger function with an enhanced dynamical constraint on the convective initiation and a suitable threshold of the trigger.By reducing the occurrence of convection during peak solar radiation hours,the revised scheme was shown to be effective in delaying the appearance of early-afternoon rainfall peaks over most land areas and accentuating the nocturnal peaks that were wrongly concealed by the more substantial afternoon peak.In addition,the revised scheme enhanced the simulation capability of the precipitation probability density function,such as increasing the extremely low-and high-intensity precipitation events and decreasing small and moderate rainfall events,which contributed to the reduction of precipitation bias over mid-latitude and tropical land areas.

The Performance of Atmospheric Component Model R42L9 of GOALS/LASG

This paper examines the performance of an atmospheric general circulation model (AGCM) developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics (LASG/IAP). It is a spectral model truncated at R42(2.8125°long×1.66°lat) resolution and with nine vertical levels, and referred to as R42L9/LASG hereafter. It is also the new version of atmospheric component model R15L9 of the global ocean-atmosphere-land system (GOALS/LASG). A 40-year simulation in which the model is forced with the climatological monthly mean sea surface temperature is compared with the 40-year (1958-97) U.S. National Center for Environmental Prediction (NGEP) global reanalysis and the 22-year (1979-2000) Xie-Arkin monthly precipitation climatology. The mean DJF and JJA geographical distributions of precipitation, sea level pressure, 500-hPa geopotential height, 850-hPa and 200-hPa zonal wind, and other fields averaged for the last 30-year integration of the R42L9 model are analyzed. Results show that the model reproduces well the observed basic patterns, particularly precipitation over the East Asian region. Comparing the new model with R15L9/LASG, the old version with coarse resolution (nearly 7.5°long×4.5°lat), shows an obvious improvement in the simulation of regional climate, especially precipitation. The weaknesses in simulation and future improvements of the model are also discussed.

Dynamical Framework of IAP Nine-Level Atmospheric General Circulation Model

The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.

Prediction of the Asian-Australian Monsoon Interannual Variations with the Grid-Point Atmospheric Model of IAP LASG(GAMIL)

Seasonal prediction of Asian-Australian monsoon （A-AM） precipitation is one of the most important and challenging tasks in climate prediction. In this paper, we evaluate the performance of Grid Atmospheric Model of IAP LASG （GAMIL） on retrospective prediction of the A-AM interannual variation （IAV）, and determine to what extent GAMIL can capture the two major observed modes of A-AM rainfall IAV for the period 1979-2003. The first mode is associated with the turnabout of warming （cooling） in the Nifio 3.4 region, whereas the second mode leads the warming/cooling by about one year, signaling precursory conditions for ENSO. We show that the GAMIL one-month lead prediction of the seasonal precipitation anomalies is primarily able to capture major features of the two observed leading modes of the IAV, with the first mode better predicted than the second. It also depicts the relationship between the first mode and ENSO rather well. On the other hand, the GAMIL has deficiencies in capturing the relationship between the second mode and ENSO. We conclude： （1） successful reproduction of the E1 Nifio-excited monsoon-ocean interaction and E1 Nifio forcing may be critical for the seasonal prediction of the A-AM rainfall IAV with the GAMIL; （2） more efforts are needed to improve the simulation not only in the Nifio 3.4 region but also in the joining area of Asia and the Indian-Pacific Ocean; （3） the selection of a one-tier system may improve the ultimate prediction of the A-AM rainfall IAV. These results offer some references for improvement of the GAMIL and associated seasonal prediction skill.

Spatial and seasonal patterns of atmospheric nitrogen deposition in North China

As one of the regions with intensive agriculture and rapid economic development in China,North China also has a high nitrogen(N)deposition.This study characterized the spatial pattern of N deposition in North China,combining the tropospheric columns from satellite measurements and the simulated profiles from an atmospheric chemistry transport model.The total N deposition fluxes ranged from 16.3 to 106.5 kg N ha−1 yr−1,with an average of 54.5±17.2 kg N ha−1 yr−1.The high values were concentrated in urban and farmland areas,while low values were found in forests and grasslands with less human activities.Of the total N deposition,36%was deposited via precipitation,12%was deposited through dry particulate deposition,and the remaining 52%was comprised of dry gaseous deposition.For the seasonal variation of dry deposition,gaseous HNO3 and particulate NO3−were higher in winter and autumn,but lower in spring and summer.In contrast,gaseous NH3 and particulate NH4+were higher in spring and summer,but lower in winter and autumn.This is possibly caused by the seasonal differences in emission intensity between NOx and NH3 emission sources.The gaseous NO2 deposition did not show strong seasonal variation.The wet deposition was mainly affected by precipitation,with high values in summer and low values in winter.This research provides an objective spatial perspective and insight into the total N deposition in North China.

Analysis and Evaluation of the Global Aerosol Optical Properties Simulated by an Online Aerosol-coupled Non-hydrostatic Icosahedral Atmospheric Model

Aerosol optical properties are simulated using the Spectral Radiation Transport Model I~）r Aerosol Species （SPRINTARS） coupled with the Non-hydrostatic ICosahedral Atmospheric Model （NICAM）. The 3-year global mean all-sky aerosol optical thickness （AOT） at 550 nm, theAngstr/Sm Exponent （AE） based on AOTs at 440 and 870 nm, and the single scattering albedo （SSA） at 550 nm are estimated at 0.123, 0.657 and 0.944, respectively. For each aerosol species, the mean AOT is within the range of the AeroCom models. Both the modeled all-sky and clear-sky results are compared with observations from the Moderate Resolution Imaging Spectroradiometer （MODIS） and the Aerosol Robotic Network （AERONET）. The simulated spatiotemporal distributions of all-sky AOTs can generally reproduce the MODIS retrievals, and the correlation and model skill can be slightly improved using the clear-sky results over most land regions. The differences between clear-sky and all-sky AOTs are larger over polluted regions. Compared with observations from AERONET, the modeled and observed all-sky AOTs and AEs are generally in reasonable agreement, whereas the SSA variation is not well captured. Although the spatiotemporal distributions of all-sky and clear-sky results are similar, the clear-sky results are generally better correlated with the observations. The clear-sky AOT and SSA are generally lower than the all-sky results, especially in those regions where the aerosol chemical composition is contributed to mostly by sulfate aerosol. The modeled clear-sky AE is larger than the all-sky AE over those regions dominated by hydrophilic aerosol, while the＇opposite is found over regions dominated by hydrophobic aerosol.

Basic Pattern in Atmospheric Turbulence Model

From the controlling equations of atmosphere motion, Prandtl's mixing length theory is used to derive the atmospheric turbulence models, such as Burgers equation model and Burgers-KdV equation model. And then the projective Riccati equations are applied to solve these atmospheric turbulence models, where much more patterns are obtained, including solitary wave pattern, singular pattern, and so on.

Generative adversarial network-based atmospheric scattering model for image dehazing

This paper presents a trainable Generative Adversarial Network(GAN)-based end-to-end system for image dehazing,which is named the DehazeGAN.DehazeGAN can be used for edge computing-based applications,such as roadside monitoring.It adopts two networks:one is generator(G),and the other is discriminator(D).The G adopts the U-Net architecture,whose layers are particularly designed to incorporate the atmospheric scattering model of image dehazing.By using a reformulated atmospheric scattering model,the weights of the generator network are initialized by the coarse transmission map,and the biases are adaptively adjusted by using the previous round's trained weights.Since the details may be blurry after the fog is removed,the contrast loss is added to enhance the visibility actively.Aside from the typical GAN adversarial loss,the pixel-wise Mean Square Error(MSE)loss,the contrast loss and the dark channel loss are introduced into the generator loss function.Extensive experiments on benchmark images,the results of which are compared with those of several state-of-the-art methods,demonstrate that the proposed DehazeGAN performs better and is more effective.

Two-moment Bulk Stratiform Cloud Microphysics in the Grid-point Atmospheric Model of IAP LASG (GAMIL)

A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LASG （GAMIL） as an effort to enhance the model＇s capability to simulate aerosol indirect effects. Unlike the previous one-moment cloud microphysics scheme, the new scheme produces a reasonable rep- resentation of cloud particle size and number concentration. This scheme captures the observed spatial variations in cloud droplet number concentrations. Simulated ice crystal number concentrations in cirrus clouds qualitatively agree with in situ observations. The longwave and shortwave cloud forcings are in better agreement with observations. Sensitivity tests show that the column cloud droplet number concentrations calculated from two different droplet activation parameterizations are similar. However, ice crystal number concentration in mixed-phased clouds is sensitive to different heterogeneous ice nucleation formulations. The simulation with high ice crystal number concentration in mixed-phase clouds has less liquid water path and weaker cloud forcing. ~rthermore, ice crystal number concentration in cirrus clouds is sensitive to different ice nucleation parameterizations. Sensitivity tests also suggest that the impact of pre-existing ice crystals on homogeneous freezing in old clouds should be taken into account.

Evaluation and Improvement of a SVD-Based Empirical Atmospheric Model

An empirical atmospheric model(EAM) based on the singular value decomposition(SVD) method is evaluated using the composite El Ni(?)o/Southern Oscillation(ENSO) patterns of sea surface temperature (SST) and wind anomalies as the target scenario.Two versions of the SVD-based EAM were presented for comparisons.The first version estimates the wind anomalies in response to SST variations based on modes that were calculated from a pair of global wind and SST fields(i.e.,conventional EAM or CEAM).The second version utilizes the same model design but is based on modes that were calculated in a region-wise manner by separating the tropical domain from the remaining extratropical regions(i.e.,region-wise EAM or REAM). Our study shows that,while CEAM has shown successful model performance over some tropical areas, such as the equatorial eastern Pacific(EEP),the western North Pacific(WNP),and the tropical Indian Ocean(TIO),its performance over the North Pacific(NP) seems poor.When REAM is used to estimate the wind anomalies instead of CEAM,a marked improvement over NP readily emerges.Analyses of coupled modes indicate that such an improvement can be attributed to a much stronger coupled variability captured by the first region-wise SVD mode at higher latitudes compared with that captured by the conventional one. The newly proposed way of constructing the EAM(i.e.,REAM) can be very useful in the coupled studies because it gives the model a wider application beyond the commonly accepted tropical domain.

The Investigation of Atmospheric Angular Momentum as a Contributor to Polar Wobble and Length of Day Change with AMIP II GCM Data

The atmospheric angular momentum (AAM) functions in terms of contribution to polar wobble and length of day change, are calculated from the output data of GSM9603 global circulation model (GCM) of Japan Meteorological Agency (JMA), from the reanalysis data of the National Centers for the Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR), and from the operational objective analysis data of JMA, respectively. The comparison shows that during the period from 1985 to 1995, the values of the pressure terms in the equatorial components of AAM functions calculated from three data sets agree with each other better along 90°E longitude than along Greenwich meridian direction. The axial component of relative AAM function estimated from GSM 9603 agrees well with those from the other two data sets in terms of seasonal variations with the moderate amplitudes, but not so well with the composite axial component of relative AAM functions estimated from 23 GCM models anticipating in the first phase of AMIP. In addition, its interannual variation from 1979 to 1996 shows the main characteristics of ENSO evolution, just as does the axial component of relative AAM function estimated from NCEP reanalysis data except for the period of anomalous ENSO from 1991 to 1993.

Study on the Total Amount Control of Atmospheric Pollutant Based on GIS

To provide effective environmental management for total amount control of atmospheric pollutants. Methods An atmospheric diffusion model of sulfur dioxide on the surface of the earth was established and tested in Shantou of Guangdong Province on the basis of an overall assessment of regional natural environment, social economic state of development, pollution sources and atmospheric environmental quality. Compared with actual monitoring results in a studied region, simulation values fell within the range of two times of error and were evenly distributed in the two sides of the monitored values. Predicted with the largest emission model method, the largest emission of sulfur dioxide would be 54 279.792 tons per year in 2010. Conclusion The mathematical model established and revised on the basis of GIS is more rational and suitable for the regional characteristics of total amount control of air pollutants.

Response of the North Pacific Oscillation to global warming in the models of the Intergovernmental Panel on Climate Change Fourth Assessment Report

Based on 22 of the climate models from phase 3 of the Coupled Model Intercomparison Project, we investigate the ability of the models to reproduce the spatiotemporal features of the wintertime North Pacific Oscillation(NPO), which is the second most important factor determining the wintertime sea level pressure field in simulations of the pre-industrial control climate, and evaluate the NPO response to the future most reasonable global warming scenario(the A1B scenario). We reveal that while most models simulate the geographic distribution and amplitude of the NPO pattern satisfactorily, only 13 models capture both features well. However, the temporal variability of the simulated NPO could not be significantly correlated with the observations. Further analysis indicates the weakened NPO intensity for a scenario of strong global warming is attributable to the reduced lower-tropospheric baroclinicity at mid-latitudes, which is anticipated to disrupt large-scale and low-frequency atmospheric variability, resulting in the diminished transfer of energy to the NPO, together with its northward shift.

The African Climate as Predicted by the IAP Grid-Point Nine-Layer Atmospheric General Circulation Model (IAP-9L-AGCM)

ＴｈｅＡｆｒｉｃａｎＣｌｉｍａｔｅａｓＰｒｅｄｉｃｔｅｄｂｙｔｈｅＩＡＰＧｒｉｄ－ＰｏｉｎｔＮｉｎｅ－ＬａｙｅｒＡｔｍｏｓｐｈｅｒｉｃＧｅｎｅｒａｌＣｉｒｃｕｌａｔｉｏｎＭｏｄｅｌ（ＩＡＰ－９Ｌ－ＡＧＣＭ）ＣｈｉｎｅｋｅＴｈｅｏＣｈｉｄｉｅｚｉｅ①，...

Aerosol Indirect Effects on Warm Clouds in the Grid-Point Atmospheric Model of IAP LASG(GAMIL)

Aerosol indirect effects on warm clouds are estimated in the Grid-point Atmospheric Model of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics,Institute of Atmospheric Physics(IAP LASG)(GAMIL) with a new two-moment cloud microphysics scheme using two different physically-based aerosol activation parameterizations:Abdul-Razzak and Ghan,and Nenes and Seinfeld.The annual global mean changes in shortwave cloud forcing from preindustrial times to present day(a measure of the aerosol indirect effects) estimated from these two parameterizations are remarkably similar:0.76 W m?2 with the Abdul-Razzak and Ghan parameterization,and 0.78 W m?2 with the Nenes and Seinfeld parameterization.Physically-based parameterizations can provide robust representations of aerosol effects on droplet nucleation,meaning that aerosol activation is no longer the most uncertain factor in modeling aerosol indirect effects.