Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 46443.
Published online 2017 April 18. doi:  10.1038/srep46443
PMCID: PMC5394461

Simulated climate effects of desert irrigation geoengineering


Geoengineering, the deliberate large-scale manipulation of earth’s energy balance to counteract global warming, is an attractive proposition for sparsely populated deserts. We use the BNU and UVic Earth system models to simulate the effects of irrigating deserts under the RCP8.5 scenario. Previous studies focused on increasing desert albedo to reduce global warming; in contrast we examine how extending afforestation and ecological projects, that successfully improve regional environments, fair for geoengineering purposes. As expected desert irrigation allows vegetation to grow, with bare soil or grass gradually becoming shrub or tree covered, with increases in terrestrial carbon storage of 90.3 Pg C (UVic-ESCM) – 143.9 Pg C (BNU-ESM). Irrigating global deserts makes the land surface temperature decrease by 0.48 °C and land precipitation increase by 100 mm yr−1. In the irrigated areas, BNU-ESM simulates significant cooling of up to 4.2 °C owing to the increases in low cloud and latent heat which counteract the warming effect due to decreased surface albedo. Large volumes of water would be required to maintain global desert irrigation, equivalent 10 mm/year of global sea level (BNU-ESM) compensate for evapotranspiration losses. Differences in climate responses between the deserts prompt research into tailored albedo-irrigation schemes.

Desertification over the past several decades has provided a significant cooling effect on the earth’s surface, equivalent to about 20% of the warming due to global anthropogenic CO2 emissions over the same period1. This strong climatic effect promoted suggestions for increasing desert albedo e.g., by covering it with reflective material, as a potential geoengineering scheme2. However, simulations suggest that increasing desert albedo or expanding deserts3 would cause significant changes in global-scale circulation, the hydrological cycle, regional-scale precipitation patterns, and a large reduction in the intensity of the Indian and African monsoons4. In addition, desert expansion would also lead to the loss of productive land and the modification of natural arid ecosystems5.

Afforestation, the conversion of croplands or marginal lands into forests, is in some locations such as China, a cost-effective way to sequester CO2 from the atmosphere and an environmentally friendly way to ameliorate sand storms and mitigate desert encroachment6. The United Nations considers afforestation to be a key governmental climate-change mitigation strategy7. Afforestation not only contributes to increased carbon storage8 but also alters local albedo and turbulent energy fluxes, and enhances evapotranspiration that feedbacks on local and regional climate9,10. An expansion of relatively dark forest area increases the absorption of solar energy and increases surface temperature, particularly in regions where the land surface is unable to compensate with latent heat flux due to water limitation11. Afforestation can therefore result in net climate warming, particularly at high latitudes7. However, in the tropics, afforestation tends to have a net cooling effect because increased evaporative cooling overwhelms the heating effect of lowered surface albedo12.

Large increases in plant cover at geoengineering scales would require the addition of water to arid regions where they cannot grow at present. The climatic impact of irrigation in agriculture and/or desert regions has been variously assessed previously8,9,13,14,15,16. Agricultural irrigation can produce significant regional cooling with strong seasonal variability13,14 and regional increases in precipitation14, but has a negligible effect on global average near-surface temperatures15,16. There have been few desert irrigation or afforestation simulation studies8,9. Ornstein et al.9 used GISS ModelE General Circulation Model to simulate Saharan and Australian deserts irrigation with prescribed vegetation and phenology. Keller et al.8 used the University of Victoria Earth System Climate Model (UVic-ESCM) to study climate effects and carbon uptakes of afforestation geoengineering in the deserts of North Africa and Australia.

In this study, we use the BNU-ESM and the intermediate complexity UVic-ESCM Earth System Models to conduct idealized deserts irrigation geoengineering simulations during the 21st century. We conducted five sets of simulation experiments: no irrigation as a baseline (GE_none), irrigation of global deserts (GE_Globe), irrigation of Northwest China deserts (GE_China), irrigation of Australian deserts (GE_Australia), and irrigation of North African deserts (GE_Africa) (Supplementary Table S1). In each scenario we simulate irrigation by forcing soil moisture to be constant at 360 kg m−2 for UVic (following Keller et al.8) and by daily restoring the moisture content in the 3.433 m deep BNU-ESM soil column to be 800 mm thick, so as to avoid excessive runoff into ocean from the irrigated regions. All simulations are under the Intergovernmental Panel on Climate Change (IPCC) representative concentration pathway, RCP8.517 CO2 concentration scenario, which can only take the biogeophysical effects from irrigation, but exclude biogeochemical responses. We perform these simulations to illustrate the likely response of carbon reservoirs and climate to irrigating global and individual deserts geoengineering during the 21st century. The UVic-ESCM uses a reduced complexity atmosphere model with no parametrization for clouds and their feedbacks, nor atmospheric transport dynamic processes18. Therefore, the UVic-simulated climate effect of irrigation might be less reliable than that of BNU-ESM, and for results from UVic-ESCM simulations, we focus on the effect of irrigation on carbon reservoirs and land surface cover. A detailed description of the models used and simulation experiments is provided in the Method section.


