Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

Stratospheric temperature and radiative forcing response to 11-year solar cycle changes in irradiance and ozone

The 11-yr solar cycle temperature response to spectrally resolved solar irradiance changes and associated ozone changes is calculated using a fixed dynamical heating (FDH) model. Imposed ozone changes are from satellite observations, in contrast to some earlier studies. A maximum of 1.6 K is found in the equatorial upper stratosphere and a secondary maximum of 0.4 K in the equatorial lower stratosphere, forming a double peak in the vertical. The upper maximum is primarily due to the irradiance changes while the lower maximum is due to the imposed ozone changes. The results compare well with analyses using the 40-yr ECMWF Re-Analysis (ERA-40) and NCEP/NCAR datasets. The equatorial lower stratospheric structure is reproduced even though, by definition, the FDH calculations exclude dynamically driven temperature changes, suggesting an important role for an indirect dynamical effect through ozone redistribution. The results also suggest that differences between the Stratospheric Sounding Unit (SSU)/Microwave Sounding Unit (MSU) and ERA-40 estimates of the solar cycle signal can be explained by the poor vertical resolution of the SSU/MSU measurements. The adjusted radiative forcing of climate change is also investigated. The forcing due to irradiance changes was 0.14 W m−2, which is only 78% of the value obtained by employing the standard method of simple scaling of the total solar irradiance (TSI) change. The difference arises because much of the change in TSI is at wavelengths where ozone absorbs strongly. The forcing due to the ozone change was only 0.004 W m−2 owing to strong compensation between negative shortwave and positive longwave forcings.

Convectively coupled equatorial waves in high resolution Hadley centre climate models

Current global atmospheric models fail to simulate well organised tropical phenomena in which convection interacts with dynamics and physics. A new methodology to identify convectively coupled equatorial waves, developed by NCAS-Climate, has been applied to output from the two latest models of the Met Office/Hadley Centre which have fundamental differences in dynamical formulation. Variability, horizontal and vertical structures, and propagation characteristics of tropical convection and equatorial waves, along with their coupled behaviour in the models are examined and evaluated against a previous comprehensive study of observations. It is shown that, in general, the models perform well for equatorial waves coupled with off-equatorial convection. However they perform poorly for waves coupled with equatorial convection. The vertical structure of the simulated wave is not conducive to energy conversion/growth and does not support the correct physical-dynamical coupling that occurs in the real world. The following figure shows an example of the Kelvin wave coupled with equatorial convection. It shows that the models fail to simulate a key feature of convectively coupled Kelvin wave in observations, namely near surface anomalous equatorial zonal winds together with intensified equatorial convection and westerly winds in phase with the convection. The models are also not able to capture the observed vertical tilt structure and the vertical propagation of the Kelvin wave into the lower stratosphere as well as the secondary peak in the mid-troposphere, particularly in HadAM3. These results can be used to provide a test-bed for experimentation to improve the coupling of physics and dynamics in climate and weather models.

Annular variability and eddy-zonal flow interactions in a simplified atmospheric GCM. Part 1: Characterization of high and low frequency behaviour

Experiments have been performed using a simplified, Newtonian forced, global circulation model to investigate how variability of the tropospheric jet can be characterized by examining the combined fluctuations of the two leading modes of annular variability. Eddy forcing of this variability is analyzed in the phase space of the leading modes using the vertically integrated momentum budget. The nature of the annular variability and eddy forcing depends on the time scale. At low frequencies the zonal flow and baroclinic eddies are in quasi equilibrium and anomalies propagate poleward. The eddies are shown primarily to reinforce the anomalous state and are closely balanced by the linear damping, leaving slow evolution as a residual. At high frequencies the flow is strongly evolving and anomalies are initiated on the poleward side of the tropospheric jet and propagate equatorward. The eddies are shown to drive this evolution strongly: eddy location and amplitude reflect the past baroclinicity, while eddy feedback on the zonal flow may be interpreted in terms of wave breaking associated with baroclinic life cycles in lateral shear.

Climate impacts of recent multidecadal changes in Atlantic Ocean Sea Surface Temperature: A multimodel comparison

