INTERDECADAL VARIABILITY IN A MODEL ATMOSPHERE

By using the simulation results of an AGCM, which had been run from 1945 to 1993 forced by COADS SST, the interdecadal variability of the model atmosphere was investigated and compared with that of NCEP reanalysis data. It was found that, interdecadal variability exists significantly in both the tropical Pacific wind fields and the mid-high latitude atmospheric circulation of the model atmosphere. The tendency of time variation and spatial distributions of the interdecadal variability of the model atmosphere are basically consistent with observation. Relative to the mid-high latitude atmospheric circulation, the simulation of tropical Pacific wind is more satisfying, which suggests that anomalous variation of SST is still the main factor for the interdecadal variability of tropical Pacific wind. It might have more significant influence on the tropical wind than on the mid-high latitude atmosphere. However, there is still obvious difference between the simulation and observation. They could be attributed to both the simulation capability of the model and absence of other factors in the model which are important for the interdecadal climate variation.

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 Application of Flux-Form Semi-Lagrangian Transport Scheme in a Spectral Atmosphere Model

A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.

Decoy State Quantum Key Distribution via Beam-Wandering Modeled Atmosphere Channel

We investigate the decoy state quantum key distribution via the atmosphere channels. We consider the efficient decoy state method with one-signal state and two-decoy states. Our results show that the decoy state method works even in the channels with fluctuating transmittance. Nevertheless, the key generation rate will be dra-matically decreased by atmosphere turbulence, which sheds more light on the characterization of atmosphere turbulence in realistic free-space based quantum key distributions.

Impacts of Vegetation on the Intraseasonal Oscillation Simulated by the Community Atmosphere Model (CAM3)

The influences of vegetation on intraseasonal oscillation （ISO） were examined using the Community Atmosphere Model version 3 （CAM3）. Two 15-year numerical experiments were completed： the first was performed with a realistic vegetation distribution （VEG run）, and the second was identical to the VEG run except without land vegetation （NOVEG run）. Generally speak- ing, CAM3 was able to reproduce the spatial distribution of the ISO, but the ISO intensity in the simulation was much weaker than that observed in nature： the 1SO has a relatively much stronger signal. A comparison of the VEG run with the NOVEG run revealed that the presence of vegetation usually produces a weak ISO. The vegetation effects on ISO intensity were significant over West Africa and South Asia, especially in the summer half-year. Vegetation also plays an important role in modulating ISO propagation. The eastward propagation of the ISO in the VEG run was clearer than that in the NOVEG run over the West African and Maritime Continent regions. The northward propagation of the ISO in the VEG run was more consistent with observation than that in the NOVEG run.

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.

Continuous Variable Quantum Teleportation in Beam-Wandering Modeled Atmosphere Channel

We investigate the continuous variable quomtum teleportation in atmosphere channels. The beam-wandering mode/is employed to analyze the teleportation of the unknown single-mode coherent state. Two methods, one is deterministic by increasing the aperture size of the detecting device and one is probabilistic by entanglement distillation, are proposed to improve the teleportation fidelity in the presence of atmosphere noises.

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.

Image Dehazing with Hybrid λ2-λ0 Penalty Mode

Due to the presence of turbid media, such as microdust and water vapor in the environment, outdoor pictures taken under hazy weather circumstances are typically degraded. To enhance the quality of such images, this work proposes a new hybrid λ2-λ0 penalty model for image dehazing. This model performs a weighted fusion of two distinct transmission maps, generated by imposing λ2 and λ0 norm penalties on the approximate regression coefficients of the transmission map. This approach effectively balances the sparsity and smoothness associated with the λ0 and λ2 norms, thereby optimizing the transmittance map. Specifically, when the λ2 norm is penalized in the model, an updated guided image is obtained after implementing λ0 penalty. The resulting optimization problem is effectively solved using the least square method and the alternating direction algorithm. The dehazing framework combines the advantages of λ2 and λ0 norms, enhancing sparse and smoothness, resulting in higher quality images with clearer details and preserved edges.

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.

A Study on Sulfate Optical Properties and Direct Radiative Forcing Using LASG-IAP General Circulation Model

