Testing a GCM land surface scheme against catchment-scale runoff data

A GCM land surface scheme was used, in off-line mode, to simulate the runoff, latent and sensible heat fluxes for two distinct Australian catchments using observed atmospheric forcing. The tropical Jardine River catchment is 2500 km² and has an annual rainfall of 1700 mm y^–1 while the Canning River catchment is 540 km², has a Mediterranean climate (annual rainfall of 800 mm y^–1) and is ephemeral for half the year. It was found that the standard version of a land surface scheme developed for a GCM, and initialised as for incorporation into a GCM, simulated similar latent and sensible heat fluxes compared to a basin-scale hydrological model (MODHYDROLOG) which was calibrated for each catchment. However, the standard version of the land surface scheme grossly overestimated the observed peak runoff in the wet Jardine River catchment at the expense of runoff later in the season. Increasing the soil water storage permitted the land surface scheme to simulate observed runoff quite well, but led to a different simulation of latent and sensible heat compared to MODHYDROLOG. It is concluded that this 2-layer land surface scheme was unable to simulate both catchments realistically. The land surface scheme was then extended to a three-layer model. In terms of runoff, the resulting control simulations with soil depths chosen as for the GCM were better than the best simulations obtained with the two-layer model. The three-layer model simulated similar latent and sensible heat for both catchments compared to MODHYDROLOG. Unfortunately, for the ephemeral Canning River catchment, the land surface scheme was unable to time the observed runoff peak correctly. A tentative conclusion would be that this GCM land surface scheme may be able to simulate the present day state of some larger and wetter catchments but not catchments with peaky hydrographs and zero flows for part of the year. This conclusion requires examination with a range of GCM land surface schemes against a range of catchments. Received: 9 June 1995 / Accepted: 4 April 1996     

The impact of new land surface physics on the GCM simulation of climate and climate sensitivity

 Recent improvements to the Hadley Centre climate model include the introduction of a new land surface scheme called “MOSES” (Met Office Surface Exchange Scheme). MOSES is built on the previous scheme, but incorporates in addition an interactive plant photosynthesis and conductance module, and a new soil thermodynamics scheme which simulates the freezing and melting of soil water, and takes account of the dependence of soil thermal characteristics on the frozen and unfrozen components. The impact of these new features is demonstrated by comparing 1×CO₂ and 2×CO₂ climate simulations carried out using the old (UKMO) and new (MOSES) land surface schemes. MOSES is found to improve the simulation of current climate. Soil water freezing tends to warm the high-latitude land in the northern Hemisphere during autumn and winter, whilst the increased soil water availability in MOSES alleviates a spurious summer drying in the mid-latitudes. The interactive canopy conductance responds directly to CO₂, supressing transpiration as the concentration increases and producing a significant enhancement of the warming due to the radiative effects of CO₂ alone. Received: 16 March 1998 / Accepted: 4 August 1998     

GCM sensitivity experiments with locally modified land surface properties over tropical South America

Ensembles of 1-year-long experiments with a relatively high-resolution ECMWF model were conducted in order to investigate the impact of modified land surface properties on local, regional and large-scale atmospheric circulations. The modifications consisted of changes to land cover and increased albedo over the northern part of South America. In many respects the experimental design resembles the setting of classical deforestation experiments. The local model response to imposed modifications, which includes a reduction in precipitation as well as in evaporation and an increase in surface temperature, was found to be stronger in dry (July–September, JAS) than in wet (January–March, JFM) season, and in the ensemble with higher albedo value. Local drying is discussed in terms of locally generated overturning that resembles a direct thermal circulation. The effects of this circulation seem to be dominant over the reduction in large-scale moisture supply from the adjacent ocean. On large scales, changes to the Pacific branch of the Walker circulation lead, through modified divergent flow, to a tropics-wide impact on precipitation. In addition to South America, the largest changes are seen in the south Pacific convergence zone in JFM, while the impact on the Atlantic inter-tropical convergence zone is stronger in JAS. In the extratropics, there is little change in precipitation. In the upper troposphere, a distinctive teleconnection wave-pattern could be seen in the Pacific/North American region during JFM. A notable feature in the upper-air model response in JAS is a wave train extending from South America, over the northern Atlantic into Europe. With regard to the interaction between the land surface response and model systematic errors, our results suggest that the erroneous shift of the downward branch of the Pacific/South American Walker circulation is likely to be a cause, rather than a consequence, of the rainfall deficit over South America in the model climatology.     