During the twentieth century sea surface temperatures in the Atlantic Ocean exhibited prominent multidecadal variations. The source of such variations has yet to be rigorously established—but the question of their impact on climate can be investigated. Here we report on a set of multimodel experiments to examine the impact of patterns of warming in the North Atlantic, and cooling in the South Atlantic, derived from observations, that is characteristic of the positive phase of the Atlantic Multidecadal Oscillation (AMO). The experiments were carried out with six atmospheric General Circulation Models (including two versions of one model), and a major goal was to assess the extent to which key climate impacts are consistent between the different models. The major climate impacts are found over North and South America, with the strongest impacts over land found over the United States and northern parts of South America. These responses appear to be driven by a combination of an off-equatorial Gill response to diabatic heating over the Caribbean due to increased rainfall within the region and a Northward shift in the Inter Tropical Convergence Zone (ITCZ) due to the anomalous cross-equatorial SST gradient. The majority of the models show warmer US land temperatures and reduced Mean Sea Level Pressure during summer (JJA) in response to a warmer North Atlantic and a cooler South Atlantic, in line with observations. However the majority of models show no significant impact on US rainfall during summer. Over northern South America, all models show reduced rainfall in southern hemisphere winter (JJA), whilst in Summer (DJF) there is a generally an increase in rainfall. However, there is a large spread amongst the models in the magnitude of the rainfall anomalies over land. Away from the Americas, there are no consistent significant modelled responses. In particular there are no significant changes in the North Atlantic Oscillation (NAO) over the North Atlantic and Europe in Winter (DJF). Additionally, the observed Sahel drying signal in African rainfall is not seen in the modelled responses. Suggesting that, in contrast to some studies, the Atlantic Multidecadal Oscillation was not the primary driver of recent reductions in Sahel rainfall.

Atmosphere feedbacks during ENSO in a coupled GCM with a modified atmospheric convection scheme

The too diverse representation of ENSO in a coupled GCM limits one’s ability to describe future change of its properties. Several studies pointed to the key role of atmosphere feedbacks in contributing to this diversity. These feedbacks are analyzed here in two simulations of a coupled GCM that differ only by the parameterization of deep atmospheric convection and the associated clouds. Using the Kerry–Emanuel (KE) scheme in the L’Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL CM4; KE simulation), ENSO has about the right amplitude, whereas it is almost suppressed when using the Tiedke (TI) scheme. Quantifying both the dynamical Bjerknes feedback and the heat flux feedback in KE, TI, and the corresponding Atmospheric Model Intercomparison Project (AMIP) atmosphere-only simulations, it is shown that the suppression of ENSO in TI is due to a doubling of the damping via heat flux feedback. Because the Bjerknes positive feedback is weak in both simulations, the KE simulation exhibits the right ENSO amplitude owing to an error compensation between a too weak heat flux feedback and a too weak Bjerknes feedback. In TI, the heat flux feedback strength is closer to estimates from observations and reanalysis, leading to ENSO suppression. The shortwave heat flux and, to a lesser extent, the latent heat flux feedbacks are the dominant contributors to the change between TI and KE. The shortwave heat flux feedback differences are traced back to a modified distribution of the large-scale regimes of deep convection (negative feedback) and subsidence (positive feedback) in the east Pacific. These are further associated with the model systematic errors. It is argued that a systematic and detailed evaluation of atmosphere feedbacks during ENSO is a necessary step to fully understand its simulation in coupled GCMs.

The role of atmosphere feedbacks during ENSO in the CMIP3 models

Several studies using ocean-atmosphere GCMs suggest that the atmospheric component plays a dominant role in the modelled ENSO. To help elucidate these findings, the two main atmosphere feedbacks relevant to ENSO, the Bjerknes positive feedback (µ) and the heat flux negative feedback (), are analysed here in 12 coupled GCMs. We find that the models generally underestimate both feedbacks, leading to an error compensation. The strength of is inversely related to the ENSO amplitude in the models and the latent heat and shortwave flux components of this feedback dominate. Furthermore, the shortwave component could help explain the model diversity in both overall and ENSO amplitude.

Understanding land-sea warming contrast in response to increasing greenhouse gases. Part I: Transient adjustment

