We describe the baseline model configuration and simulation characteristics of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Land Model version 4.1 (LM4.1), which builds on component and coupled model developments over 2013–2019 for the coupled carbon-chemistry-climate Earth System Model Version 4.1 (ESM4.1) simulation as part of the sixth phase of the Coupled Model Intercomparison Project. Analysis of ESM4.1/LM4.1 is focused on biophysical and biogeochemical processes and interactions with climate. Key features include advanced vegetation dynamics and multi-layer canopy energy and moisture exchanges, daily fire, land use representation, and dynamic atmospheric dust coupling. We compare LM4.1 performance in the GFDL Earth System Model (ESM) configuration ESM4.1 to the previous generation component LM3.0 in the ESM2G configuration. ESM4.1/LM4.1 provides significant improvement in the treatment of ecological processes from GFDL's previous generation models. However, ESM4.1/LM4.1 likely overestimates the influence of land use and land cover change on vegetation characteristics, particularly on pasturelands, as it overestimates the competitiveness of grasses versus trees in the tropics, and as a result, underestimates present-day biomass and carbon uptake in comparison to observations.
The local climatic impacts of historical expansion of irrigation are substantial, but the distant impacts are poorly understood, and their governing mechanisms generally have not been rigorously analyzed. Our experiments with an earth-system model suggest that irrigation in the Middle East and South Asia may enhance rainfall in a large portion of the Sahel-Sudan Savanna (SSS) to an extent comparable and opposite to its suppression by other anthropogenic climate drivers during the last several decades. The enhancement arises through a reduction in the meridional gradient of moist static energy from the Sahara Desert to the tropical rainforests. An implication of this study is that remote irrigation is a possible factor affecting the risk of drought and famine and, thus, future water security in the SSS region.
We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR (Seamless System for Prediction and EArth System Research) is built from component models recently developed at GFDL ‐ the AM4 atmosphere model, MOM6 ocean code, LM4 land model and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0o (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1o to 0.25o. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM 4 models, but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time‐mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM 4 models.
We describe the baseline coupled model configuration and simulation characteristics of GFDL's Earth System Model Version 4.1 (ESM4.1), which builds on component and coupled model developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation contributing to the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's CM4.0 development effort that focuses on ocean resolution for physical climate, ESM4.1 focuses on comprehensiveness of Earth system interactions. ESM4.1 features doubled horizontal resolution of both atmosphere (2° to 1°) and ocean (1° to 0.5°) relative to GFDL's previous‐generation coupled ESM2‐carbon and CM3‐chemistry models. ESM4.1 brings together key representational advances in CM4.0 dynamics and physics along with those in aerosols and their precursor emissions, land ecosystem vegetation and canopy competition, and multiday fire; ocean ecological and biogeochemical interactions, comprehensive land‐atmosphere‐ocean cycling of CO2, dust and iron, and interactive ocean‐atmosphere nitrogen cycling are described in detail across this volume of JAMES and presented here in terms of the overall coupling and resulting fidelity. ESM4.1 provides much improved fidelity in CO2 and chemistry over ESM2 and CM3, captures most of CM4.0's baseline simulations characteristics, and notably improves on CM4.0 in (1) Southern Ocean mode and intermediate water ventilation, (2) Southern Ocean aerosols, and (3) reduced spurious ocean heat uptake. ESM4.1 has reduced transient and equilibrium climate sensitivity compared to CM4.0. Fidelity concerns include (1) moderate degradation in sea surface temperature biases, (2) degradation in aerosols in some regions, and (3) strong centennial scale climate modulation by Southern Ocean convection.
The sensitivity of river discharge to climate-system warming is highly uncertain, and the processes that govern river discharge are poorly understood, which impedes climate-change adaptation. A prominent exemplar is the Colorado River, where meteorological drought and warming are shrinking a water resource that supports more than 1 trillion dollars of economic activity per year. A Monte Carlo simulation with a radiation-aware hydrologic model resolves the longstanding, wide disparity in sensitivity estimates and reveals the controlling physical processes. We estimate that annual mean discharge has been decreasing by 9.3% per degree Celsius of warming because of increased evapotranspiration, mainly driven by snow loss and a consequent decrease in reflection of solar radiation. Projected precipitation increases likely will not suffice to fully counter the robust, thermodynamically induced drying. Thus, an increasing risk of severe water shortages is expected.
We describe GFDL's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the pre‐industrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasi‐periodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.
Although interannual streamflow variability is primarily a result of precipitation variability, temperature also plays a role. The relative weakness of the temperature effect at the annual time scale hinders understanding, but may belie substantial importance on climatic time scales. Here we develop and evaluate a simple theory relating variations of streamflow and evapotranspiration (E) to those of precipitation (P) and temperature. The theory is based on extensions of the Budyko water‐balance hypothesis, the Priestley‐Taylor theory for potential evapotranspiration (Ep), and a linear model of interannual basin storage. The theory implies that the temperature affects streamflow by modifying evapotranspiration through a Clausius‐Clapeyron‐like relation and through the sensitivity of net radiation to temperature. We apply and test (1) a previously introduced “strong” extension of the Budyko hypothesis, which requires that the function linking temporal variations of the evapotranspiration ratio (E/P) and the index of dryness (Ep/P) at an annual time scale is identical to that linking inter‐basin variations of the corresponding long‐term means, and (2) a “weak” extension, which requires only that the annual evapotranspiration ratio depends uniquely on the annual index of dryness, and that the form of that dependence need not be known a priori nor be identical across basins. In application of the weak extension, the readily observed sensitivity of streamflow to precipitation contains crucial information about the sensitivity to potential evapotranspiration and, thence, to temperature. Implementation of the strong extension is problematic, whereas the weak extension appears to capture essential controls of the temperature effect efficiently.
Human communities, economies, and natural ecosystems require a reliable supply of water, and this often is provided by rivers. It has at times been noted that rivers deliver less water in warm years than in cool years, even after adjustment for variations in precipitation. This dependence on temperature raises concerns about the effect of heat waves or climatic warming on water supply. Just why and how much the flow of rivers depends on temperature has not been well understood, but the answers to these questions are relevant for ensuring future water security. Here we present and evaluate a process‐based theory that attempts to answer both questions. In those river basins where long‐term observations of river flow, temperature, and precipitation data are available, the theory is consistent overall with observed sensitivities of river flows to temperature. This success implies potential applicability of the theory also where such observations are not available.
In this two-part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part I, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode – with prescribed sea surface temperatures (SSTs) and sea ice distribution – is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top-of-atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part II, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.
In Part II of this two-part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part I. Part II provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.
Milly, P C., and Krista A Dunne, August 2017: A Hydrologic Drying Bias in Water-Resource Impact Analyses of Anthropogenic Climate Change. Journal of the American Water Resources Association, 53(4), DOI:10.1111/1752-1688.12538. Abstract
For water-resource planning, sensitivity of freshwater availability to anthropogenic climate change (ACC) often is analyzed with “offline” hydrologic models that use precipitation and potential evapotranspiration (Ep) as inputs. Because Ep is not a climate-model output, an intermediary model of Ep must be introduced to connect the climate model to the hydrologic model. Several Ep methods are used. The suitability of each can be assessed by noting a credible Ep method for offline analyses should be able to reproduce climate models’ ACC-driven changes in actual evapotranspiration in regions and seasons of negligible water stress (Ew). We quantified this ability for seven commonly used Ep methods and for a simple proportionality with available energy (“energy-only” method). With the exception of the energy-only method, all methods tend to overestimate substantially the increase in Ep associated with ACC. In an offline hydrologic model, the Ep-change biases produce excessive increases in actual evapotranspiration (E), whether the system experiences water stress or not, and thence strong negative biases in runoff change, as compared to hydrologic fluxes in the driving climate models. The runoff biases are comparable in magnitude to the ACC-induced runoff changes themselves. These results suggest future hydrologic drying (wetting) trends likely are being systematically and substantially overestimated (underestimated) in many water-resource impact analyses.
By various measures (drought area1 and intensity2, climatic aridity index3, and climatic water deficits4), some observational analyses have suggested that much of the Earth’s land has been drying during recent decades, but such drying seems inconsistent with observations of dryland greening and decreasing pan evaporation5. ‘Offline’ analyses of climate-model outputs from anthropogenic climate change (ACC) experiments portend continuation of putative drying through the twenty-first century3, 6, 7, 8, 9, 10, despite an expected increase in global land precipitation9. A ubiquitous increase in estimates of potential evapotranspiration (PET), driven by atmospheric warming11, underlies the drying trends4, 8, 9, 12, but may be a methodological artefact5. Here we show that the PET estimator commonly used (the Penman–Monteith PET13 for either an open-water surface1, 2, 6, 7, 12 or a reference crop3, 4, 8, 9, 11) severely overpredicts the changes in non-water-stressed evapotranspiration computed in the climate models themselves in ACC experiments. This overprediction is partially due to neglect of stomatal conductance reductions commonly induced by increasing atmospheric CO2 concentrations in climate models5. Our findings imply that historical and future tendencies towards continental drying, as characterized by offline-computed runoff, as well as other PET-dependent metrics, may be considerably weaker and less extensive than previously thought.
“LM3” is a new model of terrestrial water, energy, and carbon, intended for use in global hydrologic analyses and as a component of earth-system and physical-climate models. It is designed to improve upon the performance and extend the scope of the predecessor Land Dynamics (LaD) and LM3V models, by quantifying better the physical controls of climate and biogeochemistry and by relating more directly to components of the global water system that touch human concerns. LM3 includes multi-layer representations of temperature, liquid-water content, and ice content of both snow pack and macroporous soil/bedrock; topography-based description of saturated area and groundwater discharge; and transport of runoff to the ocean via a global river and lake network. Sensible heat transport by water mass is accounted throughout for a complete energy balance. Carbon and vegetation dynamics and biophysics are represented as in the model LM3V. In numerical experiments, LM3 avoids some of the limitations of the LaD model and provides qualitatively (though not always quantitatively) reasonable estimates, from a global perspective, of observed spatial and/or temporal variations of vegetation density, albedo, streamflow, water-table depth, permafrost, and lake levels. Amplitude and phase of annual cycle of total water storage are simulated well. Realism of modeled lake levels varies widely. The water table tends to be consistently too shallow in humid regions. Biophysical properties have an artificial step-wise spatial structure, and equilibrium vegetation is sensitive to initial conditions. Explicit resolution of thick (>100 m) unsaturated zones and permafrost is possible, but only at the cost of long (>>300 y) model spin-up times.
We describe carbon system formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate good climate fidelity as described in Part I while incorporating explicit and consistent carbon dynamics. The two models differ almost exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. On land, both ESMs include a revised land model to simulate competitive vegetation distributions and functioning, including carbon cycling among vegetation, soil and atmosphere. In the ocean, both models include new biogeochemical algorithms including phytoplankton functional group dynamics with flexible stoichiometry. Preindustrial simulations are spun up to give stable, realistic carbon cycle means and variability. Significant differences in simulation characteristics of these two models are described. Due to differences in oceanic ventilation rates (Part I) ESM2M has a stronger biological carbon pump but weaker northward implied atmospheric CO2 transport than ESM2G. The major advantages of ESM2G over ESM2M are: improved representation of surface chlorophyll in the Atlantic and Indian Oceans and thermocline nutrients and oxygen in the North Pacific. Improved tree mortality parameters in ESM2G produced more realistic carbon accumulation in vegetation pools. The major advantages of ESM2M over ESM2G are reduced nutrient and oxygen biases in the Southern and Tropical Oceans.
We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.
Milly, P C., and Krista A Dunne, January 2011: On the hydrologic adjustment of climate-model projections: The potential pitfall of potential evapotranspiration. Earth Interactions, 15(1), DOI:10.1175/2010EI363.1. Abstract
Hydrologic models often are applied to adjust projections of
hydroclimatic change that come from climate models. Such adjustment includes
climate-bias correction, spatial refinement (‘‘downscaling’’), and consideration
of the roles of hydrologic processes that were neglected in the climate model.
Described herein is a quantitative analysis of the effects of hydrologic adjustment
on the projections of runoff change associated with projected twenty-first-century
climate change. In a case study including three climatemodels and 10 river basins
in the contiguous United States, the authors find that relative (i.e., fractional or
percentage) runoff change computed with hydrologic adjustment more often
than not was less positive (or, equivalently, more negative) than what was projected
by the climate models. The dominant contributor to this decrease in runoff
was a ubiquitous change in runoff (median 211%) caused by the hydrologic
model’s apparent amplification of the climate-model-implied growth in potential
evapotranspiration. Analysis suggests that the hydrologic model, on the basis of
the empirical, temperature-based modified Jensen–Haise formula, calculates a
change in potential evapotranspiration that is typically 3 times the change implied
by the climate models, which explicitly track surface energy budgets. In comparison
with the amplification of potential evapotranspiration, central tendencies
of other contributions from hydrologic adjustment (spatial refinement, climate-bias adjustment, and process refinement) were relatively small. The authors’ findings
highlight the need for caution when projecting changes in potential evapotranspiration
for use in hydrologic models or drought indices to evaluate climate change
impacts on water.
The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.
Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.
The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic.
Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).
Manuscript received 8 December 2004, in final form 18 March 2005
Water availability on the continents is important for human health, economic activity, ecosystem function and geophysical processes. Because the saturation vapour pressure of water in air is highly sensitive to temperature, perturbations in the global water cycle are expected to accompany climate warming6. Regional patterns of warming-induced changes in surface hydroclimate are complex and less certain than those in temperature, however, with both regional increases and decreases expected in precipitation and runoff. Here we show that an ensemble of 12 climate models exhibits qualitative and statistically significant skill in simulating observed regional patterns of twentieth-century multidecadal changes in streamflow. These models project 10–40% increases in runoff in eastern equatorial Africa, the La Plata basin and high-latitude North America and Eurasia, and 10–30% decreases in runoff in southern Africa, southern Europe, the Middle East and mid-latitude western North America by the year 2050. Such changes in sustainable water availability would have considerable regional-scale consequences for economies as well as ecosystems.
for climate research developed at the Geophysical Fluid Dynamics Laboratory (GFDL) are presented. The atmosphere model, known as AM2, includes a new gridpoint dynamical core, a prognostic cloud scheme, and a multispecies aerosol climatology, as well as components from previous models used at GFDL. The land model, known as LM2, includes soil sensible and latent heat storage, groundwater storage, and stomatal resistance. The performance of the coupled model AM2–LM2 is evaluated with a series of prescribed sea surface temperature (SST) simulations. Particular focus is given to the model's climatology and the characteristics of interannual variability related to E1 Niño– Southern Oscillation (ENSO).
One AM2–LM2 integration was performed according to the prescriptions of the second Atmospheric Model Intercomparison Project (AMIP II) and data were submitted to the Program for Climate Model Diagnosis and Intercomparison (PCMDI). Particular strengths of AM2–LM2, as judged by comparison to other models participating in AMIP II, include its circulation and distributions of precipitation. Prominent problems of AM2– LM2 include a cold bias to surface and tropospheric temperatures, weak tropical cyclone activity, and weak tropical intraseasonal activity associated with the Madden–Julian oscillation.
An ensemble of 10 AM2–LM2 integrations with observed SSTs for the second half of the twentieth century permits a statistically reliable assessment of the model's response to ENSO. In general, AM2–LM2 produces a realistic simulation of the anomalies in tropical precipitation and extratropical circulation that are associated with ENSO.
Developments in analysis and modeling of continental water and energy balances are hindered by the limited availability and quality of observational data. The lack of information on error characteristics of basin water supply is an especially serious limitation. Here we describe the development and testing of methods for quantifying several errors in basin mean precipitation, both in the long-term mean and in the monthly and annual anomalies. To quantify errors in long-term mean, two error indices are developed and tested with positive results. The first provides an estimate of the variance of the spatial sampling error of long-term basin mean precipitation obtained from a gauge network, in the absence of orographic effects; this estimate is obtained by use only of the gauge records. The second gives a simple estimate of the basin mean orographic bias as a function of the topographic structure of the basin and the locations of gauges therein. Neither index requires restrictive statistical assumptions (such as spatial homogeneity) about the precipitation process. Adjustments of precipitation for gauge bias and estimates of the adjustment errors are made by applying results of a previous study. Additionally, standard correlation-based methods are applied for the quantification of spatial sampling errors in the estimation of monthly and annual values of basin mean precipitation. These methods also perform well, as indicated by network subsampling tests in densely gauged basins. The methods are developed and applied with data for 175 large (median area of 51,000 km2 ) river basins of the world for which contemporaneous, continuous (missing fewer than 2% of data values), long-term (median record length of 54 years) river discharge records are also available. Spatial coverage of the resulting river basin data set is greatest in the middle latitudes, though many basins are located in the tropics and the high latitudes, and the data set spans the major climatic and vegetation zones of the world. This new data set can be applied in diagnostic and theoretical studies of water balance of large basins and in the evaluation of performance of global models of land water balance.
Radiative effects of anthropogenic changes in atmospheric composition are expected to cause climate changes, in particular an intensification of the global water cycle with a consequent increase in flood risk. But the detection of anthropogenically forced changes in flooding is difficult because of the substantial natural variability; the dependence of streamflow trends on flow regime further complicates the issue. Here we investigate the changes in risk of great floods- that is, floods with discharges exceeding 100-year levels from basins larger than 200,000 km2- using both streamflow measurements and numerical simulations of the anthropogenic climate change associated with greenhouse gases and direct radiative effects of sulphate aerosols. We find that the frequency of great floods increased substantially during the twentieth century. The recent emergence of a statistically significant positive trend in risk of great floods is consistent with results from the climate model, and the model suggests that the trend will continue.
Shmakin, A B., P C D Milly, and Krista A Dunne, 2002: Global modeling of land water and energy balances. Part III: Interannual variability. Journal of Hydrometeorology, 3(3), 311-321. Abstract PDF
The Land Dynamics (LaD) model is tested by comparison with observations of interannual variations in discharge from 44 large river basins for which relatively accurate time series of monthly precipitation (a primary model input) have recently been computed. When results are pooled across all basins, the model explains 67% of the interannual variance of annual runoff ratio anomalies (i.e., anomalies of annual discharge volume, normalized by long-term mean precipitation volume). The new estimates of basin precipitation appear to offer an improvement over those from a state-of-the-art analysis of global precipitation (the Climate Prediction Center Merged Analysis of Precipitation, CMAP), judging from comparisons of parallel model runs and of analyses of precipitation-discharge correlations. When the new precipitation estimates are used, the performance of the LaD model is comparable to, but not significantly better than, that of a simple, semiempirical water-balance relation that uses only annual totals of surface net radiation and precipitation. This implies that the LaD simulations of interannual runoff variability do not benefit substantially from information on geographical variability of land parameters or seasonal structure of interannual variability of precipitation.
The aforementioned analyses necessitated the development of a method for downscaling of long-term monthly precipitation data to the relatively short timescales necessary for running the model. The method merges the long-term data with a reference dataset of 1-yr duration, having high temporal resolution. The success of the method, for the model and data considered here, was demonstrated in a series of model-model comparisons and in the comparisons of modeled and observed interannual variations of basin discharge.
Milly, P C., and Krista A Dunne, 2001: Trends in evaporation and surface cooling in the Mississippi River basin. Geophysical Research Letters, 28(7), 1219-1222. Abstract PDF
A synthesis of available data for the Mississippi River basin (area 3x106 km2) reveals an upward trend in evaporation during recent decades, driven primarily by increases in precipitation and secondarily by human water use. A cloud-related decrease in surface net radiation appears to have accompanied the precipitation trend. Resultant evaporative and radiative cooling of the land and lower atmosphere quantitatively explains downward trends in observed pan evaporation. These cooling tendencies also reconcile the observed regional atmospheric cooling with the anticipated regional "greenhouse warming." If recent high levels of precipitation (which correlate with the North Atlantic Oscillation) are mainly caused by an internal climatic fluctuation, an eventual return to normal precipitation could reveal heretofore-unrealized warming in the basin. If, instead, they are caused by some unidentified forcing that will continue to grow in the future, then continued intensification of water cycling and suppression of warming in the basin could result.
Milly, P C., and Krista A Dunne, 1998: Non-detectability of 20th-Century trends in river discharge from large basins: Observational and model-based results In Ninth Symposium on Global Change Studies and Namias Symposium on the Status and Prospects for Climate Prediction, Boston, MA, American Meteorological Society, 162-163.
Dunne, Krista A., and C J Willmott, 1996: Global distribution of plant-extractable water capacity of soil. International Journal of Climatology, 16, 841-859. Abstract
Plant-extractable water capacity of soil is the amount of water that can be extracted from the soil to fulfill evapotranspiration demands. It is often assumed to be spatially invariant in large-scale computations of the soil-water balance. Empirical evidence, however, suggests that this assumption is incorrect. In this paper, we estimate the global distribution of the plant-extractable water capacity of soil.
A representative soil profile, characterized by horizon (layer) particle size data and thickness, was created for each soil unit mapped by FAO (Food and Agriculture Organization of the United Nations)/Unesco. Soil organic matter was estimated empirically from climate data. Plant rooting depths and ground coverages were obtained from a vegetation characteristic data set. At each 0.5° x 0.5° grid cell where vegetation is present, unit available water capacity (cm water per cm soil) was estimated from the sand, clay, and organic content of each profile horizon, and integrated over horizon thickness. Summation of the integrated values over the lesser of profile depth and root depth produced an estimate of the plant-extractable water capacity of soil.
The global average of the estimated plant-extractable water capacities of soil is 8.6 cm (Greenland, Antarctica and bare soil areas excluded). Estimates are less than 5, 10 and 15 cm--over approximately 30, 60, and 89 per cent of the area, respectively. Estimates reflect the combined effects of soil texture, soil organic content, and plant root depth or profile depth. The most influential and uncertain parameter is the depth over which the plant-extractable water capacity of soil is computed., which is usually limited by root depth. Soil texture exerts a lesser, but still substantial, influence. Organic content, except where concentrations are very high, has relatively little effect.
Milly, P C., and Krista A Dunne, 1994: Sensitivity of the global water cycle to the water-holding capacity of land. Journal of Climate, 7(4), 506-526. Abstract PDF
The sensitivity of the global water cycle to the water-holding capacity of the plant-root zone of continental soils is estimated by simulations using a mathematical model of the general circulation of the atmosphere, with prescribed ocean surface temperatures and prescribed cloud. With an increase of the globally constant storage capacity, evaporation from the continents rises and runoff falls, because a high storage capacity enhances the ability of the soil to store water from periods of excess for later evaporation during periods of shortage. In addition to this direct effect, atmospheric feedbacks associated with the resulting higher precipitation and lower potential evaporation drive further changes in evaporation and runoff. Most of the changes in evaporation and runoff occur in the tropics and in the northern middle- latitude rain belts. Global evaporation from land increases by about 7 cm for each doubling of storage capacity in the range from less than 1 cm to almost 60 cm. Sensitivity is negligiblefor capacity above 60 cm. In the tropics and in the extratropics, the increased continental evaporation is split, in approximately equal parts, between increased continental precipitation and decreased convergence of atmospheric water vapor from ocean to land. In the tropics, this partitioning is strongly affected by induced circulation changes, which are themselves forced by changes in latent heating. The increased availability of water at the continental surface leads to an intensification of the Hadley circulation and a weakening of the monsoonal circulations. In the northern middle and high latitudes, the increased continental evaporation moistens the atmosphere. This change in humidity of the atmosphere is greater above the continents than above the oceans, and the resulting reduction in the sea-land humidity gradient causes a decreased onshore transport of water vapor by transient eddies. Results established here may have implications for certain problems in global hydrology and climate dynamics, including the effects of water resource development on global precipitation, climatic control of plant rooting characteristics, climatic effects of tropical deforestation, and climate-model errors induced by errors in land-surface hydrologic parameterizations.