Response of terrestrial carbon cycle to irrigation

In the each of the irrigated regions, desert irrigation lead to decreases in the area of bare soil and C4 grasses (Supplementary Figs S1 and S2) which are less adapted to wetter environments19. At the same time, the fractional area of C3 grasses increases initially, but then declines as shrubs take over (Supplementary Fig. S2). Eventually irrigation increases the fractional area of broadleaf trees (Supplementary Figs S1 and S2). By the end of simulation, the modeled predominant plant functional type in the irrigated areas is broadleaf tress for BNU-ESM and shrubs for UVic-ESCM, owing to the different vegetation dynamics used: the Lund–Potsdam–Jena (LPJ) model20 for BNU-ESM; and the Lotka-Volterra approach21 in UVic-ESCM. However, with BNU-ESM, irrigation could not make vegetation grow in Australia, likely due to simulated limitation in soil organic matter content. These changes in vegetation cover have a remarkable effect on terrestrial carbon cycle.

Increased shrub and tree cover in the irrigated areas increases both vegetation and soil terrestrial carbon reservoirs. By the end of simulation, BNU-ESM suggests increases in global vegetation and soil carbon caused by global desert irrigation of 73.1 and 70.8 Pg C respectively (Supplementary Fig. S3a,c). For UVic-ESCM, corresponding increases are 45.4 and 44.9 Pg C (Supplementary Fig. S3b,d). The greatest simulated increase is in vegetation carbon, reflecting the delay between production of additional living biomass and transfer to the litter and soil pools. GE_Globe increases total terrestrial carbon by 143.9 Pg C (7.6%) in BNU-ESM and 90.3 Pg C (4%) in UVic-ESCM compared with GE_none (Fig. 1a,b), and the increase in terrestrial carbon under GE_China, GE_Australia and GE_Africa are 7.5% (6.8 Pg C), 30.3% (27.3 Pg C) and 62.7% (56.6 Pg C) of GE_Globe using UVic-ESCM. These models predict increases in global net primary productivity (NPP), soil respiration (RES) and net ecosystem productivity (NEP, NPP minus RES) of 10.9, 9.2 and 1.7 Pg C yr−1 individually for BNU-ESM and 8.1, 7 and 1.1 Pg C yr−1 individually for UVic-ESCM with GE_Globe compared with GE_none averaged over the whole simulations. The increase of NPP in the irrigated area is due to vegetation growth, however, the increase of RES is affected by the increase of soil moisture, and the accumulation of litter and organic matter in the irrigated area. During 2020 to 2030, the average NEP anomaly under GE_Globe scenarios compared with GE_none is about 3.3 Pg C yr−1 for BNU-ESM and 6.3 Pg C yr−1 for UVic-ESCM.

Figure 1
Simulated temporal evolution of carbon changes anomalies, relative to GE_none, due to irrigation desert geoengineering.

