The GEFS-Based Daily Reference Evapotranspiration (ETo) Forecast and Its Implication for Water Management in the Southeastern United States

Di Tian Department of Agricultural and Biological Engineering, University of Florida, Gainesville, Florida

Search for other papers by Di Tian in
Current site
Google Scholar
PubMed
Close
and
Christopher J. Martinez Department of Agricultural and Biological Engineering, University of Florida, Gainesville, Florida

Search for other papers by Christopher J. Martinez in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

NOAA’s second-generation retrospective forecast (reforecast) dataset was created using the currently operational Global Ensemble Forecast System (GEFS). It has the potential to accurately forecast daily reference evapotranspiration ETo and can be useful for water management. This study was conducted to evaluate daily ETo forecasts using the GEFS reforecasts in the southeastern United States (SEUS) and to incorporate the ETo forecasts into irrigation scheduling to explore the usefulness of the forecasts for water management. ETo was estimated using the Penman–Monteith equation, and ensemble forecasts were downscaled and bias corrected using a forecast analog approach. The overall forecast skill was evaluated using the linear error in probability space skill score, and the forecast in five categories (terciles and 10th and 90th percentiles) was evaluated using the Brier skill score, relative operating characteristic, and reliability diagrams. Irrigation scheduling was evaluated by water deficit WD forecasts, which were determined based on the agricultural reference index for drought (ARID) model driven by the GEFS-based ETo forecasts. All forecast skill was generally positive up to lead day 7 throughout the year, with higher skill in cooler months compared to warmer months. The GEFS reforecast improved ETo forecast skill for all lead days over the SEUS compared to the first-generation reforecast. The WD forecasts driven by the ETo forecasts showed higher accuracy and less uncertainty than the forecasts driven by climatology, indicating their usefulness for irrigation scheduling, hydrological forecasting, and water demand forecasting in the SEUS.

Corresponding author address: Di Tian, Department of Agricultural and Biological Engineering, University of Florida, P.O. Box 110570, Gainesville, FL 32611-0570. E-mail: tiandi@ufl.edu

Abstract

NOAA’s second-generation retrospective forecast (reforecast) dataset was created using the currently operational Global Ensemble Forecast System (GEFS). It has the potential to accurately forecast daily reference evapotranspiration ETo and can be useful for water management. This study was conducted to evaluate daily ETo forecasts using the GEFS reforecasts in the southeastern United States (SEUS) and to incorporate the ETo forecasts into irrigation scheduling to explore the usefulness of the forecasts for water management. ETo was estimated using the Penman–Monteith equation, and ensemble forecasts were downscaled and bias corrected using a forecast analog approach. The overall forecast skill was evaluated using the linear error in probability space skill score, and the forecast in five categories (terciles and 10th and 90th percentiles) was evaluated using the Brier skill score, relative operating characteristic, and reliability diagrams. Irrigation scheduling was evaluated by water deficit WD forecasts, which were determined based on the agricultural reference index for drought (ARID) model driven by the GEFS-based ETo forecasts. All forecast skill was generally positive up to lead day 7 throughout the year, with higher skill in cooler months compared to warmer months. The GEFS reforecast improved ETo forecast skill for all lead days over the SEUS compared to the first-generation reforecast. The WD forecasts driven by the ETo forecasts showed higher accuracy and less uncertainty than the forecasts driven by climatology, indicating their usefulness for irrigation scheduling, hydrological forecasting, and water demand forecasting in the SEUS.

Corresponding author address: Di Tian, Department of Agricultural and Biological Engineering, University of Florida, P.O. Box 110570, Gainesville, FL 32611-0570. E-mail: tiandi@ufl.edu

1. Introduction

Medium-range (1–2 weeks ahead) reference evapotranspiration ETo forecasts can provide valuable information for water resource and irrigation management. Such ETo forecasts can be incorporated into hydrological models, water demand models, and irrigation scheduling models to reduce risk and improve reliability of decision making (e.g., Srivastava et al. 2014). When adequate meteorological data input are provided, ETo can be estimated using the Food and Agriculture Organization (FAO) Irrigation and Drainage Paper 56 (FAO-56) Penman–Monteith (PM) equation, which is considered to be the standard method to compute ETo (Allen et al. 1998).

Daily medium-range ETo forecasts can be made using statistical models or numerical weather prediction (NWP) models. A typical statistical model may use prior ETo values calculated from ground-based observations to make a forecast (e.g., Trajkovic et al. 2003). The main drawbacks of such statistical models are that they are black-box models and cannot be implemented in other locations without local training (e.g., Kumar et al. 2011), and the ground-based observations are often not available (e.g., Srivastava et al. 2013). In recent years, several works have been conducted to predict daily ETo using weather forecast data as inputs to the PM equation. The operational forecast products include the Experimental Forecast Reference Evapotranspiration (FRET) produced by the National Weather Service using outputs of the NWP models (e.g., www.wrh.noaa.gov/forecast/evap/FRET/FRET.php?wfo=sto). Many studies mainly included public weather forecast messages (Cai et al. 2007, 2009) and forecasts from mesoscale NWP models (Ishak et al. 2010; Silva et al. 2010; Srivastava et al. 2013) to predict ETo. However, the availability of the complete public weather forecast messages for the PM equation could not always be guaranteed over a large area at fine resolution (Cai et al. 2007), and while mesoscale NWP models could produce sufficient meteorological data for the PM equation over a large region at a high resolution, it has been noted to be very time consuming and the output suffered from systematic errors (Ishak et al. 2010; Silva et al. 2010).

As an alternative to the above methods, a retrospective forecast (reforecast) archive of a global NWP could be used to forecast daily ETo and correct systematic errors. Archives of the first-generation reforecast of the Global Forecast System (GFS; hereafter reforecast1; Hamill et al. 2006) with a forecast analog technique (Hamill and Whitaker 2006) have been applied to produce skillful ensemble forecasts of daily ETo at high resolution (Tian and Martinez 2012a,b). However, because the reforecast1 did not include sufficient meteorological variables for the PM equation, the forecasting skill was largely impaired because of the need to approximate several variables. The second-generation reforecast dataset (hereafter reforecast2) was recently developed through archiving the forecast outputs of the operational 2012 Global Ensemble Forecast System (GEFS; Hamill et al. 2013). Compared to reforecast1, reforecast2 used a newer version of the weather forecast model and archived a larger set of output data at a higher spatial and temporal resolution (Hamill et al. 2013). Thus, using the forecast analog technique with reforecast2 could potentially improve the daily ETo forecast skill.