Climate model simulations consistently show that surface temperature over land increases more rapidly than over sea in response to greenhouse gas forcing. The enhanced warming over land is not simply a transient effect caused by the land–sea contrast in heat capacities, since it is also present in equilibrium conditions. This paper elucidates the transient adjustment processes over time scales of days to weeks of the surface and tropospheric climate in response to a doubling of CO2 and to changes in sea surface temperature (SST), imposed separately and together, using ensembles of experiments with an atmospheric general circulation model. These adjustment processes can be grouped into three stages: immediate response of the troposphere and surface processes (day 1), fast adjustment of surface processes (days 2–5), and adjustment of the whole troposphere (days 6–20). Some land surface warming in response to doubled CO2 (with unchanged SSTs) occurs immediately because of increased downward longwave radiation. Increased CO2 also leads to reduced plant stomatal resistance and hence restricted evaporation, which increases land surface warming in the first day. Rapid reductions in cloud amount lead in the next few days to increased downward shortwave radiation and further warming, which spreads upward from the surface, and by day 5 the surface and tropospheric response is statistically consistent with the equilibrium value. Land surface warming in response to imposed SST change (with unchanged CO2) is slower. Tropospheric warming is advected inland from the sea, and over land it occurs at all levels together rather than spreading upward from the surface. The atmospheric response to prescribed SST change in about 20 days is statistically consistent with the equilibrium value, and the warming is largest in the upper troposphere over both land and sea. The land surface warming involves reduction of cloud cover and increased downward shortwave radiation, as in the experiment with CO2 change, but in this case it is due to the restriction of moisture supply to the land (indicated by reduced soil moisture), whereas in the CO2 forcing experiment it is due to restricted evaporation despite increased moisture supply (indicated by increased soil moisture). The warming over land in response to SST change is greater than over the sea and is the dominant contribution to the land–sea warming contrast under enhanced CO2 forcing.

Life cycle of the QBO-modulated 11-year solar cycle signals in the Northern Hemispheric winter

This paper provides some insights on the quasi-biennial oscillation (QBO) modulated 11-year solar cycle (11-yr SC) signals in Northern Hemisphere (NH) winter temperature and zonal wind. Daily ERA-40 Reanalysis and ECMWF Operational data for the period of 1958-2006 were used to examine the seasonal evolution of the QBO-solar cycle relationship at various pressure levels up to the stratopause. The results show that the solar signals in the NH winter extratropics are indeed QBO-phase dependent, moving poleward and downward as winter progresses with a faster descent rate under westerly QBO than under easterly QBO. In the stratosphere, the signals are highly significant in late January to early March and have a life span of 30-50 days. Under westerly QBO, the stratospheric solar signals clearly lead and connect to those in the troposphere in late March and early April where they have a life span of 10 days. As the structure changes considerably from the upper stratosphere to the lower troposphere, the exact month when the maximum solar signals occur depends largely on the altitude chosen. For the low-latitude stratosphere, our analysis supports a vertical double-peaked structure of positive signature of the 11-yr SC in temperature, and demonstrates that this structure is further modulated by the QBO. These solar signals have a longer life span (3-4 months) in comparison to those in the extratropics. The solar signals in the lower stratosphere are stronger in early winter but weaker in late winter, while the reverse holds in the upper stratosphere.

High Resolution Dynamics Limb Sounder measurements of gravity wave activity in the 2006 Arctic stratosphere

Observations from the High Resolution Dynamics Limb Sounder (HIRDLS) instrument on NASA's Aura satellite are used to quantify gravity wave momentum fluxes in the middle atmosphere. The period around the 2006 Arctic sudden stratospheric warming (SSW) is investigated, during which a substantial elevation of the stratopause occurred. Analysis of the HIRDLS results, together with analysis of European Centre for Medium-Range Weather Forecasting zonal winds, provide direct evidence of wind filtering of the gravity wave spectrum during this period. This confirms previous hypotheses from model studies and further contributes to our understanding of the effects of gravity wave driving on the winter polar stratopause.

Carbon dioxide induced stomatal closure increases radiative forcing of climate via a rapid reduction in low cloud

