On the tropospheric response to anomalous stratospheric wave drag and radiative heating

Observational and numerical evidence suggest that variability in the extratropical stratospheric circulation has a demonstrable impact on tropospheric variability on intraseasonal time scales. In this study, it is demonstrated that the amplitude of the observed tropospheric response to vacillations in the stratospheric flow is quantitatively similar to the zonal-mean balanced response to the anomalous wave forcing at stratospheric levels. It is further demonstrated that the persistence of the tropospheric response is consistent with the impact of anomalous diabatic heating in the polar stratosphere as stratospheric temperatures relax to climatology. The results contradict previous studies that suggest that variations in stratospheric wave drag are too weak to account for the attendant changes in the tropospheric flow. However, the results also reveal that stratospheric processes alone cannot account for the observed meridional redistribution of momentum within the troposphere.

Radiative-dynamical climatology of the first-generation Canadian middle atmosphere model

The Canadian Middle Atmosphere Modelling (MAM) project is a collaboration between thé Atmospheric Environment Service (AES) of Environment Canada and several Canadian universities. Its goal is thé development of a comprehensive General Circulation Model of the troposphere-stratosphere-mesosphere System, starting from the AES/CCCma third-generation atmospheric General Circulation Model. This paper describes the basic features of the first-generation Canadian MAM and some aspects of its radiative-dynamical climatology. Standard first-order mean diagnostics are presented for monthly means and for the annual cycle of zonal-mean winds and temperatures. The mean meridional circulation is examined, and comparison is made between thé steady diabatic, downward controlled, and residual stream functions. It is found that downward control holds quite well in the monthly mean through most of the middle atmosphere, even during equinoctal periods. The relative roles of different drag processes in determining the mean downwelling over the wintertime polar middle stratosphere is examined, and the vertical structure of the drag is quantified.

Improved evaporative flux partitioning and carbon flux in the land surface model JULES: impact on the simulation of land surface processes in temperate Europe

The primary role of land surface models embedded in climate models is to partition surface available energy into upwards, radiative, sensible and latent heat fluxes. Partitioning of evapotranspiration, ET, is of fundamental importance: as a major component of the total surface latent heat flux, ET affects the simulated surface water balance, and related energy balance, and consequently the feedbacks with the atmosphere. In this context it is also crucial to credibly represent the CO2 exchange between ecosystems and their environment. In this study, JULES, the land surface model used in UK weather and climate models, has been evaluated for temperate Europe. Compared to eddy covariance flux measurements, the CO2 uptake by the ecosystem is underestimated and the ET overestimated. In addition, the contribution to ET from soil and intercepted water evaporation far outweighs the contribution of plant transpiration. To alleviate these biases, adaptations have been implemented in JULES, based on key literature references. These adaptations have improved the simulation of the spatio-temporal variability of the fluxes and the accuracy of the simulated GPP and ET, including its partitioning. This resulted in a shift of the seasonal soil moisture cycle. These adaptations are expected to increase the fidelity of climate simulations over Europe. Finally, the extreme summer of 2003 was used as evaluation benchmark for the use of the model in climate change studies. The improved model captures the impact of the 2003 drought on the carbon assimilation and the water use efficiency of the plants. It, however, underestimates the 2003 GPP anomalies. The simulations showed that a reduction of evaporation from the interception and soil reservoirs, albeit not of transpiration, largely explained the good correlation between the carbon and the water fluxes anomalies that was observed during 2003. This demonstrates the importance of being able to discriminate the response of individual component of the ET flux to environmental forcing.

A multimodel assessment of the influence of regional anthropogenic emission reductions on aerosol direct radiative forcing and the role of intercontinental transport