In this study, we downscaled and bias corrected medium-range probabilistic forecasts of daily ETo at high spatial resolution using the forecast analog approach with the reforecast2 dataset in the southeastern United States [SEUS; Alabama, Georgia (GA), Florida (FL), North Carolina (NC), and South Carolina]. The ETo was calculated using the PM equation and the probabilistic forecasts of daily ETo were evaluated compared to the climatology of the forcing dataset of phase 2 of the North American Land Data Assimilation System (NLDAS-2) for all months and leads over the SEUS. A water deficit forecast driven by the ETo forecasts was evaluated to explore its usefulness for water management.

2. Data

More than 28 years (from December 1984 to July 2013) of reforecast2 data are currently archived (Hamill et al. 2013), and five variables are particularly useful for ETo estimation. This work used the ensemble mean of the 11-member forecasts of 2-m mean temperature Tmean, maximum temperature Tmax, minimum temperature Tmin, incoming solar radiation Rs, and 10-m wind speed u10 (or Wind) for all 16 lead days at 1.0° × 1.0° (~100 km) over the SEUS (24°–28°N, 90°–74°W; www.esrl.noaa.gov/psd/forecasts/reforecast2/download.html; Fig. 1). The forecasts are archived every 3 h for the first 8 lead days and 6 h from lead day 9 to 16. From lead day 1 to 16, Tmean and Rs were averaged into daily values, daily Tmax and Tmin were obtained by comparing 3- or 6-h Tmax and Tmin on each day, and u10 was estimated from the magnitude of the orthogonal U and V components of 10-m wind speed and then averaged into daily values. Wind speed at 2-m height u2 was estimated from u10 (Allen et al. 1998):
e1
where z is the measurement height (10 m). For this work, we used reforecast2 data from 1 October 1985 to 31 September 2012.
Fig. 1.
Fig. 1.

Study area and an example of a subset of 16 GEFS grid points and corresponding NLDAS-2 grid points in the Tampa Bay area.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

Because long-term continuous daily observed data for ETo estimation were generally not available over the SEUS, we used the forcing dataset of the NLDAS-2 (Xia et al. 2012a,b) as a surrogate for observations. The forcing data of NLDAS-2 has a resolution of 0.125° × 0.125° (~12 km; Fig. 1) over North America. We obtained relevant variables at an hourly time step (temperature, Rs, and u10) over the SEUS (25°–37°N, 89°–75°W). Daily Tmax and Tmin were obtained by comparing 1-h values on each day, daily Tmean was the average of daily Tmax and Tmin, and daily Rs and u10 were converted by averaging hourly values into daily values. The u10 was converted to u2 using Eq. (1).

Observations from two weather stations in the SEUS (Fig. 2) were selected to compare with the NLDAS-2 ETo, Tmax, Tmin, Rs, and Wind extracted at the same locations. These two stations included Lake Alfred, FL (28.10°N, 81.71°W), and Raleigh, NC (38.88°N, 78.79°W). The Lake Alfred station data was from the Florida Automated Weather Network, and the Clayton station (identification CLAY) data was from the North Carolina Climate Retrieval and Observations Network of the Southeast Database (NC CRONOS). Both stations are daily observations covering the period from 1 January 1998 to 31 December 2012. Both of the Lake Alfred and Clayton stations included Tmax, Tmin, dewpoint temperature Tdew, Rs, Wind, and ETo, where ETo was calculated by the FAO-56 PM equation.

Fig. 2.
Fig. 2.

Locations of the observed stations and selected locations for WD forecast.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

3. Methods

a. ETo calculation

The FAO-56 PM equation (Allen et al. 1998) was used to estimate daily ETo (mm day−1):
e2
where Rn is the net radiation at the crop surface (MJ m−2 day−1), G is soil heat flux density (MJ m−2 day−1) and was considered to be negligible for the calculations, Tmean is the mean air temperature at 2-m height (°C), u2 is wind speed at 2-m height (m s−1), es is saturation vapor pressure (kPa), ea is actual vapor pressure (kPa), Δ is the slope of the saturation vapor pressure curve (kPa °C−1), and γ is the psychrometric constant (kPa °C−1). Details of the calculation of each of the terms in Eq. (2) can be found in Allen et al. (1998). Allen et al. (1998) recommended using the following equation to approximated Tdew using Tmin:
e3
where Ko = 0°C for humid and subhumid climates and Ko = 2°C for arid and semiarid climates. Because the reforecast2 and forcing data of NLDAS-2 did not include dewpoint temperature, which are required to calculate ea, we approximate Tdew using Eq. (3) with Ko = 0°C. It is worthwhile to note that by assuming the mixing ratio is approximately equal to the specific humidity (SH; kg kg−1), ea can also be estimated with SH and surface pressure P (kPa) fields available in the reforecast2 and forcing data of NLDAS-2 using Eq. (4):
e4
To check the approximation of ea using Tmin, we compared the daily ea(SH) (ea calculated using NLDAS-2 daily SH and P) with the daily ea(Tmin) (ea approximated using NLDAS-2 daily Tmin) from 1985 to 2012 at two locations, Saint Leo, FL, and Smithfield, NC, in 12 months (Fig. 2). Figure 3 shows that ea was slightly over estimated by Tmin in June, July, and August, but ea(Tmin) matched quite well with ea(SH) in other months, which indicated that the approximation using Tmin was generally valid.
Fig. 3.
Fig. 3.

Box-and-whisker plots of the NLDAS-2 ea(Tmin) (red) and the NLDAS-2 ea(SH) (blue) in 12 months at (a) Saint Leo, FL, and (b) Smithfield, NC.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

b. The analog method

The fine-resolution ETo forecasts were produced using a moving spatial window forecast analog approach (Hamill and Whitaker 2006; Hamill et al. 2006). This method is briefly described below. The first step was to find the closest analog dates for one lead day forecast within a domain of 16 grid points of GEFS (Fig. 1). That is, the root-mean-square error (RMSE) over the 16 grid points of GEFS was computed by comparing a given day’s ETo forecast with the historical ETo forecasts within a ±45-day window of the same day of the year at the same lead. For example, a 15 February 2012 forecast at lead day 5 over central Florida (the 16 GEFS grid points) was compared against the reforecasts from 1 January to 1 April between 1985 and 2011 (i.e., 2457 forecasts) at the same lead day. The RMSE between the current forecast and each reforecast was computed and averaged over the 16 grid points in Fig. 1. The analog dates for the day’s ETo forecast were selected as the 75 dates with the lowest RMSE. In the second step, the analogs for the forecast of the day were selected on the dates of the 75 analogs at each of the grid points of the NLDAS-2 dataset in Fig. 1. Thus, 75 analogs were chosen for this forecast to generate an ensemble forecast for each grid point of the NLDAS-2. This process was repeated for the other forecast lead days and the other locations across the SEUS, and the forecast over the entire SEUS were generated by tiling together the local analog forecasts. Leave-one-out cross validation was conducted by excluding the current year from the list of potential analogs.

