Precipitation Amount and Intensity Trends Across Southwest Saudi Arabia

This study examines precipitation accumulation and intensity trends across a region in southwest Saudi Arabia characterized by distinct seasonal weather patterns and mountainous terrain. The region is an example of an arid/semiarid area faced with maintaining sustainable water resources with a growing population. Annual and seasonal trends in precipitation amount were examined from 29 rain gages divided among four geographically unique regions from 1945/1946 to 2009. Two of the regions displayed significantly declining annual trends over the time series using a Mann‐Kendall test modified for autocorrelation (α < 0.05). Seasonal analysis revealed insignificant declining trends in at least two of the regions during each season. The largest and most consistent declining trends occurred during wintertime where all regions experienced negative trends. Several intensity metrics were examined in the study area from four additional stations containing daily data from 1985 to 2011. Intensity metrics included total precipitation, wet day count, simple intensity index, maximum daily annual rainfall, and upper/lower precipitation distribution changes. In general, no coherent trends were found among the daily stations suggesting precipitation is intensifying across the study area. The work represents the first of its size in the study area, and one of few in the region due to the lack of available long‐term data needed to properly examine precipitation changes.     

EVALUATION OF PROBABILITY DENSITY FUNCTIONS IN PRECIPITATION MODELS FOR THE PACIFIC NORTHWEST1

ABSTRACT: Recent work has found that a one-parameter Weibull model of wet day precipitation amount based on the Weibull distribution provides a better fit to historical daily precipitation data for eastern U.S. sites than other one-parameter models. The general two-parameter Weibull distribution was compared in this study to other widely used distributions for describing the distribution of daily precipitation event sizes at 99 sites from the U.S. Pacific Northwest. Surprisingly little performance was sacrificed by reducing the two-parameter Weibull to a single-parameter distribution. Advantages of the single-parameter model included requiring only the mean wet day precipitation amount for calibration, invertibility for simulation purposes, and ease of analytical manipulation. The fit of the single-parameter Weibull to the 99 stations included in this study was significantly better than other single-parameter models tested, and performed as well as the widely endorsed, more cumbersome, two-parameter gamma model. Both the one-and two-parameter Weibull distributions are shown to have b-moments that are consistent with historical precipitation data, while the ratio of b-skew and b-variance in the gamma model is inconsistent with the historical recerd by this measure. In addition, it was found that the two-parameter gamma distribution was better fit using the method of moments estimators than maximum likelihood estimates. These findings suggested that the distribution in precipitation among sites in the Pacific Northwest with dramatically different settings are nearly identical if expressed in proportion to the mean site event size.     

DISTRIBUTION AND STOCHASTIC GENERATION OF ANNUAL AND MONTHLY PRECIPITATION ON A MOUNTAINOUS WATERSHED IN SOUTHWEST IDAHO1

ABSTRACT: The 18-year precipitation record from the dense gage network on the Reynolds Creek Experimental Watershed located in southwest Idaho was used to determine the spatial distribution of annual and monthly precipitation on a mountainous watershed. Analyses of these data showed a linear relationship between annual amounts and elevation. This relationship was best when the gages were grouped into downwind and upwind sites. This grouping was appropriate because most of the winter storms moved over the watershed from the west and southwest, and the heaviest precipitation was on the west (downwind) side of the watershed. Gage sites along the western and southern watershed borders were most representative of the upwind gages on the east side, because they measured the precipitation from the air moving upwind onto the watershed. The maximum annual precipitation on the watershed was just leeward of the western watershed boundary. The monthly precipitation and elevation relationship was also best represented by grouping the gage sites into upwind and downwind sites. However, during the summer when there are only small amounts of pre cipitation and thunderstorms are the source of most precipitation, one equation can be used to represent the elevation relationship. This study also showed that the log-normal distribution could be used to generate the annual synthetic series, and the cube-root-normal distribution could be used to generate monthly synthetic series for all locations on the watershed.     

Probabilistic Streamflow Generation Model for Data Sparse Arid Watersheds1

