Comparison Between GAMIL, and CAM2 on Interannual Variability Simulation

Recently, a new atmospheric general circulation model （GAMIL： Grid-point Atmospheric Model of IAP LASG） has been developed at the Institute of Atmospheric Physics （IAP）, Chinese Academy of Sciences （CAS）, which is based on the Community Atmospheric Model Version 2 （CAM2） of the National Center for Atmospheric Research （NCAR）. Since the two models have the same physical processes but different dynamical cores, the interannual variability simulation performances of the two models are compared. The ensemble approach is used to reduce model internal variability. In general, the simulation performances of the two models are similar. Both models have good per- formance in simulating total space-time variability and the Southern Oscillation Index. GAMIL performs better in the Eastern Asian winter circulation simulation than CAM2, and the model internal variability of GAMIL has a better response to external forcing than that of CAM2. These indicate that the improvement of the dynamic core is very important. It is also verified that there is less predictability in the middle and high latitudes than in the low latitudes.     

Parameterization of the Absorption of the H2O Continuum, CO2, O2, and Other Trace Gases in the Fu-Liou Solar Radiation Program

The absorption properties of the water vapor continuum and a number of weak bands for H2O, O2, CO2, CO, N2O, CH4, and O3 in the solar spectrum are incorporated into the Fu-Liou radiation parameterization program by using the correlated k-distribution method (CKD) for the sorting of absorption lines. The overlap absorption of the H2O lines and the H2O continuum (2500-14500 cm^-1) are treated by taking the two gases as a single-mixture gas in transmittance calculations. Furthermore, in order to optimize the computation efforts, CO2 and CH4 in the spectral region 2850-5250 cm^-1 are taken as a new singlemixture gas as well. For overlap involving other absorption lines in the Fu-Liou spectral bands, the authors adopt the multiplication rule for transmittance computations under which the absorption spectra for two gases are assumed to be uncorrelated. Compared to the line-by-line (LBL) computation, it is shown that the errors in fluxes introduced by these two approaches within the context of the CKD method are small and less than 0.48% for the H20 line and continuum in the 2500-14500 cm^-1 solar spectral region, -1% for H2O (line) H2O (continuum) CO2 CH4 in the spectral region 2850-5250 cm^-1, and -1.5% for H2O (line) H2O (continuum) O2 in the 7700-14500 cm^-1 spectral region. Analysis also demonstrates that the multiplication rule over a spectral interval as wide as 6800 cm^-1 can produce acceptable errors with a maximum percentage value of about 2% in reference to the LBL calculation. Addition of the preceding gases increases the absorption of solar radiation under all sky conditions. For clear sky, the increase in instantaneous solar absorption is about 9%-13% (-12 W m^-2) among which the H20 continuum produces the largest increase, while the contributions from O2 and CO2 rank second and third, respectively. In cloudy sky, the addition of absorption amounts to about 6-9 W m^-2. The new, improved program with the incorporation of the preceding gases produces a smaller solar absorption in clouds due to the reduced solar flux reaching the cloud top.     

Parameterization of the Absorption of the H_2O Continuum,CO_2,O_2,and Other Trace Gases in the Fu-Liou Solar Radiation Program

The absorption properties of the water vapor continuum and a number of weak bands for H2O, O2, CO2, CO, N2O, CH4, and O3 in the solar spectrum are incorporated into the Fu-Liou radiation parameterization program by using the correlated k-distribution method (CKD) for the sorting of absorption lines. The overlap absorption of the H2O lines and the H2O continuum (2500-14500 cm-1) are treated by taking the two gases as a single-mixture gas in transmittance calculations. Furthermore, in order to optimize the computation efforts, CO2 and CH4 in the spectral region 2850-5250 cm-1 are taken as a new single-mixture gas as well. For overlap involving other absorption lines in the Fu-Liou spectral bands, the authors adopt the multiplication rule for transmittance computations under which the absorption spectra for two gases are assumed to be uncorrelated. Compared to the line-by-line (LBL) computation, it is shown that the errors in fluxes introduced by these two approaches within the context of the CKD method are small and less than 0.48% for the H2O line and continuum in the 2500-14500 cm-1 solar spectral region, -1% for H2O (line) H2O (continuum) CO2 CH4 in the spectral region 2850-5250 cm-1, and -1.5% for H2O (line) H2O (continuum) O2 in the 7700-14500 cm-1 spectral region. Analysis also demonstrates that the multiplication rule over a spectral interval as wide as 6800 cm-1 can produce acceptable errors with a maximum percentage value of about 2% in reference to the LBL calculation. Addition of the preceding gases increases the absorption of solar radiation under all sky conditions. For clear sky, the increase in instantaneous solar absorption is about 9%-13% (~12 W m~2) among which the H2O continuum produces the largest increase, while the contributions from O2 and CO2 rank second and third, respectively. In cloudy sky, the addition of absorption amounts to about 6-9 W m-2. The new, improved program with the incorporation of the preceding gases produces a smaller solar absorption in clouds due to the reduced solar flux reaching the cloud top.     

