An annual cycle of vegetation in a GCM. Part I: implementation and impact on evaporation

The sensitivity of evaporation to a prescribed vegetation annual cycle is examined globally in the Met Office Hadley Centre Unified Model (HadAM3) which incorporates the Met Office Surface Exchange Scheme (MOSES2) as the land surface scheme. A vegetation annual cycle for each plant functional type in each grid box is derived based on satellite estimates of Leaf Area Index (LAI) obtained from the nine-year International Satellite Land Surface Climatology Project II dataset. The prescribed model vegetation seasonality consists of annual cycles of a number of structural vegetation characteristics including LAI as well as canopy height, surface roughness, canopy water capacity, and canopy heat capacity, which themselves are based on empirical relationships with LAI. An annual cycle of surface albedo, which in the model is a function of soil albedo, surface soil moisture, and LAI, is also modelled and agrees reasonably with observed estimates of the surface albedo annual cycle. Two 25-year numerical experiments are completed and compared: the first with vegetation characteristics held at annual mean values, the second with prescribed realistic seasonally varying vegetation. Initial analysis uncovered an unrealistically weak relationship between evaporation and vegetation state that is primarily due to the insensitivity of evapotranspiration to LAI. This weak relationship is strengthened by the adjustment of two MOSES2 parameters that together improve the models LAI-surface conductance relationship by comparison with observed and theoretical estimates. The extinction coefficient for photosynthetically active radiation, k _par , is adjusted downwards from 0.5 to 0.3, thereby enhancing the LAI-canopy conductance relationship. A canopy shading extinction coefficient, k _sh , that controls what fraction of the soil surface beneath a canopy is directly exposed to the overlying atmosphere is increased from 0.5 to 1.0, which effectively reduces soil evaporation under a dense canopy. When the experiments are repeated with the adjusted parameters, the relationship between evaporation and vegetation state is strengthened and is more spatially consistent. At nearly all locations, the annual cycle of evaporation is enhanced in the seasonally varying vegetation experiment. Evaporation is stronger during the peak of the growing season and, in the tropics, reduced transpiration during the dry season when LAI is small leads to significantly lower total evaporation.     

Regional climate simulation with a high resolution GCM: surface hydrology

Aspects of the surface hydrology of high resolution (T106) versions of the ECHAM3 and ECHAM4 general circulation models are analysed over the European region and compared with available observations. The focus is on evaporation, and surface measurements are shown to be useful for the identification of systematic deficiencies in the regional-scale performance of climate models on an annual and seasonal basis, such as the excessive summer dryness over continents. The annual mean evaporation at the available European observation sites is overestimated by 4 mm/month by the ECHAM3 T106, quantitatively consistent with an overestimated surface net radiation of 4 Wm^–2 over Europe. In winter, ECHAM3 shows an overestimated evaporation which compensates for an overestimated downward sensible heat flux. This is primarily related to a too strong zonalisation of the large-scale flow and associated overestimated warm air advection and windspeed. Inaccurate local land surface parameters (e.g. leaf area index, roughness length) are minor contributors to the overestimation. In early summer, the excessive solar radiation at the surface calculated with the ECHAM3 radiation scheme generates a too large evaporation and an excessive depletion of the soil moisture reservoirs. This favours the subsequent excessive summer dryness over Europe with too low values of evaporation, convective precipitation and soil moisture content, leading to a too high surface temperature. In the ECHAM4 T106 simulation, the problem of the European summer dryness is largely reduced, and the simulated evaporation as well as convective precipitation, cloud amount and soil moisture content during summer are substantially improved. The new ECHAM4 radiation scheme appears to be an important factor for this improvement, since it calculates smaller insolation values in better agreement with observations and subsequently may avoid an excessive drying of the soil. Received: 20 September 1995 / Accepted: 10 May 1996     

Empirical Estimates of Global Climate Sensitivity： An Assessment of Strategies Using a Coupled GCM