We performed an ensemble of twelve five-year experiments using a coupled climate-carbon-cycle model with scenarios of prescribed atmospheric carbon dioxide concentration; CO2 was instantaneously doubled or quadrupled at the start of the experiments. Within these five years, climate feedback is not significantly influenced by the effects of climate change on the carbon system. However, rapid changes take place, within much less than a year, due to the physiological effect of CO2 on plant stomatal conductance, leading to adjustment in the shortwave cloud radiative effect over land, due to a reduction in low cloud cover. This causes a 10% enhancement to the radiative forcing due to CO2, which leads to an increase in the equilibrium warming of 0.4 and 0.7 K for doubling and quadrupling. The implications for calibration of energy-balance models are discussed.

Climate predictability in the second year

In this paper, the predictability of climate arising from ocean heat content (OHC) anomalies is investigated in the HadCM3 coupled atmosphere–ocean model. An ensemble of simulations of the twentieth century are used to provide initial conditions for a case study. The case study consists of two ensembles started from initial conditions with large differences in regional OHC in the North Atlantic, the Southern Ocean and parts of the West Pacific. Surface temperatures and precipitation are on average not predictable beyond seasonal time scales, but for certain initial conditions there may be longer predictability. It is shown that, for the case study examined here, some aspects of tropical precipitation, European surface temperatures and North Atlantic sea-level pressure are potentially predictable 2 years ahead. Predictability also exists in the other case studies, but the climate variables and regions, which are potentially predictable, differ. This work was done as part of the Grid for Coupled Ensemble Prediction (GCEP) eScience project.

Using idealized snow forcing to test teleconnections with the Indian summer monsoon in the Hadley Centre GCM

Anomalous heavy snow during winter or spring has long been regarded as a possible precursor of deficient Indian monsoon rainfall during the subsequent summer. However previous work in this field is inconclusive, in terms of the mechanism that communicates snow anomalies to the monsoon summer, and even the region from which snow has the most impact. In this study we explore these issues in coupled and atmosphere-only versions of the Hadley Centre model. A 1050-year control integration of the HadCM3 coupled model, which well represents the seasonal cycle of snow cover over the Eurasian continent, is analysed and shows evidence for weakened monsoons being preceded by strong snow forcing (in the absence of ENSO) over either the Himalaya/Tibetan Plateau or north/west Eurasia regions. However, empirical orthogonal function (EOF) analysis of springtime interannual variability in snow depth shows the leading mode to have opposite signs between these two regions, suggesting that competing mechanisms may be possible. To determine the dominant region, ensemble integrations are carried out using HadAM3, the atmospheric component of HadCM3, and a variety of anomalous snow forcing initial conditions obtained from the control integration of the coupled model. Forcings are applied during spring in separate experiments over the Himalaya/Tibetan Plateau and north/west Eurasia regions, in conjunction with climatological SSTs in order to avoid the direct effects of ENSO. With the aid of idealized forcing conditions in sensitivity tests, we demonstrate that forcing from the Himalaya region is dominant in this model via a Blanford-type mechanism involving reduced surface sensible heat and longwave fluxes, reduced heating of the troposphere over the Tibetan Plateau and consequently a reduced meridional tropospheric temperature gradient which weakens the monsoon during early summer. Snow albedo is shown to be key to the mechanism, explaining around 50% of the perturbation in sensible heating over the Tibetan Plateau, and accounting for the majority of cooling through the troposphere.

Tracer correlations in the northern high latitude lowermost stratosphere: Influence of cross‐tropopause mass exchange

We present an analysis of trace gas correlations in the lowermost stratosphere. In‐situ aircraft measurements of CO, N2O, NOy and O3, obtained during the STREAM 1997 winter campaign, have been used to investigate the role of cross‐tropopause mass exchange on tracer‐tracer relations. At altitudes several kilometers above the local tropopause, undisturbed stratospheric air was found with NOy/NOy * ratios close to unity, NOy/O3 about 0.003–0.006 and CO mixing ratios as low as 20 ppbv (NOy * is a proxy for total reactive nitrogen derived from NOy–N2O relations measured in the stratosphere). Mixing of tropospheric air into the lowermost stratosphere has been identified by enhanced ratios of NOy/NOy * and NOy/O3, and from scatter plots of CO versus O3. The enhanced NOy/O3 ratio in the lowermost stratospheric mixing zone points to a reduced efficiency of O3 formation from aircraft NOx emissions.