Agricultural and Forest Meteorology 223 (2016) 151–167
Contents lists available at ScienceDirect
Agricultural and Forest Meteorology
journal homepage: www.elsevier.com/locate/agrformet
Assessment and simulation of global terrestrial latent heat flux by
synthesis of CMIP5 climate models and surface eddy covariance observations
a,∗ a b a c
Yunjun Yao , Shunlin Liang , Xianglan Li , Shaomin Liu , Jiquan Chen ,
a a a a d a a
Xiaotong Zhang , Kun Jia , Bo Jiang , Xianhong Xie , Simon Munier , Meng Liu , Jian Yu
e f g h i
, Anders Lindroth , Andrej Varlagin , Antonio Raschi , Asko Noormets , Casimiro Pio ,
j,k l m,n o,p
Georg Wohlfahrt , Ge Sun , Jean-Christophe Domec , Leonardo Montagnani ,
q m,n r s t
Magnus Lund , Moors Eddy , Peter D. Blanken , Thomas Grünwald , Sebastian Wolf , u
Vincenzo Magliulo
a
State Key Laboratory of Remote Sensing Science, School of Geography, Beijing Normal University, Beijing, 100875, China
b
College of Global Change and Earth System Science, Beijing Normal University, Beijing, 100875, China
c
CGCEO/Geography, Michigan State University, East Lansing, MI 48824, USA
d
Laboratoire d’études en géophysique et océanographie spatiales, LEGOS/CNES/CN RS/IRD/UPS-UMR5566, Toulouse, France
e
GeoBiosphere Science Centre, Lund University, Sölvegatan 12, 223 62 Lund, Sweden
f
A.N.Severtsov Institute of Ecology and Evolution RAS 123103, Leninsky pr.33, Moscow, Russia
g
CNR-IBIMET, National Research Council, Via Caproni 8, 50145 Firenze, Italy
h
Dept. Forestry and Environmental Resources, North Carolina State University, 920 Main Campus Drive, Ste 300, Raleigh, NC 27695, USA
i
CESAM & Departamento de Ambiente e Ordenamento, Universidade de Aveiro, 3810-193 Aveiro, Portugal
j
Institute for Ecology, University of Innsbruck, Sternwartestr. 15, A-6020 Innsbruck, Austria
k
European Academy Bolzano, Drususallee 1, 39100 Bolzano, Italy
l
Eastern Forest Environmental Threat Assessment Center, Southern Research Station, United States Department of Agriculture Forest Services, Raleigh, NC
27606, USA
m
Climate Change and Adaptive Land and Water Management, Alterra Wageningen UR, PO Box 47, 6700 AA Wageningen, The Netherlands
n
VU University Amsterdam, Boelelaan 1085, 1081HV Amsterdam, The Netherlands
o
Faculty of Science and Technology, Free University of Bolzano, Piazza Universita 5, 39100 Bolzano, Italy
p
Forest Services, Autonomous Province of Bolzano, 39100 Bolzano, Italy
q
Department of Bioscience, Aarhus University, Frederiksborgvej 399, DK-4000 Roskilde, Denmark
r
Department of Geography, University of Colorado at Boulder, 260 UCB, Boulder, CO 80309-0260, USA
s
Technische Universität Dresden, Institute of Hydrology and Meteorology, Chair of Meteorology, D-01062 Dresden, Germany
t
ETH Zurich, Institute of Agricultural Sciences, Universitaetsstr. 2, 8092 Zurich, Switzerland
u
CNR-ISAFOM, National Research Council, Via Patacca 85, 80040 Ercolano (Napoli), Italy
a r t i c l e i n f o a b s t r a c t
Article history: The latent heat flux (LE) between the terrestrial biosphere and atmosphere is a major driver of the global
Received 19 August 2015
hydrological cycle. In this study, we evaluated LE simulations by 45 general circulation models (GCMs)
Received in revised form 24 March 2016
in the Coupled Model Intercomparison Project Phase 5 (CMIP5) by a comparison with eddy covariance
Accepted 29 March 2016
(EC) observations from 240 globally distributed sites from 2000 to 2009. In addition, we improved global
terrestrial LE estimates for different land cover types by synthesis of seven best CMIP5 models and EC
Keywords:
observations based on a Bayesian model averaging (BMA) method. The comparison results showed sub-
Global terrestrial LE
stantial differences in monthly LE among all GCMs. The model CESM1-CAM5 has the best performance with
CMIP5
GCMs the highest predictive skill and a Taylor skill score (S) from 0.51–0.75 for different land cover types. The
BMA cross-validation results illustrate that the BMA method has improved the accuracy of the CMIP5 GCM’s LE
2
Taylor skill score simulation with a decrease in the averaged root-mean-square error (RMSE) by more than 3 W/m when
compared to the simple model averaging (SMA) method and individual GCMs. We found an increas-
2 −1
ing trend in the BMA-based global terrestrial LE (slope of 0.018 W/m yr , p < 0.05) during the period
1970–2005. This variation may be attributed directly to the inter-annual variations in air temperature
(Ta), surface incident solar radiation (Rs) and precipitation (P). However, our study highlights a large
difference from previous studies in a continuous increasing trend after 1998, which may be caused by
∗
Corresponding author.
E-mail address: [email protected] (Y. Yao).
http://dx.doi.org/10.1016/j.agrformet.2016.03.016
0168-1923/© 2016 Elsevier B.V. All rights reserved.
152 Y. Yao et al. / Agricultural and Forest Meteorology 223 (2016) 151–167
the combined effects of the variations of Rs, Ta, and P on LE for different models on these time scales. This
study provides corrected-modeling evidence for an accelerated global water cycle with climate change.
© 2016 Elsevier B.V. All rights reserved.
1. Introduction potential to be used as a reference dataset to assess the accuracy of
CMIP5 LE results. Yet a detailed comparison between CMIP5 mod-
Latent heat flux (LE) is the flux of heat from the earth’s surface to eled versus global EC observed LE among different land cover types
the atmosphere for soil evaporation, plant transpiration, and evap- has not been performed.
oration from intercepted precipitation by vegetation canopies. LE To maximize the value of GCMs or other multiple datasets, sev-
is a fundamental quantity for understanding ecosystem processes eral merging algorithms have been effectively used to estimate
and functions (Sun et al., 2011) and developing general circulation global terrestrial climatic and hydrologic variables (e.g., air tem-
models (GCMs) and global climatic forecasting and land surface perature, P and LE) with high accuracy (Duan et al., 2007; Miao
models (LSMs) (Liang et al., 2010; Wang and Dickinson, 2012; et al., 2013; Miao et al., 2014; Yang et al., 2012; Yao et al., 2014a;
Wild et al., 2015; Yao et al., 2013). During the past two decades, Yao et al., 2014b). Recent studies have demonstrated that even
eddy covariance (EC) measurement system (e.g. FLUXNET) has been a simple multi-model ensemble, such as simple model averaging
established to measure LE and sensible heat flux (H) exchanges (SMA), is superior to an individual model (Buser et al., 2009; Duan
between the atmosphere and land surface (Baldocchi et al., 2001; and Phillips, 2010; Lambert and Boer, 2001; Weigel et al., 2008).
Liu et al., 2013; Twine et al., 2000; Yao et al., 2015). However, these Sophisticated multi-model ensemble approaches have acquired the
short-term point-based LE measurements by EC are limited due to weights of single-model contributions to improve performance by
the spare coverage. Remote sensing technology has a large spa- training ground-measured observations. Among these complicated
tial coverage but satellites also do not directly measure LE, which ensemble methods, Bayesian model averaging (BMA) is one of the
hampers accurately understanding the long-term variations of ter- most promising methods that combines simulations from multi-
restrial LE due to the substantial uncertainties of the individual ple datasets and a unanimous probability density function (PDF)
datasets. (Duan and Phillips, 2010; Raftery et al., 2005; Wu et al., 2012). Some
The Coupled Model Intercomparison Project phase 5 (CMIP5) studies have highlighted the application of GCM data and the BMA
from the latest 5th Intergovernmental Panel on Climate Change method in global P and air temperature (Ta) simulations (Miao et al.,
(IPCC) assessment report (IPCC-AR5) provides an opportunity to 2013; Miao et al., 2014). However, there is a lack of similar studies
assess global terrestrial LE variations and attributions by coupling that simulate global terrestrial LE by using the BMA method driven
land-atmosphere interaction processes (Dirmeyer et al., 2013; by CMIP5 models and surface EC observations. As a result, little is
Taylor et al., 2012). Relative to the previous beta version, CMIP5 accurately understood regarding spatiotemporal characterization
includes more than 45 GCMs from different modeling groups with of the response of global terrestrial LE to climate change over long
higher spatial and temporal resolution and multiple models for a periods.
single experiment (Taylor et al., 2012).The state-of-the-art GCMs In this study, we evaluated LE simulations of 45 CMIP5 mod-
that are available through CMIP5 are now widely used to investigate els and improved global terrestrial LE simulations among different
the theoretical mechanisms of climatic changes (Covey et al., 2003; land cover types by synthesis of seven best CMIP5 models and
Miao et al., 2013). Compared to the Coupled Model Intercomparison FLUXNET EC observations based on the BMA method. Our study
Project phase 3 (CMIP3) in the 4th IPCC assessment report (IPCC- has three specific objectives: (1) evaluate LE simulations from the
AR4), the GCMs in IPCC-AR5 have improved many more model types, state-of-the-art GCMs of 45 CMIP5 models with a comprehensive
including the Earth System Models (ESMs), with more interactive ground-measured LE flux data set; (2) use the BMA method to merge
components, including aerosols, dynamic vegetation, atmospheric the seven best CMIP5 models to generate a global terrestrial long-
physics and carbon and hydrological cycles (Liu et al., 2013; Miao term (1970–2005) monthly LE; and (3) analyze the spatiotemporal
et al., 2014). Most dynamic, physical and chemical algorithms were variability in the global terrestrial LE and its attributions by com-
also improved in the IPCC AR5 models (Moss et al., 2010; Wild et al., paring the changes of the relevant climatic variables.
2013; Wild et al., 2015). These improvements will help the long-
term climate forecasts, including global terrestrial inter-annual LE prediction.
2. Data
Recently, a much broader comparison of CMIP5 (or CMIP3) data
and other gridded datasets had been attempted to assess GCM
2.1. CMIP5 GCM latent heat flux simulation
values of the global terrestrial latent and sensible heat fluxes. As
reported by Mueller et al. (2011) and Wang and Dickinson (2012),
IPCC-AR5 used more than 45 state-of-the-art GCM results as part
the standard deviation (SD) of the IPCC AR4 simulations within each
2 of CMIP5 for the World Climate Research Programme (WCRP). In
category (4.6 W/m ) is lower than those of the reference datasets,
CMIP5, the same numerical experiments were performed by differ-
including satellite, reanalysis and LSMs datasets (SD varying from
2 ent models of the same protocols, which are accepted for a direct
4.9 to 5.6 W/m ). Wild et al. (2015) reported that the CMIP5 mod-
2 2 comparison of these models (Guilyardi et al., 2013; Ma et al., 2014).
els varied greatly (32–46 W/m for LE, 16–43 W/m for H) in their
The LE simulations of the historical experiments were used for the
calculation of the land mean LE and H, with a global land (including
2 2 GCMs in this study. All of the LE outputs were simulated using
Antarctica) mean LE and H of 38 W/m and 32 W/m , respectively.
the same initial approach, initial time, and rattled physics with an
These inter-comparison studies, however, focus mainly on evalu-
ensemble member set at r1i1p1 (Taylor et al., 2012). Monthly CMIP5
ating the global annual mean and the errors of the surface LE and
◦
GCM LE simulations with 0.56–3.75 spatial resolution were used
H based on gridded datasets, with a few using eddy covariance
in this study. Most CMIP5 GCM LE simulations spanned from 1850 to
observations (Jiménez et al., 2011; Wild et al., 2013). Meanwhile, a
2005, including 45 GCMs used in this study. When combining some
large number of EC observations from the FLUXNET project has the
gridded GCM datasets with different spatial resolutions, they were
Y. Yao et al. / Agricultural and Forest Meteorology 223 (2016) 151–167 153
Table 1
Description of the 45 CMIP5 GCMs used in this study, mean annual global terrestrial LE (excluding Antarctica) and trend of the global terrestrial mean of LE (excluding
Antarctica) from 1970 to 2005.
No. Model name Source Spatial Mean annual global Trend of the global
resolution terrestrial LE terrestrial mean of LE
(excluding (excluding Antarctica)
2 2
Antarctica) (W/m ) (W/m per decade)
◦ ◦
×
1 ACCESS-1-0 Commonwealth Scientific and Industrial Research Organization 1.88 1.24 43.7 0.14
2 ACCESS-1-3 (CSIRO) and Bureau of Meteorology (BOM), Australia 46.9 0.12
◦ ◦
×
3 BCC-CSM-1-1-m Beijing Climate Center, China Meteorological 1.13 1.13 36.0 0.28
◦ ◦
×
4 BCC-CSM-1-1 Administration 2.81 2.81 38.0 0.30
◦ ◦
5 BNU-ESM College of Global Change and Earth System 2.81 × 2.81 44.1 0.43
Science, Beijing Normal University
◦ ◦
×
6 CanCM4 Canadian Centre for Climate Modelling and 2.81 2.81 39.3 0.23
7 CanESM2 Analysis 37.5 0.23
◦ ◦
8 CCSM4 National Center for Atmospheric Research 1.25 × 0.94 43.5 0.21
◦ ◦
×
9 CESM1-BGC Community Earth System 1.25 0.94 43.4 0.12
10 CESM1-CAM5 Model Contributors 40.4 0.07
11 CESM1-FASTCHEM 43.5 0.20
◦ ◦
12 CESM1-WACCM 2.50 × 1.88 43.0 0.17
◦ ◦
×
13 CMCC-CESM Centro Euro-Mediterraneo per I Cambiamenti 3.75 3.75 42.8 0.41 Climatici
◦ ◦
14 CMCC-CM 0.75 × 0.75 36.0 0.25
◦ ◦
15 CMCC-CMS 1.88 × 1.88 39.2 0.34
◦ ◦
×
16 CNRM-CM5-2 Centre National de Recherches Meteorologiques/Centre 1.41 1.41 40.2 0.26
17 CNRM-CM5 Europeen de Recherche et Formation Avancees en Calcul 40.6 0.23 Scientifique
◦ ◦
×
18 CSIRO-Mk-3-6-0 Commonwealth Scientific and Industrial Research 1.88 1.88 39.4 0.14
Organization in collaboration with Queensland Climate
◦ ◦
×
19 CSIRO-Mk3L-1-2 Change Centre of Excellence 5.63 3.21 37.8 0.01
◦ ◦
20 EC-EARTH EC-EARTH consortium 1.13 × 1.00 40.1 0.17
◦ ◦
21 FGOALS-g2 LASG, Institute of Atmospheric Physics, 2.81 × 3.00 43.8 0.23
Chinese Academy of Sciences and CESS,
Tsinghua University
◦ ◦
22 FIO-ESM The First Institute of Oceanography, SOA, China 2.81 × 2.81 44.6 0.18
◦ ◦
23 GFDL-CM3 NOAA Geophysical Fluid 2.50 × 2.00 43.9 0.15
24 GFDL-ESM2G Dynamics Laboratory 43.0 0.31
25 GFDL-ESM2M 42.7 0.35
◦ ◦
26 GISS-E2-H-CC NASA Goddard Institute for Space Studies 2.50 × 2.00 47.2 0.08
27 GISS-E2-H 47.2 0.13
28 GISS-E2-R-CC 45.7 0.15
29 GISS-E2-R 45.8 0.25
◦ ◦
×
30 HadGEM2-AO Met Office Hadley Centre (additional HadGEM2-ES 1.88 1.24 43.8 0.13
realizations contributed by Instituto Nacional de Pesquisas
31 HadGEM2-ES Espaciais) 44.8 0.24
◦ ◦
32 inmcm4 Institute for Numerical Mathematics 2.00 × 1.50 43.6 0.10
◦ ◦
33 IPSL-CM5A-LR Institut Pierre-Simon Laplace 3.75 × 1.88 37.8 0.24
◦ ◦
34 IPSL-CM5A-MR 2.50 ×1.26 37.7 0.17
◦ ◦
35 IPSL-CM5B-LR 3.75 ×1.88 36.0 0.25
◦ ◦
×
36 MIROC-ESM-CHEM Japan Agency for Marine-Earth Science and Technology, 2.81 2.81 47.5 0.43
Atmosphere and Ocean Research Institute (The University
37 MIROC-ESM of Tokyo), and National Institute for Environmental Studies 48.1 0.26
◦ ◦
38 MIROC4h Atmosphere and Ocean Research Institute (The University of 0.56 × 0.56 41.1 0.21
Tokyo), National Institute for Environmental Studies, and
Japan Agency for Marine-Earth Science and Technology
◦ ◦
39 MPI-ESM-LR Max Planck Institute for Meteorology 1.88 × 1.88 41.5 0.15
40 MPI-ESM-MR 42.5 0.36
41 MPI-ESM-P 41.4 0.18
◦ ◦
×
42 MRI-CGCM3 Meteorological Research Institute 1.13 1.13 38.0 0.08
43 MRI-ESM1 38.2 0.22
◦ ◦
44 NorESM1-M Norwegian Climate Centre 2.50 × 1.88 43.0 0.21
45 NorESM1-ME 42.7 0.19
◦ ◦
46 SMA simple model averaging in this study 1.00 × 1.00 41.9 0.19
◦ ◦
47 BMA Bayesian model averaging in this study 1.00 × 1.00 39.7 0.18
interpolated to one degree using the bilinear interpolation method. degree for 1984–2005, which were derived from the Global Energy
Detailed information on the CMIP5 GCMs is summarized in Table 1. and Water Cycle Experiment (GEWEX) Surface Radiation Budget
(SRB) products because of its higher accuracy when compared to
other satellite, reanalysis and GCM R products (Liang et al., 2010;
2.2. Reanalysis and satellite datasets s
Zhang et al., 2010a). To identify the land cover types, before 1999,
◦
the UMD Global Land Cover Classification (AVHRR: UMD CLCF: 1 )
To investigate the impact of the downward surface solar
product (Hansen et al., 1998), which was generated from Advanced
shortwave radiation (Rs), near surface air temperature (Ta) and
Very High Resolution Radiometer (AVHRR) satellites between 1981
precipitation (P) on the variations in surface LE across different
and 1994, was used and after 1999, the Collection 4 MODIS land
regions, reanalysis and satellite datasets are needed. We used the
◦
cover (MOD12C1: CMG, 0.05 ) product (Friedl et al., 2002) was
global monthly surface Rs products at a spatial resolution of one
154 Y. Yao et al. / Agricultural and Forest Meteorology 223 (2016) 151–167
used. The MOD12C1 product was interpolated into one degree using the average of xo. ıf is the normalized standard deviation of the
nearest neighbor interpolation method. Considering that the Euro- LE simulations over the standard deviation of the corresponding LE
pean Centre for Medium Range Weather Forecasts (ECMWF) has observations:
the most accurate T among the available datasets (NASA Data a N
1 2
Assimilation Office, DAO, and National Centers for Environmental − N (xs x¯ s)
ıs s=1
Prediction Reanalysis, NCEP), the monthly Ta data for 1984–2005 = =
ıf (3)
were extracted from the ERA-Interim reanalysis product with a ıo N 1 2