Understanding Land–Atmosphere–Climate Coupling from the Canadian Prairie Dataset

Analysis of the hourly Canadian Prairie data for the past 60 years has transformed our quantitative understanding of land–atmosphere–cloud coupling. The key reason is that trained observers made hourly estimates of the opaque cloud fraction that obscures the sun, moon, or stars, following the same protocol for 60 years at all stations. These 24 daily estimates of opaque cloud data are of sufficient quality such that they can be calibrated against Baseline Surface Radiation Network data to yield the climatology of the daily short-wave, long-wave, and total cloud forcing (SWCF, LWCF and CF, respectively). This key radiative forcing has not been available previously for climate datasets. Net cloud radiative forcing changes sign from negative in the warm season, to positive in the cold season, when reflective snow reduces the negative SWCF below the positive LWCF. This in turn leads to a large climate discontinuity with snow cover, with a systematic cooling of 10 ◦C or more with snow cover. In addition, snow cover transforms the coupling between cloud cover and the diurnal range of temperature. In the warm season, maximum temperature increases with decreasing cloud, while minimum temperature barely changes; while in the cold season with snow cover, maximum temperature decreases with decreasing cloud, and minimum temperature decreases even more. In the warm season, the diurnal ranges of temperature, relative humidity, equivalent potential temperature, and the pressure height of the lifting condensation level are all tightly coupled to the opaque cloud cover. Given over 600 station-years of hourly data, we are able to extract, perhaps for the first time, the coupling between the cloud forcing and the warm season imbalance of the diurnal cycle, which changes monotonically from a warming and drying under clear skies to a cooling and moistening under cloudy skies with precipitation. Because we have the daily cloud radiative forcing, which is large, we are able to show that the memory of water storage anomalies, from precipitation and the snowpack, goes back many months. The spring climatology shows the memory of snowfall back through the entire winter, and the memory in summer, goes back to the months of snowmelt. Lagged precipitation anomalies modify the thermodynamic coupling of the diurnal cycle to the cloud forcing, and shift the diurnal cycle of the mixing ratio, which has a double peak. The seasonal extraction of the surface total water storage is a large damping of the interannual variability of precipitation anomalies in the growing season. The large land-use change from summer fallow to intensive cropping, which peaked in the early 1990s, has led to a coupled climate response that has cooled and moistened the growing season, lowering cloud-base, increasing equivalent potential temperature, and increasing precipitation. We show a simplified energy balance of the Prairies during the growing season, and its dependence on reflective cloud.


Introduction
Understanding land-atmosphere-climate coupling is challenging, because so many coupled processes are involved: soil temperature and moisture, vegetation types, properties and coverage, near-surface temperature and humidity, the atmospheric boundary layer, the shallow and deep cloud fields which determine the surface radiation balance and surface precipitation, and the soil hydraulic properties that determine the surface and deep runoff, to name only the local components. In the cold season, precipitation falls as snow, and the surface accumulation increases the albedo, and stores water until snowpack melt.
The coupling between the energy and water cycles at the land surface is central to hydrometeorology, and important to weather forecasts on timescales from days to seasons. Earlier reviews [1,2] looked at hydrometeorology from the global modeling perspective using model reanalysis data, which showed how net long-wave and short-wave radiation, cloud cover, surface fluxes, diurnal temperature range, soil moisture, and cloud-base height were coupled on daily timescales over river basins [3]. On daily timescales, the land-atmosphere system is fully coupled, so that errors in the model representation of processes in the soil, vegetation, boundary layer, and cloud fields can rapidly bias a model forecast. Nonetheless, this model perspective was a strong motivation for our analyses of the Canadian Prairie data, and the search for a quantitative description of the fully coupled observed system.
Historically, many climate and hydrometeorology studies have been largely based on precipitation, temperature, and humidity, for which long-term records are available [4][5][6]. However, the diurnal cycle is driven primarily by the surface radiation balance, which depends critically on the daily cloud fields, which are generally unknown in climate records, until satellite-based estimates became available. We cannot study the fully coupled nature of the land-atmosphere-climate system without the surface radiation budget. www.videleaf.com The Canadian Prairie data is, however, an exception, because observers, typically at most major airports, were trained to estimate hourly the opaque cloud fraction in tenths, by cloud level and in total. The definition of opaque cloud is -opaque to the sun, moon, or stars‖; and this protocol has been followed by trained observers hourly for 60 years across the Prairies. With 24 observations per day (almost none are missing), we have representative estimates of the fraction of the daytime short-wave clear-sky (SWCS) flux reaching the surface, and the fraction of the sky that is opaque to outgoing long-wave (LW) radiation for over 600 station-years of data. Because there are 17 years of Baseline Surface Radiation Network (BSRN) data just 25 km south of Regina, SK, we were able to calibrate the opaque cloud data in terms of the LW and SW cloud forcing (Section 3.3). This is transformative, as it meant that we were able to determine quantitatively the climate coupling between the cloud radiative forcing, and the diurnal and seasonal cycle. In addition, simply because we can separate the large radiative impact of clouds from the impact of precipitation, we can better quantify the hydrometeorological processes that couple the energy and water cycle, and observe the long-term memory of precipitation anomalies. In recent years, data from the Gravity Recovery and Climate Experiment (GRACE) [7,8] give estimates of the seasonal drawdown of total water storage. Canadian archives also record agricultural crops grown on the Prairies back to 1955, so we could assess the large climate impact of the shift away from summer fallowing to continuous cropping. This paper is not a conventional review of the literature on landatmosphere-climate coupling. Instead, it is a synthesis of our key conclusions from a series of Prairie data analysis papers [9][10][11][12][13][14][15][16]. Readers interested in more details, or in the evolution of our thinking, can refer back to these original papers. It is remarkable that the long-term Prairie climate dataset, with better cloud observations, have taken our understanding of land-atmosphereclimate coupling to a new level. In retrospect, much of our analysis could have been done two decades ago, but the data was not widely accessible. www.videleaf.com Section 2 discusses the Prairie data and our analysis methods. Section 3 outlines how the climate is coupled to opaque cloud and snow cover on daily timescales, and shows the difference in cloud forcing between warm and cold season with snow. Section 4 looks at the long-term memory of precipitation anomalies, both using multiple regression for the cold and warm season memory, and the dependence of the diurnal coupling on opaque cloud and precipitation anomalies. Section 5 looks at how the seasonal extraction of the surface total water storage dampens the interannual variability of precipitation anomalies in the growing season, and how the large land-use change from summer fallowing to intensive cropping has led to a coupled climate response. Finally, we return to reanalysis data to show how the growing season surface and top-of-atmosphere (TOA) budgets change with cloud cover. Section 6 summarizes our conclusions.  Table 1 lists the station locations and elevation, and the two letter code is used to identify stations in the figures and text. These have an hourly pressure, temperature, relative humidity, wind speed and direction, opaque cloud, and derived radiation, starting in 1953, for all stations, except RG, and MJ, which start in 1954, and ED, which starts in 1961. We accessed the data through June 2011. The hours of missing data are remarkably small. For key stations, such as Calgary, Regina and Winnipeg, more than 99.9% of the days have no missing hours in the first 40 years. In more recent years since 1994, the number of days with less than 23 h of data is typically less than 1%. A few stations (PS in 1992;MJ in 1998;LE and MH in 2006) shifted to daytime-only observation in recent years, because of reduced staffing. The stations also have daily precipitation and snow depth (except for PO), although the last year with complete precipitation data was 1994 for SW, 2005 for LE and MH, 2007  for WI, 2008 for RG, and 2009 for SK. The snow depth data  begins in 1955, and ends in 1994 for SW, 1997 for MJ, 2002 for  LE, 2003 for WI, 2005 for RG and SK, 2006 for ES, GP, MH,  PA, RD, RG, and TP, and 2010 for ED. This synthesis paper extracts significant results from many analyses [9][10][11][12][13][14][15][16], which use different subsets of the data, ranging from all station-years with snow depth (e.g., Section 3.1) to selected representative stations, which we will identify in the text. Figure 1: Climate station locations, Canadian ecozones, regional zones, agricultural regions, and boreal forest (adapted from [14]). www.videleaf.com

