1. Introduction
It is now well established that the chaotic nature of the atmosphere severely limits the predictability of weather, while the slowly varying sea surface temperature (SST) and land surface states can enhance the predictability of longer-time-scale atmospheric variations (e.g., Shukla 1993, 1998; Shukla et al. 2000; Graham et al. 1994; Koster et al. 2000; Dirmeyer et al. 2003; Quan et al. 2004). The Global Land–Atmosphere Coupling Experiment (GLACE; see Koster et al. 2004, 2006) builds a framework to objectively estimate the potential contribution of land states to atmospheric predictability (called the land–atmosphere coupling strength) in numerical weather and climate models. By averaging the estimated land–atmosphere coupling strength from 12 models participating in GLACE, an ensemble average coupling strength is obtained. However, the coupling strength varies widely among models. The discrepancy is certainly related to differences in the parameterization of processes and their complex interactions, from soil hydrology and vegetation physiology to boundary layer, cloud, and precipitation processes. It is difficult to determine what causes the relatively strong or weak coupling strengths seen in the individual models.
Guo et al. (2006) separated the total impact of soil moisture on precipitation into the ability of soil moisture to affect evaporation and the ability of evaporation to affect precipitation. They find that most of the differences between the models are associated with the first process; however, the second process also plays a role. This water cycle pathway for land–atmosphere coupling has been established as the major means for affecting precipitation predictability (e.g., Dirmeyer 2006; Wang et al. 2007). Lawrence and Slingo (2005) investigated the weak land–atmosphere coupling strength in one GLACE participating model, the third Hadley Centre Atmosphere Model (HadAM3), and found that the lack of coupling between the soil moisture and boundary layer moist static energy, together with the overly frequent rainfall, may be the reason for the weak coupling. These studies do not pinpoint the exact source (within either the land model or the atmospheric model) that determines the unusual coupling strength. In our companion paper (Wei et al. 2010, hereafter Part I), we coupled the Center for Ocean–Land–Atmosphere Studies (COLA) atmospheric general circulation model (AGCM) to three different land models. Here we compare the land–atmosphere coupling strength of the three coupled configurations, and locate the approximate source of the variations in coupling strength.
Guo et al. (2006) and Lawrence and Slingo (2005) analyzed the coupling strength from a local point of view, without relating it to large-scale processes. Is there a connection between the pattern of land–atmosphere coupling and global climate? It is already known that some climate regimes (e.g., transitional zones between arid and humid climates) tend to have stronger land–atmosphere coupling (Koster et al. 2000, 2004; Wang et al. 2007), but this may not be the whole story. A theoretical study found that the coupling strength is also related to climate variability, and the strength of the external forcing can affect the coupling strength and the location of coupling hot spots (Wei et al. 2006). Even less is known about how the land–atmosphere coupling is related to the different time scales of climate variability. The intraseasonal variability of precipitation has a strong influence on the soil moisture variability (Wei et al. 2008), but little study has been performed on the connection between this variability and land–atmosphere coupling. By comparison, the connection between intraseasonal variability and air–sea interaction has been studied much better (e.g., Pegion and Kirtman 2008a,b).
In this paper, GLACE-type experiments with the COLA AGCM coupled to three land models are performed. We examine the large-scale connections between precipitation predictability, land–atmosphere coupling strength, and climate variability. Based on the connections found, the model-estimated land–atmosphere coupling strength can be calibrated to account for errors in the simulation of precipitation variability, a quantity that is observable in the large scale and found to be closely related to the coupling strength. The next section introduces the experiments and observational datasets, followed by results in section 3 and conclusions and discussion in section 4.
2. Experiments and datasets
In Part I, we have shown results from long-term integrations with the COLA AGCM coupled to three land models: the Simplified Simple Biosphere Model (SSiB), the Community Land Model (CLM) 3.5, and the Noah land model. In this study, GLACE-type experiments are performed individually with these three coupled configurations. Our experiments are performed exactly the same as in GLACE. A brief description of GLACE and its experimental design are given in the appendix.
Three observational precipitation datasets are used to evaluate the model results. A brief description is given below.
a. GPCP-1DD
The Global Precipitation Climatology Project one-degree daily (GPCP-1DD) precipitation dataset provides daily, global 1° × 1° gridded fields of precipitation totals from October 1996 to the present (available online at http://precip.gsfc.nasa.gov/gpcp_daily_comb.html; Huffman et al. 2001). It is a combined observation-only dataset, which is based on satellite estimates of rainfall and is scaled to match the monthly accumulation provided by the GPCP satellite–gauge (SG) product, which combines monthly satellite and gauge observations on a 2.5° × 2.5° grid (Huffman et al. 1997; Adler et al. 2003). In this study, we use data from 1997 to 2007.
b. TMPA
The Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA) 3B42 (available online at http://trmm.gsfc.nasa.gov/3b42.html; Huffman et al. 2007) provides a calibration-based scheme for combining precipitation estimates from multiple satellites and gauge analyses where feasible, at fine scales (0.25° × 0.25° and 3 hourly) and over 50°N–50°S. We use the research product, which uses a better calibrator than the real-time product and rescales the monthly sums of the 3-hourly fields to a monthly gauge analysis. The data we used here is from 1998 to 2007.
c. CMORPH
The Climate Prediction Center (CPC) morphing technique (CMORPH; Joyce et al. 2004) produces global precipitation analyses at very high spatial and temporal resolution (half-hourly and 8 km) by using satellite microwave observations only. Instead of using IR brightness temperature algorithms to estimate rainfall rate, as in GPCP-1DD and TMPA, CMORPH uses IR-derived advection vectors to propagate the relatively high-quality precipitation estimates derived from passive microwave data. The shape and intensity of the precipitation features are modified (morphed) during the time between the microwave sensor scans by performing a time-weighted linear interpolation. This process yields spatially and temporally complete microwave-derived precipitation analyses, independent of the IR temperature field. The data are at 0.25° × 0.25° and 3-hourly resolution, ranging from 60°N to 60°S. Data from 2003 to 2008 are used in this study.
Although these datasets are all based on satellites observations, they used different techniques to retrieve precipitation, so intercomparison among these datasets is a kind of verification. Verifications of these datasets have been done by comparing them with gauge- or radar-based estimates (e.g., McPhee and Margulis 2005; Skomorowski et al. 2001; Gebremichael et al. 2005; Su et al. 2008; Zeweldi and Gebremichael 2009). There is generally a good agreement between these datasets and ground-based estimates, especially over large spatial and long temporal scales.
3. Results
a. Results from GLACE-type experiments
Figure 1 shows the Ω values of total precipitation for ensembles W [16-member control experiment for June–August (JJA)] and S (subsurface soil wetness specified in all ensemble members from an arbitrarily chosen member of W) and their difference Ω(S) − Ω(W) (see the appendix for complete definitions). No matter with which land model the AGCM is coupled, Ω shows similar patterns, with the largest values in the tropical rain belt where the SST forcing has the strongest influence (Shukla 1998). The patterns of Ω(W) and Ω(S) are very similar, with large differences [Ω(S) − Ω(W)] mainly over the regions with high Ω values. This indicates that the land–atmosphere coupling strength may be strongly influenced by the external forcing. By “external,” we mean that the forcing is from outside of the land–atmosphere system, such as that from SST. The patterns of Ω(S) − Ω(W) for different model configurations have much lower similarity than for Ω [spatial correlations are 0.18–0.23 for Ω(S) − Ω(W) and 0.6–0.72 for Ω]. Globally, COLA–SSiB has the strongest land–atmosphere coupling strength, while COLA–Noah has the weakest. The potential contribution of subsurface soil moisture to the total potential predictability {[Ω(S) − Ω(W)]/Ω(S)} is approximately 30%, 28%, and 18% for COLA–SSiB, COLA–CLM, and COLA–Noah, respectively. The differences seen mainly should be from the land models’ different connections between soil moisture and surface fluxes, because they are coupled to the same AGCM.
b. A decomposed view of land–atmosphere coupling strength
As discussed above, the slowly varying boundary forcing may play an important role in the similarity of the precipitation time series in different ensemble members (magnitude of Ω). It is very likely that the “fingerprints” of these slow forcings also exist in the precipitation time series. An effective way to examine this is to decompose the time series by frequency bands. After ignoring the first 8 days of integration of each JJA range to avoid possible problems associated with the initial shock to the model atmosphere, as in calculating Ω, there are 84 days that are left for analysis. We performed a discrete Fourier transform [discussed in detail in Ruane and Roads (2007)] to decompose the daily time series into three frequency bands: fast synoptic (2–6 days), slow synoptic (6–20 days), and intraseasonal (20–84 days). The choice of these frequency bands is arbitrary; other comparable choices give similar results. Note that the time series may contain a portion of the seasonal cycle, but because of the length of the time series we refer to the 20–84-day variation as intraseasonal.
Figure 2 shows the variance percentages of precipitation in these three bands for model simulations and three observational datasets. The three model simulations are highly consistent in their variance distributions, as are the three observational datasets (note that for TMPA and CMORPH the latitudinal extent is limited). However, all of the model simulations underestimate the high-frequency (fast synoptic) variance and overestimate the low-frequency (intraseasonal) variance, especially over the tropics and subtropics. Multiyear simulations of these models, as described Part I, have similar variance percentage distributions as these GLACE-type simulations (not shown). This overestimation of low-frequency variance is also consistent with the overestimation of precipitation persistence shown in Part I. Although these observational datasets have uncertainties and errors, Sun et al. (2006) found that no matter what observational dataset is used, this model bias is relatively large compared to the uncertainties between the observations. This bias of the models may be related to a well-known problem in AGCM parameterizations: premature triggering of convection so that precipitation falls too frequently but too lightly in intensity (Trenberth et al. 2003; Sun et al. 2006; Ruane and Roads 2007).
For theoretical white noise, the variance at each frequency is the same, so the variance percentages are determined by the widths of the frequency bands. Therefore, the variance percentages of the three bands (from fast to slow) for white noise are 69%, 21%, and 10%, respectively. Overall, both the model results and observations follow a red spectrum, with variance percentages less than white-noise values at high frequencies and greater than white-noise values at low frequencies.
In Fig. 2, the spatial correlations between Ω and the percentage of intraseasonal variance (IV) are very high (right column), but the correlations of Ω with the other two frequency bands are negative (left two columns). Ensemble S (not shown) has similar results as ensemble W. This demonstrates that regions with a larger percentage of IV tend to have a higher value of Ω, no matter whether soil moisture is interactive (W) or not (S). This is not unexpected because, as we discussed above, most of the precipitation predictability (or Ω) is from the slowly varying boundary forcing. Regions with stronger boundary forcing may be constrained to show more low-frequency variation, and the precipitation time series will be more similar in an ensemble (larger Ω). For ensemble S, the prescribed subsurface soil moisture is also one of the slow boundary forcings. However, compared to the ensemble without the constraint of this slow forcing (W), ensemble S does not show significant change in the global pattern of variance distribution. These results show that different land models or land states do not matter much for the global pattern of precipitation variance distribution, which may be determined by other factors, such as global climate (SST, radiation, etc.) and the convection scheme. Ruane and Roads (2008) obtained similar results from a global assimilation system. They found that two different land models did not produce a noticeable difference in variance distribution of precipitation, but two different convection schemes can have significantly different effects. Wilcox and Donner (2007) also showed that the convection parameterization of a GCM can greatly impact the frequency distribution of rain rate, and their model with a relaxed Arakawa–Schubert formulation of cumulus convection (also used in the COLA AGCM) exhibits a strong bias toward excessive light rain events and too-few heavy rain events.
The above shows that neither the land model nor the soil moisture has a great impact on the global pattern of precipitation variability and predictability. However, their impact may be strong at regional scales. The difference Ω(S) − Ω(W) shows the impact of soil moisture. It tries to remove the effects of the same strong external forcing on both S and W and highlights the role of soil moisture, although we understand that the effects of those forcings cannot be completely removed in a nonlinear system (more discussion on this aspect follows in section 3c). The spatial correlations between percentage of IV and Ω(S) − Ω(W) are also shown in Fig. 2 (as the last number in right column). They are generally weaker than the correlation with Ω, indicating that something other than the low-frequency external forcing is playing an important role in Ω(S) − Ω(W). This should be the impact of subsurface soil moisture.
To demonstrate the role of each frequency band in Ω and Ω(S) − Ω(W), we recalculate the variables in Fig. 1 using the decomposed precipitation time series in each ensemble member. The result calculated with only the intraseasonal precipitation (higher frequencies are filtered out) is shown in Fig. 3; the other two frequency bands result in very weak values and are not shown. It can be seen in Fig. 3 that the patterns of Ω and Ω(S) − Ω(W) are very similar to those in Fig. 1. This is because the intraseasonal component of precipitation, mostly caused by the same low-frequency external forcing, has high consistency among the ensemble members, while the high-frequency component of precipitation, mostly from chaotic atmospheric dynamics, is generally incoherent among the ensemble members. (We do not compare the amplitudes of the two figures because their time series have different degrees of freedom.) It is interesting that the goal of GLACE, as described in Koster et al. (2006), is to examine land surface’s control on precipitation variability on time scales of ∼1 week (synoptic scale), but our results here show that the main impact of the land surface in the models is on intraseasonal time scales.
c. Conceptual relationships
In our three model simulations, the atmospheric model is the same, and we have shown that Ω (or F) is similar for all three of the models, but Ω(S) − Ω(W) has larger spatial differences (Fig. 1). From Eq. (2), the differences between the models should be mainly from differences in α(S) − α(W), which appears to be generally smaller when the COLA AGCM is coupled to the Noah model than to the other two land models. However, we must recognize that the impacts of external forcing and local land states may not be explicitly separable as in Eqs. (1) and (2). This simplification is for the conceptual understanding of the relationships.
d. Calibration of the estimated GLACE land–atmosphere coupling strength
To examine whether our results are consistent with those from other models, we look at the GLACE dataset (Koster et al. 2006). Figure 4 shows the percentage of IV for 12 models participating in GLACE, and their respective spatial correlations with Ω and Ω(S) − Ω(W). As for the COLA AGCM above, ensemble S (not shown) shows very similar results to ensemble W. Also, all of the models here have overestimated the percentage of IV. The global average percentages for the three observational datasets are consistently 15% (Fig. 2), while the 12 models show values of 21%–28%. (The COLA AGCM used in the original GLACE experiment is an older version than the one used here for coupling to multiple land schemes, and the versions of SSiB used are also different.) The spatial distributions of the variance percentage differ widely, and the overestimations are not limited to the tropics and subtropics. The spatial correlations of variance percentage and Ω are always high; even the lowest value [0.32 from the University of California, Los Angeles (UCLA)] is well over the 99% confidence level (assuming that the grid points are independent). The correlations with Ω(S) − Ω(W) are always lower than with Ω, and some models even show weak negative correlations, indicating the spatially heterogeneous influence of soil moisture [α(S) − α(W) in Eq. (2)]. Therefore, the GLACE models and our models show similar relationships between Ω, Ω(S) − Ω(W), and the percentage of IV.
Although the spatial correlations between Ω and the percentage of IV are very significant for all of the models, the correlation coefficients are not necessarily high (e.g., UCLA). This means that the strong connection between Ω and F described in Eq. (1) may not happen over all of the regions. The noise may dampen the connection in some cases. We then calculated the correlation between Ω and the percentage of IV across the 12 GLACE models (a correlation with a sample size of 12; see Fig. 5). It indicates a general connection between Ω and IV for all of the 12 models at each grid point. The regions with strong positive correlations are where Ω and IV have a strong connection for almost all of the models. Over other regions where the correlation is positive but not so strong, their relationship may not be so consistent for the models, but larger F can generally lead to larger Ω. We can see in Fig. 5 that over 90% of the land areas show positive correlations, which supports our assumptions of their relationship.
Based on the above analysis, we perform two separate calibrations to the original Ω(S) − Ω(W)—a strong one to calibrate all of the regions where the correction in Fig. 5 is positive, and a weak one to calibrate only the regions where the correlation in Fig. 5 is over 95% confidence level (0.576); both use Eq. (3). The results are shown in Fig. 6. It can be seen that the coupling strength reduces by approximately half after strong calibration, and only a few hot spots survive. The result after weak calibration is similar to the original, but shows weakened coupling strength in the U.S. Great Plains and Mexico, northern India and Pakistan, and Mauritania. Generally, the calibrated coupling strengths are weaker than that from the GLACE, but the pattern does not change much.
In this study, the impact of the model overestimation of low-frequency variance on land–atmosphere coupling is discussed from the perspective of atmosphere, but the effect of this overestimation also shows over land. More low-frequency variance usually means more sustained dry and wet periods and higher evaporation variability, which can lead to a more robust precipitation response (higher predictability). Also, as mentioned, this model bias is manifested as overly frequent rainfall at reduced intensity. This increases infiltration and reduces runoff, enhancing subsequent evapotranspiration and increasing soil moisture memory, which then leads to the overestimation of land–atmosphere coupling.
Because of the changing and heterogeneous nature of the relationship between Ω and F, our calibration method is not flawless. However, the results illustrate how the amplitude and distribution of coupling strength may change after some rectification of the model bias. The unique design of GLACE makes its results difficult to evaluate by directly comparing them with observational variables, even if these large-scale observations exist, because the GLACE metric is based on ensemble statistics, and observations present us with only one “ensemble member.” Some recent studies using observationally based data have cast doubt on the strong coupling strength in the Great Plains (Ruiz-Barradas and Nigam 2005, 2006; Zhang et al. 2008), but whether their results are comparable to the GLACE results need further study. On the other hand, Wang et al. (2007) have shown that less restrictive metrics than measuring ensemble coherence, such as the change in overall precipitation variance between S and W cases, reveals even more areas of apparently strong coupling strength.
4. Conclusions and discussion
Coupling one AGCM to three different land models gives us an unprecedented opportunity to study the role of different components in land–atmosphere interaction. The behavior of the coupled models and their land–atmosphere interaction are investigated by a set of GLACE-type experiments. The connections found from our model results are consistent with that from the GLACE dataset. Our findings lead to a calibration method for the land–atmosphere coupling strength based on evident systematic error in AGCMs. The main findings are summarized below:
Different land models or subsurface soil moisture have little influence on the global pattern of precipitation predictability (Ω) and variance distribution because of the stronger control of other factors. While the regional effect of soil moisture can be highlighted by the difference of Ω from two ensembles (one with a specified identical soil moisture time series across all ensemble members and one that is free running), which shows contrasting patterns for the three models. This indicates that the segment of the water cycle pathway that connects soil moisture to evaporation is a key component of the coupling strength, in agreement with Guo et al. (2006).
Most of the precipitation predictability and land–atmosphere coupling strength [Ω(S) − Ω(W)] are associated with the intraseasonal component of precipitation in the models, although they only account for a small percentage (∼⅕) of the total variance.
The GLACE coupling strength can be conceptually decomposed into the impact of the slow external forcing F and the local impact of soil moisture α. Both F and the potential soil moisture influence (α difference between two ensembles) combine to determine the pattern of the coupling strength.
The models participating in GLACE have overestimated the low-frequency variance percentage and underestimated the high-frequency variance percentage of precipitation. Based on the findings in this paper, we calibrate the land–atmosphere coupling strength estimated by GLACE. It is found that the calibrated coupling strengths are generally weaker than that from GLACE but the patterns are nearly the same.
The integration period of these GLACE-type experiments is JJA, so the longest time scale in this study is intraseasonal. The effect of land may be different for longer time scales. Our judgment on the connections between external forcing, low-frequency precipitation variability, and Ω and Ω(S) − Ω(W) is based on model results, so they may not be absolutely true in reality. Some intraseasonal variations of precipitation may not come from external forcing, but from the flow instabilities (especially in the midlatitudes; e.g., Charney and Shukla 1981) that are not fully understood nor properly simulated. The observational datasets used in this study are all based on satellite observations, and they may have consistent biases. In spite of these limitations, this study qualitatively separates the role of external forcing and local soil moisture on precipitation variability and predictability, and increases our understanding of land-atmosphere interaction.
Acknowledgments
We thank Alex C. Ruane for help with the variance decomposition, and all of the model groups participating in GLACE for providing their experimental results. Comments from three reviewers have enhanced the quality of the paper. This research was supported by National Oceanic and Atmospheric Administration Award NA06OAR4310067. The computing was completed on NCAR supercomputers.
REFERENCES
Adler, R. F., and Coauthors, 2003: The Version-2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979–present). J. Hydrometeor., 4 , 1147–1167.
Charney, J. G., and J. Shukla, 1981: Predictability of monsoons. Monsoon Dynamics, Sir J. Lighthill and R. P. Pearce, Eds., Cambridge University Press, 99–109.
Dirmeyer, P. A., 2006: The hydrologic feedback pathway for land–climate coupling. J. Hydrometeor., 7 , 857–867.
Dirmeyer, P. A., M. J. Fennessy, and L. Marx, 2003: Low skill in dynamical prediction of boreal summer climate: Grounds for looking beyond sea surface temperature. J. Climate, 16 , 995–1002.
Gebremichael, M., W. F. Krajewski, M. L. Morrissey, G. J. Huffman, and R. F. Adler, 2005: A detailed evaluation of GPCP 1° daily rainfall estimates over the Mississippi River basin. J. Appl. Meteor., 44 , 665–681.
Graham, N. E., T. P. Barnett, R. Wilde, M. Ponater, and S. Schubert, 1994: On the roles of tropical and midlatitude SSTs in forcing interannual to interdecadal variability in the winter Northern Hemisphere circulation. J. Climate, 7 , 1416–1441.
Guo, Z., and Coauthors, 2006: GLACE: The Global Land–Atmosphere Coupling Experiment. Part II: Analysis. J. Hydrometeor., 7 , 611–625.
Huffman, G. J., and Coauthors, 1997: The Global Precipitation Climatology Project (GPCP) combined precipitation dataset. Bull. Amer. Meteor. Soc., 78 , 5–20.
Huffman, G. J., R. F. Adler, M. Morrissey, D. T. Bolvin, S. Curtis, R. Joyce, B. McGavock, and J. Susskind, 2001: Global precipitation at one-degree daily resolution from multisatellite observations. J. Hydrometeor., 2 , 36–50.
Huffman, G. J., and Coauthors, 2007: The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeor., 8 , 38–55.
Joyce, R. J., J. E. Janowiak, P. A. Arkin, and P. Xie, 2004: CMORPH: A method that produces global precipitation estimates from passive microwave and infrared data at high spatial and temporal resolution. J. Hydrometeor., 5 , 487–503.
Koster, R. D., M. J. Suarez, and M. Heiser, 2000: Variance and predictability of precipitation at seasonal-to-interannual time scales. J. Hydrometeor., 1 , 26–46.
Koster, R. D., and Coauthors, 2004: Regions of strong coupling between soil moisture and precipitation. Science, 305 , 1138–1140.
Koster, R. D., and Coauthors, 2006: GLACE: The Global Land–Atmosphere Coupling Experiment. Part I: Overview. J. Hydrometeor., 7 , 590–610.
Lawrence, D. M., and J. M. Slingo, 2005: Weak land–atmosphere coupling strength in HadAM3: The role of soil moisture variability. J. Hydrometeor., 6 , 670–680.
McPhee, J., and S. A. Margulis, 2005: Validation and error characterization of the GPCP-1DD precipitation product over the contiguous United States. J. Hydrometeor., 6 , 441–459.
Pegion, K., and B. P. Kirtman, 2008a: The impact of air–sea interactions on the predictability of the tropical intraseasonal oscillation. J. Climate, 21 , 5870–5886.
Pegion, K., and B. P. Kirtman, 2008b: The impact of air–sea interactions on the simulation of tropical intraseasonal variability. J. Climate, 21 , 6616–6635.
Quan, X. W., P. J. Webster, A. M. Moore, and H. R. Chang, 2004: Seasonality in SST-forced atmospheric short-term climate predictability. J. Climate, 17 , 3090–3108.
Ruane, A. C., and J. O. Roads, 2007: 6-hour to 1-year variance of five global precipitation sets. Earth Interactions, 11 .[Available online at http://EarthInteractions.org].
Ruane, A. C., and J. O. Roads, 2008: Diurnal to annual precipitation sensitivity to convective and land surface schemes. Earth Interactions, 12 .[Available online at http://EarthInteractions.org].
Ruiz-Barradas, A., and S. Nigam, 2005: Warm season rainfall variability over the U.S. Great Plains in observations, NCEP and ERA-40 reanalyses, and NCAR and NASA atmospheric model simulations. J. Climate, 18 , 1808–1830.
Ruiz-Barradas, A., and S. Nigam, 2006: Great Plains hydroclimate variability: The view from North American Regional Reanalysis. J. Climate, 19 , 3004–3010.
Shukla, J., 1993: Predictability of short-term climate variations. Prediction of Interannual Climate Variations, J. Shukla, Ed., NATO ASI Series I: Global Environmental Change, Vol. 6, Springer-Verlag, 217–232.
Shukla, J., 1998: Predictability in the midst of chaos: A scientific basis for climate forecasting. Science, 282 , 728–731.
Shukla, J., and Coauthors, 2000: Dynamical seasonal prediction. Bull. Amer. Meteor. Soc., 81 , 2593–2606.
Skomorowski, P., F. Rubel, and B. Rudolf, 2001: Verification of GPCP-1DD global satellite precipitation products using MAP surface observations. Phys. Chem. Earth, 26 , 403–409.
Su, F., Y. Hong, and D. P. Lettenmaier, 2008: Evaluation of TRMM Multisatellite Precipitation Analysis (TMPA) and its utility in hydrologic prediction in the La Plata basin. J. Hydrometeor., 9 , 622–640.
Sun, Y., S. Solomon, A. Dai, and R. W. Portmann, 2006: How often does it rain? J. Climate, 19 , 916–934.
Trenberth, K. E., A. Dai, R. M. Rasmussen, and D. B. Parsons, 2003: The changing character of precipitation. Bull. Amer. Meteor. Soc., 84 , 1205–1217.
Wang, G., Y. Kim, and D. Wang, 2007: Quantifying the strength of soil moisture–precipitation coupling and its sensitivity to changes in surface water budget. J. Hydrometeor., 8 , 551–570.
Wei, J., R. E. Dickinson, and N. Zeng, 2006: Climate variability in a simple model of warm climate land–atmosphere interaction. J. Geophys. Res., 111 , G03009. doi:10.1029/2005JG000096.
Wei, J., R. E. Dickinson, and H. Chen, 2008: A negative soil moisture precipitation relationship and its causes. J. Hydrometeor., 9 , 1364–1376.
Wei, J., P. A. Dirmeyer, Z. Guo, L. Zhang, and V. Misra, 2010: How much do different land models matter for climate simulation? Part I: Climatology and variability. J. Climate, 23 , 3120–3134.
Wilcox, E. M., and L. J. Donner, 2007: The frequency of extreme rain events in satellite rain-rate estimates and an atmospheric general circulation model. J. Climate, 20 , 53–69.
Yamada, T. J., R. D. Koster, S. Kanae, and T. Oki, 2007: Estimation of predictability with a newly derived index to quantify similarity among ensemble members. Mon. Wea. Rev., 135 , 2674–2687.
Zeweldi, D. A., and M. Gebremichael, 2009: Evaluation of CMORPH precipitation products at fine space–time scales. J. Hydrometeor., 10 , 300–307.
Zhang, J., W-C. Wang, and J. Wei, 2008: Assessing land–atmosphere coupling using soil moisture from the Global Land Data Assimilation System and observational precipitation. J. Geophys. Res., 113 , D17119. doi:10.1029/2008JD009807.
APPENDIX
Global Land–Atmosphere Coupling Experiment


Theoretically, if the 16 members of an ensemble have identical time series of X,
The GLACE parameter Ω for precipitation from ensembles (left) W and (middle) S, and (right) their difference. (top) COLA–SSiB, (middle) COLA–CLM, and (bottom) COLA–Noah. The global mean (land only) value of each panel is shown at the left corner.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1
The average variance (%) of JJA daily precipitation time series in three frequency bands: (left) fast synoptic (2–6 days), (middle) slow synoptic (6–20 days), and (right) intraseasonal (20–84 days). (top three) Model simulations: COLA–SSiB, COLA–CLM, and COLA–Noah (all from ensemble W; ensemble S has similar results), and (bottom three) observations: GPCP-1DD, TMPA, and CMORPH. TMPA is available in 50°S–50°N, and CMORPH is available in 60°S–60°N. The value (or three values) at the left corner of each panel is the global mean percentage [the spatial correlations of the variance percentage with Ω(W), and with Ω(S) − Ω(W)]. The three observational datasets are shown at 1° × 1° grid; interpolating them to model grid does not affect the results.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1
As in Fig. 1, but computed only with the intraseasonal component of precipitation time series.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1
The variance (%) of 20–84-day precipitation (IV) during JJA for 12 models participating in GLACE (all from ensemble W; ensemble S has very similar results). The definitions of the three values at the left corner of each panel are the same as in Fig. 2.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1
The correlation between Ω and the percentage of IV across the 12 models participating in GLACE. Data are from ensemble W.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1
The estimated land–atmosphere coupling strength [Ω(S) − Ω(W)]. (a) From the original GLACE dataset; (b) all the regions with positive correlation in Fig. 5 are calibrated; and (c) only the regions where correlation in Fig. 5 exceeds 0.576 (95% confidence level) are calibrated. The global average value is shown at the left corner of each panel.
Citation: Journal of Climate 23, 11; 10.1175/2010JCLI3178.1