Abstract: The authors present a model that generates streamflow for ephemeral arid streams. The model consists of a stochastic hourly precipitation point process model and a conceptual model that transforms precipitation into flow. It was applied to the Santa Cruz River at the border crossing from Mexico into Southern Arizona. The model was constructed for four different seasons and three categories of inter‐annual variability for the wet seasons of summer and winter. The drainage area is ungauged and precipitation information was inferred from a precipitation gauge downstream. The precipitation gauge record was evaluated against simulated precipitation from a mesoscale numerical weather prediction model, and was found to be the representative of the regional precipitation variability. The flow generation was found to reproduce the variability in the observed record at the daily, seasonal and annual time scales, and it is suitable for use in planning studies for the study site.     

ESTIMATION OF SOIL MOISTURE WITH API ALGORITHMS AND MICROWAVE EMISSION1

Large area soil moisture estimations are required to describe input to cloud prediction models, rainfall distribution models, and global crop yield models. Satellite mounted microwave sensor systems that as yet can only detect moisture at the surface have been suggested as a means of acquiring large area estimates. Relations previously discovered between microwave emission at the 1.55 cm wavelength and surface moisture as represented by an antecedent precipitation index were used to provide a pseudo infiltration estimation. Infiltration estimates based on surface wetness on a daily basis were then used to calculate the soil moisture in the surface 0–23 cm of the soil by use of a modified antecedent precipitation index. Reasonably good results were obtained (R²= 0.7162) when predicted soil moisture for the surface 23 cm was compared to measured moisture. Where the technique was modified to use only an estimate of surface moisture each three days an R² value of 0.7116 resulted for the same data set. Correlations between predicted and actual soil moisture fall off rapidly for repeat observations more than three days apart. The algorithms developed in this study may be used over relatively flat agricultural lands to provide improved estimates of soil moisture to a depth greater than the depth of penetration for the sensor.     

SPATIAL VARIATION OF PRECIPITATION ON THE TEXAS HIGH PLAINS1

ABSTRACT: Space autocorrelation techniques have been used to reveal the nature and spatial distribution of precipitation in the Texas High Plains. Correlation in precipitation amounts varies with both distance and direction, dropping off rapidly with distance, particularly during the warm season. The analyses can be used to estimate storm characteristics in conjunction with a wide variety of problems dealing with irrigation, crop yields, drainage and water supply, and evaluation of artificial weather modification efforts.     

CHARACTERISTICS OF EXTREME PRECIPITATION1

ABSTRACT: Records of extreme precipitation were investigated using the Discrete Autoregressive Moving Average (DABMA) process, which can explain long persistences of wet and dry spells that exist in daily precipitation data. The results show that the daily precipitation with strong autocorrelation is inclined to be better fit by a Discrete Autoregressive (DAB) model. On the other hand, those data with weak autocorrelations tend to be best fit by a Discrete Moving Average (DMA) model. It can also be concluded that based on the records from extremely wet and dry regions there is no geographic preference regarding the selection of the best model.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.     

Impacts of Changes in Precipitation Amount and Distribution on Water Resources Studied Using a Model Rainwater Harvesting System

Water supply reliability is expected to be affected by both precipitation amount and distribution changes under recent and future climate change. We compare historical (1951‐2010) changes in annual‐mean and annual‐maximum daily precipitation in the global set of station observations from Global Historical Climatology Network and climate models from the Inter‐Sectoral Impact Model Intercomparison Project (ISI‐MIP), and develop the study to 2011‐2099 for model projections under high radiative forcing scenario (RCP8.5). We develop a simple rainwater harvesting system (RWHS) model and drive it with observational and modeled precipitation. We study the changes in mean and maximum precipitation along with changes in the reliability of the model RWHS as tools to assess the impact of changes in precipitation amount and distribution on reliability of precipitation‐fed water supplies. Results show faster increase in observed maximum precipitation (10.14% per K global warming) than mean precipitation (7.64% per K), and increased reliability of the model RWHS driven by observed precipitation by an average of 0.2% per decade. The ISI‐MIP models show even faster increase in maximum precipitation compared to mean precipitation. However, they imply decreases in mean reliability, for an average 0.15% per decade. Compared to observations, climate models underestimate the increasing trends in mean and maximum precipitation and show the opposite direction of change in reliability of a model water supply system.     