The Influence of O3, Relative Humidity, NO and NO2 on the Oxidation of α-Pinene and Δ3-Carene

Upto 13% of -pinene and ³-carene had reacted after 213 s in this dark experimental set-up, where O₃, NO and NO₂ were mixed with terpenes at different relative humidities (RHs). The different experiments were planned according to an experimental design, where O₃, NO₂, NO, RH and reaction time were varied between high and low settings (25 and 75 ppb, 15 and 42%, 44 and 213 s). An increased amount of -pinene and ³-carene reacted in the chamber was observed, when the level of O₃, NO and reaction time was increased and RH was decreased. In the study, it was found that different interactions affected the amount of terpene reacted as well. These interactions were between O₃ and NO, O₃ and reaction time, NO and RH, and between NO and reaction time.     

Tropical Cyclones and Polar Lows: Velocity, Size, and Energy Scales, and Relation to the 26℃ Cyclone Origin Criteria

The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones (TC) and polar lows (PL) by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms: (1) as expressed through the buoyancy flux b and the Coriolis parameter lc for rotating fluids convection, and (2) as expresse...     

Diffuse radiation,twilight, and photochemistry — I

A photochemical scheme which includes a detailed treatment of multiple scattering up to solar zenith angles of 96° (developed for use in a GCM) has been used to study partitioning within chemical families. Attention is drawn to the different zenith angle dependence of diffuse radiation for the two spectral regions <310 nm and >310 nm. The effect that this has on the so-called 40 km ozone problem is discussed. The importance of correctly including multiple scattering for polar ozone studies is emphasised.     

Tropical cyclones and polar lows: Velocity,size, and energy scales,and relation to the 26°C cyclone origin criteria

The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones (TC) and polar lows (PL) by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms: (1) as expressed through the buoyancy flux b and the Coriolis parameter l_c for rotating fluids convection, and (2) as expressed with the cube of velocity times the drag coefficient through the formula for total kinetic energy dissipation in the atmospheric boundary layer. In the quasistationary case the dissipation equals the generation of the energy. In both cases the velocity scale can be expressed through temperature and humidity differences between the ocean and the atmosphere in terms of the reduced gravity, and both forms produce quite comparable velocity scales. Using parameters b and l_c, we can form scales of the area and, by adding the mass of a unit air column, a scale of the total kinetic energy as well. These scales nicely explain the much smaller size of a PL, as compared to a TC, and the total kinetic energy of a TC is of the order 10¹⁸-10¹⁹ J. It will be shown that wind of 33 m s^-1 is produced when the total enthalpy fluxes between the ocean and the atmosphere are about 700 W m^-2 for a TC and 1700 W m^-2 for a PL, in association with the much larger role of the latent heat in the first case and the stricter geostrophic constraints and larger static stability in the second case. This replaces the mystical role of 26^oC as a criterion for TC origin. The buoyancy flux, a product of the reduced gravity and the wind speed, together with the atmospheric static stability, determines the rate of the penetrating convection. It is known from the observations that the formation time for a PL reaching an altitude of 5--6 km can be only a few hours, and a day, or even half a day, for a TC reaching 15--18 km. These two facts allow us to construct curves on the plane of T^s and ΔT=T_s-T_a to determine possibilities for forming an intense vortex. Here, T_a is the atmospheric temperature at the height z=10 m. A PL should have ΔT>20^oC in accordance with the observations and numerical simulations. The conditions for a TC are not so straightforward but our diagram shows that the temperature difference of a few degrees, or possibly even a fraction of a degree, might be sufficient for TC development for a range of static stabilities and development times.     