The direct radiative forcing （DRF） of sulfate aerosols depends highly on the atmospheric sulfate loading and the meteorology, both of which undergo strong regional and seasonal variations. Because the optical properties of sulfate aerosols are also sensitive to atmospheric relative humidity, in this study we first examine the scheme for optical properties that considers hydroscopic growth. Next, we investigate the seasonal and regional distributions of sulfate DRF using the sulfate loading simulated from NCAR CAM-Chem together with the meteorology modeled from a spectral atmospheric general circulation model （AGCM） developed by LASG-IAP. The global annual-mean sulfate loading of 3.44 mg m-2 is calculated to yield the DRF of -1.03 and -0.57 W m-2 for clear-sky and all-sky conditions, respectively. However, much larger values occur on regional bases. For example, the maximum all-sky sulfate DRF over Europe, East Asia, and North America can be up to -4.0 W m-2. The strongest all-sky sulfate DRF occurs in the Northern Hemispheric July, with a hemispheric average of -1.26 W m-2. The study results also indicate that the regional DRF are strongly affected by cloud and relative humidity, which vary considerably among the regions during different seasons. This certainly raises the issue that the biases in model-sinmlated regional meteorology can introduce biases into the sulfate DRF. Hence, the model processes associated with atmospheric humidity and cloud physics should be modified in great depth to improve the simulations of the LASG-IAP AGCM and to reduce the uncertainty of sulfate direct effects on global and regional climate in these simulations.

An Introduction to the Coupled Model FGOALS1.1-s and Its Performance in East Asia

The spectral version 1.1 of the Flexible Global Ocean–atmosphere–land System （FGOALS1.1-s） model was developed in the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophys- ical Fluid Dynamics at the Institute of Atmospheric Physics （LASG/IAP）. This paper reports the major modifications to the physical parameterization package in its atmospheric component, including the radiation scheme, convection scheme, and cloud scheme. Furthermore, the simulation of the East Asian Summer Monsoon （EASM） by FGOALS1.1-s is examined, both in terms of climatological mean state and interannual variability. The results indicate that FGOALS1.1-s exhibits significant improvements in the simulation of the balance of energy at the top of the atmosphere： the net radiative energy flux at the top was 0.003 W m-2 in the 40 years fully coupled integration. The distribution of simulated sea surface temperature was also quite reasonable, without obvious climate drift. FGOALS1.1-s is also capable of capturing the major features of the climatological mean state of the EASM： major rainfall maximum centers, the annual cycle of precipitation, and the lower-level monsoon circulation flow were highly consistent with observations in the EASM region. Regarding interannual variability, simulation of the EASM leading patterns and their relationship with sea surface temperature was examined. The results show that FGOALS1.1-s can reproduce the first leading pattern of the EASM and its close relationship with the decaying phase of the ENSO. However, the model lacked the ability to capture either the second major mode of the EASM or its relationship with the developing phase of the ENSO.

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.

Effects of Nested Area Size upon Regional Climate Model Simulations

This paper presents a numerical study on the 1998 summer rainfall over the Yangtze River valley in central and eastern China, addressing effect of a nested area size on simulations in terms of the technique of nesting a regional climate model (RCM) upon a general circulation model (GCM). Evidence suggests that the size exerts greater impacts upon regional climate of the country, revealing that a larger nested size is superior to a small one for simulation in mitigating errors of GCM-provided lateral boundary forcing. Also, simulations show that the RCM should incorporate regions of climate systems of great importance into study and a low-resolution GCM yields more pronounced errors as a rule when used in the research of the Tibetan Plateau, and, in contrast, our PσRCM can do a good job in describing the plateau’s role in a more realistic and accurate way. It is for this reason that the tableland should be included in the nested area when the RCM is employed to investigate the regional climate. Our PσRCM nesting upon a GCM reaches more realistic results compared to a single GCM used.

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.

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.

Modelling the January and July Climate of 9000 Years before Present

The January and July climate of 9000 years before present was simulated through a set of perpetual experimentsby means of the newly designed 9 Level Atmospheric General Circulation Model (AGCM) in the institute of Atmospheric Physics (IAP). The results were analysed and compared with previous results simulated by the IAP 2Level AGCM. There exists good agreement between them. It is found that the temperature is higher in July and lower inJanuary in 9000 yBP than that at present. The temperature difference is more obvious in the Northern Hemispherethan in the Southern Hemisphere and greater in July than in January. These results prove the potential abilities of theg-L AGCM in the climate simulation and climate prediction.