Effects of land cover change on the tropical circulation in a GCM

Multivariate statistics are used to investigate sensitivity of the tropical atmospheric circulation to scenario-based global land cover change (LCC), with the largest changes occurring in the tropics. Three simulations performed with the fully coupled Parallel Climate Model (PCM) are compared: (1) a present day control run; (2) a simulation with present day land cover and Intergovernmental Panel on Climate Change (IPCC) Special Report on Emission Scenarios (SRES) A2 greenhouse gas (GHG) projections; and (3) a simulation with SRES A2 land cover and GHG projections. Dimensionality of PCM data is reduced by projection onto a priori specified eigenvectors, consisting of Rossby and Kelvin waves produced by a linearized, reduced gravity model of the tropical circulation. A Hotelling T ² test is performed on projection amplitudes. Effects of LCC evaluated by this method are limited to diabatic heating. A statistically significant and recurrent signal is detected for 33% of all tests performed for various combinations of parameters. Taking into account uncertainties and limitations of the present methodology, this signal can be interpreted as a Rossby wave response to prescribed LCC. The Rossby waves are shallow, large-scale motions, trapped at the equator and most pronounced in boreal summer. Differences in mass and flow fields indicate a shift of the tropical Walker circulation patterns with an anomalous subsidence over tropical South America.     

Evaluation of an urban land surface scheme over a tropical suburban neighborhood

A two-step statistical downscaling method has been reviewed and adapted to simulate twenty-first-century climate projections for the Gulf of Fonseca (Central America, Pacific Coast) using Coupled Model Intercomparison Project (CMIP5) climate models. The downscaling methodology is adjusted after looking for good predictor fields for this area (where the geostrophic approximation fails and the real wind fields are the most applicable). The method’s performance for daily precipitation and maximum and minimum temperature is analysed and revealed suitable results for all variables. For instance, the method is able to simulate the characteristic cycle of the wet season for this area, which includes a mid-summer drought between two peaks. Future projections show a gradual temperature increase throughout the twenty-first century and a change in the features of the wet season (the first peak and mid-summer rainfall being reduced relative to the second peak, earlier onset of the wet season and a broader second peak).     

The impact of climatic and non-climatic factors on land surface temperature in southwestern Romania

Land surface temperature is one of the most important parameters related to global warming. It depends mainly on soil type, discontinuous vegetation cover, or lack of precipitation. The main purpose of this paper is to investigate the relationship between high LST, synoptic conditions and air masses trajectories, vegetation cover, and soil type in one of the driest region in Romania. In order to calculate the land surface temperature and normalized difference vegetation index, five satellite images of LANDSAT missions 5 and 7, covering a period of 26 years (1986–2011), were selected, all of them collected in the month of June. The areas with low vegetation density were derived from normalized difference vegetation index, while soil types have been extracted from Corine Land Cover database. HYSPLIT application was employed to identify the air masses origin based on their backward trajectories for each of the five study cases. Pearson, logarithmic, and quadratic correlations were used to detect the relationships between land surface temperature and observed ground temperatures, as well as between land surface temperature and normalized difference vegetation index. The most important findings are: strong correlation between land surface temperature derived from satellite images and maximum ground temperature recorded in a weather station located in the area, as well as between areas with land surface temperature equal to or higher than 40.0 °C and those with lack of vegetation; the sandy soils are the most prone to high land surface temperature and lack of vegetation, followed by the chernozems and brown soils; extremely severe drought events may occur in the region.     

The surface heat flux density of a bare soil

Abstract

Half‐hourly measurements of soil surface heat flux density (G₀ ), solar irradiance (S), and the surface energy balance components were made at Agassiz, b.c., in the spring and early summer of 1978 at two adjacent bare‐soil sites, one of which was culti‐packed while the other was disc‐harrowed. G₀ was calculated using the null‐alignment procedure from half‐hourly measurements of soil temperature at 30 depths down to 1 m, and volumetric soil heat capacity calculated from measurements of bulk density, organic matter fraction, and moisture content. The latent and sensible heat flux densities were measured using the energy balance/Bowen ratio technique.