Tropical Cyclones and Polar Lows: Velocity, Size,and Energy Scales, and Relation to the 26°C Cyclone Origin Criteria

The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones (TC) and polar lows (PL) by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms: (1) as expressed through the buoyancy flux b and the Coriolis parameter lc for rotating fluids convection, and (2) as expressed with the cube of velocity times the drag coefficient through the formula for total kinetic energy dissipation in the atmospheric boundary layer. In the quasistationary case the dissipation equals the generation of the energy. In both cases the velocity scale can be expressed through temperature and humidity differences between the ocean and the atmosphere in terms of the reduced gravity, and both forms produce quite comparable velocity scales. Using parameters b and lc., we can form scales of the area and, by adding the mass of a unit air column, a scale of the total kinetic energy as well. These scales nicely explain the much smaller size of a PL, as compared to a TC, and the total kinetic energy of a TC is of the order 1018 - 1019 J. It will be shown that wind of 33 m s-1 is produced when the total enthalpy fluxes between the ocean and the atmosphere are about 700 W m-2 for a TC and 1700 W m-2 for a PL, in association with the much larger role of the latent heat in the first case and the stricter geostrophic constraints and larger static stability in the second case. This replaces the mystical role of 26~C as a criterion for TC origin. The buoyancy flux, a product of the reduced gravity and the wind speed, together with the atmospheric static stability, determines the rate of the penetrating convection. It is known from the observations that the formation time for a PL reaching an altitude of 5-6 km can be only a few hours, and a day, or even half a day, for a TC reaching 15-18 km. These two facts allow us to construct curves on the plane of Ts and △T= Ts - Ta to determine possibilities for forming an intense vortex. Here, Ta is the atmospheric temperature at the height z = 10 m. A PL should have △T > 20℃ in accordance with the observations and numerical simulations. The conditions for a TC are not so straightforward but our diagram shows that the temperature difference of a few degrees, or possibly even a fraction of a degree, might be sufficient for TC development for a range of static stabilities and development times.     

Comparison of Products from ERA-40, NCEP-2, and CRU with Station Data for Summer Precipitation over China

Summer precipitation products from the 45-Year European Centre for Medium-Range Weather Forecast (ECMWF) Reanalysis (ERA-40), and NCEP-Department of Energy (DOE) Atmospheric Model Intercomparison Project (AMIP-II) Reanalysis (NCEP-2), and Climatic Research Unit (CRU) TS 2.1 dataset are compared with the corresponding observations over China in order to understand the quality and utility of the reanalysis datasets for the period 1979–2001. The results reveal that although the magnitude and location of the rainfall belts differ among the reanalysis, CRU, and station data over South and West China, the spatial distributions show good agreement over most areas of China. In comparison with the observations in most areas of China, CRU best matches the observed summer precipitation, while ERA-40 reports less precipitation and NCEP-2 reports more precipitation than the observations. With regard to the amplitude of the interannual variations, CRU is better than either of the reanalyses in representing the corresponding observations. The amplitude in NCEP-2 is stronger but that of ERA-40 is weaker than the observations in most study domains. NCEP-2 has a more obvious interannual variability than ERA-40 or CRU in most areas of East China. Through an Empirical orthogonal function (EOF) analysis, the main features of the rainfall belts produced by CRU agree better with the observations than with those produced by the reanalyses in the Yangtze-Huaihe River valley. In East of China, particularly in the Yangtze-Huaihe River valley, CRU can reveal the quasi-biennial oscillation of summer precipitation represented by the observations, but the signal of ERA-40 is comparatively weak and not very obvious, whereas that of NCEP-2 is also weak before 1990 but very strong after 1990. The results also suggest that the magnitude of the precipitation difference between ERA-40 and the observations is smaller than that between NCEP-2 and the observations, but the variations represented by NCEP-2 are more reasonable than those given by ERA-40 in most areas of East China to some extent.     