c. Forecast verification

The overall skill for the probabilistic forecasts was evaluated using the modified linear error in probability space (LEPS) skill score (Potts et al. 1996). The modified LEPS measures the distance between the forecast and the observation in the cumulative distribution space and the LEPS score S is defined as
e5
where Pf and Po are the probabilities of the forecast and observation from their respective empirical cumulative distributions. The value of S can be expressed as a skill score with a range from −100% to 100%, where a skill score of 100% indicates a perfect forecast and a skill of 0% indicates a skill score equivalent to climatology (Piechota et al. 2001; Potts et al. 1996), where the climatology was an ensemble of a running ±30-day window. The LEPS skill score was evaluated at each grid point for each month and lead day.
The Brier skill score (BSS) was used to evaluate the skill of probabilistic forecasts for five categories, including terciles (upper, middle, and lower) and extremes (10th and 90th percentiles). The BSS is written as (Wilks 2011)
e6
where BSf is the Brier score of the forecast and BSc is the Brier score of climatology, which was computed from daily values ±30 days of the forecast date. The BSS ranges between −∞ and 1.0, and values of the BSS of 1.0 indicate a perfect forecast and values of 0 indicate that the skill of the forecast is equivalent to climatology. The BSS was evaluated at each grid point for each month and lead day.
Reliability diagrams (Sobash et al. 2011; Wilks 2011) and relative operating characteristic (ROC) diagrams (Wilks 2011) were used to evaluate the resolution and reliability across the study area and period. The reliability diagram indicates the degree that forecast probabilities match observed frequencies. An overall measure of the reliability of the forecasts can be assessed by the deviation of the reliability curve from the diagonal. For a perfectly reliable forecast system, the reliability curve is aligned along the diagonal. Curves below (above) the diagonal indicate overforecasting (underforecasting), and the nearer the curve to horizontal, the less reliable the forecast. The ROC diagram compares hit rates H to false alarm rates F at different forecast probability levels and is a measure of how well the probabilistic forecast discriminates between events and nonevents. The H and F are defined as follows:
e7
e8
where CE is the correctly forecasted events, FE is the falsely forecasted events, FN is the falsely forecasted nonevents, and CN is the correctly forecasted nonevents. An ROC curve that lies along the 1:1 line indicates no resolution, and a curve that is far toward the upper-left corner indicates high resolution. The ROC and reliability diagrams were evaluated by pooling all grid points and months for each lead day.

d. Irrigation scheduling model

A simple model based on the agricultural reference index for drought (ARID; Woli et al. 2012, 2013; Khare et al. 2013) was used for irrigation scheduling. The ARID is a generic model that uses a budget-based simple soil water balance for a reference grass having a 400-mm soil layer with evenly distributed roots. The irrigation scheduling was made in terms of the water deficit WD (mm) component of the ARID model, which is defined as
e9
where ET is actual evapotranspiration (mm) and ETo is the reference evapotranspiration (mm). The water balance components include precipitation, surface runoff, deep drainage, and evapotranspiration. The soil water content on the ith day is calculated from the soil water content on the previous day and the day’s inputs and outputs:
e10
where Wi is soil water (mm) on the ith day, Wi−1 is soil water (mm) on the i − 1th day, Pi is precipitation (mm) on the ith day, Di is deep drainage (mm) on the ith day, and Ri is surface runoff (mm) on the ith day. The ETo is estimated using the PM equation. The ET is estimated based on the water uptake model developed by Meinke et al. (1993):
e11
where θ is the volumetric soil moisture content (m3 m−3), θwp is the volumetric soil moisture content at the wilting point (m3 m−3), MUF is the maximum uptake factor (day−1), and Z is the root zone depth. The surface runoff was calculated using the U.S. Department of Agriculture (USDA) Natural Resources Conservation Service (NRCS) curve number method (NRCS 1986):
e12
where RCN is the runoff curve number and is dependent on soil type and land cover, I is the initial abstraction (mm), and S is the potential maximum retention (mm). In ARID, deep drainage is calculated as
e13
where θfc is the volumetric soil content at field capacity (m3 m−3) and DC is a drainage coefficient (m3 m−3 day−1).

Because ARID assumes a reference grass having a 400-mm soil layer with evenly distributed roots, it only requires site-specific weather data as a site-specific input. Following Woli et al. (2012), default values for the other parameters were used in this study (DC = 0.55 m3 m−3 day−1, Z = 400 mm, RCN = 65, MUF = 0.096 day−1, θwp = 0.06 m3 m−3, and θfc = 0.12 m3 m−3). The 27 years of daily WDs (from 1 October 1985 to 31 September 2012) at five locations in the SEUS (Fig. 2), including Smithfield, NC (35.52°N, 78.35°W); Gainesville, GA (34.30°N, 83.86°W); Tifton, GA (31.45°N, 83.48°W); Saint Leo, FL (28.34°N, 82.26°W); and Perrine, FL (25.58°N, 80.44°W), were calculated using the ARID model with the input from the daily NLDAS-2 precipitation (aggregated from hourly precipitation; mm day−1) and ETo (mm day−1).

e. Evaluation of analog forecast of WD

The ETo forecast analogs and NLDAS-2 precipitation were used as inputs into the ARID model to generate the WD analog forecast. As a comparison, the WD climatology was generated using the NLDAS-2 ETo climatology (ensemble of a running ±30-day window) and NLDAS-2 precipitation. The calculated NLDAS-2 WDs were used to verify the WD forecasts. The ensemble forecast of both analog and climatology were evaluated using the p factor and r factor (Yang et al. 2008). The p factor is the percent of observations within the 2.5% and 97.5% percentiles of the ensemble forecast. The r factor is the relative width of the 2.5% and 97.5% percentiles of the ensemble forecast. It is calculated as
e14
where n is the number of days, and represent the 2.5% and 97.5% percentiles of the ensemble forecast at day ti, and σobs is the standard deviation of the observations. The goodness of the ensemble forecast is judged on the basis of the closeness of the p factor to 100% (i.e., all observations bracketed by the ensemble forecast) and the r factor to 1 (i.e., achievement of rather small ensemble forecast range).
The deterministic WD forecasts were defined as the mean of the analog forecast or climatology ensemble, respectively. The Nash–Sutcliffe coefficient E and root-mean-square error were employed to evaluate both the deterministic climatology forecast and the deterministic analog forecast. The Nash–Sutcliffe coefficient is written as (Nash and Sutcliffe 1970)
e15
where n is the number of days; and represent the observation and forecast at day ti, respectively; and is the average value of the observations over the entire period. The E ranges from −∞ to 1, with a value of 1 corresponding to a perfect match of the forecast to the observations and a value of less than 0 representing that the forecast is worse than the mean of observations.