Diurnal Range Definition
The diurnal range of temperature, DTR, is defined as the difference between the maximum temperature, T x , and the minimum temperature, T n : Similarly for relative humidity, RH, (and other variables), we define the diurnal range, DRH, as the difference between the maximum, RH x , and the minimum, RH n : www.videleaf.com DRH = RH x − RH n (1b) In our early papers [9][10][11][12][13] we generally reduced the hourly data to daily means, T m , RH m , and recorded T x , T n , and DTR. The difference in relative humidity (RH), DRH, between T n and T x was used as an approximation of the diurnal range. However, there has been considerable discussion in recent years about the difference between DTR, T x , and T n derived from the monthly means of hourly data, and the conventional monthly mean of daily values of DTR, T x , and T n [12,[17][18][19].
We explored this issue [14], using stratifications by month and by opaque cloud cover, and found systematic biases, especially in winter, and even in summer under cloudy conditions. We concluded that the radiatively-forced diurnal cycle, that is, the lagged response to the diurnally varying radiation field, which is dependent on opaque cloud cover, is represented best by first binning the hourly data for groups of many days, and then by determining the diurnal ranges from the composites. Specifically, we found that this radiatively-forced diurnal cycle has a smaller amplitude than the corresponding average of the daily ranges. The reason is transparent. Without an advection of temperature, T n is near sunrise and T x is in the mid-afternoon, but advection can shift the daily minimum temperature away from the time of sunrise to a lower value than the temperature at sunrise, and similarly, advection can shift the daily maximum temperature away from the mid-afternoon to a higher value than the midafternoon temperature. Either will give a larger diurnal range.
Our dataset has around 240,000 days, so coarse stratifications may have 2000 days in each bin, and detailed sub-stratifications typically have >200 days in each bin. This means that the radiatively-forced diurnal cycle emerges from composites of the hourly data, since the advection of temperature and humidity varies from day to day. This leads to a fundamental quantitative improvement in our understanding of the coupling between the diurnal cycle and the opaque cloud cover that determines the cloud radiative forcing. www.videleaf.com We also derived from T, RH, and surface pressure, PS, the other thermodynamic variables: the mixing ratio (Q), the potential temperature (θ), the equivalent potential temperature (θ E ), and the saturation pressure (p*) at the lifting condensation level (LCL). We defined the pressure height to the LCL, P LCL = PS − p* [2], which in the warm season, is often an indicator of the height of cloud base [9]. We calculated the diurnal ranges that are related to moist convective processes:

Opaque Cloud Bins
Since opaque cloud reflects the solar flux and traps the outgoing long-wave, we used the daily mean of the hourly opaque cloud measurements to stratify the daily mean data, and the diurnal ranges of temperature and humidity and derived variables (see Section 3). In fact, we computed two daily averages from the allsky opaque cloud cover estimates to use for stratification. The first is the simple mean of the 24 hourly values, OPAQ m . The second, OPAQSW, is a mean of the hourly opaque cloud values during daylight hours, weighted by a fit to the downward clear sky flux derived from the reanalysis known as ERA-Interim (details in [13]).

Cloud Radiative Forcing
In the short-wave radiation budget, we can define an effective cloud albedo (ECA) and the short-wave cloud forcing (SWCF) in terms of a downwelling SW clear-sky flux, SWCS dn , based on a fit to the clear-sky fluxes from the nearest grid-point of the reanalysis ERA-Interim [13,20]: The dimensionless ECA, with a range from 0 to 1, is a useful measure of the impact of the reflective cloud field on the surface shortwave radiation budget [2,3]. SWCF becomes increasingly negative as ECA increases, while SWCS dn has a large increase from the winter to the summer solstice.
Similarly, we can define a long-wave cloud forcing (LWCF) in terms of a downwelling clear-sky flux LWCS dn , also from ERA-Interim, as: LWCS dn is the smaller term, and LW dn increases with increasing cloud cover, so that LWCF is positive. The total cloud forcing (CF) of the downwelling radiative fluxes is the sum: In the warm season, the SWCF dominates, and CF is negative. The net cloud forcing can be defined as: where the mean surface albedo, defined as: ranges for Saskatchewan from about 0.18 in summer to 0.73 in winter with snow cover [11,12]. When there is a snow cover, the positive LWCF dominates, because the lower solar elevation and larger surface albedo greatly reduce the net SWCF.
We computed the net LW flux: LW n = LW dn − LW up (8) using observations for LW dn , and estimating LW up from the daily mean air temperature, T m (°C), from: with T k (K) = T m + 273.15, ζ = 5.67 × 10 −8 (W m −2 K −4 ) and the emissivity ε set to 1.

Climate Coupling to Opaque Cloud and Snow Cover
This section will present several topics: the monthly diurnal cycle with and without snow cover, the relationship between snow cover, opaque cloud, and cloud radiative forcing, the climate impact of snow cover, the coupling between opaque cloud and warm season diurnal thermodynamic ranges, and the dependence of the 24 hr imbalances of the diurnal cycle on opaque cloud cover.

Forcing of Diurnal Cycle by Cloud and Snow Cover
We start with the dependence of the monthly diurnal cycle of temperature on cloud and snow cover [14]. Taking the data from all stations-years in Table 1 that have snow depth data, we first stratified by temperature and snow cover: selecting the warm group of days with T m > 0 °C and no snow cover (141,160 days), and the cold group of days with T m < 0°C with surface snow cover (74,260 days). Here, we exclude the much smaller mixed group of days, above freezing with snow cover and below freezing without snow (see [14]). Figure 2 shows the mean diurnal cycle of temperature by month, stratified into 10 bins of daily mean opaque cloud, OPAQ m . In the warm season from May to October, we see a steep increase of maximum temperature T x and diurnal temperature range DTR with decreasing opaque cloud, and a rather small fall of minimum temperature in summer. The changing day-length is clearly visible by September and October. In sharp contrast, in the cold season with snow, from December to February, T x decreases with decreasing opaque cloud, and T n decreases even more steeply to its lowest minimum at sunrise under clear skies.
Beside the September and December plots, we show OPAQ m legends in ascending and descending order to illustrate this reversal of the diurnal cycle coupling to opaque cloud between warm and cold seasons.
For the transition months, November, March, and April between warm and cold seasons, both regimes with and without snow cover are well-represented; it is clear that the distributions are non-overlapping. Note that the temperature range shown for the transition months is broader (32 K) than for the single months with a single regime (21 K). This is a large dataset with about 20,000 days per month, so that each cloud bin has typically about 2000 days in summer and winter. For the transition months, where the data is also split unevenly, the number in each bin varies from about 200 to 1500.
It is clear that snow cover has two large climate impacts. First, it cools the mean climate, represented by T m , by about 10 °C; and second, it reverses the sign of the coupling to opaque cloud cover. Snow cover acts as a climate switch between nonoverlapping regimes [11,14]. We will explore the climate impact of snow cover further in Section 3.4, but first we will show the seasonal impact on the cloud radiative forcing.

Change of Cloud Forcing between Warn and Cold Season
The dramatic differences in the diurnal cycles of temperature shown in Figure 2 are related to the reversal of the sign of the net cloud forcing between the warm season and the cold season with snow cover. We computed this using Equations (3) to (7), and data from the Baseline Surface Radiation Network (BSRN) Prairie site at Bratt's Lake, Saskatchewan at 50.204° N, 104.713° W, elevation 588 m [13]. We have 17 years of the downwelling fluxes, SW dn and LW dn , at Bratt's Lake, which we first averaged from 1-min data to hourly means, and then to daily means. Figure 3 shows that CF net from Equation (6b) reverses the sign from increasing negative with cloud cover in the warm season, to increasing positive in the cold season with cloud cover. This is consistent with the daily mean temperature response seen in Figure 2 to the changing opaque cloud cover. www.videleaf.com

Relationship between Opaque Cloud and Cloud Radiative Forcing
We then binned the BSRN data from Bratt's Lake for the downward SW and LW fluxes using the opaque cloud measurements at Regina, 25 km to the north, simply defining the warm season as days with T m > 0°C and the cold season as days with T m < 0°C , because we have no snow cover data for Bratt's Lake. For the SW comparison, we compared the daytime weighted opaque cloud, OPAQSW (see Section 2.3) with ECA from Equation (3). For the LW comparison, we compared the 24 hr mean OPAQ m with LW n computed from Equation (8).
Figure 4 (left) shows the relationship between ECA and OPAQSW for the warm season above freezing, and the cold season below freezing. ECA increases more steeply with increasing opaque cloud in the warm season than in the cold season. We show the mean and standard error of the binned data, and quadratic regression fits to the daily data, which could be used to convert opaque cloud to ECA. For the warm season, the fit is (R 2 = 0. For the cold season, the fit is (R 2 = 0.71): The uncertainty in ECA on a daily basis is of the order of ±0.08 in the warm season and ±0.11 in the cold season. The standard errors (SE) shown for the climatological fits are much smaller, because they are reduced by the large number of days. Figure 4 (middle) shows the dependence of LW n on opaque cloud for days above freezing (3245 days) for three bins of daily mean RH m (<60, 60-75, >75%). The outgoing LW n flux for the same cloud cover increases as RH falls. The temperature dependence is very small when T m > 0°C (not shown). The right panel shows the dependence of LW n on opaque cloud for temperatures below freezing (2198 days) for three bins of daily mean T m (<−20, −20 to −10, −10 to 0 °C). The outgoing LW n flux now decreases with colder temperatures, probably because the surface cools under a stable BL in the cold season [13].
In the warm season, multiple regression of the daily values of LW n on quadratic opaque cloud and RH m gives (R 2 = 0.91): www.videleaf.com  Figure 2 shows that the impact of snow cover on the Prairies on the diurnal cycle of temperature is very large. The transition months show that the cooling with snow cover is large, and show a reversal of the response to cloud cover, consistent with the reversal of the net cloud forcing between cold and warm seasons shown in Figure 3. This section addresses the resulting mean climate impact of snow cover. Figure 5 shows four different analyses of the climate impact of snow cover. The top-left is a composite of the six climate stations in Saskatchewan for eight days before and after fresh snowfall in November, showing a mean of about 270 snowfall events, with a mean date of November 15 (adapted from [11]). We see the fall of daily mean temperature across the snow event, from near 0 °C a week before, to −9.4 ± 0.7 °C for days 2 to 8 afterwards. The climate transition from fall to winter often comes abruptly with these snow events [11], as the snowpack may not melt till spring. Similar composites for individual stations and the means for other provinces are shown in [11]. All of these suggest that as the albedo of the Prairies changes from about 0.2 with no snow cover to above 0.7 with snow cover [11,12], there is a fall of temperature of nearly 10 °C, and the reverse change occurs in spring with snow melt (see [11]).

Climate Impact of Snow Cover
The top-right (adapted from [14]) shows the fall of mean daily temperature, δT m , with snow cover, derived from Figure 2 by calculating the difference of the diurnal composites with and without snow for the transition months, November and March, for each opaque cloud cover bin. We made a correction of about www.videleaf.com 2 °C, based on the mean seasonal cycle [14], to allow for the fact that the mean date of the snow-free composite is about 15 days earlier in November, and later in March than the composite with snow. The curves are a little noisy, because the independent sampling in opaque cloud bins, with and without snow, is far from homogeneous, and in these transition months, the number of days in each bin ranges widely from 184 to 1869 (not shown). Nonetheless, we see a larger degree of cooling as the opaque cloud decreases. The climate cooling with snow, averaged across all cloud bins (open circles), is −11.8 °C (−10.7 °C) for November (March). We also show quadratic fits (dashed) as a useful smooth reference for the impact of cloud cover. We note that the radiative forcing is stronger in March than November, but we cannot assess whether the small difference between the November and March curves is significant, given the inhomogeneity across the cloud bins.
The bottom-left plot shows the monthly mean temperature across the cold season (black line) and the partition into days with snow cover (blue) and days with no snow cover (red line) for a single station (Lethbridge, Alberta), together with the mean snow depth. The difference between the blue and red curves (the magenta curve) shows the monthly climate cooling of snow cover with a mean value of ΔT = −10.4 ± 0.4 °C. The standard errors shown are small because of the large number of days in the 49-year record. Other stations show similar plots [15], suggesting that the cold season climatology with and without snow (red and blue curves) are distinct and non-overlapping. Conventionally, they are merged to the black curve, so this can be misleading.
The bottom-right panel shows the mean temperature T M , for October to April against the fraction of days with snow cover (FDS) for five stations in Alberta, Lethbridge, Medicine Hat, Calgary, Red Deer, and Grande Prairie, listed in order of increasing latitude (adapted from [11,15]). The line fit shown is for 326 years of data, and we show the station means (black circles) that lie close to this line fit. Since it is clear that the southern three stations (red points) have warmer temperatures and lower FDS than the northern two stations, we also computed the linear regression slopes for these two groups. www.videleaf.com These agree within the uncertainty, which increases for fewer stations. The corresponding plot for Saskatchewan is similar [11]. We conclude that the climate coupling between the fraction of days with snow cover and the mean cool season temperature is a robust feature of the Prairie landscape. The shift of the station means with increasing latitude suggests that reduced insolation is also playing a tightly coupled role. Figure 5 confirms that snow cover has a large cooling impact on the mean temperatures in the cold season: snow cover acts as a climate switch between the two non-overlapping regimes. On daily timescales, the cooling is about −10 °C for the Prairies, where the surface albedo with snow cover is in the order of 0.7. The larger slope of −14.6 °C in fit (12a) for the change of mean cold season temperature with the fraction of days with snow cover suggests that there may be coupling to larger scales that enhance the regional cooling with snow cover. www.videleaf.com

Coupling of Warm Season Diurnal Ranges and 24-h Imbalances to Opaque Cloud
This very large hourly dataset allowed us for the first time to extract the radiatively forced diurnal ranges, shown in Equations (1) and (2), for the key thermodynamic variables [14]. Here, we will just show the warm season; the cold season can be found in [14]. From Figure 2, we extracted DTR as a function of opaque cloud and month, and we extracted DRH, Dθ E , and DP LCL from similar diurnal composites (not shown). Close examination of Figure 2 shows that there is a discontinuity across local midnight that changes with opaque cloud cover. So we calculated, also for the first time, this 24 hr imbalance of the diurnal cycle as a function of opaque cloud and month. These are key conceptual www.videleaf.com improvements in our understanding of the diurnal cycle over land in the warm season, and our results are robust as there are about 20,000 days per month. Figure 6 (top left panel) shows the mean diurnal ranges of temperature, DTR, relative humidity, DRH, and mean daily precipitation for the warm season months April to September with no snow. Remarkably, the diurnal ranges are tightly clustered [9,14], so we also show the 6-month warm season mean. The quadratic regression fits for the dependence of the 6month mean DTR and DRH on OPAQ m are: The leading coefficient is the clear-sky diurnal range, which is a rise of 16.7 °C to the afternoon maximum, coupled to a fall of 47.5% in RH from the morning maximum at sunrise. The cloudy limit for OPAQ m = 1, given by these fits, are the small values (DTR, DRH) = (1.4 °C, 6.0%).
Monthly mean precipitation is very low for OPAQ m < 0.4, and the increase of precipitation with OPAQ m is largest in summer, peaking in July when T and the mixing ratio Q also peak. However, because June has substantially greater opaque cloud cover [12], mean June precipitation (2.28 mm d −1 ) is greater than July (1.91 mm d −1 ). www.videleaf.com   24 , which we calculated from the discontinuities across local midnight [14]. We see that over the range of OPAQ m from 0.05 to 0.95 (nearly clear to nearly opaque cloud cover), the mean (ΔT 24 , ΔRH 24 ) change monotonically from (+2 °C, −6%) to (−1.5 °C, +6%). Under nearly clear skies, the warming, and drying over the diurnal cycle is slightly larger in April, May, and June when the mean temperature is increasing seasonally, and slightly smaller in August and September. Under cloudy skies, there is a larger increase in ΔRH 24 in April and May. The SE of the hourly binned data from which Figure 6 is derived as ≈0.1 K for T, ≤0.5% for RH.
The warming of +2 °C and a drying of −6% over the diurnal cycle under nearly clear skies is about 12% of both the DTR and DRH. The cooling of −1.5 °C and a moistening of +6% under very cloudy skies may be coupled to both the evaporation of rain and downdraft transports. The uniform progression of the diurnal imbalance with increasing cloud is not surprising. However, this means that a steady state diurnal cycle only exists under partly cloudy conditions: for the 6-month mean, ΔT 24 , and ΔRH 24 cross zero for OPAQ m = 0.45. This presents a conceptual challenge for equilibrium models for the non-precipitating convective BL over land [21]. Figure 6 (lower panels) are the corresponding warm season diurnal ranges and 24 hr imbalances for θ E and P LCL . The spread in the diurnal ranges and the diurnal imbalances is again small from April to September. For the 6-month means, the quadratic regression fits for the OPAQ m dependence are: Again, the leading coefficient is the clear-sky diurnal range, which is a rise of (θ E , P LCL ) = (19.7 K, 181.4 hPa) from the morning sunrise minimums. an asymmetric structure that is consistent with ΔRH 24 , since a lower RH is tightly coupled to a higher P LCL . Figure 6 shows the remarkably tight climatological coupling from April to September which links opaque cloud cover, the diurnal ranges and the 24 hr diurnal imbalances, despite substantial differences in the solar zenith angle. Over the diurnal cycle, under nearly clear skies, we see a warming and drying, and a rise of θ E and P LCL . At the cloudy extreme with rain, we see 24 hr imbalances of opposite sign that are generally slightly smaller. These rather precise warm-season patterns across opaque cloud cover, and therefore cloud radiative forcing, set a clear target for modeling the partially cloudy boundary layer over land.

Hydrometeorological Memory on Monthly Timescales
The close coupling between the energy and water cycles at the land surface is central to hydrometeorology, and important to weather forecasts on timescales from days to seasons. An earlier review looked at hydrometeorology using global model reanalysis data [1], which showed how net long-wave and shortwave radiation, cloud cover, surface fluxes, diurnal temperature range, soil moisture, and cloud-base height were coupled on daily timescales over river basins. Reanalysis data contain all the key variables, but historically, the observed near-surface climate variables were temperature and precipitation, along with pressure, wind, relative humidity, and snow-depth. Section 3.5 shows that the warm season diurnal cycle is dominated by the radiative forcing of the opaque cloud cover. But on monthly and longer timescales, soil moisture anomalies are linked to precipitation anomalies, both for the current month and several preceding months.
Here, we summarize some key results from [16], who merged the 12 stations in Table 1 in Alberta and Saskatchewan for the years when precipitation is available. For this monthly analysis, the hourly data were processed as intact monthly mean diurnal cycles for each station for each year. As noted in Section 2.1, the hourly dataset is remarkably complete. Days were omitted if <20 hr of www.videleaf.com data were available. Months were omitted if they had fewer than 28 days remaining, except for February, where this threshold was reduced to 25 days. From the monthly diurnal cycles of T, RH, and PS, we computed the derived thermodynamic variables, Q, θ E , and P LCL , and the diurnal ranges defined in Equations (1) and (2).
For each variable, Y, we extracted from the monthly mean diurnal cycles, the daily mean, Y m , the maximum and minimum, Y x and Y n , and the times of the maximum and minimum [16]. We then computed the long-term station monthly mean, and used these to compute monthly anomalies, δY. For the daily precipitation and snow-depth, we also computed monthly means, the long-term station monthly means, and used these to compute monthly anomalies for each station. The monthly anomalies of opaque cloud, precipitation, snow depth, and snow cover frequency were then standardized by their monthly standard deviation (SD). For the temperature anomalies, δT m , δT x , δT n , and the diurnal temperature range, δDTR, we standardized by the monthly SD of δT m . Similarly for the variables, δRH m , δRH x , δRH n and the diurnal RH range δDRH, we standardized by the monthly SD of δRH m . The corresponding set of anomalies for equivalent potential temperature, δθ E , and pressure-height to the LCL, δP LCL , were standardized by the monthly SD of δθ Em , and δP LCLm respectively.
We used multiple linear regression to explore the correlation between variables. Following [12,16], our starting format was to regress a standardized thermodynamic anomaly, δY, on opaque cloud anomalies (δOPAQ m ) for the current month, and lagged precipitation anomalies for the current month (δPR0) and preceding months (δPR1, δPR2, δPR3, δPR4, δPR5) in the form: δY = A * δOPAQ + B * δPR0 + C * δPR1 + D * δPR2 + E * δPR3 + F * δPR4 + G * δPR5 (15) Multiple regression shows no memory of cloud for previous months. Since we are using anomalies, the leading coefficient is of order zero, so it is not shown. After standardization, all variables are dimensionless. www.videleaf.com

Memory of Cold Season Precipitation in April Climatology
On the Prairies, precipitation memory lasts through winter, as water is stored until the snowpack melts in late March or April. The reflective snow cover on the Prairies, with an albedo ≈0.7, acts as a climate switch that reduces T m by 10 °C (Figures 2 and  5). April is the month when the snowpack finally melts and the ground thaws. The upper group in Table 2 shows selected April standardized anomalies regressed on standardized anomalies of opaque cloud for April; and precipitation from April back to November (coefficients A to G in Equation (15)). We see that the April monthly anomalies show memories of the anomalies of precipitation 5 months back through the entire cold season to November, when typically the ground begins to freeze, and the first lasting snow occurs ( Figure 5). Some of this memory remains in the March snowpack depth (not shown here, see [16]).
For the first row, δOPAQ-Apr, the large negative coefficients for the monthly anomalies δDTR, δT x , and δP LCLx , mean that these variables decrease with increasing opaque cloud cover, while the positive sign for the δRH n and δRH m means that they increase with opaque cloud. For δT x and δDTR (and δT m , not shown), the negative coefficients B to G, for the months March back to November, mean that the positive cold season precipitation anomalies are correlated with cold April temperatures. For δRH n , δRH m (and δRH x , not shown), the positive coefficients, B to G, mean that positive cold season precipitation anomalies are correlated with higher RH in April. Most coefficients for δDTR, δRH n , δRH m , and δP LCLx (representative of afternoon cloud-base) have a 99% confidence (p < 0.01).
There are several physical processes that are probably involved. The precipitation over the cold season is mostly stored in the snowpack till spring, when the melt absorbs energy and cools the surface; the melt also provides water for evaporation, which also cools and increases RH. In addition, the freeze-up of the soil in November may similarly preserve November precipitation anomalies as soil ice through the cold season until spring melt. www.videleaf.com In April, the high albedo of the remaining snowpack, as well as fresh snow, also play a direct climate role, as discussed in Section 3.4 and shown in Figure 2, because snow cover acts as a climate switch. Thus, we computed the standardized April snow cover frequency anomaly from the fraction of days in April with snow depth >0, and added this to the multiple regression (15) to give: δY-Apr = A * δOPAQ m -Apr + B * δPR-Apr + C * δPR-Mar + D * δPR-Feb + E * δPR-Jan + F * δPR-Dec + G * δPR-Nov + S * δSnowCover-Apr The lower group in Table 2 shows the coefficients from Equation (16). There is an increase in R 2 for all variables, and especially for T x , with the addition of snow cover. For maximum temperature, T x , snow cover frequency anomalies have as large an impact as opaque cloud anomalies. Note that the coefficients G for δPR-Nov for δRH n , δRH m , and δP LCLx are not significant, but the coefficients G for δDTR and δT x have a confidence >99% in Table 2. It is possible that this is the cooling impact in April coming from the melt of soil-ice that was frozen back in November. www.videleaf.com  Table 2 shows that the climate in April, when the snow pack finally melts, has memory of precipitation through the entire previous winter up to November. Most of the variability in the April climate is explained by anomalies of winter precipitation and the fraction of days in April with residual snow cover.

Growing Season Memory of Precipitation
After snowmelt on the Prairies, the transition into the growing season months May to August (MJJA) is rapid, and typically, the memory of precipitation for the months May to August only goes back to March or April [16]. Merging the 2466 MJJA growing season months gives a unified description for the growing season correlation of the thermodynamic anomalies with opaque cloud and lagged precipitation, as shown in Table 3 and adapted from [16]. We retain the precipitation anomalies for four months. Table 3 shows that only some anomalies, such as δDTR, δRH n , δRH m , δP LCLx with high R 2 values, are correlated with precipitation anomalies going back three months. As in Table 2, the OPAQ coefficients A are typically the largest, except notably for δQ m . www.videleaf.com The first groups are the regression coefficients for the temperature anomalies, δT x , δT m , δT n , and δDTR, which were all standardized by the SD of δT m . The fit represented by R 2 is largest for DTR, and it decreases from δT x to δT n . All the temperature variable anomalies show a strong inverse correlation with opaque cloud anomalies that reflect the downward SW radiation. The warm season is dominated by negative SWCF as shown in Figure 3. The negative values of A decrease from δT x to δT n . δDTR has a negative correlation to both cloud anomalies, and to the precipitation anomalies going back three months. Note that because all the temperatures were standardized by the SD of δT m , the coefficients for the diurnal range are the difference of the corresponding coefficients for the maximum and minimum. . We see that the coefficients B change sign in the sequence from δT x to δT m to δT n . We also see that T m falls strongly with cloud, but its coupling to precipitation is weak because the coefficients B and C have opposite signs. This regression analysis clearly shows that mean temperature anomalies, δT m , are strongly coupled to cloud, and therefore solar forcing, but rather weakly to precipitation, while δDTR (and δT x ) decrease with both cloud and precipitation. We cannot infer causality from multiple regressions, but negative B for δT x is consistent with evaporation from moist soils reducing T x , and the positive B for δT n is consistent with the fact that under wetter conditions, the fall of T n at night is limited by saturation.
The next group are the four RH anomalies, δRH x , δRH m , δRH n , and δDRH. For the first three, the regression coefficients show that positive RH anomalies are correlated with positive cloud and precipitation anomalies, and the coefficients are significant for both the present and three past months. The coefficients for δDRH are negative because δRH n increases faster with cloud and precipitation than δRH x , and the coefficients are significant for only one past month. The R 2 fit decreases monotonically from the afternoon minimum δRH n to δRH m to the sunrise maximum δRH x to δDRH. The diurnal cycle of T and RH have an inverse dependence on opaque cloud, reaching T x and RH n in the afternoon at the same time [16]. This is related to the fact that www.videleaf.com mixing ratio Q is tightly constrained by BL transports, which we will discuss in Section 4.3. However, over land, near-surface RH is constrained by the availability of soil moisture for evaporation from bare soil and transpiration (which is often modeled as a stomatal resistance to evaporation [22,23]. Soil moisture anomalies are related in turn to precipitation anomalies. We see that afternoon RH n and mean RH m anomalies have a strong positive correlation to precipitation anomalies, and a large R 2 . However, RH x , which increases with precipitation, is limited if the surface saturation is reached and dew forms before sunrise. Because the latent heat release slows the temperature fall, it is consistent that RH x and T n anomalies are both positively coupled to wetter precipitation anomalies for the current month (coefficient B).
The third group in Table 3 is the four P LCL anomalies: P LCLx is generally representative of afternoon cloud-base [9]. P LCL has a strong dependence on RH and a weak dependence on T, and we see that negative P LCL anomalies are coupled to positive cloud and precipitation anomalies. The coefficients are largest for afternoon δP LCLx , for which R 2 is high. The coefficients for δP LCLx , δP LCLm , and δP LCLn are all 99% significant for both the present and three past months, showing that cloud-based anomalies have a long memory of precipitation anomalies in the growing season.
The fourth group in Table 3 shows the coefficients for δθ Ex , δθ Em , δθ En , and δDθ E . The first three show the decrease with increased cloud, but an increase with precipitation. The R 2 values are small, even though the coefficients have 99% confidence. The diurnal range of θ E is dominated by the dependence of DTR on opaque cloud. The two afternoon anomalies, δθ Ex and δP LCLx , are related to moist convective instability, which is favored by a higher θ Ex and a lower cloud base.
The diurnal variation of the mixing ratio, Q, has a double maxima and minima, which we will show in Section 4.3, and so Table 3 shows only the coefficients for the mean anomaly δQ m . The R 2 fit is much smaller for Q than for RH. The negative correlation to opaque cloud is small, because T and RH have an www.videleaf.com inverse diurnal dependence on cloud. The positive correlation to precipitation anomalies goes back two months, consistent with positive precipitation anomalies increasing evapotranspiration. Table 3 summarizes the multiple regression correlation coefficients between warm season near-surface variables and opaque cloud and lagged precipitation, and gives a quantitatively useful target for the evaluation of the coupled processes in models. Two important conceptual results emerge for the monthly mean climate on the Canadian Prairies. Afternoon anomalies of δT x , δRH n , δP LCLx are strongly correlated to opaque cloud anomalies. Correlation with precipitation anomalies are weaker, but stretch back for three past months for these key variables. Anomalies of Q m are coupled to precipitation anomalies with memory of two months past, but they have weak correlations to opaque cloud. Figure 6 showed the very tight coupling in the warm season between opaque cloud and the diurnal range of key thermodynamic variables on daily timescales. Table 3 used multiple linear regression to show the correlation of the monthly anomalies of thermodynamic variables to anomalies of opaque cloud and precipitation. Table 3 confirms the strong correlation with opaque cloud, but shows that the coefficients for the lagged precipitation anomalies differ considerably for different variables.

Growing Season Coupling of the Diurnal Cycle to Precipitation and Cloud
For a graphical representation [16] we approximate by defining a weighted precipitation anomaly δPRwt, based on precipitation for just the current month and the past month: δPRwt = 0.6 * δPR0 + 0.4 * δPR1 (17) This simplification, with this choice of coefficients in the ratio of 1.5, captures much of the precipitation dependence for the variables that have the highest R 2 , such as DTR, RH n , and P LCLx , www.videleaf.com because these have the ratio of the coefficients B/C ≈ 1.5 in Table 3.
The x-axis of Figure 7 is 0.1 bins of OPAQ m = δOPAQ m + 0.46, where 0.46 is the mean opaque cloud over all the months. For each MJJA month (total 2466 months), we computed the weighted anomaly δPRwt from (17), and added the MJJA mean precipitation rate of 1.8 mm d −1 to give PRwt = δPRwt + 1.8. We then stratified the data into three ranges of PRwt of <1.2 mm d −1 ; 1.2 to 2 mm d −1 , and >2 mm d −1 , which have mean values of 0.9, 1.6, and 2.6 mm d −1 . There are (531, 1103, 832) months in these three PRwt bins. To generate Figure 7, we compute for each variable bin, the mean and standard error (SE) of the anomalies, and add back the MJJA variable means. Figure 7 (top-left) shows DTR and its components, T x and T n , the top-right shows DRH, RH x and RH n , the bottom-left shows Dθ E , θ Ex , and θ En and the bottom-right is DP LCL , P LCLx , and P LCLn . The strong dependence on opaque cloud, seen in Figure 6, clearly dominates most of these climate variables, since T falls and RH increases with increasing cloud. This is turn is connected to the weak dependence of Q on cloud ( Table 3). The color scheme is red and blue, respectively, for the dry and wet precipitation bins. As PRwt falls, DTR increases faster than T x . Figure 7 (top-right) shows that RH x and RH n (and RH m , not shown) increase with both cloud and PRwt, but because afternoon RH n increases faster than RH x , DRH decreases with increasing PRwt. Note the rise of RH x with PRwt towards saturation. If RH x reaches saturation at the surface on individual days, condensation of dew and the release of latent heat limit the fall of T n . www.videleaf.com  Figure 7 (bottom panels) show the variables that determine the BL coupling to clouds and precipitation. Afternoon P LCLx and θ Ex determine the cloud-base height and moist adiabat. Both θ Ex and θ En increase with PRwt, but the diurnal range Dθ E depends primarily on cloud instead of precipitation, as shown in Table 3.
All of the P LCL variables decrease with increasing PRwt. The sunrise minimum of P LCLn falls with PRwt, as the surface moves towards saturation. Thus, higher precipitation, which we can loosely associate with increased daytime evapotranspiration (ET), corresponds with a lower monthly mean cloud base and a higher θ E in the afternoon, which would both favor increased convective instability.
In the warm season on the Prairies, the diurnal cycle of mixing ratio Q has two maxima and minima, except under cloudy conditions [9,14,16]. We can graph this dependence on www.videleaf.com anomalies of opaque cloud cover, δOPAQ m , and weighted precipitation anomalies, δPRwt, from Equation (17) in mm d −1 , from the MJJA growing season merge of 2466 months. Figure 8 (left panel) shows the stratification by δOPAQ m into four ranges: δOPAQ m < −0.08; −0.08 to 0; 0 to 0.08, and >0.08, based on the SD of δOPAQ m ≈ 0.08. There are (371, 839, 909, 347) months in these respective bins. We averaged in bins the diurnal cycle of the anomalies from the station monthly means, calculate the SE, and added back the 12-station MJJA mean of Q. The legend shows the mean value for each δOPAQ m bin, and in parentheses the corresponding mean of δPRwt. As the mean δOPAQ m increases from −0.12 to +0.12, the mean δPRwt increases from −0.43 to +0.42 mm d −1 . We have binned by δOPAQ m , but mean OPAQm and precipitation increase together ( Figure 6). The small increase in Q m with δOPAQ m is consistent with Table 3. The sunrise minimum of Q occurs at the minimum temperature, when the night-time BL is shallow with a strong temperature inversion. As the surface net radiation turns positive after sunrise, it drives increasing surface sensible and latent heat fluxes. This warms and moistens a shallow ML trapped beneath the stable nocturnal inversion, and there is a steep rise of Q. When the surface potential temperature reaches that of the top of the capping inversion in mid-morning, the ML deepens more rapidly, www.videleaf.com typically mixing into a deep drier residual ML from the previous day, so that Q falls to the afternoon minimum. With less cloud and more solar forcing, the ML can grow deeper, and mix with more dry air from above, so both the morning rise and mid-day fall of Q are larger in Figure 8. Finally Q rises again to an evening maximum as the surface layer cools and decouples from the deep BL, while ET continues. There are (387, 961, 745, 373) months in these respective bins. The legend shows the mean value for each δPRwt bin, and again in parentheses the corresponding mean of δOPAQ m . With increasing δPRwt, there is a large upward shift of the mean diurnal cycle of Q, as Q m increases with precipitation anomalies, which we can associate with increased soil moisture and ET. As the mean δPRwt increases from −0.99 to +1.23 mm d −1 , mean δOPAQ m increases from −0.05 to +0.05, and the fall of Q from mid-morning maximum to afternoon minimum is reduced as in the left panel.
Clearly, we are dealing with a fully coupled system, but Figure 8 shows that climatologically, while the amplitude of the diurnal cycle of Q increases a little with reduced cloud cover (presumably increased solar forcing and vertical mixing), there is a large upward shift in the diurnal cycle with increased weighted precipitation, presumably from increased ET.

Seasonal Extraction of Surface Total Water Storage
We made a simplified estimate of the growing season water and energy budgets of the Prairies in [12], using fits between opaque cloud and cloud forcing (see Section 3.3) to estimate the surface radiation budget. We calculated the seasonal change in total water storage from the GRACE, using the gridded 1 × 1 degree monthly land mass grids [7,8], version JPL-RL05.DSTvSCS1401, for liquid water equivalent thicknesses for the decade of 2002-2012. This allowed us to link the draw-down www.videleaf.com of total water storage (TWS) [24,25] to precipitation anomalies during the growing season. We will present a brief overview, referring to [12] for discussion of some of the uncertainties. Figure 9 shows the mean annual cycle of the monthly anomalies of TWS from the annual mean, with the small 2002-2012 trend removed. Although the data have been interpolated to a 1 × 1 degree grid for user convenience, the effective spatial resolution of the GRACE data is about 300 km, so we averaged the data over the Prairie regions of the three provinces. The amplitude of the mean annual cycle is 90 mm. We computed ΔTWS:MJJA, the growing season drawdown of TWS for MJJA, as the difference between 31 August (mean of August and September) and 1 May (mean of April and May) for each year and province. Mean drawdown is ΔTWS:MJJA = −79 mm. We regressed the anomalies of ΔTWS:MJJA for the years 2002-2012 onto the corresponding anomalies of MJJA precipitation, www.videleaf.com derived from the monthly archive (see [26]: http://ec.gc.ca/dcchaahccd/) of the second generation adjusted precipitation dataset. The linear regression fit [12] with R 2 = 0.56, is: ΔTWS:MJJA = −0.59(±0.08) + 0.56(±0.09) δPrecip(MJJA) (18) The mean value of Precip(MJJA) for 2002-2012 (details in [12]) is 2.32 mm d −1 , while the mean drawdown of TWS is 0.59 mm d −1 for δPrecip(MJJA) = 0; or 25% of the mean precipitation.
Thus Equation (18) shows that, as δPrecip(MJJA) decreases from +1 to −1 mm d −1 , ΔTWS:MJJA increases from near zero to −1.15 mm d −1 , which corresponds to −141 mm over the 123 day growing season. The coupling coefficient of 0.56(±0.09) in Equation (18) is effectively a 56 ± 9% damping coefficient for precipitation anomalies in the growing season by changes in the drawdown of TWS. This drawdown of stored water in the growing season means that ET > precipitation, and this difference increases in dry summers.

Impact of Land-Use Change on Growing Season Climate
In recent decades, there has been a major change in land use across the Canadian Prairies, as more than five million hectares of summer fallow have been converted to continuous cropping [10]. This large increase in the area of cropland has increased summer transpiration, which in turn have reduced the maximum temperatures in the growing season over the Prairies [27]. Other analyses of US Midwest summer temperature maxima also show a cooling from land-use change to cropland [28] and cropland intensification [29]. Figure 10 summarizes the long-term climate impact of the reduction of summer fallow in Saskatchewan [10,15]. The left panel shows the land-use trends in total cropland, pasture, and summer fallow around five climate stations in Saskatchewan. We generated local averages of the ecodistrict crop data [10] within the 50 km radius circles around each station, as shown in Figure  1. The climate station time-series was split into two periods: a longer historic period, 1954-1991, when summer fallow cover was large (although slowly decreasing), and a recent 20-year period, 1992-2011, when summer fallow has fallen rapidly to its present low value [30]. Figure 10: Long term trends in total cropland, pasture, and summer fallow around five climate stations in Saskatchewan (left); (center) RH m , Q tx and mean precipitation in southern Saskatchewan, and (right) mean changes in annual cycle of T x , P LCLtx , and θ Etx for Saskatoon, Regina, and Estevan (adapted from [10,15]).
The center panel shows the changes between the two time periods. For RH m and mixing ratio Q tx at the time of the afternoon T x . We averaged the 10-day means from Saskatoon, Regina, and Estevan, the three southern stations in Saskatchewan with complete datasets [10]. We show the standard errors (SE) of the difference between the two mean time series as an indication of significance. For precipitation, which has much more variability than temperature and humidity, we used monthly precipitation for the 21 stations in Saskatchewan south of 53.22° N, from the second generation-adjusted precipitation dataset [26]. Error bars are the SE of each monthly mean. The right panel shows the corresponding changes in the annual cycle of T x , P LCLtx , and θ Etx between the two time periods.
It is clear that there are significant changes in the growing season climate between the historic period, 1953-1991, and the more recent period since 1992. The vertical dashed lines mark the period 140 ≤ DOY < 240 (20 May-27 August) considered in [10] to be representative of the crop growing season. Over this time window, the growing season is cooler since 1992, with a drop of T x of −0.93 ± 0.09 K, and significantly moister with a rise of (RH m , Q tx ) of (6.9 ± 0.2%, 0.70 ± 0.04 gkg −1 ). There is a corresponding fall of the P LCLx of 22.3 ± 1.1 hPa, and a small rise of θ Etx of 1.1 ± 0.2 K, both at the time of afternoon T x . There is www.videleaf.com also an increase of summer (June, July, and August) precipitation of 25.9 ± 4.6 mm.
It is reasonable to conclude that we are seeing the fully coupled response to the large shift from summer fallow to intensive cropping. It seems that more intensive agriculture has increased transpiration, which has cooled and moistened the growing season climate, lowering the cloud-base and increasing the equivalent potential temperature. The coupled increase of Q with precipitation is consistent with Figure 8; and the increase of summer precipitation is consistent with the increase in moist instability, represented by lower afternoon P LCLx and higher θ Etx . Not shown here is a distribution shift in cloud frequency: with 6% fewer days with 2-4 tenths, and 7% more with 7-10 tenths cloud cover [10].

Warm Season Atmospheric and Surface Energy Budgets
The Prairie climate data have no surface nor TOA fluxes, so we will use flux data from ERA-Interim to put our analyses in the context of the atmospheric and surface energy budgets for the MJJA warm season. We used gridpoint data co-located with the four climate stations Estevan, Regina, Saskatoon, and Prince Albert in Saskatchewan that were used in [31] to analyze the near-surface biases of temperature in ERA-Interim. We computed the terms in the atmospheric TOA and surface energy budget, first on a daily basis, and then for each MJJA stationyear, dropping a few days in May with surface snow cover. Figure 11 is a simplified graphical representation of the key flux terms for this region of the central Prairies from 49.2 to 53.2° N. www.videleaf.com On the left is the SW budget terms. Part of the incoming TOA solar flux is absorbed by the atmosphere, part reflected by clouds and atmospheric aerosols, and part is reflected by the surface; so that surface SW n = 212 W m −2 is about half the TOA downward flux. On the right is the LW budget terms, where the upward surface LW emission of 400 W m −2 is mostly absorbed by the atmospheric greenhouse gases, primarily water vapor and CO 2 , as well as by clouds. About 10% escapes to space through the atmospheric infrared window. Clouds and atmosphere re-emit to space, and back to the surface, giving a surface LW n = −73 Wm −2 . The resulting surface R n = 139 W m −2 balances the upward sensible and latent heat fluxes (40 and 82 W m −2 respectively), and the growing season warming of the ground of 17 W m −2 , which is probably too large [31]. For the MJJA surface water budget, precipitation is 2 mm d −1 and ET is 2.8 mm d −1 . ET exceeds precipitation in the growing season, because of the substantial drawdown of surface water storage, as discussed in Section 5.1. The TOA SW n > outgoing LW for the MJJA warm season. Table 4 stratifies the 140 station years of data into ECA bins, showing how reflective cloud cover changes the TOA and surface energy budgets. The mean values across all the data in bold are those that are shown as rounded values in Figure 11. The model sign convention is downward fluxes are positive and upward fluxes are negative. The number of station-years in each bin, denoted by K, is not uniform, nor is the distribution across the four stations. Nonetheless, there is a uniform progression as reflective cloud (ECA) increases: TOA SW n , surface SW dn , net radiation, R n , the sensible heat flux, SH, Bowen ratio, BR, and the 2 m temperatures all decrease, while the latent heat flux, LH, and precipitation increase. There is compensation between the surface SW n and LW n fluxes as reflective cloud changes, which reduces the change in R n . If there is less reflection by clouds, SW dn increases, and this warms the surface, and LW up increases. Because the decrease of cloud cover and precipitation are coupled, and ET changes much less than precipitation, we see that the surface BR increases by a factor of two as cloud cover and precipitation decrease. www.videleaf.com Table 4. Surface and TOA fluxes for the MJJA warm season from ERA-Interim for 1979-2013, stratified by ECA, derived from the gridpoints for Estevan, Regina, Saskatoon, and Prince Albert, SK. Mean data in bold are shown in Figure 11. Despite this analysis correction, we found in [31] that on daily timescales, ERA-Interim in the warm season has a cold bias in T x and a warm bias in T n , so that DTR is biased low; and these biases increase under clear skies. Table 4 stratifies four-month composites by ECA, so that it is not directly comparable to the daily analysis in [31], but it is likely that T x and T n in Table 4 have respectively cold and warm biases of about 1 °C.

Conclusions
We have reviewed progress in the quantitative understanding of the coupling between the land surface, clouds, precipitation, snow cover and the climate system that came from analysis of the long-term hourly Canadian Prairie dataset. This progress was possible because, along with conventional hourly measurements, trained observers have recorded opaque cloud fraction hourly across the Prairies for the past 60 years. These 24 daily estimates of opaque cloud are of sufficient quality such that they can be calibrated against BSRN data to obtain the climatology of the daily short-wave, long-wave, and total cloud forcing. This key radiative forcing of the surface energy budget has previously not been available for long-term climate datasets. We found that net cloud radiative forcing changes sign from negative in the warm season to positive in the cold season with snow cover. This transforms the coupling between cloud cover and the diurnal cycle between warm and cold seasons. In the warm season, www.videleaf.com maximum temperature increases with decreasing opaque cloud, while minimum temperature barely changes. In contrast in the cold season with snow cover, maximum temperature falls with decreasing cloud while minimum temperature falls even more steeply. Although our results are stratified by the observed opaque cloud cover, we show the regression Equations (10) and (11) that convert opaque cloud into cloud radiative forcing.
We have shown the many ways in which snow cover acts as a climate switch between two non-overlapping climate regimes, producing a systematic cooling of 10 °C or more. With fresh snow cover in November, the temperature falls 10 °C during the transition. Simply separating cold season days into those with and without snow cover shows that these are two distinct climates, again separated by 10 °C. The traditional merging of these two climates is unrepresentative, since typically the nearsurface boundary layer also changes from unstable to stable with snow cover. As cloud cover falls, the climate cooling with snow increases. Mean cold season temperatures fall by almost 1.5 °C for each 10% increase in the fraction of days with snow cover.
In the warm season with no snow cover, the diurnal ranges of temperature, relative humidity, equivalent potential temperature, and the pressure height of the lifting condensation level all fall steeply with increasing opaque cloud cover. Remarkably, this tight coupling between diurnal range and cloud cover is almost unchanged from April to September, so we provide quadratic fits to the mean profiles. Given 600 station-years of hourly data, we are able to extract, perhaps for the first time, the relationship between cloud forcing and the warm season imbalance of the diurnal cycle. This imbalance over 24 hr changes monotonically from a +2 °C warming and a −6% drying under clear skies, to a −1.5 °C cooling and 6% moistening under cloudy skies with precipitation. Correspondingly, there is a 24 hr rise of θ E and LCL under clear skies, and a fall of both under cloudy skies.
On daily timescales, the radiative forcing dominates, and the changing cloud forcing drives the changing diurnal cycle of the thermodynamic variables. On longer timescales, such as monthly and seasonal, precipitation anomalies, stored as soil moisture anomalies, change the ET and modify the diurnal response. Since we know the cloud radiative forcing which is large on daily timescales, we can show statistically that the memory of water storage anomalies, from precipitation and the snowpack, goes back many months. The spring climatology for April shows the memory of snowfall back through the entire winter, and the memory in summer months goes back to the months of snowmelt.
In addition, we showed how the thermodynamic coupling of the diurnal cycle to the cloud forcing is modified by lagged precipitation anomalies. With reduced precipitation, the diurnal ranges of T, RH, and P LCL , but not θ E , increase. With increased precipitation, afternoon T x falls a little, while RH n and θ Ex increase and P LCLx , representing cloud-base, falls.
Climatologically, the diurnal cycle of Q has a sunrise minimum, a rise to a mid-morning maximum while evaporation is trapped beneath the nocturnal inversion, then a fall to an afternoon minimum, as water vapor is rapidly transported upward into a deep daytime BL; followed by a second rise to an evening minimum as the surface cools and starts to uncouple from the mixed BL. We showed that while the amplitude of the diurnal cycle of Q increases a little with reduced cloud cover (presumably from increased solar forcing and vertical mixing), there is a large upward shift in the diurnal cycle of Q as weighted precipitation increases, presumably from increased evaporation.
The GRACE satellite data from a recent period show that the seasonal extraction of the surface total water storage is a large damping of the interannual variability of precipitation anomalies in the growing season. Over a range of precipitation anomalies of ±1 mm d −1 , the seasonal extraction of ground water increases from near-zero with high precipitation, to 1.15 mm d −1 for low precipitation.
Also on seasonal timescales, the large land-use change on the Prairies from summer fallowing to intensive cropping, with the most rapid transition in the early 1990s, has led to a coupled climate response that has cooled and moistened the growing www.videleaf.com season, lowering cloud-base, increasing equivalent potential temperature, and increasing precipitation.
Finally, we compute the atmospheric energy and water budgets for four stations in Saskatchewan for the 35-y period of the ERA-Interim reanalysis, to show graphically the mean TOA and surface budgets that are representative of the Canadian Prairies for the MJJA growing season. We also show the dependence of these budgets on reflective cloud cover for this fully coupled model system.
More broadly, these long-term climate datasets from the Prairies can be used to improve the representation of clouds and landsurface processes in numerical forecast models. A preliminary step was the assessment of the biases in the surface diurnal cycle of temperature in ERA-Interim [31], but much more is possible. For example, knowledge of the biases is providing meteorological information for regions where no data were previously available. We plan to investigate the impact of these biases on model estimates of crop yields and greenhouse gas emission estimates using the DNDC model [33]. With the improved knowledge of the energy budget for the Prairies, we can now investigate which management practices could significantly change the energy budget and allow farmers opportunities to minimize their impact on climate change.