Magnitude,Scale, and Dynamics of the 2020 Mei-yu Rains and Floods over China※

Large parts of East and South Asia were affected by heavy precipitation and flooding during early summer 2020. This study provides both a statistical and dynamical characterization of rains and floods affecting the Yangtze River Basin (YRB). By aggregating daily and monthly precipitation over river basins across Asia, it is shown that the YRB is one of the areas that was particularly affected. June and July 2020 rainfall was higher than in the previous 20 years, and the YRB experienced anomalously high rainfall across most of its sub-basins. YRB discharge also attained levels not seen since 1998/1999. An automated method detecting the daily position of the East Asian Summer Monsoon Front (EASMF) is applied to show that the anomalously high YRB precipitation was associated with a halted northward progression of the EASMF and prolonged mei-yu conditions over the YRB lasting more than one month. Two 5-day heavy-precipitation episodes (12?16 June and 4?8 July 2020) are selected from this period for dynamical characterization, including Lagrangian trajectory analysis. Particular attention is devoted to the dynamics of the airstreams converging at the EASMF. Both episodes display heavy precipitation and convergence of monsoonal and subtropical air masses. However, clear differences are identified in the upper-level flow pattern, substantially affecting the balance of airmass advection towards the EASMF. This study contextualizes heavy precipitation in Asia in summer 2020 and showcases several analysis tools developed by the authors for the study of such events.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Drought in Northeast Brazil—past,present, and future

This study provides an overview of the drought situation in Northeast Brazil for the past, present, and future. Droughts affect more people than any other natural hazard owing to their large scale and long-lasting nature. They are recurrent in the region and while some measures have been taken by the governments to mitigate their impacts, there is still a perception that residents, mainly in rural areas, are not yet adapted to these hazards. The drought affecting the Northeast from 2012 to 2015, however, has had an intensity and impact not seen in several decades and has already destroyed large swaths of cropland, affecting hundreds of cities and towns across the region, and leaving ranchers struggling to feed and water cattle. Future climate projections for the area show large temperature increases and rainfall reductions, which, together with a tendency for longer periods with consecutive dry days, suggest the occurrence of more frequent/intense dry spells and droughts and a tendency toward aridification in the region. All these conditions lead to an increase in evaporation from reservoirs and lakes, affecting irrigation and agriculture as well as key water uses including hydropower and industry, and thus, the welfare of the residents. Integrating drought monitoring and seasonal forecasting provides efficient means of assessing impacts of climate variability and change, identifying vulnerabilities, and allowing for better adaptation measures not only for medium- and long-term climate change but also for extremes of the interannual climate variability, particularly droughts.

     

Air,Water and Soil Environmental Pollution and Regulation Mechanisms in Beijing and Its Ambient Area

1 Physical and Chemical Structural Characteristics of Urban ABLObtained preliminarily the atmospheric dynamic structure of the urban ABL: Urban heat island and urban wind, the zero-plane displacement and roughness of urban building blocks, the distribution and dynamic and thermal impacts of urban fog. In the aspect of city-scale simulation research, the effect of the topography around the outskirts of Beijing on the urban wind field and the effect of buildings on the urban atmospheric circulation are studied. The possible influencing mechanisms of the U-shaped mountains to the west and north of Beijing on the daily variation of the urban wind field, and the formation, transport and cumulation of urban pollutants are analyzed and some distribution charts are obtained.     

Fourier transform infrared studies of the reaction of Cl atoms with PAN,PPN, CH3OOH,HCOOH, CH3COCH3 and CH3COC2H5 at 295±2 K