4. Results

a. Comparison between observations and NLDAS-2

Figure 4 shows a comparison of the observed and NLDAS-2 daily ETo, Tmax, Tmin, Rs, and Wind for 12 months at Lake Alfred, FL, and Clayton, NC, in the SEUS. While the NLDAS-2 ETo showed greater positive bias at the Lake Alfred station than the Clayton station throughout the year, the Tmax and Tmin matched quite well with the observations at both stations. The Rs showed higher bias at the Lake Alfred station than the Clayton station throughout the year. The wind also showed slightly smaller positive bias at the Clayton station than the Lake Alfred station across all months. Because ETo was estimated by the FAO-56 PM equation with the inputs of Tmax, Tmin, Rs, and Wind, any bias of ETo was due to one or more of these variables. At the Lake Alfred station in Florida, the high bias of ETo was likely caused by the high bias of both Rs and Wind. The ETo bias was relatively lower in the Clayton station than the bias at the Lake Alfred station, which was likely due to only Wind showing positive bias at this station.

Fig. 4.
Fig. 4.

Box-and-whisker plots of NLDAS-2 (red) and observations (blue) for ETo, Tmax, Tmin, Rs, and Wind in 12 months at (a) Alfred Lake, FL, and (b) Clayton, NC.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

b. Skill of the analog ETo forecast

Figure 5 shows the skill scores as a function of month and lead day. For lead day 1, the lower and upper tercile forecast skills (BSS) were above 0.4 in cool months and above 0.2 in warm months (Figs. 5b,d); the lower and upper extreme forecast skills (BSS) were between 0.2 and 0.4 throughout the year, with the skill higher in cool months than in warm months (Figs. 5a,e); the middle tercile forecast skills (BSS) were the lowest among the five categories, with the skill scores above 0 and below 0.2 (Fig. 5c); and the overall forecast skill (LEPS skill score) was above 40% in cool months and above 20% in warm months. The forecast skills for tercile, extreme, and overall forecasts were generally positive up to lead day 7, with the skill scores higher in cool months than in warm months and the skill decreased with the increase in lead day.

Fig. 5.
Fig. 5.

BSS and LEPS skill score (%) as a function of month and lead day. (a)–(e) The BSS reflected the extreme and tercile forecasting skill. (f) The LEPS skill score reflected the overall forecasting skill.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

Figure 6 shows the reliability and ROC diagrams for tercile and extreme forecasts for lead days 1 and 5. While the reliability curves showed slight overforecasting bias, the curves were along the diagonal for all five categorical forecasts at lead days 1 and 5 (Figs. 6a,b), which indicated the forecast probability matched very well with the observed relative frequency. All ROC curves with the exception of the middle tercile forecasts were far toward the upper-left corner for both lead days 1 and 5 (Figs. 6c,d), indicating the upper and lower tercile and extreme forecasts had very good resolution and could discriminate between events and nonevents very well. While the reliability and ROC curves all showed high skill at lead day 5, the BSS only showed modest skill. Our explanation is that the reliability and ROC curves were created by pooling the forecasts and observations of all grid points and months instead of evaluating at each grid point and month before averaging, so that the variability of ETo across locations and seasons produced artificially high skill (Hamill and Juras 2006).

Fig. 6.
Fig. 6.

(a),(b) Reliability and (c),(d) ROC diagrams for tercile and extreme forecasts at lead days 1 and 5.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

The skill scores for each grid point were calculated by averaging over all months for different lead days. Figure 7 shows the overall, tercile, and extreme forecast skill for each grid point over the SEUS for four lead days. The skill decreased with the increase of lead days. All forecasts indicated the Florida peninsula showed the lowest skill, and the northern portion of the study area generally showed higher skill than the southern.

Fig. 7.
Fig. 7.

(top) LEPS skill score (%) and (below) BSS over the SEUS. The LEPS skill score reflected the overall forecasting skill. The BSS reflected the tercile and extreme forecasting skill.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

c. Evaluation of WD forecast

Table 1 shows the evaluation of the analog forecasts and climatology of WDs (mm) at lead days 1 and 5 at five locations in the SEUS. The analog forecast showed better E, RMSE, and r factor than climatology at all the five locations, with the stations in Georgia and North Carolina showing better performance than the two stations in Florida for both lead days 1 and 5. In terms of the E and RMSE, which evaluates the deterministic forecast, the analog-based WD forecast is more accurate than the climatology-based WD forecast. There is a trade-off between the p factor and r factor, that is, the percent of observations falling within the upper and lower boundaries of the ensemble and the width of the boundaries of the ensemble. While the climatology ensemble only shows a slightly higher p factor than the analog ensemble forecast, the analog ensemble forecast shows much lower r factor than the climatology ensemble. This indicates the analog ensemble forecast has a much narrower boundary than the climatology ensemble but includes a similar percent of observations as the climatology ensemble, that is, the analog ensemble forecast achieved a better trade-off than the climatology ensemble.

Table 1.

Evaluation of ensemble and deterministic forecast of water deficits (mm) at five locations in the SEUS.

Table 1.

5. Concluding remarks

In this work, daily medium-range probabilistic ETo forecasts at high resolution were produced using the analog approach with the reforecast2 dataset. The daily ETo probabilistic forecast showed skill up to lead day 7. The WD forecasts were generated using the daily ETo analog forecast at five locations in the SEUS. The analog-based WD forecast showed high accuracy and small uncertainty at lead days 1 and 5. Because the WD is considered to be the amount of water needed for irrigation, this WD forecast based on the ETo analog forecast will be useful for irrigation scheduling.

The daily ETo forecast showed higher skill in cool months than in warm months. The lower skill in warm months than in cool months was likely due to convective heating. Convective heating is more frequent in summertime than in wintertime. Such convection could cause different weather conditions such as cloudiness and rainfall at small scales, which may not have been well captured by the GEFS model resolution. Over space, the skill was higher in the north than in the south of the SEUS, with the Florida peninsula showing the lowest skill. This is probably due to the influence of the sea breeze over the Florida peninsula (Marshall et al. 2004; Misra et al. 2011).