In the North African irrigated areas, the simulated land carbon rises to more than 7.1 and 10.2 kg C m−2, and NEP (NPP minus RES) to 0.065 and 0.095 kg C m−2 yr−1 (Fig. 2a,b,g,h) during 2071 to 2100 for UVic-ESCM and BNU-ESM, individually, because of the different vegetation dynamics between the two models we used (Supplementary Figs S1c and S2c). However, in Australian irrigated areas, there are no increases in land carbon and NEP using BNU-ESM (Fig. 2a,g) as irrigation fails to make vegetation grow (Supplementary Fig. S1c). In the Northwest China irrigated areas UVic-ESCM simulated NEP anomaly is about −0.031 kg C m−2 yr−1 over the same period ((Fig. 2h) as the carbon uptake from shrubs is offset by decreases in C3 grasses (Supplementary Fig. S2b).

Figure 2
Maps showing model-simulated land carbon changes for GE_Globe compared with GE_none for 2071–2100.

Response of surface albedo to irrigation

The models’ terrestrial carbon cycle modules determine the exchange of CO2 between the land and the atmosphere. They are coupled to the physical climate through changes in vegetation growth (Supplementary Fig. S1), vegetation distribution and leaf area index, which affect the surface albedo, the evapotranspiration flux and so on. Relative to the GE_none simulation, BNU-ESM simulated global mean surface albedo under GE_Globe and GE_Africa abruptly decrease in the first 7 years (Fig. 3a) as the desert is covered with vegetation that has a lower albedo, particularly in the North African irrigated areas (Fig. 4a), then albedo continues to gradually decrease with inter-annual variability until 2053 (Fig. 3a). In western North America increased snow cover (Supplementary Fig. S4) leads to an albedo increase of over 0.1 under global desert irrigation scenario compared with the GE_none during 2071–2100 (Fig. 4a). The snow increases are due to greater moisture supply from the northern Pacific driven by sea level pressure rises there (Supplementary Figs S5 and S6). Subsequent positive albedo feedback cools the surface locally. In contrast, under GE_China simulation albedo changes are small because the bare soil fraction in Northwest China is much less before irrigation starts than in North Africa (Supplementary Fig. 1b,d; Table 1). Under GE_Australia simulation, no increased vegetation occurs, however, the increase in the soil moisture leads to the small decrease of albedo locally (Fig. 4a). During 2071 to 2100, the average surface albedo changes in irrigated areas and non-irrigated area on land under GE_Globe scenarios compared with GE_none are −0.152 and −0.004 (Supplementary Fig. S7a; Table 1), with the local responses to irrigation dominating changes in remote regions.

Figure 3
BNU-ESM simulated temporal evolution of climate changes anomalies, relative to GE_none, due to irrigation desert geoengineering.
Figure 4
Maps showing BNU-ESM simulated climate impacts for GE_Globe compared with GE_none for 2071–2100.
Table 1
Desert irrigation geoengineering areas, water required and induced changes in carbon and climate using BNU-ESM during 2071–2100.

Response of surface temperature to irrigation

In our simulations, irrigation affects surface temperature mainly through its effect on surface albedo, latent heat and cloud effect. Decreases in surface albedo associated with vegetation growth due to irrigation act to warm the surface, but the increases in evapotranspiration and low cloud cover due to irrigation and vegetation growth also cool the surface. Compared to CO2 induced warming, irrigation has a tiny effect on surface temperatures. Simulated global mean land surface temperature with global desert irrigation for the period 2020 and 2100 cools by 0.31 °C (Fig. 3b), which is, of course, much less than the 4.8 °C rise in global mean land temperature caused by increasing atmospheric CO2 concentrations under the GE_none scenario driven by RCP8.5 concentrations. However, the global mean surface temperature caused by global desert irrigation warms by 0.08 °C due to a rise of 0.24 °C in mean ocean temperature (Supplementary Fig. S8b,d). Under the GE_Globe and GE_Africa scenarios, latent heat flux abruptly increases in the first 5 years due to the increase of evapotranspiration in irrigated areas (Supplementary Fig. S7d), producing a net global land cooling of 1 °C (Fig. 3b).

Averaged over the 2071 to 2100 period, GE_Globe and GE_Africa produce global land average temperature anomalies of about −0.48 and −0.46 °C respectively relative to GE_none (Fig. 3b). However, global mean land temperature anomalies are near zero (Fig. 3b and Supplementary Table S2) under the GE_China and GE_Australia for the same period. The cooling effects are concentrated in the irrigated and adjacent areas (Fig. 4b), where decrease of up to −4.2 °C (Supplementary Fig. S7b). The cooling effect is due to TOA cloud net radiative forcing (longwave and shortwave) caused by the increases in low cloud cover fraction (Fig. 4e and Supplementary Fig. S9) and latent heat flux due to increases in evapotranspiration (Fig. 4f), which overwhelm the warming from decreased surface albedo (Fig. 4a). Modeled surface air temperatures cool over western North America by up to 2 °C due to increased snow cover (Fig. 4a and Supplementary Fig. S4). Over ocean, the temperature changes are consistent with cloud fraction, TOA cloud net radiative forcing and sea surface albedo changes. For example the positive temperature anomaly over the north Pacific is linked to positive TOA net cloud radiative forcing and decreased low cloud fraction (Fig. 4e and Supplementary Fig. S9), and the small positive temperature anomaly over the Southern Ocean to albedo changes probably caused by sea ice loss and changes in Antarctic Circumpolar Current (ACC) (Fig. 4a and Supplementary Fig. S10a,b).

Response of precipitation and evapotranspiration to irrigation

As a result of increased evapotranspiration associated with vegetation growth, desert irrigation increases global average precipitation, especially in and around the irrigated areas. Compared with GE_none in 2071 to 2100, increases in global mean precipitation are about 30.4, 2.3, 1.0 and 26.5 mm yr−1 for GE_Globe, GE_China, GE_Australia and GE_Africa respectively (Fig. 3c). Over land areas precipitation increases by up to 99.6 mm yr−1 (Supplementary Fig. S7c), but by only 1.7 mm yr−1 over the oceans during the 2071 to 2100 period (Supplementary Fig. S8h). Positive precipitation anomalies equator-ward and westward of the irrigation in North Africa and Australia, and westward of the irrigation in Northwest China (Fig. 4c) are driven by climatological winds, geostrophic adjustment and diffusion carrying moisture from the irrigated regions (Supplementary Fig. S5). Increases in tropical easterlies is consistent with enhancement of the Walker circulation due to increased sea-level pressure differences between eastern Pacific Ocean (130–80° W, 10° S–10° N) and Indonesia (110–160° E, 10° S–10° N) (Supplementary Figs S6 and S11), and in-line with model simulations under greenhouse gas forcing22 and recent observations23.

Evapotranspiration increases in irrigated areas because of increases in both vegetation cover and soil moisture (Fig. 4f). Global mean precipitation minus evapotranspiration (P-E) (Fig. 3d) decreases for all irrigation cases (Table 1) relative to GE_none. Averaged over years 2071–2100 regional irrigation in North Africa, Australia Northwest China causes decreases in P-E of up to −276.2 mm yr−1 (Fig. 4d, Supplementary Fig. S7d). For the same period, P-E over land increases by 4.5 mm yr−1 while it decreases by 1.8 mm yr−1 over oceans relative to GE_none (Supplementary Fig. S7d).


We may validate our simulations by comparison with some previous published work on the regional climatic impacts of irrigation13,14,15,16; and with irrigation impacts on carbon storage6,24,25. The models we use have also been thoroughly validated against observations (see Methods), and via model intercomparisons of future greenhouse and geoengineered climates26,27,28. Kueppers et al.13 used a regional climate model (RegCM3) to deduce that a local annual mean cooling of 1.6 °C (1981–2000) was caused by irrigated agriculture in California, which may be compared with the BNU-ESM simulated local cooling of 4.2 °C in desert irrigated areas. This cooling is more comparable with atmospheric general circulation model (GISS ModelE) simulations of 20th century irrigation14, which suggested peak regional cooling of about 3 °C in the northwestern portions of the Indian subcontinent. While global agricultural irrigation has altered climate significantly in some regions during the late 20th century, there has been a negligible effect on average near-surface temperatures14,15,16, in agreement with our very small (0.08 °C) mean temperature difference simulated by desert irrigation over 2071–2100. Our simulated changes in carbon stocks are comparable with field measurements in deserts. The UVic-ESCM and BNU-ESM simulated average NEP (NPP minus RES) in the North African irrigated areas rises to 0.065 and 0.095 kg C m−2 yr−1 (Fig. 2a,b,g,h), which are both consistent with net uptake of carbon ranging from 0.033 to 0.127 kg C m−2 yr−1 in field experiments in the Tengger Desert by Yang et al.6, northern Chihuahuan Desert by Petrie et al.24 and Mojave Desert by Jasoni et al.25.

There are only a few existing studies on climate effect of desert irrigation or afforestation8,9. Ornstein et al.9, using the GISS ModelE coupled climate model, simulated the effect of afforestation by replacing the deserts of Sahara and Australia with evergreen tropical rain forest and changing the corresponding desert soil to a composition typical for rain forest. They found that afforestation has very localized effects with temperature decreases of 2 to 3 °C and precipitation increases of 2 to 3 mm day−1 over North African and Australian deserts. In contrast, Keller et al.8, using the UVic-ESCM, simulated the climate effect of afforestation by forcing the soil moisture to be 360 kg m−2 in the deserts of North Africa and the Australia. They found that the albedo change caused by simulated afforestation resulted in a global mean temperature warming of 0.1 °C and an 18 mm yr−1 increase in precipitation in 2100 relative to no climate engineering. We, using the BNU-ESM, find that irrigating global and North African deserts produces significant cooling in the irrigated and adjacent areas, where decreases of up to −4.2 °C occur by 2071–2100, mainly as a result of increases in low cloud cover and latent heat flux, which counteracts the warming effect associated with reduced surface albedo. However, there is a warming of 0.08 °C in global mean surface temperature, which the cooling of 0.31 °C in global mean land surface temperature is overwhelmed by the warming of 0.24 °C in global mean ocean surface temperature, from irrigation of global deserts during the simulations. Global average precipitation by 2071–2100 under GE_Globe increases by 30.4 mm yr−1 (Fig. 3j) compared with GE_none. The increases in precipitation under GE_China, GE_Australia and GE_Africa are 7.6% (2.3 mm yr−1), 3.3% (1.0 mm yr−1) and 87.2% (26.5 mm yr−1) of GE_Globe. Puma et al.14 and Sacks et al.15 estimated that agriculture irrigation increased global land average precipitation by 9.5–4.3 mm yr−1, during the late 20th century, both of which are much lower than the increase in precipitation under GE_Globe, but this may be explained by the different irrigation rates16. Evapotranspiration is high in the irrigated areas because of increases in both vegetation cover and soil moisture, hence average precipitation minus evapotranspiration (P-E) anomalies in the irrigated area are negative (Fig. 4d and Table 1). In contrast, Ornstein et al.9 simulations produced considerable extra rainfall over irrigated new forests, almost offsetting the 2 to 3 mm day−1 increases in evaporation. GISS ModelE as used by Ornstein et al.9 included a rather simple vegetation module compared with that in the BNU-ESM29.

In our simulations using the BNU-ESM and UVic-ESCM model we rely on natural vegetation succession, though actually a more likely scenario would be to irrigate and plant. In practice irrigation would be done with a blend of drought-resistant trees and shrub species which use water efficiently and can tolerate the arid environment30. In our simulations, global desert irrigation increased the terrestrial carbon uptake by 6.4 Pg C yr−1 by the fifth year of irrigation for BNU-ESM and by 9.4 Pg C yr−1 in only the third year of irrigation for UVic-ESCM, which has a faster vegetation succession than the BNU-ESM. GE_China, GE_Australia and GE_Africa peak sequestration rates are 34% (11%), 14% (30%) and 72% (60%) of GE_Globe using BNU-ESM (UVic-ESCM) (Fig. 1g,h). Keller et al.8 finds a similar global peak removal rate, while Ornstein et al.9 estimated the terrestrial carbon reservoir flux from desert afforestation at about 8 Pg C yr−1 in the first 10 years. These rates are more than the carbon emission gap between Paris declaration of INDCs (Intended Nationally Determined Contributions) scenarios and emissions needed to meet a 2 °C peak warming31, which are about 2 Pg C yr−1.

Desert irrigation geoengineering needs huge amounts of fresh water. A first year rise in soil water content to daily restoring the moisture content in the 3.433 m deep BNU-ESM soil column to be 800 mm thick across the global deserts based on present-day soil water contents is equivalent to 7.6 mm of global mean sea level (2722 billion tons) (Supplementary Fig. S12f and Table 1). Ornstein et al.9 estimated an initial requirement of 4900 billion tons fresh water per year for Sahara desert irrigation (which they prescribe as requiring 500 mm yr−1 of precipitation), though in their simulation, this becomes self-sustaining due to increased rainfall over newly forested deserts. The required initial amount of fresh water estimated by Ornstein et al.9 is about twice our initial irrigation requirement, due to the different methods of irrigation or afforestation. More significantly in our simulations, global desert irrigation would require at least 3766 billion tons (10 mm of global sea level) per year to supplement net water lost by evapotranspiration (Fig. 4d and Supplementary Fig. S12h).

Keller et al.8 estimated that the higher soil moisture and precipitation caused by irrigation would raise the global mean sea level by about 13 cm by 2100 unless the water used for irrigation is desalinated seawater. However, economical irrigation limits water supply to maintain under-saturated soil, preventing runoff of excess water. The water required may be desalinated from sea water using reverse osmosis or permeable membrane desalination or extracted from deep aquifers9. This process would involve additional energy cost which will bring additional carbon emissions. Runoff and rainfall from desert irrigation would also regionally increase freshwater input to oceans thereby reducing coastal salinity. Lowered salinity ocean water in the upper ocean of the North Atlantic (Supplementary Fig. S13) under GE_Globe and GE_Africa may weaken Atlantic Meridional Overturning Circulation (AMOC)32 (Supplementary Fig. S10, Table S2).

Uncertainty in irrigation results appear quite large given the differences between our BNU-ESM results and those from the GISS ModelE, while UVic-ESCM seems reliable for carbon budgets, the lack of realistic atmospheric processes limits its usefulness in predicting regional climatic effects. The vegetation model exerts a large impact on the carbon budget. Furthermore, our simulations were driven by CO2 concentrations, and it only takes into account the biogeophysical effects from irrigation, but excludes biogeochemical responses. The increase of vegetation by irrigation or afforestation in coastal areas leads to increased sea-surface temperature owing to greater roughness and weaker winds over the adjacent coastal ocean, and then greater humidity, which is transported into the monsoon region by anomalous winds33. Furthermore, mineral dust aerosols affect atmosphere radiation budget by scattering and partly absorbing shortwave and longwave radiation34,35. Increasing soil moisture by irrigation and promoting vegetation growth means that less mineral dust aerosols will be mobilized and transported in the atmosphere36,37,38. Dust also has important fertilization effects in the ocean, and reductions may then lead to changes in ocean carbon budget. Last but not least, vegetation growth by irrigation or afforestation may induce more frequent occurrences of extreme rainfall events (flooding) over the coastal region and more frequent occurrence of heat waves and droughts over the semi-arid region39.

For these desert irrigation simulations, the changes in latent heat flux and low cloud are two dominant factors causing cooling (Fig. 4e,f). So we raise the idea that a combination of locally balancing artificial increases in albedo with irrigation driven increases in evapotranspiration may have much better prospects as a tool for desert climate engineering. This could be realized for example, by having both irrigated and highly reflective areas in the desert regions. An investigation of optimal desert geoengineering requires additional simulations in the future with more sophisticated models that contain high-resolution land cover data for reducing modeling uncertainties40.


Model description

The Beijing Normal University—Earth System Model (BNU-ESM)29 is a fully coupled Earth system model comprising elements representing the atmosphere (NCAR-CAM3.5), and land surface (BNU-CoLM3) at a T42 horizontal spectral resolution (approximately 2.81° × 2.81° transform grid), ocean (GFDL-MOM4p1) and sea-ice (LANL-CICE4.1) at a nominal latitude-longitude resolution of 1° (down to 1/3° within 10o of the equatorial tropics) with 360 longitudinal grids and 200 latitudinal grids; and a coupler (improved NCAR-CPL6). It has an interactive carbon cycle model in the land component (BNU-DGVM (C/N)) based on LPJ and an ecosystem-biogeochemical module in the ocean component (IBGC). It has participated in the Coupled Model Intercomparison Project Phase 5 (CMIP5)26 and Geoengineering Model Intercomparison Project (GeoMIP)27.

As a member of CMIP5, BNU-ESM is based on several widely evaluated climate model components41,42,43,44. Ji et al.29 have validated the BNU-ESM’s performance in terms of the mean model state and the internal variability by comparisons between CMIP5 piControl and historical simulations and the observed datasets. BNU-ESM can in general, simulate observed features of the earth climate system as well as any other Earth System Model, including the climatological annual cycle of surface-air temperature and precipitation, annual cycle of cloud fraction and forcing, terrestrial gross primary productions (GPP) and soil organic carbon stocks (see Figs 4–9, refs 23 and 24 in Ji et al.29). BNU-ESM has proven to be a useful modeling tool and is being actively used by many researchers in prognostic simulations for both anthropogenic and geoengineering forcing scenarios45,46,47,48,49,50,51.

The University of Victoria Earth System Climate Model (UVic-ESCM) is a coupled climate-carbon cycle model with a horizontal resolution of 3.6° longitude × 1.8° latitude. It includes a vertically integrated energy–moisture balance atmospheric model, a three-dimensional general circulation model of the ocean and a terrestrial and ocean carbon cycle model18. The terrestrial component is a dynamic global vegetation model, including a land surface scheme52 based on the Hadley Center model TRIFFID (Top-down Representation of Interactive Foliage and Flora Including Dynamics)53. UVic-ESCM has been widely used in a variety of research topics such as climate feedbacks of land cover change52,54,55; climate change and the carbon cycle8,56,57.

UVic-ESCM has also been used for multi-century climate projections in support of the IPCC Fifth Assessment Report58. Eby et al.28 have compared the climate and carbon cycle responses of Earth System Models of Intermediate Complexity (EMICs), including UVic-ESCM, over the historical period against observational estimates. Overall 20th century trends in carbon uptake are reasonably well simulated when compared to observed trends (see Fig. 5 in Eby et al.28). Weaver et al.18 showed that present-day atmospheric surface air temperature simulated by UVic-ESCM shows good agreement with annually-averaged NCEP reanalysis climatology59, especially on the zonal mean (see Fig. 10 in Weaver et al.18). However, the UVic-ESCM uses a reduced complexity atmosphere model with no parametrization for clouds and their feedbacks, nor atmospheric transport dynamic processes18. Therefore, the UVic-simulated climate effect of irrigation may be less reliable than that of BNU-ESM, and therefore we focus on the effect of irrigation on carbon reservoirs and land surface cover results from UVic-ESCM.

Simulation experiments

All simulations are forced by the IPCC high-end CO2 concentration scenario, RCP8.517. We use the BNU-ESM and UVic-ESCM to conduct idealized desert irrigation geoengineering simulations from year 2020 to 2100, which include simulations of irrigating deserts in Northwest China, Australia and North Africa, individually, and all together (Fig. 2a and Supplementary Table S1). In these simulations, we consider only CO2 radiative forcing. We neglect radiative forcing from aerosols and other greenhouse gases, which roughly cancel each other, will have little effect on our simulated impact of irrigation on carbon reservoir and climate change. In each scenario we simulate irrigation by forcing soil moisture to be constant at 360 kg m−2 for UVic (following Keller et al.8) and by daily restoring the moisture content in the 3.433 m deep BNU-ESM soil column to be 800 mm thick, so as to avoid excessive runoff into ocean from the irrigated regions.

Statistical tests

We use the Student’s t-test since the model output are quite Normal applied to each grid point on the maps in the Figs 2,,44 and Supplementary Figs S4–6,S9 and S13. Hatched areas in these figures are regions where changes are not statistically significant at the 5% level. The Null hypotheses are “GE_Globe = GE_none” over the period of 2071 to 2100 (n = 30) for GE_Globe results, “GE_China = GE_none” for GE_China results, “GE_Australia = GE_none” for GE_Australia results and “GE_Africa = GE_none” for GE_Africa results. There are 8192 grid points of temperature, albedo, precipitation and precipitation minus evapotranspiration (P–E) annual mean BNU-ESM data, and 3527 land grid points for BNU-ESM and 2543 land grid points for UVic-ESCM of land carbon, net primary productivity (NPP), soil respiration (RES) and net ecosystem productivity (NEP, NPP minus RES) annual mean data. The returned probability for each grid is two-tailed. The Student’s t-test and plot of these maps in Figs 2,,44 and Supplementary Figs S4–6,S9 and S13 were produced using NCAR Command Language (NCL)60 version 6.1.2.

Additional Information

How to cite this article: Cheng, W. et al. Simulated climate effects of desert irrigation geoengineering. Sci. Rep. 7, 46443; doi: 10.1038/srep46443 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:


Research was funded by the National Basic Research Program of China grant number 2015CB953600.


The authors declare no competing financial interests.

Author Contributions W.C., J.C.M., L.C., D.J. and L.Z. conceived, interpreted the results and wrote the paper. W.C. conducted the analysis.


  • Rotenberg E. & Yakir D. Contribution of Semi-Arid Forests to the Climate System. Science 327, 451–454 (2010). [PubMed]
  • Gaskill A. Summary of meeting with US DOE to discuss geoengineering options to prevent long-term climate change (2004).
  • Dirmeyer P. A. & Shukla J. The effect on regional and global climate of expansion of the world’s deserts. Q. J. R. Meteorol. Soc. 122, 451–482 (1996).
  • Irvine P. J., Ridgwell A. & Lunt D. J. Climatic effects of surface albedo geoengineering. J. Geophys. Res. 116, D24112 (2011).
  • Yang X., Jia Z. & Ci L. Assessing effects of afforestation projects in China. Nature 466, 315 (2010). [PubMed]
  • Yang H. et al. . Carbon sequestration capacity of shifting sand dune after establishing new vegetation in the Tengger Desert, northern China. Sci. Total Environ. 478, 1–11 (2014). [PubMed]
  • Arora V. K. & Montenegro A. Small temperature benefits provided by realistic afforestation efforts. Nat. Geosci. 4, 514–518 (2011).
  • Keller D. P., Feng E. Y. & Oschlies A. Potential climate engineering effectiveness and side effects during a high carbon dioxide-emission scenario. Nat. Commun. 5, 3304 (2014). [PMC free article] [PubMed]
  • Ornstein L., Aleinov I. & Rind D. Irrigated afforestation of the Sahara and Australian Outback to end global warming. Clim. Change 97, 409–437 (2009).
  • Peng S. et al. . Afforestation in China cools local land surface temperature. Proc. Natl. Acad. Sci. USA 111, 2915–2919 (2014). [PubMed]
  • Swann A. L. S., Fung I. Y. & Chiang J. C. H. Mid-latitude afforestation shifts general circulation and tropical precipitation. Proc. Natl. Acad. Sci. USA 109, 712–716 (2012). [PubMed]
  • Betts R. A. Climate science: Afforestation cools more or less. Nat. Geosci. 4, 504–505 (2011).
  • Kueppers L. M., Snyder M. A. & Sloan L. C. Irrigation cooling effect: Regional climate forcing by land-use change. Geophys. Res. Lett. 34, L03703 (2007).
  • Puma M. J. & Cook B. I. Effects of irrigation on global climate during the 20th century. J. Geophys. Res. 115, D16120 (2010).
  • Sacks W. J., Cook B. I., Buenning N., Levis S. & Helkowski J. H. Effects of global irrigation on the near-surface climate. Clim. Dyn. 33, 159–175 (2009).
  • Cook B. I., Puma M. J. & Krakauer N. Y. Irrigation induced surface cooling in the context of modern and increased greenhouse gas forcing. Clim. Dyn. 37, 1587–1600 (2011).
  • Meinshausen M. et al. . The RCP greenhouse gas concentrations and their extensions from 1765 to 2300. Clim. Change 109, 213–241 (2011).
  • Weaver A. J. et al. . The UVic Earth System Climate Model: Model description, climatology, and applications to past, present and future climates. Atmos.-Ocean 39, 361–428 (2001).
  • Ward J. K., Tissue D. T., Thomas R. B. & Strain B. R. Comparative responses of model C3 and C4 plants to drought in low and elevated CO2. Glob. Change Biol. 5, 857–867 (1999).
  • Sitch S. et al. . Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ dynamic global vegetation model. Glob. Change Biol. 9, 161–185 (2003).
  • Silvertown J. & Charlesworth D. Introduction to plant population ecology (second edition) 122–152 (Longman Scientic and Technical, 1987).
  • Bayr T., Dommenget D., Martin T. & Power S. B. The eastward shift of the Walker Circulation in response to global warming and its relationship to ENSO variability. Clim. Dyn. 43, 2747–2763 (2014).
  • L’Heureux M. L., Lee S. & Lyon B. Recent multidecadal strengthening of the Walker circulation across the tropical pacific. Nat. Clim. Chang. 3, 571–576 (2013).
  • Petrie M. D., Collins S. L., Swann A. M., Ford P. L. & Litvak M. E. Grassland to shrubland state transitions enhance carbon sequestration in the northern Chihuahuan Desert. Glob. Change Biol. 21, 1226–1235 (2015). [PubMed]
  • Jasoni R. L., Smith S. D. & Arnone III J. A. Net ecosystem CO2 exchange in Mojave Desert shrublands during the eighth year of exposure to elevated CO2. Glob. Change Biol. 11, 749–756 (2005).
  • Flato G. et al. . Evaluation of Climate Models In: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (ed. Stocker T. F. et al. .) 741–866 (Cambridge University Press, 2013).
  • Kravitz B. et al. . Climate model response from the Geoengineering Model Intercomparison Project (GeoMIP). J. Geophys. Res. 118, 8320–8332 (2013).
  • Eby M. et al. . Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity. Clim. Past 9, 1111–1140 (2013).
  • Ji D. et al. . Description and basic evaluation of Beijing Normal University Earth System Model (BNU-ESM) version 1. Geosci. Model Dev. 7, 2039–2064 (2014).
  • McDowell N. et al. . Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytol. 178, 719–739 (2008). [PubMed]
  • Conference of the Parties United Nations Framework Convention on Climate Change. Aggregate effect of intended nationally determined contributions: an update (2016).
  • Behrens E., Biastoch A. & Boning C. W. Spurious AMOC trends in global ocean sea-ice models related to subarctic freshwater forcing. Ocean Model. 69, 39–49 (2013).
  • Ma D., Notaro M., Liu Z., Chen G. & Liu Y. Simulated impacts of afforestation in East China monsoon region as modulated by ocean variability. Clim. Dyn. 41, 2439–2450 (2013).
  • Choobari O. A., Zawar-Reza P. & Sturman A. The global distribution of mineral dust and its impacts on the climate system: A review. Atmos. Res. 138, 152–165 (2014).
  • Spyrou C. et al. . Modeling the radiative effects of desert dust on weather and regional climate. Atmos. Chem. Phys. 13, 5489–5504 (2013).
  • Cowie S. M., Knippertz P. & Marsham J. H. Are vegetation-related roughness changes the cause of the recent decrease in dust emission from the Sahel? Geophys. Res. Lett. 40, 1868–1872 (2013). [PMC free article] [PubMed]
  • Li X. & Zhang H. Soil moisture effects on sand saltation and dust emission observed over the Horqin sandy land area in China. J. Meteorol. Res. 28, 444–452 (2014).
  • Pierre C. et al. . Impact of vegetation and soil moisture seasonal dynamics on dust emissions over the Sahel. J. Geophys. Res. 117, D06114 (2012).
  • Abiodun B. J., Salami A. T., Matthew O. J. & Odedokun S. Potential impacts of afforestation on climate change and extreme events in Nigeria. Clim. Dyn. 41, 277–293 (2013).
  • Di Vittorio A. V. et al. . From land use to land cover: restoring the afforestation signal in a coupled integrated assessment-earth system model and the implications for CMIP5 RCP simulations. Biogeosciences 11, 6435–6450 (2014).
  • Neale R. B. et al. . The Mean Climate of the Community Atmosphere Model (CAM4) in Forced SST and Fully Coupled Experiments. J. Clim. 26, 5150–5168 (2013).
  • Griffies S. M. Elements of MOM4p1, GFDL Ocean Group Technical Report No. 6 (2010).
  • Hunke E. C. & Lipscomb W. H. CICE: The Los Alamos sea ice model user’s manual, version 4.1’s_manual.pdf (2010).
  • Dai Y. et al. . The Common Land Model (CLM). Bull. Amer. Meteorol. Soc. 84, 1013–1023 (2003).
  • Wu R. G., Chen J. P. & Wen Z. P. Precipitation-surface temperature relationship in the IPCC CMIP5 Models. Adv. Atmos. Sci. 30, 766–778 (2013).
  • Bracegirdle T. J. et al. . Assessment of surface winds over the Atlantic, Indian, and Pacific Ocean sectors of the Southern Ocean in CMIP5 models: historical bias, forcing response, and state dependence. J. Geophys. Res. 118, 547–562 (2013).
  • Chen L., Yu Y. & Sun D. Z. Cloud and Water Vapor Feedbacks to the El Niño Warming: Are They Still Biased in CMIP5 Models? J. Clim. 26, 4947–4961 (2013).
  • Vial J., Dufresne J. L. & Bony S. On the interpretation of intermodal spread in CMIP5 climate sensitivity estimates. Clim. Dyn. 41, 3339–3362 (2013).
  • Wei T. et al. . Developed and developing world responsibilities for historical climate change and CO2 mitigation. Proc. Natl. Acad. Sci. USA 109, 12911–12915 (2012). [PubMed]
  • Tian D. et al. . Climate System Responses to a Common Emission Budget of Carbon Dioxide. J. Clim. 29, 2433–2442 (2016).
  • Lo Y. T. E., Charlton-Perez A. J., Lott F. C. & Highwood E. J. Detecting sulphate aerosol geoengineering with different methods. Sci. Rep. 6, 39169 (2016). [PMC free article] [PubMed]
  • Meissner K. J., Weaver A. J., Matthews H. D. & Cox P. M. The role of land surface dynamics in glacial inception: a study with the UVic Earth System Model. Clim. Dyn. 21, 515–537 (2003).
  • Cox P. M. Description of the “TRRIFFID” dynamic global vegetation model. Hadley centre technical report technical note 24 (2001).
  • Matthews H. D., Weaver A. J., Eby M. & Meissner K. J. Radiative forcing of climate by historical land cover change. Geophys. Res. Lett. 30, 1055 (2003).
  • Matthews H. D., Weaver A. J., Meissner K. J., Gillett N. P. & Eby M. Natural and anthropogenic climate change: incorporating historical land cover change, vegetation dynamics and the global carbon cycle. Clim. Dyn. 22, 461–479 (2004).
  • Schmittner A., Oschlies A., Matthews H. D. & Galbraith E. D. Future changes in climate, ocean circulation, ecosystems, and biogeochemical cycling simulated for a business-as-usual CO2 emission scenario until year 4000 AD. Glob. Biogeochem. Cycles 22, GB1013 (2008).
  • Wania R. et al. . Carbon-nitrogen feedbacks in the UVic-ESCM. Geosci. Model Dev. 5, 1137–1160 (2012).
  • Zickfeld K. et al. . Long-term climate change commitment and reversibility: an EMIC Intercomparison. J. Clim. 26, 5782–5809 (2013).
  • Kalnay E. et al. . The NCEP/NCAR 40-year reanalysis project. Bull. Amer. Meteorol. Soc. 77, 437–471 (1996).
  • Boulder, Colorado: UCAR/NCAR/CISL/TDD. The NCAR Command Language (Version 6.1.2) [Software]. (2015).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group