In this study, we assess changes of aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four major pollution regions in the Northern Hemisphere by using results from nine global models in the framework of the Hemispheric Transport of Air Pollution (HTAP). DRF at top of atmosphere (TOA) and surface is estimated based on AOD results from the HTAP models and AOD-normalized DRF (NDRF) from a chemical transport model. The multimodel results show that, on average, a 20% reduction of anthropogenic emissions in North America, Europe, East Asia, and South Asia lowers the global mean AOD (all-sky TOA DRF) by 9.2% (9.0%), 3.5% (3.0%), and 9.4% (10.0%) for sulfate, particulate organic matter (POM), and black carbon (BC), respectively. Global annual average TOA all-sky forcing efficiency relative to particle or gaseous precursor emissions from the four regions (expressed as multimodel mean ± one standard deviation) is  ±3.5 ±0.8,  ±4.0 ±1.7, and 29.5 ±18.1mWm ±2 per Tg for sulfate (relative to SO2), POM, and BC, respectively. The impacts of the regional emission reductions on AOD and DRF extend well beyond the source regions because of intercontinental transport (ICT). On an annual basis, ICT accounts for 11 ±5% to 31 ±9% of AOD and DRF in a receptor region at continental or subcontinental scale, with domestic emissions accounting for the remainder, depending on regions and species. For sulfate AOD, the largest ICT contribution of 31 ±9% occurs in South Asia, which is dominated by the emissions from Europe. For BC AOD, the largest ICT contribution of 28 ±18% occurs in North America, which is dominated by the emissions from East Asia. The large spreads among models highlight the need to improve aerosol processes in models, and evaluate and constrain models with observations.

Quantifying the importance of galactic cosmic rays in cloud microphysical processes

Galactic Cosmic Rays are one of the major sources of ion production in the troposphere and stratosphere. Recent studies have shown that ions form electrically charged clusters which may grow to become cloud droplets. Aerosol particles charge by the attachment of ions and electrons. The collision efficiency between a particle and a water droplet increases, if the particle is electrically charged, and thus aerosol-cloud interactions can be enhanced. Because these microphysical processes may change radiative properties of cloud and impact Earth's climate it is important to evaluate these processes' quantitative effects. Five different models developed independently have been coupled to investigate this. The first model estimates cloud height from dew point temperature and the temperature profile. The second model simulates the cloud droplet growth from aerosol particles using the cloud parcel concept. In the third model, the scavenging rate of the aerosol particles is calculated using the collision efficiency between charged particles and droplets. The fourth model calculates electric field and charge distribution on water droplets and aerosols within cloud. The fifth model simulates the global electric circuit (GEC), which computes the conductivity and ionic concentration in the atmosphere in altitude range 0–45 km. The first four models are initially coupled to calculate the height of cloud, boundary condition of cloud, followed by growth of droplets, charge distribution calculation on aerosols and cloud droplets and finally scavenging. These models are incorporated with the GEC model. The simulations are verified with experimental data of charged aerosol for various altitudes. Our calculations showed an effect of aerosol charging on the CCN concentration within the cloud, due to charging of aerosols increase the scavenging of particles in the size range 0.1 µm to 1 µm.

An assessment of the impact of local processes on dust lifting in martian climate models

Simulation of the lifting of dust from the planetary surface is of substantially greater importance on Mars than on Earth, due to the fundamental role that atmospheric dust plays in the former’s climate, yet the dust emission parameterisations used to date in martian global climate models (MGCMs) lag, understandably, behind their terrestrial counterparts in terms of sophistication. Recent developments in estimating surface roughness length over all martian terrains and in modelling atmospheric circulations at regional to local scales (less than O(100 km)) presents an opportunity to formulate an improved wind stress lifting parameterisation. We have upgraded the conventional scheme by including the spatially varying roughness length in the lifting parameterisation in a fully consistent manner (thereby correcting a possible underestimation of the true threshold level for wind stress lifting), and used a modification to account for deviations from neutral stability in the surface layer. Following these improvements, it is found that wind speeds at typical MGCM resolution never reach the lifting threshold at most gridpoints: winds fall particularly short in the southern midlatitudes, where mean roughness is large. Sub-grid scale variability, manifested in both the near-surface wind field and the surface roughness, is then considered, and is found to be a crucial means of bridging the gap between model winds and thresholds. Both forms of small-scale variability contribute to the formation of dust emission ‘hotspots’: areas within the model gridbox with particularly favourable conditions for lifting, namely a smooth surface combined with strong near-surface gusts. Such small-scale emission could in fact be particularly influential on Mars, due both to the intense positive radiative feedbacks that can drive storm growth and a strong hysteresis effect on saltation. By modelling this variability, dust lifting is predicted at the locations at which dust storms are frequently observed, including the flushing storm sources of Chryse and Utopia, and southern midlatitude areas from which larger storms tend to initiate, such as Hellas and Solis Planum. The seasonal cycle of emission, which includes a double-peaked structure in northern autumn and winter, also appears realistic. Significant increases to lifting rates are produced for any sensible choices of parameters controlling the sub-grid distributions used, but results are sensitive to the smallest scale of variability considered, which high-resolution modelling suggests should be O(1 km) or less. Use of such models in future will permit the use of a diagnosed (rather than prescribed) variable gustiness intensity, which should further enhance dust lifting in the southern hemisphere in particular.

