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.

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.

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.

A Local Ensemble Transform Kalman Filter Data Assimilation System for the Global FSU Atmospheric Model

Data assimilation is the process by which measurements and model predictions are combined to obtain an accurate representation of the state of the modeled system. We implemented a data assimilation scheme called LETKF （local ensemble transform Kalman filter） with FSUGSM （Florida State University Global Spectral Model） and made an experiment to evaluate the initial condition generated to numerical weather prediction to FSUGSM model. The LETKF analysis carries out independently at each grid point with the use of ＂local＂ observations. An ensemble of estimates in state space represents uncertainty. The FSUGSM is a multilevel （27 vertical levels） spectral primitive equation model, where the variables are expanded horizontally in a truncated series of spherical harmonic functions （at resolution T63） and a transform technique is applied to calculate the physical processes in real space The assimilation cycle runs on the period 01/01/2001 to 31/01/2001 at （00, 06, 12 and 18 GMT） for each day. We examined the atmospheric fields during the period and the OMF （observation-minus-forecast） and the OMA （observation-minus-analysis） statistics to verify the analysis quality comparing with forecasts and observations. The analyses present stability and show suitable to initiate the weather predictions.

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.

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.

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 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.

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.

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.

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.

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.

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.

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

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

Estimate Totel Number of the Earth Atmospheric Particle with Standard Atmosphere model

The total number of atmospheric particle (AP) is an important datum for planetary science and geoscience. Estimating entire AP number is also a familiar question in general physics. With standard atmosphere model, considering the number difference of AP caused by rough and uneven in the earth surface below, the sum of dry clean atmosphere particle is . So the whole number of AP including water vapor is . The rough estimation for the total number of AP on other planets (or satellites) in condensed state is also discussed on the base of it.

A tool model for predicting atmospheric chemical kinetics with sensitivity analysis

A package(a tool model) for program of predicting atmospheric chemical kinetics with sensitivity analysis is presented. The new direct method of calculating the first order sensitivity coefficients using sparse matrix technology to chemical kinetics is included in the tool model, it is only necessary to triangularize the matrix related to the Jacobian matrix of the model equation. The Gear type procedure is used to integrate a model equation and its coupled auxiliary sensitivity coefficient equations. The FORTRAN subroutines of the model equation, the sensitivity coefficient equations, and their Jacobian analytical expressions are generated automatically from a chemical mechanism. The kinetic representation for the model equation and its sensitivity coefficient equations, and their Jacobian matrix is presented. Various FORTRAN subroutines in packages, such as SLODE, modified MA28, Gear package, with which the program runs in conjunction are recommended. The photo\|oxidation of dimethyl disulfide is used for illustration.

Comparison of some atmospheric chemical modeling schemes

The comparison of five atmospheric chemical modeling schemes was presented from accuracy,efficiency and progamming in term of Carbon bond Mechanism IV of atmospheric chemical reactions. The major results were as follows. The classical Gear scheme can provide an accurate solution and it is easy for programming by a computer automatically. The sparse matrix Gear type scheme can also provide an accurate solution but much faster than classical Gear scheme. QSSA,EBI and hybrid schemes can run with much longer time step without sacrificing of accuracy, therefore, much efficiently. If analytical solution is obtained by EBI scheme the accuracy and efficiency are much better.