Compared to the reforecast1 ETo using either natural analogs of the ensemble-mean forecast (reforecast1a; Tian and Martinez 2012a) or constructed analogs of all forecast members (reforecast1b; Tian and Martinez 2012b), the skills of reforecast2 were greatly improved throughout the year (Fig. 8). For example, the overall forecast skill (LEPS skill score) for lead day 1 in January increased from 28% (reforecast1a) and 37% (reforecast1b) to 43%. Moreover, the ETo forecasts using reforecast1a could only produce skillful forecasts up to lead day 5. In contrast to the reforecast1a and reforecast1b ETo, the reforecast2 ETo did not show the lowest skill in the mountainous regions at the confluence of North Carolina, South Carolina, and Georgia, and the skill of reforecast2 was higher than reforecast1 at all locations over the SEUS (Fig. 8). In moving from the reforecast1 data used in Tian and Martinez (2012a) to the reforecast2 data here, the improvement of ETo forecast skill is due to various reasons. First, the reforecast2 dataset was generated by the operational version GEFS (Hamill et al. 2011), which has more advanced configurations with higher horizontal and vertical resolution and more advanced methods for data assimilation and ensemble initialization than the 1998 version GFS used by reforecast1. Second, reforecast2 provided more of the required variables for the PM equation, including Tmax, Tmin, and Rs, so that only one approximation (using Tmin to replace Tdew) was made in this study. By contrast, reforecast1 had a number of variables (Tmax, Tmin, Rs, Tdew, and RH) that needed to be approximated in order to use the PM equation (Tian and Martinez 2012a). Third, reforecast2 (~100 km) has a much higher spatial resolution than reforecast1 (~250 km). Therefore, the analogs could be selected from a smaller spatial domain (reforecast2), which was more likely to find appropriate analogs compared to a larger spatial domain (reforecast1).

Fig. 8.
Fig. 8.

Comparisons of LEPS skill scores (%) for reforecast2, reforecast1a, and reforecast2b at lead days (a) 1 and (b) 5 throughout the year. (c) Comparisons of LEPS skill scores (%) for (left) reforecast2, (middle) reforecast1a, and (right) reforecast1b at lead days (top) 1 and (bottom) 5 over the SEUS.

Citation: Journal of Hydrometeorology 15, 3; 10.1175/JHM-D-13-0119.1

For the forecast verification, the use of the forcing dataset of NLDAS-2 instead of the North American Regional Reanalysis (NARR) as used by Tian and Martinez (2012a,b) could also affect the forecast. The NLDAS-2 forcing data were based on the interpolation of the NARR fields. The Rs was bias corrected from the interpolated NARR using satellite-derived Rs (Xia et al. 2012a) over each grid cell using the ratio of their monthly average diurnal cycle. The details of the method for bias correction can be found in Berg et al. (2003) and Pinker et al. (2003). The NLDAS-2 2-m temperature was directly interpolated from NARR with adjustments to account for the vertical difference between the NARR and NLDAS-2 fields of terrain height (Xia et al. 2012b). Nevertheless, the downscaled and bias-corrected ETo forecasts produced by this work are still affected by the biases in the NLDAS-2 fields. By comparing the NLDAS-2 fields with station-based observations, we can see biases for ETo, Rs, and Wind in Florida. Because the ETo was calculated by the PM equation, the biases of Rs and Wind would be the major contribution to the biases of ETo. The ETo has a smaller bias in North Carolina than in Florida, which was due to the smaller bias of Rs and Wind. While the forcing dataset of NLDAS-2 showed bias, we should note that the skill of the forecast is not related to the bias of the NLDAS-2 data, because the NLDAS-2 data were used for both validation and the selection of forecast analogs.

It is worth noting that further research could be conducted to find the optimum number of analogs and size of the search window to improve the ETo forecast skill (e.g., Hamill et al. 2006). Also, ETo forecasts could be verified using ground-based observations rather than the forcing data of NLDAS-2. However, such ground-based observations would need to be of a similar length of the reforecast2 dataset in order for forecast analog methods to be effective. In this work, Tmin was used to approximate Tdew, which was found to be generally valid for the humid SEUS. This approximation should be validated before being applied to other, less humid locations (Allen et al. 1998).

Acknowledgments

This research was supported by the NOAA/Climate Program Office SARP and RISA program Grants NA10OAR4310171 and NA12OAR4310130. The second-generation NOAA Global Ensemble Forecast System Reforecast was provided by the NOAA/ESRL/PSD, Boulder, Colorado (www.esrl.noaa.gov/psd/forecasts/reforecast2/). The forcing data of NLDAS-2 used in this effort were acquired as part of the activities of NASA’s Science Mission Directorate and are archived and distributed by the Goddard Earth Sciences (GES) Data and Information Services Center (DISC). The authors thank three anonymous reviewers for their valuable comments.