Vertical structure and physical processes of the Madden-Julian oscillation: exploring key model physics in climate simulations

Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one-fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models, but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high and low rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.

Vertical structure and physical processes of the Madden-Julian oscillation: Biases and uncertainties at short range

An analysis of diabatic heating and moistening processes from 12-36 hour lead time forecasts from 12 Global Circulation Models are presented as part of the "Vertical structure and physical processes of the Madden-Julian Oscillation (MJO)" project. A lead time of 12-36 hours is chosen to constrain the large scale dynamics and thermodynamics to be close to observations while avoiding being too close to the initial spin-up for the models as they adjust to being driven from the YOTC analysis. A comparison of the vertical velocity and rainfall with the observations and YOTC analysis suggests that the phases of convection associated with the MJO are constrained in most models at this lead time although the rainfall in the suppressed phase is typically overestimated. Although the large scale dynamics is reasonably constrained, moistening and heating profiles have large inter-model spread. In particular, there are large spreads in convective heating and moistening at mid-levels during the transition to active convection. Radiative heating and cloud parameters have the largest relative spread across models at upper levels during the active phase. A detailed analysis of time step behaviour shows that some models show strong intermittency in rainfall and differences in the precipitation and dynamics relationship between models. The wealth of model outputs archived during this project is a very valuable resource for model developers beyond the study of the MJO. In addition, the findings of this study can inform the design of process model experiments, and inform the priorities for field experiments and future observing systems.

Radiative capture of nucleons at astrophysical energies with single-particle states

Radiative capture of nucleons at energies of astrophysical interest is one of the most important processes for nucleosynthesis. The nucleon capture can occur either by a compound nucleus reaction or by a direct process. The compound reaction cross sections are usually very small, especially for light nuclei. The direct capture proceeds either via the formation of a single-particle resonance or a non-resonant capture process. In this work we calculate radiative capture cross sections and astrophysical S-factors for nuclei in the mass region A < 20 using single-particle states. We carefully discuss the parameter fitting procedure adopted in the simplified two-body treatment of the capture process. Then we produce a detailed list of cases for which the model works well. Useful quantities, such as spectroscopic factors and asymptotic normalization coefficients, are obtained and compared to published data. (C) 2010 Elsevier Inc. All rights reserved.

Energy transfer processes in Yb(3+)-Tm(3+) co-doped sodium alumino-phosphate glasses with improved 1.8 mu m emission

Sodium alumino-phosphate glasses co-doped with Yb(3+) and Tm(3+) ions have been prepared with notably low OH(-) content, and characterized from the viewpoint of their spectroscopic properties. In these glasses, Yb(3+) acts as an efficient sensitizer of excitation energy at 0.98 mu m - which can be provided by high power and low cost diode lasers, and subsequently undergoes non-resonant energy transfer to Tm(3+) ions ((2)F(5/2), (3)H(6) --> (2)F(7/2), (3)H(5)). Through this process, the emitting level (3)F(4) is rapidly populated, generating improved emission at 1.8 mu m ((3)F(4) --> (3)H(6)). In order to guarantee the efficiency of such favorable energy transfer, energy losses via multiphonon decay, Yb-Yb radiative trapping, and non- radiative transfer to OH(-) groups were evaluated, and minimized when possible. The dipole - dipole energy transfer microscopic parameters corresponding to Yb(3+) --> Tm(3+), Yb(3+) --> Yb(3+) and Tm(3+) --> Tm(3+) transfers, calculated by the Forster-Dexter model, are C(Yb-Tm) = 2.9 x 10(-40) cm(6) s(-1), C(Yb-Yb) = 42 x 10(-40) cm(6) s(-1) and C(Tm-Tm) = 43 x 10(-40) cm(6) s(-1), respectively.