The relative rate technique has been used to measure rate constants for the reaction of chlorine atoms with peroxyacetylnitrate (PAN), peroxypropionylnitrate (PPN), methylhydroperoxide, formic acid, acetone and butanone. Decay rates of these organic species were measured relative to one or more of the following reference compounds; ethene, ethane, chloroethane, chloromethane, and methane. Using rate constants of 9.29×10^–11, 5.7×10^–11, 8.04×10^–12, 4.9×10^–13, and 1.0×10^–13 cm³ molecule^–1 sec^–1 for the reaction of Cl atoms with ethene, ethane, chloroethane, chloromethane, and methane respectively, the following rate constants were derived, in units of cm³ molecule^–1 s^–1: PAN, <7×10^–15; PPN, (1.14±0.12)×10^–12; HCOOH, (2.00±0.25)×10^–13; CH₃OOH, (5.70±0.23)×10^–11; CH₃COCH₃, (2.37±0.12)×10^–12; and CH₃COC₂H₅, (4.13±0.57)×10^–11. Quoted errors represent 2 and do not include possible systematic errors due to errors in the reference rate constants. Experiments were performed at 295±2 K and 700 torr total pressure of nitrogen or synthetic air. The results are discussed with respect to the previous literature data and to the modelling of nonmethane hydrocarbon oxidation in the atmosphere.In recent discussions with Dr. R. A. Cox of Harwell Laboratory, UKAEA, we learnt of a preliminary value for the rate constant of the reaction of Cl with acetone of (2.5±1.0)×10^–12 cm³ molecule^–1 sec^–1 measured by R. A. Cox, M. E. Jenkin, and G. D. Hayman using molecular modulation techniques. This value is in good agreement with our results.     

Tropical cyclones,climate change,and scientific uncertainty: what do we know,what does it mean,and what should be done?

The question of whether and to what extent global warming may be changing tropical cyclone (TC) activity is of great interest to decision makers. The presence of a possible climate change signal in TC activity is difficult to detect because interannual variability necessitates analysis over longer time periods than available data allow. Projections of future TC activity are hindered by computational limitations and uncertainties about changes in regional climate, large scale patterns, and TC response. This review discusses the state of the field in terms of theory, modeling studies and data. While Atlantic TCs have recently become more intense, evidence for changes in other basins is not persuasive, and changes in the Atlantic cannot be clearly attributed to either natural variability or climate change. However, whatever the actual role of climatic change, these concerns have opened a “policy window” that, if used appropriately, could lead to improved protection against TCs.     

Existing infrastructure and the 2°C target

To clarify the link between existing infrastructure legacy and the 2°C target, we extend the work of Davis et al. (Science 329:1330–1333, 2010) by introducing non-CO₂ greenhouse gases and the inertia in transportation-needs drivers. We conclude that climate policies able to maintain climate change below 2°C cannot disregard existing infrastructure.     

Infra‐red FTS measurements of CH4, N2O,O3, HNO3, HCl,CFC‐11 and CFC‐12 from the MANTRA balloon campaign

Abstract

The Middle Atmosphere Nitrogen TRend Assessment (MANTRA) campaign is intended to address the question of whether possible changes in the mid‐latitude nitrogen budget can account for discrepancies between predicted and observed ozone loss at these latitudes. In addition, MANTRA seeks to establish consistency between old and new measurement techniques. We report here the results of infra‐red spectral measurements made with a high resolution Fourier Transform Spectrometer (FTS) during the August 1998 MANTRA flight. Vertical mixing ratio profiles of methane, nitrous oxide, ozone, nitric acid, hydrogen chloride, CFC‐11 and CFC‐12 are presented. The error estimates on the retrieved values are larger than typical for this type of measurement, due mostly to instrument problems that occurred during flight. In addition, the CFC values are larger than expected. Possible reasons for this are discussed.     

Statistical Analyses of Runoff,River Stage and Peak Discharge,Tao River,China

Empirical equations relating runoff,maximum river stageand peak flow to rainfall were developed for the Tao riverdrainage basin at Li Jia Cun.A basin recession curve wasderived and a frequency analysis of extreme flood events wasconducted.Results showed that storm runoff was highlycorrelated with rainfall.However,maximum river stage andpeak flow were not as well correlated with rainfall.It wasfound necessary to separate rainfall into stratiform andconvective types.The magnitudes of the once-in-50 year,once-in-100 year and once-in-200 year events,calculated by thelog Pearson type Ⅲ distribution were closely equal to thoseestimated from a log-probability graph.Recommendations forimproving the predictive ability of the empirical equations arepresented.