A control integration with the normal solar constant and one with it increased by 2.5% in the National Center for Atmospheric Research （NCAR） coupled atmosphere-ocean Climate System Model were conducted to see how well the actual realized global warming could be predicted just by analysis of the control results. This is a test, within a model context, of proposals that have been advanced to use knowledge of the present day climate to make ＂empirical＂ estimates of global climate sensitivity. The scaling of the top-of-the-atmosphere infrared flux and the planetary albedo as functions of surface temperature was inferred by examining four different temporal and geographical variations of the control simulations. Each of these inferences greatly overestimates the climate sensitivity of the model, largely because of the behavior of the cloud albedo. In each inference the control results suggest that cloudiness and albedo decrease with increasing surface temperature. However, the experiment with the increased solar constant actually has higher albedo and more cloudiness at most latitudes. The increased albedo is a strong negative feedback, and this helps account for the rather weak sensitivity of the climate in the NCAR model. To the extent that these model results apply to the real world, they suggest empirical evaluation of the scaling of global-mean radiative properties with surface temperature in the present day climate provides little useful guidance for estimates of the actual climate sensitivity to global changes.     

Sensitivity of the hydrological cycle to the parametrization of soil hydrology in a GCM

 The sensitivity of the hydrological cycle to soil hydrology is investigated with the LMD GCM. The reference simulation includes the land-surface scheme SECHIBA, with a two-reservoir scheme for soil water storage and runoff at saturation. We studied a non-linear drainage parametrization, and a distributed surface runoff parametrization, accounting for the subgrid scale variability (SSV) of soil moisture capacity, through a distribution where the shape parameter was b. GCM results show that the drainage parametrization induces significant reductions in soil moisture and evaporation rate compared to the reference simulation. They are related to changes in moisture convergence in the tropics, and to a precipitation decrease in the extratropics. When drainage is implemented, the effect of the SSV parametrization (b=0.2) is also to reduce soil moisture and evaporation rates compared to the simulation with drainage only. These changes are much smaller than the former, but the sensitivity of the hydrological cycle to the SSV parametrization is shown to be larger in dry periods, and to be enhanced by an increase of the shape parameter b. The comparison of simulated total runoffs with observed data shows that the soil hydrological parametrizations does not reduce the GCM systematic errors in the annual water balance, but that they can improve the representation of the total runoff’s annual cycle.     

Comparison of vertical mixing schemes embedded in a global ocean GCM. Part II: Large scale circulation and water mass formation

Summary The sensitivity of the circulation and water mass properties in a global ocean circulation model (OGCM) to the stability dependent vertical mixing parameterisations of Pacanowski and Philander (1981) and Henderson-Sellers (1985) is investigated. The work extends a previous study which examined upper ocean charateristics and mixed layer evolution resulting from these schemes incorporated in the OGCM and made a recommendation as to the appropriateness of the latter scheme for global models. Under the assumption of constant vertical eddy coefficients (the control case), the model climatology displays acceptable values of North Atlantic Deep Water formation, Antarctic Circumpolar Current strength, and Indonesian throughflow but an excessively deep and diffuse pycnocline structure with weak stratification in the deep ocean. It is found that these circulation and water mass properties are sensitive to the choice of parametrisation of vertical mixing and that the two stability dependent schemes are unable to perform satisfactorily over the global domain, instead being better suited to the tropics. Under conditions optimal for representing the tropical current and temperature structure, these schemes result in significant weakening of major currents (particularly, the ACC) and reductions in the rates of deep water formation and poleward heat transports. These deficiencies can only be remedied at the expense of the improvements to the simulation in the tropical part of the domain. The results presented indicate that the suggestions made in the previous study do not extend to situations where the deep ocean, and particularly, the global thermohaline circulation is important.With 8 Figures     

Impact of global warming on the interannual and interdecadal climate modes in a coupled GCM