REFERENCES

  • Allen, R. G., Pereira L. S. , Raes D. , and Smith M. , 1998: Crop evapotranspiration: Guidelines for computing crop water requirements. FAO Irrigation and Drainage Paper 56, 300 pp. [Available online at www.fao.org/docrep/X0490E/X0490E00.htm.]

  • Berg, A. A., Famiglietti J. S. , Walker J. P. , and Houser P. R. , 2003: Impact of bias correction to reanalysis products on simulations of North American soil moisture and hydrological fluxes. J. Geophys. Res., 108, 4490, doi:10.1029/2002JD003334.

    • Search Google Scholar
    • Export Citation
  • Cai, J., Liu Y. , Lei T. , and Pereira L. S. , 2007: Estimating reference evapotranspiration with the FAO Penman–Monteith equation using daily weather forecast messages. Agric. For. Meteor., 145, 2235, doi:10.1016/j.agrformet.2007.04.012.

    • Search Google Scholar
    • Export Citation
  • Cai, J., Liu Y. , Xu D. , Paredes P. , and Pereira L. S. , 2009: Simulation of the soil water balance of wheat using daily weather forecast messages to estimate the reference evapotranspiration. Hydrol. Earth Syst. Sci., 13, 10451059, doi:10.5194/hess-13-1045-2009.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., and Juras J. , 2006: Measuring forecast skill: Is it real skill or is it the varying climatology? Quart. J. Roy. Meteor. Soc., 132, 29052923, doi:10.1256/qj.06.25.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., and Whitaker J. S. , 2006: Probabilistic quantitative precipitation forecasts based on reforecast analogs: Theory and application. Mon. Wea. Rev., 134, 32093229, doi:10.1175/MWR3237.1.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Whitaker J. S. , and Mullen S. L. , 2006: Reforecasts: An important dataset for improving weather predictions. Bull. Amer. Meteor. Soc., 87, 3346, doi:10.1175/BAMS-87-1-33.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Whitaker J. S. , Fiorino M. , and Benjamin S. G. , 2011: Global ensemble predictions of 2009’s tropical cyclones initialized with an ensemble Kalman filter. Mon. Wea. Rev., 139, 668688, doi:10.1175/2010MWR3456.1.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Bates G. T. , Whitaker J. S. , Murray D. R. , Fiorino M. , Galarneau T. J. , Zhu Y. , and Lapenta W. , 2013: NOAA’s second-generation global medium-range ensemble reforecast dataset. Bull. Amer. Meteor. Soc., 94, 15531565, doi:10.1175/BAMS-D-12-00014.1.

    • Search Google Scholar
    • Export Citation
  • Ishak, A. M., Bray M. , Remesan R. , and Han D. W. , 2010: Estimating reference evapotranspiration using numerical weather modelling. Hydrol. Processes, 24, 34903509, doi:10.1002/hyp.7770.

    • Search Google Scholar
    • Export Citation
  • Khare, Y. P., Martinez C. J. , and Muñoz-Carpena R. , 2013: Parameter variability and drought models: A study using the Agricultural Reference Index for Drought (ARID). Agron. J., 105, 14171432, doi:10.2134/agronj2013.0167.

    • Search Google Scholar
    • Export Citation
  • Kumar, M., Raghuwanshi N. S. , and Singh R. , 2011: Artificial neural networks approach in evapotranspiration modeling: A review. Irrig. Sci., 29, 1125, doi:10.1007/s00271-010-0230-8.

    • Search Google Scholar
    • Export Citation
  • Marshall, C. H., Pielke R. A. Sr., Steyaert L. T. , and Willard D. A. , 2004: The impact of anthropogenic land-cover change on the Florida peninsula sea breezes and warm season sensible weather. Mon. Wea. Rev., 132, 2852, doi:10.1175/1520-0493(2004)132<0028:TIOALC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Meinke, H., Hammer G. L. , and Want P. , 1993: Potential soil water extraction by sunflower on a range of soils. Field Crops Res., 32, 5981, doi:10.1016/0378-4290(93)90021-E.

    • Search Google Scholar
    • Export Citation
  • Misra, V., Moeller L. , Stefanova L. , Chan S. , O’Brien J. J. , Smith T. J. III, and Plant N. , 2011: The influence of the Atlantic Warm Pool on the Florida panhandle sea breeze. J. Geophys. Res., 116, D00Q06, doi:10.1029/2010JD015367.

    • Search Google Scholar
    • Export Citation
  • Nash, J. E., and Sutcliffe J. V. , 1970: River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol., 10, 282290, doi:10.1016/0022-1694(70)90255-6.

    • Search Google Scholar
    • Export Citation
  • NRCS, 1986: Urban hydrology for small watersheds. Tech. Release 55, 2nd ed., 164 pp. [Available online at ftp://ftp.wcc.nrcs.usda.gov/wntsc/H&H/other/TR55documentation.pdf.]

  • Piechota, T., Chiew F. , Dracup J. , and McMahon T. , 2001: Development of exceedance probability streamflow forecast. J. Hydrol. Eng., 6, 2028, doi:10.1061/(ASCE)1084-0699(2001)6:1(20).

    • Search Google Scholar
    • Export Citation
  • Pinker, R. T., and Coauthors, 2003: Surface radiation budgets in support of the GEWEX Continental-Scale International Project (GCIP) and the GEWEX Americas Prediction Project (GAPP), including the North American Land Data Assimilation System (NLDAS) project. J. Geophys. Res., 108, 8844, doi:10.1029/2002JD003301.

    • Search Google Scholar
    • Export Citation
  • Potts, J. M., Folland C. K. , Jolliffe I. T. , and Sexton D. , 1996: Revised “LEPS” scores for assessing climate model simulations and long-range forecasts. J. Climate, 9, 3453, doi:10.1175/1520-0442(1996)009<0034:RSFACM>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Silva, D., Meza F. J. , and Varas E. , 2010: Estimating reference evapotranspiration (ETo) using numerical weather forecast data in central Chile. J. Hydrol., 382, 6471, doi:10.1016/j.jhydrol.2009.12.018.

    • Search Google Scholar
    • Export Citation
  • Sobash, R. A., Kain J. S. , Bright D. R. , Dean A. R. , Coniglio M. C. , and Weiss S. J. , 2011: Probabilistic forecast guidance for severe thunderstorms based on the identification of extreme phenomena in convection-allowing model forecasts. Wea. Forecasting, 26, 714728, doi:10.1175/WAF-D-10-05046.1.

    • Search Google Scholar
    • Export Citation
  • Srivastava, P. K., Han D. , Ramirez M. A. R. , and Islam T. , 2013: Comparative assessment of evapotranspiration derived from NCEP and ECMWF global datasets through Weather Research and Forecasting model. Atmos. Sci. Lett., 14, 118125, doi:10.1002/asl2.427.

    • Search Google Scholar
    • Export Citation
  • Srivastava, P. K., Han D. , Ramirez M. A. R. , and Islam T. , 2014: Sensitivity and uncertainty analysis of mesoscale model downscaled hydro-meteorological variables for discharge prediction. Hydrol. Processes, in press.

    • Search Google Scholar
    • Export Citation
  • Tian, D., and Martinez C. J. , 2012a: Forecasting reference evapotranspiration using retrospective forecast analogs in the southeastern United States. J. Hydrometeor., 13, 18741892, doi:10.1175/JHM-D-12-037.1.

    • Search Google Scholar
    • Export Citation
  • Tian, D., and Martinez C. J. , 2012b: Comparison of two analog-based downscaling methods for regional reference evapotranspiration forecasts. J. Hydrol., 475, 350364, doi:10.1016/j.jhydrol.2012.10.009.

    • Search Google Scholar
    • Export Citation
  • Trajkovic, S., Todorovic B. , and Stankovic M. , 2003: Forecasting of reference evapotranspiration by artificial neural networks. J. Irrig. Drain. Eng., 129, 454457, doi:10.1061/(ASCE)0733-9437(2003)129:6(454).

    • Search Google Scholar
    • Export Citation
  • Wilks, D. S., 2011: Statistical Methods in the Atmospheric Sciences.Academic Press, 476 pp.

  • Woli, P., Jones J. W. , Ingram K. T. , and Fraisse C. W. , 2012: An agricultural reference index for drought. Agron. J., 104, 287300, doi:10.2134/agronj2011.0286.

    • Search Google Scholar
    • Export Citation
  • Woli, P., Jones J. W. , and Ingram K. T. , 2013: Assessing the Agricultural Reference Index (ARID) using uncertainty and sensitivity analyses. Agron. J., 105, 150160, doi:10.2134/agronj2012.0033.

    • Search Google Scholar
    • Export Citation
  • Xia, Y., and Coauthors, 2012a: Continental-scale water and energy flux analysis and validation for the North American Land Data Assimilation System project phase 2 (NLDAS-2): 1. Intercomparison and application of model products. J. Geophys. Res., 117, D03109, doi:10.1029/2011JD016048.

    • Search Google Scholar
    • Export Citation
  • Xia, Y., and Coauthors, 2012b: Continental-scale water and energy flux analysis and validation for North American Land Data Assimilation System project phase 2 (NLDAS-2): 2. Validation of model-simulated streamflow. J. Geophys. Res., 117, D03110, doi:10.1029/2011JD016051.

    • Search Google Scholar
    • Export Citation
  • Yang, J., Reichert P. , Abbaspour K. C. , Xia J. , and Yang H. , 2008: Comparing uncertainty analysis techniques for a SWAT application to the Chaohe Basin in China. J. Hydrol., 358, 123, doi:10.1016/j.jhydrol.2008.05.012.

    • Search Google Scholar
    • Export Citation
Save
  • Allen, R. G., Pereira L. S. , Raes D. , and Smith M. , 1998: Crop evapotranspiration: Guidelines for computing crop water requirements. FAO Irrigation and Drainage Paper 56, 300 pp. [Available online at www.fao.org/docrep/X0490E/X0490E00.htm.]

  • Berg, A. A., Famiglietti J. S. , Walker J. P. , and Houser P. R. , 2003: Impact of bias correction to reanalysis products on simulations of North American soil moisture and hydrological fluxes. J. Geophys. Res., 108, 4490, doi:10.1029/2002JD003334.

    • Search Google Scholar
    • Export Citation
  • Cai, J., Liu Y. , Lei T. , and Pereira L. S. , 2007: Estimating reference evapotranspiration with the FAO Penman–Monteith equation using daily weather forecast messages. Agric. For. Meteor., 145, 2235, doi:10.1016/j.agrformet.2007.04.012.

    • Search Google Scholar
    • Export Citation
  • Cai, J., Liu Y. , Xu D. , Paredes P. , and Pereira L. S. , 2009: Simulation of the soil water balance of wheat using daily weather forecast messages to estimate the reference evapotranspiration. Hydrol. Earth Syst. Sci., 13, 10451059, doi:10.5194/hess-13-1045-2009.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., and Juras J. , 2006: Measuring forecast skill: Is it real skill or is it the varying climatology? Quart. J. Roy. Meteor. Soc., 132, 29052923, doi:10.1256/qj.06.25.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., and Whitaker J. S. , 2006: Probabilistic quantitative precipitation forecasts based on reforecast analogs: Theory and application. Mon. Wea. Rev., 134, 32093229, doi:10.1175/MWR3237.1.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Whitaker J. S. , and Mullen S. L. , 2006: Reforecasts: An important dataset for improving weather predictions. Bull. Amer. Meteor. Soc., 87, 3346, doi:10.1175/BAMS-87-1-33.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Whitaker J. S. , Fiorino M. , and Benjamin S. G. , 2011: Global ensemble predictions of 2009’s tropical cyclones initialized with an ensemble Kalman filter. Mon. Wea. Rev., 139, 668688, doi:10.1175/2010MWR3456.1.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., Bates G. T. , Whitaker J. S. , Murray D. R. , Fiorino M. , Galarneau T. J. , Zhu Y. , and Lapenta W. , 2013: NOAA’s second-generation global medium-range ensemble reforecast dataset. Bull. Amer. Meteor. Soc., 94, 15531565, doi:10.1175/BAMS-D-12-00014.1.

    • Search Google Scholar
    • Export Citation
  • Ishak, A. M., Bray M. , Remesan R. , and Han D. W. , 2010: Estimating reference evapotranspiration using numerical weather modelling. Hydrol. Processes, 24, 34903509, doi:10.1002/hyp.7770.

    • Search Google Scholar
    • Export Citation
  • Khare, Y. P., Martinez C. J. , and Muñoz-Carpena R. , 2013: Parameter variability and drought models: A study using the Agricultural Reference Index for Drought (ARID). Agron. J., 105, 14171432, doi:10.2134/agronj2013.0167.

    • Search Google Scholar
    • Export Citation
  • Kumar, M., Raghuwanshi N. S. , and Singh R. , 2011: Artificial neural networks approach in evapotranspiration modeling: A review. Irrig. Sci., 29, 1125, doi:10.1007/s00271-010-0230-8.

    • Search Google Scholar
    • Export Citation
  • Marshall, C. H., Pielke R. A. Sr., Steyaert L. T. , and Willard D. A. , 2004: The impact of anthropogenic land-cover change on the Florida peninsula sea breezes and warm season sensible weather. Mon. Wea. Rev., 132, 2852, doi:10.1175/1520-0493(2004)132<0028:TIOALC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Meinke, H., Hammer G. L. , and Want P. , 1993: Potential soil water extraction by sunflower on a range of soils. Field Crops Res., 32, 5981, doi:10.1016/0378-4290(93)90021-E.

    • Search Google Scholar
    • Export Citation
  • Misra, V., Moeller L. , Stefanova L. , Chan S. , O’Brien J. J. , Smith T. J. III, and Plant N. , 2011: The influence of the Atlantic Warm Pool on the Florida panhandle sea breeze. J. Geophys. Res., 116, D00Q06, doi:10.1029/2010JD015367.

    • Search Google Scholar
    • Export Citation
  • Nash, J. E., and Sutcliffe J. V. , 1970: River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol., 10, 282290, doi:10.1016/0022-1694(70)90255-6.

    • Search Google Scholar
    • Export Citation
  • NRCS, 1986: Urban hydrology for small watersheds. Tech. Release 55, 2nd ed., 164 pp. [Available online at ftp://ftp.wcc.nrcs.usda.gov/wntsc/H&H/other/TR55documentation.pdf.]

  • Piechota, T., Chiew F. , Dracup J. , and McMahon T. , 2001: Development of exceedance probability streamflow forecast. J. Hydrol. Eng., 6, 2028, doi:10.1061/(ASCE)1084-0699(2001)6:1(20).

    • Search Google Scholar
    • Export Citation
  • Pinker, R. T., and Coauthors, 2003: Surface radiation budgets in support of the GEWEX Continental-Scale International Project (GCIP) and the GEWEX Americas Prediction Project (GAPP), including the North American Land Data Assimilation System (NLDAS) project. J. Geophys. Res., 108, 8844, doi:10.1029/2002JD003301.

    • Search Google Scholar
    • Export Citation
  • Potts, J. M., Folland C. K. , Jolliffe I. T. , and Sexton D. , 1996: Revised “LEPS” scores for assessing climate model simulations and long-range forecasts. J. Climate, 9, 3453, doi:10.1175/1520-0442(1996)009<0034:RSFACM>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Silva, D., Meza F. J. , and Varas E. , 2010: Estimating reference evapotranspiration (ETo) using numerical weather forecast data in central Chile. J. Hydrol., 382, 6471, doi:10.1016/j.jhydrol.2009.12.018.

    • Search Google Scholar
    • Export Citation
  • Sobash, R. A., Kain J. S. , Bright D. R. , Dean A. R. , Coniglio M. C. , and Weiss S. J. , 2011: Probabilistic forecast guidance for severe thunderstorms based on the identification of extreme phenomena in convection-allowing model forecasts. Wea. Forecasting, 26, 714728, doi:10.1175/WAF-D-10-05046.1.

    • Search Google Scholar
    • Export Citation
  • Srivastava, P. K., Han D. , Ramirez M. A. R. , and Islam T. , 2013: Comparative assessment of evapotranspiration derived from NCEP and ECMWF global datasets through Weather Research and Forecasting model. Atmos. Sci. Lett., 14, 118125, doi:10.1002/asl2.427.

    • Search Google Scholar
    • Export Citation
  • Srivastava, P. K., Han D. , Ramirez M. A. R. , and Islam T. , 2014: Sensitivity and uncertainty analysis of mesoscale model downscaled hydro-meteorological variables for discharge prediction. Hydrol. Processes, in press.

    • Search Google Scholar
    • Export Citation
  • Tian, D., and Martinez C. J. , 2012a: Forecasting reference evapotranspiration using retrospective forecast analogs in the southeastern United States. J. Hydrometeor., 13, 18741892, doi:10.1175/JHM-D-12-037.1.

    • Search Google Scholar
    • Export Citation
  • Tian, D., and Martinez C. J. , 2012b: Comparison of two analog-based downscaling methods for regional reference evapotranspiration forecasts. J. Hydrol., 475, 350364, doi:10.1016/j.jhydrol.2012.10.009.

    • Search Google Scholar
    • Export Citation
  • Trajkovic, S., Todorovic B. , and Stankovic M. , 2003: Forecasting of reference evapotranspiration by artificial neural networks. J. Irrig. Drain. Eng., 129, 454457, doi:10.1061/(ASCE)0733-9437(2003)129:6(454).

    • Search Google Scholar
    • Export Citation
  • Wilks, D. S., 2011: Statistical Methods in the Atmospheric Sciences.Academic Press, 476 pp.

  • Woli, P., Jones J. W. , Ingram K. T. , and Fraisse C. W. , 2012: An agricultural reference index for drought. Agron. J., 104, 287300, doi:10.2134/agronj2011.0286.

    • Search Google Scholar
    • Export Citation
  • Woli, P., Jones J. W. , and Ingram K. T. , 2013: Assessing the Agricultural Reference Index (ARID) using uncertainty and sensitivity analyses. Agron. J., 105, 150160, doi:10.2134/agronj2012.0033.

    • Search Google Scholar
    • Export Citation
  • Xia, Y., and Coauthors, 2012a: Continental-scale water and energy flux analysis and validation for the North American Land Data Assimilation System project phase 2 (NLDAS-2): 1. Intercomparison and application of model products. J. Geophys. Res., 117, D03109, doi:10.1029/2011JD016048.

    • Search Google Scholar
    • Export Citation
  • Xia, Y., and Coauthors, 2012b: Continental-scale water and energy flux analysis and validation for North American Land Data Assimilation System project phase 2 (NLDAS-2): 2. Validation of model-simulated streamflow. J. Geophys. Res., 117, D03110, doi:10.1029/2011JD016051.

    • Search Google Scholar
    • Export Citation
  • Yang, J., Reichert P. , Abbaspour K. C. , Xia J. , and Yang H. , 2008: Comparing uncertainty analysis techniques for a SWAT application to the Chaohe Basin in China. J. Hydrol., 358, 123, doi:10.1016/j.jhydrol.2008.05.012.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Study area and an example of a subset of 16 GEFS grid points and corresponding NLDAS-2 grid points in the Tampa Bay area.

  • Fig. 2.

    Locations of the observed stations and selected locations for WD forecast.

  • Fig. 3.

    Box-and-whisker plots of the NLDAS-2 ea(Tmin) (red) and the NLDAS-2 ea(SH) (blue) in 12 months at (a) Saint Leo, FL, and (b) Smithfield, NC.

  • Fig. 4.

    Box-and-whisker plots of NLDAS-2 (red) and observations (blue) for ETo, Tmax, Tmin, Rs, and Wind in 12 months at (a) Alfred Lake, FL, and (b) Clayton, NC.

  • Fig. 5.

    BSS and LEPS skill score (%) as a function of month and lead day. (a)–(e) The BSS reflected the extreme and tercile forecasting skill. (f) The LEPS skill score reflected the overall forecasting skill.

  • Fig. 6.

    (a),(b) Reliability and (c),(d) ROC diagrams for tercile and extreme forecasts at lead days 1 and 5.

  • Fig. 7.

    (top) LEPS skill score (%) and (below) BSS over the SEUS. The LEPS skill score reflected the overall forecasting skill. The BSS reflected the tercile and extreme forecasting skill.

  • Fig. 8.

    Comparisons of LEPS skill scores (%) for reforecast2, reforecast1a, and reforecast2b at lead days (a) 1 and (b) 5 throughout the year. (c) Comparisons of LEPS skill scores (%) for (left) reforecast2, (middle) reforecast1a, and (right) reforecast1b at lead days (top) 1 and (bottom) 5 over the SEUS.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 895 477 92
PDF Downloads 275 61 8