It was found that both the daily averages and diurnal variations of Go at each site were not affected as the soil surface dried, despite reductions in evaporation rate of as much as 50% at the culti‐packed site and 75% at the disc‐harrowed site on the clear dry‐soil days. Diurnal variations of G₀ at the disc‐harrowed site were about 25% less than at the culti‐packed site, although daily averages were similar at both sites. Daily and daytime averages of G₀ at each site were linear functions of S alone, or functions of net radiation and some measure of near‐surface soil water content. Night‐time averages of G₀ at each site were linear functions of a cloudiness ratio equal to the fraction received of the clear‐day S.     

Effect of surface layer thickness on simultaneous transport of heat and water in a bare soil and its implications for land surface schemes

Abstract

A physically‐based multi‐layer numerical model is developed to determine the coupled transport of heat and water in the soil and in the soil‐atmosphere boundary layer. Using inputs of standard weather data and initial soil conditions the model is capable of predicting the surface energy balance components as well as water content and temperature profiles in the soil. It is used to predict these variables for a bare silt loam soil under two tillage treatments, viz. culti‐packed and left loose after disc‐harrowing, and the predicted results are compared with measurements. Very good agreement between the model predictions and measured evaporation and heat fluxes and soil water and temperatures for a ten‐day period shows that the model is capable of simulating the coupled transport of soil heat and soil water and their transfer across the soil surface‐atmosphere interface adequately.

Model predictions were compared with those of CLASS (Canadian Land Surface Scheme). It is shown that CLASS, version 2.6, provides good estimates of evaporation and hence the latent heat flux density, Q_E, under wetter soil conditions, but overestimates Q_E at moderately wet soil conditions and underestimates it under dry soil conditions. Under dry to moderately wet soil conditions the calculation of evaporation from bare soil is very sensitive to the thickness of the top layer particularly as the thickness approaches 10 cm.     

Probabilistic simulations of the impact of increasing leaf-level atmospheric carbon dioxide on the global land surface

Using a climate model with a sophisticated land surface scheme, simulations were conducted to explore the impact of increases in leaf-level carbon dioxide (CO₂) on evaporation, temperature and other land surface quantities. Fifty-one realizations were run, for each of four Januarys and four Julys for CO₂ concentrations at leaf-level of 280, 375, 500, 650, 840 and 1,000 ppmv. Atmospheric CO₂ concentration was held constant at 375 ppmv in all experiments. Statistically significant decreases in evaporation and increases in temperature occur in specific regions as leaf-level CO₂ is increased from 280 to 375 ppmv. These same areas expand geographically, and the magnitude of the changes increase as leaf-level CO₂ is increased further suggesting that changes are caused by the increase in leaf-level CO₂ and are not internal model variability. As leaf-level CO₂ is increased further, larger areas of the continental surface are affected by increasing amounts and a statistically significant change in precipitation is seen. The increase in leaf-level CO₂ from 280 ppmv to 375 ppmv causes statistically significant changes in the evaporation over 12% of continental surfaces in July. This increases to 25% at 500 ppmv, 35% at 650 ppmv, 41% at 840 ppmv and 47% at 1,000 ppmv. This affects temperature and rainfall by similar amounts, generally in coincident regions. An analysis of these results over key regions shows that the probability density functions of the latent heat flux and temperature are affected non-uniformly. There is a shift in the latent heat flux probability density function to lower values, mainly through the reduction in the upper tail of the distribution. The temperature probability density function shifts to higher values, mainly through an increase in the upper tail of the distribution indicating that the impact is focussed on extremes. Given that there are a suite of well evaluated land surface models that include the biogeochemical effects of increasing CO₂ we suggest that the inclusion of such a model should be a recommended component of climate models used in future assessment reports by the Intergovernmental Panel on Climate Change.     

A statistical study of the regional impact of deforestation on climate in the LMD GCM

The present study uses the general circulation model of the Laboratoire de Météorologie Dynamique (LMD-GCM) coupled to the land-surface, vegetation model SECHIBA. The impact of deforestation on climate is discussed. Replacing tropical forests by degraded pastures changes albedo, the roughness length and the hydrological properties of the surface. The experiment was carried out over eleven years using the observed sea surface temperature from 1978 to 1988, which includes two major El Niño events. The discussion of the results in this study is limited to the regional impact of deforestation. The changes found for the surface fluxes in Amazonia, Africa and Indonesia are examined in detail and compared in order to understand the impact on temperature. Special attention is paid to feedback mechanisms which compensate for the surface changes and to the statistical significance of these results within the tropical variability of climate. It is shown that the relatively small regional impact of deforestation in this study is statistically significant and largely independent of the El Niño-Southern Oscillation phenomenon.This paper was presented at the Second International Conference on Modelling of Global Climate Change and Variability (Hamburg, September 1992)     