GROUND-WATER RECHARGE THROUGH ACTWE SAND DUNES IN NORTHWESTERN NEVADA1

ABSTRACT: Most water-resouree investigations in semiarid basins of the Great Basin in western North America conclude that ground-water recharge from direct precipitation on the valley floor is negligible. However, many of these basins contain large areas covered by unvegetated, active sand dunes that may act as conduits for ground-water recharge. The potential for this previously undocumented recharge was investigated in an area covered by sand dunes in Desert Valley, northwestern Nevada, using a deep percolation model. The model uses daily measurements of precipitation and temperature th determine energy and moisture balances, from which estimates of long-term mean annual recharge are made. For the study area, the model calculated a mean annual recharge rate of as much as 1.3 inches per year, or 17 percent of the long-term mean precipitation. Model simulations also indicate that recharge would be virtually zero if the study area were covered by vegetation rather than dunes.     

Evaluating and Extending CLIGEN Precipitation Generation for the Loess Plateau of China1

Abstract: Climate generator (CLIGEN) is widely used in the United States to generate long‐term climate scenarios for use with agricultural systems models. Its applicability needs to be evaluated for use in a new region or climate. The objectives were to: (1) evaluate the reproducibility of the latest version of CLIGEN v5.22564 in generating daily, monthly, and yearly precipitation depths at 12 stations, as well as storm patterns including storm duration (D), relative peak intensity (ip), and peak intensity (rp) at 10 stations dispersed across the Loess Plateau and (2) test whether an exponential distribution for generating D and a distribution‐free approach for inducing desired rank correlation between precipitation depth and D can improve storm pattern generations. Mean absolute relative errors (MAREs) for simulating daily, monthly, annual, and annual maximum daily precipitation depth across all 12 stations were 3.5, 1.7, 1.7, and 5.0% for the mean and 5.0, 4.5, 13.0, and 13.6% for the standard deviations (SD), respectively. The model reproduced the distributions of monthly and annual precipitation depths well (p > 0.3), but the distribution of daily precipitation depth was less well produced. The first‐order, two‐state Markov chain algorithm was adequate for generating precipitation occurrence for the Loess Plateau of China; however, it underpredicted the longest dry periods. The CLIGEN‐generated storm patterns poorly. It underpredicted mean and SD of D for storms ≥10 mm by ?60.4 and ?72.6%, respectively. Compared with D, ip, and rp were slightly better reproduced. The MAREs of mean and SD were 21.0 and 52.1% for ip, and 31.2 and 55.2% for rp, respectively. When an exponential distribution was used to generate D, MAREs were reduced to 2.6% for the mean and 7.8% for the SD. However, ip estimation became much worse with MAREs being 128.9% for the mean and 241.1% for the SD. Overall, storm pattern generation needs improvement. For better storm pattern generation for the region, precipitation depth, D, and rp may be generated correlatively using Copula methods.     

Utility of Gridded Rainfall for IHACRES Daily River Flow Predictions in Southern California Watersheds1

Abstract: Studies to regionalize conceptual hydrologic models generally require rainfall and river flow data from multiple watersheds. Besides the considerable time (cost) to assemble and process rainfall data for many watersheds, investigators often need to choose from a number of candidate gauges, subjectively weighing the relative importance of proximity and elevation to select a representative rainfall dataset. The Unified Raingauge Dataset (URD) is a gridded daily rainfall dataset that covers the conterminous United States at 0.25 × 0.25 degrees spatial resolution and is available from 1948 to present. The objective of this study was to determine whether uncertainty in daily river flow predictions using the conceptual hydrologic model IHACRES in small to moderate size watersheds (50‐400 km²) in southern California would increase if URD gridded rainfall data were used in place of single rain gauge data to calibrate the model. Rain gauge data were obtained from the gauge nearest the watershed centroid and the gauge closest in elevation to the watershed mean elevation. Results from 20 randomly selected watersheds indicated that IHACRES calibration performance was similar using rainfall data from the URD grids and rain gauge data. There was some evidence of greater uncertainties associated with the URD calibrations in areas where topography may affect rainfall amounts. In contrast to the URD data, monthly gridded data produced by the Parameter‐Elevation Regressions on Independent Slopes Model (PRISM) includes adjustments for elevation and produces gridded values at a finer spatial resolution (4 km²). A limited test on two watersheds demonstrated that scaling the URD daily rainfall estimates to match the PRISM monthly values may improve rainfall estimates and model simulation performance.     

Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model1