In this study, we investigated the impact of global warming on the variabilities of large-scale interannual and interdecadal climate modes and teleconnection patterns with two long-term integrations of the coupled general circulation model of ECHAM4/OPYC3 at the Max-Planck-Institute for Meteorology, Hamburg. One is the control (CTRL) run with fixed present-day concentrations of greenhouse gases. The other experiment is a simulation of transient greenhouse warming, named GHG run. In the GHG run the averaged geopotential height at 500?hPa is increased significantly, and a negative phase of the Pacific/North American (PNA) teleconnection-like distribution pattern is intensified. The standard deviation over the tropics (high latitudes) is enhanced (reduced) on the interdecadal time scales and reduced (enhanced) on the interannual time scales in the GHG run. Except for an interdecadal mode related to the Southern Oscillation (SO) in the GHG run, the spatial variation patterns are similar for different (interannual?+?interdecadal, interannual, and interdecadal) time scales in the GHG and CTRL runs. Spatial distributions of the teleconnection patterns on the interannual and interdecadal time scales in the GHG run are also similar to those in the CTRL run. But some teleconnection patterns show linear trends and changes of variances and frequencies in the GHG run. Apart from the positive linear trend of the SO, the interdecadal modulation to the El Niño/SO cycle is enhanced during the GHG 2040?～?2099. This is the result of an enhancement of the Walker circulation during that period. La Niña events intensify and El Niño events relatively weaken during the GHG 2070?～?2090. It is interesting to note that with increasing greenhouse gas concentrations the relation between the SO and the PNA pattern is reversed significantly from a negative to a positive correlation on the interdecadal time scales and weakened on the interannual time scales. This suggests that the increase of the greenhouse gas concentrations will trigger the nonstationary correlation between the SO and the PNA pattern both on the interdecadal and interannual time scales.     

Validation of an annual cycle simulation with a T42-L20 GCM

An annual cycle of an atmospheric general circulation model (AGCM) is presented. The winter and summer zonal averages of the atmospheric fields are compared with an observed climatology. The main features of the observed seasonal means are well reproduced by the model. One of the main discrepancies is that the simulated atmosphere is too cold, particularly in its upper part. Some other discrepancies might be explained by the interannual variability. The AGCM surface fluxes are directly compared to climatological estimates. On the other hand, the calculation of meridional heat transport by the ocean, inferred from the simulated energy budget, can be compared to transport induced from climatologies. The main result of this double comparison is that AGCM fluxes generally are within the range of climatological estimates. The main deficiency of the model is poor partitioning between solar and non-solar heat fluxes in the tropical belt. The meridional heat transport also reveals a significant energy-loss by the Northern Hemisphere ocean north of 45° N. The possible implications of model surface flux deficiencies on coupling with an oceanic model are discussed.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

GCM intercomparison of global cloud regimes: present-day evaluation and climate change response

The radiative feedback from clouds remains the largest source of variation in climate sensitivity amongst general circulation models (GCMs). A cloud clustering methodology is applied to six contemporary GCMs in order to provide a detailed intercomparison and evaluation of the simulated cloud regimes. By analysing GCMs in the context of cloud regimes, processes related to particular cloud types are more likely to be evaluated. In this paper, the mean properties of the global cloud regimes are evaluated, and the cloud response to climate change is analysed in the cloud-regime framework. Most of the GCMs are able to simulate the principal cloud regimes, however none of the models analysed have a good representation of trade cumulus in the tropics. The models also share a difficulty in simulating those regimes with cloud tops at mid-levels, with only ECHAM5 producing a regime of tropical cumulus congestus. Optically thick, high top cloud in the extra-tropics, typically associated with the passage of frontal systems, is simulated considerably too frequently in the ECHAM5 model. This appears to be a result of the cloud type persisting in the model after the meteorological conditions associated with frontal systems have ceased. The simulation of stratocumulus in the MIROC GCMs is too extensive, resulting in the tropics being too reflective. Most of the global-mean cloud response to doubled CO₂ in the GCMs is found to be a result of changes in the cloud radiative properties of the regimes, rather than changes in the relative frequency of occurrence (RFO) of the regimes. Most of the variance in the global cloud response between the GCMs arises from differences in the radiative response of frontal cloud in the extra-tropics and from stratocumulus cloud in the tropics. This variance is largely the result of excessively high RFOs of specific regimes in particular GCMs. It is shown here that evaluation and subsequent improvement in the simulation of the present-day regime properties has the potential to reduce the variance of the global cloud response, and hence climate sensitivity, amongst GCMs. For the ensemble of models considered in this study, the use of observations of the mean present-day cloud regimes suggests a potential reduction in the range of climate sensitivity of almost a third. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Vegetation and climate variability: a GCM modelling study