GCM simulations of the Last Glacial Maximum surface climate of Greenland and Antarctica

 The LMDz variable grid GCM was used to simulate the Last Glacial Maximum (LGM, 21 ky Bp.) climate of Greenland and Antarctica at a spatial resolution of about 100 km.The high spatial resolution allows to investigate the spatial variability of surface climate change signals, and thus to address the question whether the sparse ice core data can be viewed as representative for the regional scale climate change. This study addresses primarily surface climate parameters because these can be checked against the, limited, ice core record. The changes are generally stronger for Greenland than for Antarctica, as the imposed changes of the forcing boundary conditions (e.g., sea surface temperatures) are more important in the vicinity of Greenland. Over Greenland, and to a limited extent also in Antarctica, the climate shows stronger changes in winter than in summer. The model suggests that the linear relationship between the surface temperature and inversion strength is modified during the LGM. The temperature dependency of the moisture holding capacity of the atmosphere alone cannot explain the strong reduction in snowfall over central Greenland; atmospheric circulation changes also play a crucial role. Changes in the high frequency variability of snowfall, atmospheric pressure and temperature are investigated and possible consequences for the interpretation of ice core records are discussed. Using an objective cyclone tracking scheme, the importance of changes of the atmospheric dynamics off the coasts of the ice sheets, especially for the high frequency variability of surface climate parameters, is illustrated. The importance of the choice of the LGM ice sheet topography is illustrated for Greenland, where two different topographies have been used, yielding results that differ quite strongly in certain nontrivial respects. This means that the paleo-topography is a significant source of uncertainty for the modelled paleoclimate. The sensitivity of the Greenland LGM climate to the prescribed sea surface conditions is examined by using two different LGM North Atlantic data sets. Received: 23 October 1997 / Accepted: 17 March 1998     

Atmospheric GCM response to an idealized anomaly of the Mediterranean sea surface temperature

The atmospheric response to an idealized 2 K cooling of the Mediterranean Sea is studied in a general circulation model (GCM). In the downstream region, from the eastern Mediterranean basin to the whole Asian continent (following the sub-tropical jet-stream), a baroclinic response is produced with high pressure near the surface and low pressure in the upper layers. It is the direct response to the Mediterranean cooling and it needs only a few days to be established. Teleconnections with strong zonal characteristics are found following the sub-tropical jet-stream to the North Pacific and North America, propagating further across the North Atlantic towards Northern Europe. Two remarkable remote features are the deepening of the Aleutian Low in the North Pacific and the weakening of the Icelandic Low in the North Atlantic. These two teleconnections form and grow in several days in the North Pacific and several tens of days in the North Atlantic. Both have a quasi-barotropic vertical structure. It is believed that they are the consequence of complex interactions between the mean flow and the transient eddies in the atmosphere.     

Sensitivity of a GCM simulation to subgrid infiltration and surface runoff

A subgrid parameterization of infiltration and surface runoff was evaluated using a land surface model coupled to an atmospheric general circulation model. Averaged over 5 year simulations, the subgrid parameterization resulted in significantly less infiltration of water into the soil compared to a simulation without subgrid hydrologic processes. As a result, the soils were drier, latent heat flux decreased, and surface air temperature increased. These results are consistent with other studies of subgrid hydrologic parameterizations, which also resulted in drier soils, decreased latent heat flux, and warmer surface temperatures. Several river basins were studied in detail. In the Amazon and Lena basins, the subgrid parameterization resulted in better annual runoff compared to observed annual river flow; surface air temperature was unchanged in the Amazon and better, compared to observations, in the Lena. In the Ob, Yenisey, and Amur basins, the subgrid parameterization resulted in too much annual runoff; July surface air temperature was unchanged or worse (Amur). Annual runoff for the Mississippi basin was better with the subgrid parameterization, but July surface air temperature was worse. These results suggest the utility of subgrid hydrologic parameterizations vary among river basins depending on the relative importance of Horton and Dunne runoff and the geologic factors affecting runoff generation.     

The effect of a regional increase in ocean surface roughness on the tropospheric circulation: a GCM experiment