Challenges in Radiative Transfer Modeling - W3 IRS5

By modeling the spectral energy distribution (SED) of the W3 IRS5 high-mass star formation region and matching this model to observed data, we can constrain the physical parameters of the basic system geometry and cloud mass distribution. From these parameters, we hope to add to the understanding of high-mass star formation processes. In particular, we hope to determine if the geometries associated with lowmass star formation carry over into the high-mass regime.

Silicon monosulphide radiative association

Rate coefficients for radiative association of silicon and sulphur atoms to form silicon monosulphide (SiS) molecule are estimated. The radiative association is due mainly to approach in the E(1)Sigma(+) and A(1)Pi states of SiS. For temperatures ranging from similar to 1000 to similar to 14 000 K, the rate coefficients are found to vary from 8.43 x 10(-17) to 2.69 x 10(-16) cm(3) s(-1). Our calculated rate coefficient is higher than the values used in modelling the chemistry of Type Ia supernovae.

The formation of CN and CN+ by direct radiative association

Rate coefficients for direct radiative association of carbon and nitrogen atoms to form CN, and of carbon ions and nitrogen atoms to form CN+ ions, are calculated for temperatures in the range of 300 to 14,700 K. For the CN molecule, the rate coefficients can be represented by the standard expression, k(CN)(T) = 7.87 x 10(-19)(T/300)(0.056) exp (-96.0/T) cm(3) s(-1) for temperatures between 300 and 2700 K and k(CN)(T) = 1.37 x 10(-18)(T/300)-0.128 exp (-520.1/T) cm(-3) s(-1) at T > 2700 K. For the CN+ ion, the corresponding expression is k(CN+)(T) = 1.08 x 10(-18)(T/300)(0.071) exp (-57.5/T) cm(-3) s(-1) for the temperature range studied. Calculated rate coefficients k(CN) are about 2 orders of magnitude lower than the canonical value used in the modeling of the chemistry of various astrophysical environments.

Formation of SO, SO+, and S-2 by radiative association

Rate coefficients for radiative association of SO, SO+, and S-2 are estimated. For temperatures ranging from 300 to 14,000 K, the direct radiative association rate coefficients are found to vary with temperature from 1.73 x 10(-19) to 7.29 x 10(-19) cm(3) s(-1) and from 1.49 x 10(-21) to 3.70 x 10(-19) cm(3) s(-1) for S-2 and SO, respectively. The rate coefficients for formation through the inverse predissociation for S-2 are found to vary from 3.59 x 10(-18) to 1.44 x 10(-20) cm(3) s(-1). For SO+, the direct rate coefficient varies rapidly with temperature from 3.62 x 10(-27) cm(3) s(-1) at 2000 K to 2.34 x 10(-20) cm(3) s(-1) at 14,000 K. The direct radiative association rate coefficients increase with the increase in temperature, but the inverse predissociation rate coefficients decrease with the increase in temperature.

Radiative association of C and P, and Si and P atoms

The rate coefficients for the formation of carbon monophosphide (CP) and silicon monophosphide (SiP) by radiative association are estimated for temperatures ranging from 300 to 14 100 K. In this temperature range, the radiative association rate coefficients are found to vary from 1.14 x 10(-18) to 1.62 x 10(-18) cm(3) s(-1) and from 3.73 x 10(-20) to 7.03 x 10(-20) cm(3) s(-1) for CP and SiP, respectively. In both cases, rate coefficients increase slowly with the increase in temperature.