Masih Ilyas, Shreedhar Maskey, Stefan Uhlenbrook, and Vladimir Smakhtin, 2011. Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model. Journal of the American Water Resources Association (JAWRA) 47(1):179‐195. DOI: 10.1111/j.1752‐1688.2010.00502.x Abstract: Reduction of input uncertainty is a challenge in hydrological modeling. The widely used model Soil Water Assessment Tool (SWAT) uses the data of a precipitation gauge nearest to the centroid of each subcatchment as an input for that subcatchment. This may not represent overall catchment precipitation conditions well. This paper suggests an alternative – using areal precipitation obtained through interpolation. The effectiveness of this alternative is evaluated by comparing its simulations with those based on the standard SWAT precipitation input procedure. The model is applied to mountainous semiarid catchments in the Karkheh River basin, Iran. The model performance is evaluated at daily, monthly, and annual scales by using a number of performance indicators at 15 streamflow gauging stations each draining an area in the range of 590‐42,620 km². The comparison suggests that the use of areal precipitation improves model performance particularly in small subcatchments in the range of 600‐1,600 km². The modified areal precipitation input results in increased reliability of simulated streamflows in the areas of low rain gauge density. Both precipitation input methods result in reasonably good simulations for larger catchments (over 5,000 km²). The use of areal precipitation input improves the accuracy of simulated streamflows with spatial resolution and density of rain gauges having significant impact on results.     

RESPONSE OF GREEN‐AMPT MEIN‐LARSEN SIMULATED RUNOFF VOLUMES TO TEMPORALLY AGGREGATED PRECIPITATION1

ABSTRACT: Data collection frequency in automated systems is user determined and can range from seconds to hours or days. Currently, there is no standard or recommended frequency interval for collecting precipitation data from automated systems for input to event‐based models such as Green‐Ampt Mein‐Larsen (GAML). Data from 47 storm events at seven locations were used to simulate the response of GAML excess rainfall to temporally aggregated precipitation data. No difference in model efficiency was recognized when comparing one‐minute interval data (R²= 1.00) to five‐minute data (R²= 1.00). Very little model efficiency was lost at a 10‐minute (R²= 0.96) interval. After 10‐minutes, decline in efficiency became more rapid with R²= 0.16 at one hour. The combined effect of time interval with respect to drainage area, hydraulic conductivity, maximum 30‐minute intensity, and total precipitation also revealed similar results.     

ESTIMATING REGIONAL SNOW WATER EQUIVALENT WITH A SIMPLE SIMULATION MODEL1

A simple simulation model designed to monitor snow-packs of the central Sierra Nevada is described. The model estimates average snow water equivalent for rectangular subregions in the area. Static subregion characteristics, daily precipitation and mean and minimum air temperatures measured at three index stations are the only needed input values. A water balance technique simulates daily snowpack changes in each subregion. Reasonable basinwide water equivalent values are produced. The procedure should be useful for estimating snow water distribution in large mountainous watersheds.     

NETWORK DESIGN FOR WATER SUPPLY FORECASTING IN THE WEST1