Vegetation is known to interact with the other components of the climate system over a wide range of timescales. Some of these interactions are now being taken into account in models for climate prediction. This study is an attempt to describe and quantify the climate–vegetation coupling at the interannual timescale, simulated with a General Circulation Model (HadSM3) coupled to a dynamic global vegetation model (TRIFFID). Vegetation variability is generally strongest in semi-arid areas, where it is driven by precipitation variability. The impact of vegetation variability on climate is analysed by using multivariate regressions of boundary layer fluxes and properties, with respect to soil moisture and vegetation fraction. Dynamic vegetation is found to significantly increase the variance in the surface sensible and latent heat fluxes. Vegetation growth always causes evapotranspiration to increase, but its impact on sensible heat is less straightforward. The feedback of vegetation on sensible heat is positive in Australia, but negative in the Sahel and in India. The sign of the feedback depends on the competing influences, at the gridpoint scale, of the turbulent heat exchange coefficient and the surface (stomatal) water conductance, which both increase with vegetation growth. The impact of vegetation variability on boundary layer potential temperature and relative humidity are shown to be small, implying that precipitation persistence is not strongly modified by vegetation dynamics in this model. We discuss how these model results may improve our knowledge of vegetation–atmosphere interactions and help us to target future model developments.     

Simple indices of global climate variability and change Part II: attribution of climate change during the twentieth century

Five simple indices of surface temperature are used to investigate the influence of anthropogenic and natural (solar irradiance and volcanic aerosol) forcing on observed climate change during the twentieth century. These indices are based on spatial fingerprints of climate change and include the global-mean surface temperature, the land-ocean temperature contrast, the magnitude of the annual cycle in surface temperature over land, the Northern Hemisphere meridional temperature gradient and the hemispheric temperature contrast. The indices contain information independent of variations in global-mean temperature for unforced climate variations and hence, considered collectively, they are more useful in an attribution study than global mean surface temperature alone. Observed linear trends over 1950–1999 in all the indices except the hemispheric temperature contrast are significantly larger than simulated changes due to internal variability or natural (solar and volcanic aerosol) forcings and are consistent with simulated changes due to anthropogenic (greenhouse gas and sulfate aerosol) forcing. The combined, relative influence of these different forcings on observed trends during the twentieth century is investigated using linear regression of the observed and simulated responses of the indices. It is found that anthropogenic forcing accounts for almost all of the observed changes in surface temperature during 1946–1995. We found that early twentieth century changes (1896–1945) in global mean temperature can be explained by a combination of anthropogenic and natural forcing, as well as internal climate variability. Estimates of scaling factors that weight the amplitude of model simulated signals to corresponding observed changes using a combined normalized index are similar to those calculated using more complex, optimal fingerprint techniques.     

大气环流的年代际变化 II.GCM数值模拟研究