The sensitivity of the atmospheric circulation to an increase in ocean surface roughness in the Southern Hemisphere storm track is investigated in a paired general circulation model experiment. Such a change in sea roughness could be induced by ocean waves generated by storms. Two extended permanent-July runs are made. One with standard sea surface roughness, the other with ten times as a large surface roughness over open sea poleward of 40° S. The regional increase in ocean surface roughness significantly modifies the tropospheric circulation in the Southern Hemisphere. The strongest effect is the reduction of tropospheric winds (by 2 m/s or 10%) above the area with increased roughness. The poleward eddy momentum flux is reduced in the upper troposphere and the meridional eddy sensible heat flux is reduced in the lower troposphere. Zonal mean and eddy kinetic energy are consistently reduced.     

The effect of land surface changes on Eemian climate

Transient experiments for the Eemian (128–113 ky BP) were performed with a complex, coupled earth system model, including atmosphere, ocean, terrestrial biosphere and marine biogeochemistry. In order to investigate the effect of land surface parameters (background albedo, vegetation and tree fraction and roughness length) on the simulated changes during the Eemian, simulations with interactive coupling between climate and vegetation were compared with additional experiments in which these feedbacks were suppressed. The experiments show that the influence of land surface on climate is mainly caused by changes in the albedo. For the northern hemisphere high latitudes, land surface albedo is changed partially due to the direct albedo effect of the conversion of grasses into forest, but the indirect effect of forests on snow albedo appears to be the major factor influencing the total absorption of solar radiation. The Western Sahara region experiences large changes in land surface albedo due to the appearance of vegetation between 128 and 120 ky BP. These local land surface albedo changes can be as much as 20%, thereby affecting the local as well as the global energy balance. On a global scale, latent heat loss over land increases more than 10% for 126 ky BP compared to present-day.     

A comparison of a GCM response to historical anthropogenic land cover change and model sensitivity to uncertainty in present-day land cover representations

This study assesses the sensitivity of the fully coupled NCAR-DOE PCM to three different representations of present-day land cover, based on IPCC SRES land cover information. We conclude that there is significant model sensitivity to current land cover characterization, with an observed average global temperature range of 0.21 K between the simulations. Much larger contrasts (up to 5 K) are found on the regional scale; however, these changes are largely offsetting on the global scale. These results show that significant biases can be introduced when outside data sources are used to conduct anthropogenic land cover change experiments in GCMs that have been calibrated to their own representation of present-day land cover. We conclude that hybrid systems that combine the natural vegetation from the native GCM datasets combined with human land cover information from other sources are best for simulating such impacts. We also performed a prehuman simulation, which had a 0.39 K ~higher average global temperature and, perhaps of greater importance, temperature changes regionally of about 2 K. In this study, the larger regional changes coincide with large-scale agricultural areas. The initial cooling from energy balance changes appear to create feedbacks that intensify mid-latitude circulation features and weaken the summer monsoon circulation over Asia, leading to further cooling. From these results, we conclude that land cover change plays a significant role in anthropogenically forced climate change. Because these changes coincide with regions of the highest human population this climate impact could have a disproportionate impact on human systems. Therefore, it is important that land cover change be included in past and future climate change simulations.     

Effects of land use change on North American climate: impact of surface datasets and model biogeophysics

This study examines the impact of historical land-cover change on North American surface climate, focusing on the robustness of the climate signal with respect to representation of sub-grid heterogeneity and land biogeophysics within a climate model. We performed four paired climate simulations with the Community Atmosphere Model using two contrasting land models and two different representations of land-cover change. One representation used a biome classification without subgrid-scale heterogeneity while the other used high-resolution satellite data to prescribe multiple vegetation types within a grid cell. Present-day and natural vegetation datasets were created for both representations. All four sets of climate simulations showed that present-day vegetation has cooled the summer climate in regions of North America compared to natural vegetation. The simulated magnitude and spatial extent of summer cooling due to land-cover change was reduced when the biome-derived land-cover change datasets were replaced by the satellite-derived datasets. The diminished cooling is partly due to reduced intensity of agriculture in the satellite-derived datasets. Comparison of the two land-surface models showed that the use of a comparatively warmer and drier land model in conjunction with satellite-derived datasets further reduced the simulated magnitude of summer cooling. These results suggest that the cooling signal associated with North American land-cover change is robust but the magnitude and therefore detection of the signal depends on the realism of the datasets used to represent land-cover change and the parametrisation of land biogeophysics.