ABSTRACT: The objective is to develop techniques to evaluate how changes in basic data networks can improve accuracy of water supply forecasts for mountainous areas. The approach used was to first quantify how additional data would improve our knowledge of winter precipitation, and second to estimate how this knowledge translates, quantitatively, into improvement in forecast accuracy. A software system called DATANET was developed to analyze each specific gage network alternative. This system sets up a fine mesh of grid points over the basin. The long-term winter mean precipitation at each grid point is estimated using a simple atmospheric model of the orographic precipitation process. The mean runoff at each grid point is computed from the long-term mean precipitation estimate. The basic runoff model is calibrated to produce the observed long-term runoff. The error analysis is accomplished by comparing the error in forecasts based on the best possible estimate of precipitation using all available data with the error in the forecasts based on the best possible estimate of winter precipitation using only the gaged data. Different data network configurations of gage sites can be compared in terms of forecast errors.     

PRECIPITATION CHANGES FROM 1956 TO 1996 ON THE WALNUT GULCH EXPERIMENTAL WATERSHED1

ABSTRACT: The climate of Southern Arizona is dominated by summer precipitation, which accounts for over 60 percent of the annual total. Summer and non‐summer precipitation data from the USDA‐ARS Walnut Gulch Experimental Watershed are analyzed to identify trends in precipitation characteristics from 1956 to 1996. During this period, annual precipitation increased. The annual precipitation increase can be attributed to an increase in precipitation during non‐summer months, and is paralleled by an increase in the proportion of annual precipitation contributed during non‐summer months. This finding is consistent with previously reported increases in non‐summer precipitation in the southwestern United States. Detailed event data were analyzed to provide insight into the characteristics of precipitation events during this time period. Precipitation event data were characterized based on the number of events, event precipitation amount, 30‐minute event intensity, and event duration. The trend in non‐summer precipitation appears to be a result of increased event frequency since the number of events increased during nonsummer months, although the average amount per event, average event intensity, and average event duration did not. During the summer “monsoon” season, the frequency of recorded precipitation events increased but the average precipitation amount per event decreased. Knowledge of precipitation trends and the characteristics of events that make up a precipitation time series is a critical first step in understanding and managing water resources in semiarid ecosystems.     

Does the Temporal Resolution of Precipitation Input Influence the Simulated Hydrological Components Employing the SWAT Model?

This study aimed to evaluate the influence of sub‐daily precipitation time steps on model performance and hydrological components by applying the Green and Ampt infiltration method using the Soil and Water Assessment Tool (SWAT). Precipitation was measured at a resolution of 0.1 mm and aggregated to 5‐, 15‐, 30‐, and 60‐min time steps. Daily discharge data over a 10‐year period were used to calibrate and validate the model. Following a global sensitivity analysis, relevant parameters were optimized through an automatic calibration procedure using SWAT‐CUP for each time step. Daily performance statistics were almost equal among all four time steps (NSE ≈ 0.47). Discharge mainly consisted of groundwater flow (55%) and tile flow (42%), in reasonable proportions for the investigated catchment. In conclusion, model outputs were almost identical, showing simulations responded nearly independently of the chosen precipitation time step. This held true for (1) the selection of sensitive parameters, (2) performance statistics, (3) the shape of the hydrographs, and (4) flow components. However, a scenario analysis revealed that the precipitation time step becomes important when saturated hydraulic conductivities are low and curve numbers are high. The study suggests that there is no need in using precipitation time steps <1 h for lowland catchments dominated by soils with a low surface runoff potential if daily flow values are being considered. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

MULTIVARIATE TREND TESTING OF LAKE WATER QUALITY1

ABSTRACT: Multivariate methods of trend analysis offer the potential for higher power in detecting gradual water quality changes as compared to multiple applications of univariate tests. Simulation experiments were used to investigate the power advantages of multivariate methods for both linear model and Mann-Kendall based approaches. The experiments focused on quarterly observations of three water quality variables with no serial correlation and with several different intervariable correlation structures. The multivariate methods were generally more powerful than the univariate methods, offering the greatest advantage in situations where water quality variables were positively correlated with trends in opposing directions. For illustration, both the univariate and multivariate versions of the Mann-Kendall based tests were applied to case study data from several lakes in Maine and New York which have been sampled as part of EPA's long term monitoring study of acid precipitation effects.