类似大气环流模式比较计划（AMIP）的数值模拟,将实际观测的海表水温（SST）资料引入模式进行40多年的数值积分,得到长时间的大气环流模拟结果。分析数值模拟结果发现,无论是大气中的主要涛动和遥相关型,还是重要大气环流系统都极为清楚地存在着年代际变化特征,包括10～20年准周期振荡和可能的30年以上的准周期振荡;而且上述主要环流系统的形势及其年代际变化大都与实际观测资料所给出的结果相一致。顺便分析中国东部气候的模拟结果,年代际变化特征（包括60年代的气候突变）也很清楚,并同大气环流变化配合十分合理。结果也表明,同研究季节和年际变化一样,大气环流模式（AGCM）数值模拟也是研究大气环流和气候年代际变化的有效方法。     

Surface energy balance complexity in GCM land surface models. Part II: coupled simulations

 Global coupled simulations with the Bureau of Meteorology Research Centre climate model and the CHAmeleon Surface Model (CHASM) are used to examine how four general extensions to the representation of the basic land surface energy balance affect simulated land-atmosphere interface variables: evaporation, precipitation, skin temperature and air temperature. The impacts of including separate surface energy balance calculations for: vegetated and non-vegetated portions of the land surface; an explicit parametrisation of canopy resistance; explicit bare ground evaporation; and explicit canopy interception are isolated and quantified. The hypothesis that these aspects of surface energy balance parametrisation do not contain substantial information at the monthly time scale (and are therefore not important to consider in a land surface model) is shown to be false. Considerable sensitivity to each of the four general surface energy balance extensions is identified in average pointwise monthly changes for important land-atmosphere interface variables. Average pointwise changes in monthly precipitation and land evaporation are equal to about 40 and 31–37% of the global-average precipitation and land evaporation respectively. Average pointwise changes for land surface skin temperature and lowest model layer air temperature are about 2 and 0.9 K respectively. The average pointwise change and average pointwise biases are statistically significant at 95% in all cases. Substantial changes to zonally average variables are also identified. We demonstrate how the globally averaged surface resistance parameter can vary from 150 to 25 s/m depending on which aspects of the surface energy balance are treated implicitly. We also show that if interception is treated implicitly, the effective surface resistance must vary geographically in order to capture the behaviour of a model which treats this process explicitly. The implication of these results for the design of land surface models is discussed. Received: 8 July 1999 / Accepted: 1 September 2000     

Natural and anthropogenic climate change: incorporating historical land cover change,vegetation dynamics and the global carbon cycle

This study explores natural and anthropogenic influences on the climate system, with an emphasis on the biogeophysical and biogeochemical effects of historical land cover change. The biogeophysical effect of land cover change is first subjected to a detailed sensitivity analysis in the context of the UVic Earth System Climate Model, a global climate model of intermediate complexity. Results show a global cooling in the range of –0.06 to –0.22 °C, though this effect is not found to be detectable in observed temperature trends. We then include the effects of natural forcings (volcanic aerosols, solar insolation variability and orbital changes) and other anthropogenic forcings (greenhouse gases and sulfate aerosols). Transient model runs from the year 1700 to 2000 are presented for each forcing individually as well as for combinations of forcings. We find that the UVic Model reproduces well the global temperature data when all forcings are included. These transient experiments are repeated using a dynamic vegetation model coupled interactively to the UVic Model. We find that dynamic vegetation acts as a positive feedback in the climate system for both the all-forcings and land cover change only model runs. Finally, the biogeochemical effect of land cover change is explored using a dynamically coupled inorganic ocean and terrestrial carbon cycle model. The carbon emissions from land cover change are found to enhance global temperatures by an amount that exceeds the biogeophysical cooling. The net effect of historical land cover change over this period is to increase global temperature by 0.15 °C.     

Impact of resolving the diurnal cycle in an ocean–atmosphere GCM. Part 1: a diurnally forced OGCM

The diurnal cycle is a fundamental time scale in the climate system, at which the upper ocean and atmosphere are routinely observed to vary. Current climate models, however, are not configured to resolve the diurnal cycle in the upper ocean or the interaction of the ocean and atmosphere on these time scales. This study examines the diurnal cycle of the tropical upper ocean and its climate impacts. In the present paper, the first of two, a high vertical resolution ocean general circulation model (OGCM), with modified physics, is developed which is able to resolve the diurnal cycle of sea surface temperature (SST) and current variability in the upper ocean. It is then validated against a satellite derived parameterization of diurnal SST variability and in-situ current observations. The model is then used to assess rectification of the intraseasonal SST response to the Madden–Julian oscillation (MJO) by the diurnal cycle of SST. Across the equatorial Indo-Pacific it is found that the diurnal cycle increases the intraseasonal SST response to the MJO by around 20%. In the Pacific, the diurnal cycle also modifies the exchange of momentum between equatorially divergent Ekman currents and the meridionally convergent geostrophic currents beneath, resulting in a 10% increase in the strength of the Ekman cells and equatorial upwelling. How the thermodynamic and dynamical impacts of the diurnal cycle effect the mean state, and variability, of the climate system cannot be fully investigated in the constrained design of ocean-only experiments presented here. The second part of this study, published separately, addresses the climate impacts of the diurnal cycle in the coupled system by coupling the OGCM developed here to an atmosphere general circulation model.     

Potential impacts of global climate change on the hydrology and ecology of ephemeral freshwater systems of the forests of the northeastern United States

Global, national, and regional assessments of the potential effects of Global Climate Change (GCC) have been recently released, but not one of these assessments has specifically addressed the critical issue of the potential impacts of GCC on ephemeral freshwater systems (EFS). I suggest that this is a major oversight as EFS occur in various forms across the globe. In the northeastern United States, these systems, whether ephemeral (“vernal”) pools or ephemeral or intermittent headwater streams are abundant and provide unique habitats critical to the maintenance of forest biodiversity. Since the hydrology of these waterbodies is strongly affected by weather patterns (in the short-term) or climate (long-term), they are especially sensitive to climate change. In this essay, I review the literature on relationships between climate and hydrology of EFS and on relationships between hydrology and ecology of these systems. I then conclude with my assessment of potential impacts of GCC on the hydrology of EFS and implications for their ecology. The focus of this essay will be on EFS of the forests of the northeastern United States, but will include literature from other regions as they relate to the general relationships between GCC and EFS.     

Continental vegetation as a dynamic component of a global climate model: a preliminary assessment

A simplified vegetation distribution prediction scheme is used in combination with the Biosphere-Atmosphere Transfer Scheme (BATS) and coupled to a version of the NCAR Community Climate Model (CCM1) which includes a mixed-layer ocean. Employed in an off-line mode as a diagnostic tool, the scheme predicts a slightly darker and slightly rougher continental surface than when BATS' prescribed vegetation classes are used. The impact of tropical deforestation on regional climates, and hence on diagnosed vegetation, differs between South America and S.E. Asia. In the Amazon, the climatic effects of removing all the tropical forest are so marked that in only one of the 18 deforested grid elements could the new climate sustain tropical forest vegetation whereas in S.E. Asia in seven of the 9 deforested elements the climate could continue to support tropical forest. Following these off-line tests, the simple vegetation scheme has been coupled to the GCM as an interactive (or two-way) submodel for a test integration lasting 5.6 yr. It is found to be a stable component of the global climate system, producing only ~ 3% (absolute) interannual changes in the predicted percentages of continental vegetation, together with globally-averaged continental temperature increases of up to + 1.5 °C and evaporation increases of 0 to 5 W m^–2 and no discernible trends over the 67 months of integration. On the other hand, this interactive land biosphere causes regional-scale temperature differences of ± 10 °C and commensurate disturbances in other climatic parameters. Tuning, similar to the q-flux schemes used for ocean models, could improve the simulation of the present-day surface climate but, in the longer term, it will be important to focus on predicting the characteristics of the continental surface rather than simple vegetation classes. The coupling scheme will also have to allow for vegetation responses occurring over longer timescales so that the coupled system is buffered from sudden shocks.