GENERATION OF SYNTHETIC AND MISSING CLIMATIC DATA FOR PUERTO RICO1

ABSTRACT: Climatic data such as temperature, solar radiation, relative humidity, and wind speed have been widely used to estimate evapotranspiration. Moat of the solar radiation data and portions of the relative humidity data are either not available or missing from the records in Puerto Rico. Depending upon the availability and data characteristics of records, three methods (including a regression technique, an averaging of historical data, and a regional average) were used to generate missing data, and a time series analysis was used to synthesize a series of climatic data. The limitations and applicability of each method are discussed. The results showed that the time series analysis method can be successfully used to synthesize a series of monthly solar radiations for several stations. The regression technique and the regional average can be successfully applied to generate missing monthly solar radiation data. The regression technique and the averaging of historical data have been satisfactorily used to interpolate missing monthly relative humidity. The explained variance (R²) varied from 0.68 to 0.88, which are both significant at the 0.05 level of significance.     

SUPPLEMENTATION OF MISSING VALUES IN WATER QUALITY DATA1

ABSTRACT The problem of estimating missing values in water quality data using linear interpolation and harmonic analysis is studied to see which one of these two methods yields better estimates for the missing values. The data used in this study consisted of midnight values of dissolved oxygen from the Ohio River collected over a period of one year at Stratton station. Various hypothetical cases of missing data are considered and the two methods of supplementing missing values are evaluated using statistical tests. The results indicate that when the percentage of missed data points exceeded ten percent of the total number in the original sample, harmonic analysis usually yielded better estimates for both the regularly and irregularly missed cases. For data that exhibit cyclic variation, examples of which are dissolved oxygen concentration and water temperature, harmonic analysis as a data generation technique appears to be superior to linear interpolation.     

Data integration across borders: a case study of the Abbotsford‐Sumas aquifer (British Columbia/Washington State)1

Abstract: Integrating spatial datasets from diverse sources is essential for cross‐border environmental investigations and decision‐making. This is a little investigated topic that has profound implications for the availability and reliability of spatial data. At present, ground‐water hydrostratigraphic models exist for both the Canadian or for the United States (U.S.) portion of the aquifer but few are integrated across the border. In this paper, we describe the challenges of integrating multiple source, large datasets for development of a ground‐water hydrostratigraphic model for the Abbotsford‐Sumas Aquifer. Growing concerns in Canada regarding excessive withdrawal south of the border and in the U.S. regarding nitrate contamination originating north of the border make this particular aquifer one of international interest. While much emphasis in GIScience is on theoretical solutions to data integration, such as current ontology research, this study addresses pragmatic ways of integrating data across borders. Numerous interoperability challenges including the availability of data, metadata, data formats and quality, database structure, semantics, policies, and cooperation are identified as inhibitors of data integration for cross‐border studies. The final section of the paper outlines two possible solutions for standardizing classification schemes for ground‐water models – once data heterogeneity has been addressed.     

TREND ASSESSMENT OF WATER QUALITY TIME SERIES1

ABSTRACT A general methodology is described for identifying and statistically modeling trends which may be contained in a water quality time series. A range of useful exploratory data analysis tools are suggested for discovering important patterns and statistical characteristics of the data such as trends caused by external interventions. To estimate the entries in an evenly spaced time series when data are available at irregular time intervals, a new procedure based upon seasonal adjustment is described. Intervention analysis is employed at the confirmatory data analysis stage to rigorously model changes in the mean levels of a series which are identified using exploratory data analysis techniques. Furthermore, intervention analysis can be utilized for estimating missing observations when they are not too numerous. The effects of cutting down a forest upon various water quality variables and also the consequences of acid rain upon the alkalinity in a stream provide illustrative applications which demonstrate the effectiveness of the methodology.     

Assessing the efficiency gains of improved spatial targeting of policy interventions; the example of an agri-environmental scheme

GIS-based spatial targeting is increasingly recognised as a potentially useful tool to design more efficient policy interventions. The use of this tool has also been advocated in the context of incentive-based agri-environmental schemes, but there has been little work to date to estimate the level of efficiency gains which it may help to achieve. This paper investigates the requirements to arrive at such estimates, using a Scottish farm woodland scheme as a case study. This agri-environmental scheme aims to provide visual amenity and biodiversity. Maps of these two benefits are used to develop improved spatial targeting scenarios that deliver significant efficiency gains in comparison to the existing scheme design. The paper discusses the nature of the spatial distribution of the relevant benefits at the landscape scale and the data requirements for the realistic estimation of efficiency gains. It concludes that although much work needs to be done, the methods available today could and should play a much greater role in improving the landscape-scale design of existing land use schemes focused on the delivery of non-market benefits.     

AN EVALUATION OF THE AVAILABLE TECHNIQUES FOR ESTIMATING MISSING FECAL COLIFORM DATA1

ABSTRACT: This paper presents the findings of a study aimed at evaluating the available techniques for estimating missing fecal coliform (FC) data on a temporal basis. The techniques investigated include: linear and nonlinear regression analysis and interpolation functions, and the use of artificial neural networks (ANNs). In all, seven interpolation, two regression, and one ANN model structures were investigated. This paper also investigates the validity of a hypothesis that estimating missing FC data by developing different models using different data corresponding to different dynamics associated with different trends in the FC data may result in a better model performance. The FC data (counts/100 ml) derived from the North Fork of the Kentucky River in Kentucky were employed to calibrate and validate various models. The performance of various models was evaluated using a wide variety of standard statistical measures. The results obtained in this study are able to demonstrate that the ANNs can be preferred over the conventional techniques in estimating missing FC data in a watershed. The regression technique was not found suitable in estimating missing FC data on a temporal basis. Further, it has been found that it is possible to achieve a better model performance by first decomposing the whole data set into different categories corresponding to different dynamics and then developing separate models for separate categories rather than developing a single model for the composite data set.     

Evaluation of Gauge‐Radar Merging Methods Using a Semi‐Distributed Hydrological Model in the Upper Thames River Basin,Canada

Gauge‐radar merging methods combine rainfall estimates from rain gauges and radar to capitalize on the strengths of the individual instruments. The performance of four well‐known gauge‐radar merging methods, including mean field bias correction, Brandes spatial adjustment, local bias correction using kriging, and conditional merging, are examined using Environment Canada radar and the Upper Thames River Basin in southwestern Ontario, Canada, as a case study. The analysis assesses the effect of gauge‐radar merging methods on: (1) the accuracy of predicted rainfall accumulations; and (2) the accuracy of predicted streamflows using a semi‐distributed hydrological model. In addition, several influencing factors (i.e., gauge density, storm type, basin type, proximity to the radar tower, and time‐step of adjustment) are analyzed to determine their effect on the performance of the rainfall estimation techniques. Confirming results of previous studies, the merging methods provide an increase in the accuracy of both rainfall accumulation estimations and predicted streamflows. The results also indicate specific factors such as gauge density, rainfall intensity, and time‐step of adjustment can reduce the accuracy of merging methods and play a key role in the examination of its use for operational purposes. Results provide guidance for hydrologists and engineers assessing how best to apply corrected radar products to improve rainfall estimation and hydrological modeling accuracy.     

Grouping Lakes for Water Quality Assessment and Monitoring: The Roles of Regionalization and Spatial Scale

Regionalization frameworks cluster geographic data to create contiguous regions of similar climate, geology and hydrology by delineating land into discrete regions, such as ecoregions or watersheds, often at several spatial scales. Although most regionalization schemes were not originally designed for aquatic ecosystem classification or management, they are often used for such purposes, with surprisingly few explicit tests of the relative ability of different regionalization frameworks to group lakes for water quality monitoring and assessment. We examined which of 11 different lake grouping schemes at two spatial scales best captures the maximum amount of variation in water quality among regions for total nutrients, water clarity, chlorophyll, overall trophic state, and alkalinity in 479 lakes in Michigan (USA). We conducted analyses on two data sets: one that included all lakes and one that included only minimally disturbed lakes. Using hierarchical linear models that partitioned total variance into within-region and among-region components, we found that ecological drainage units and 8-digit hydrologic units most consistently captured among-region heterogeneity at their respective spatial scales using all lakes (variation among lake groups = 3% to 50% and 12% to 52%, respectively). However, regionalization schemes capture less among-region variance for minimally disturbed lakes. Diagnostics of spatial autocorrelation provided insight into the relative performance of regionalization frameworks but also demonstrated that region size is only partly responsible for capturing variation among lakes. These results suggest that regionalization schemes can provide useful frameworks for lake water quality assessment and monitoring but that we must identify the appropriate spatial scale for the questions being asked, the type of management applied, and the metrics being assessed.     

COMBINATION OF RADAR AND GAUGE DATA IN A RAINFALL ARCHIVE SYSTEM1

ABSTRACT: Near real time daily rainfall estimates for the UK are available from three sources: a sparse network of gauges, radar data, or radar data adjusted by the sparse gauges. The PARAGON rainfall archive system, which has been developed by the UK Meteorological Office, is able to produce these estimates in near real time on a 5 km grid. The ability of these estimates to reproduce the 5 km grid point field derived later from a dense network of gauges is compared using case studies. Five techniques have been used to assess the relative quality of the various estimates. There is general agreement between the results of the various techniques. For the London radar there are examples of days when the rainfall estimate was improved by incorporating radar data; conversely, there are days when the radar data make it worse. Overall little evidence was found to suggest that adjusted radar data are consistently markedly better than gauge estimates. Discriminate use of radar data is recommended.     

QUALITY OF 1995 SPRING TOTAL DISSOLVED GAS DATA: COLUMBIA AND SNAKE RWERS1

ABSTRACT: The quality of the U.S. Army Corps of Engineers' (Corps) total dissolved gas (TDG) data base for the 1995 spring spill season was reviewed to determine the value of this information in real-time management decisions regarding river operations. We concluded that problems in transmitting, archiving, correcting and interpreting the records constitute significant sources of data anomalies that affect the accuracy and reliability of information necessary to manage spill and TDG in the Columbia and Snake rivers. The data base that was reviewed covers 25 selected Columbia and Snake river stations, and includes real-time TDG data needed to regulate spill operations to maintain gas levels within state water quality standards and to monitor effects on fish and aquatic life during the salmon migration season. A wide range of anomalies (daily averages missing or in error or based on incomplete records) was detected in more than one-third (37 percent) of the Corps' gas data base. Extreme anomalies (daily averages including errors and discontinuities for more than eight hours in a day) were found in 16 percent of the data base. The Fish Passage Center, also reviewed the Corps' data and reported an overall 33 percent incidence of anomalous days. Despite arriving at similar findings about the Corps' data base, we detected a 28 percent discrepancy in the type of data anomalies between our analyses. Real. time improvements in the quality of the dissolved gas data base are necessary to provide managers with a reliable product from this monitoring effort.     

GIS‐Based Spatial Precipitation Estimation: A Comparison of Geostatistical Approaches1

Abstract: As one of the primary inputs that drive watershed dynamics, the estimation of spatial variability of precipitation has been shown to be crucial for accurate distributed hydrologic modeling. In this study, a Geographic Information System program, which incorporates Nearest Neighborhood (NN), Inverse Distance Weighted (IDW), Simple Kriging (SK), Ordinary Kriging (OK), Simple Kriging with Local Means (SKlm), and Kriging with External Drift (KED), was developed to facilitate automatic spatial precipitation estimation. Elevation and spatial coordinate information were used as auxiliary variables in SKlm and KED methods. The above spatial interpolation methods were applied in the Luohe watershed with an area of 5,239 km², which is located downstream of the Yellow River basin, for estimating 10 years’ (1991‐2000) daily spatial precipitation using 41 rain gauges. The results obtained in this study show that the spatial precipitation maps estimated by different interpolation methods have similar areal mean precipitation depth, but significantly different values of maximum precipitation, minimum precipitation, and coefficient of variation. The accuracy of the spatial precipitation estimated by different interpolation methods was evaluated using a correlation coefficient, Nash‐Sutcliffe efficiency, and relative mean absolute error. Compared with NN and IDW methods that are widely used in distributed hydrologic modeling systems, the geostatistical methods incorporated in this GIS program can provide more accurate spatial precipitation estimation. Overall, the SKlm_EL_X and KED_EL_X, which incorporate both elevation and spatial coordinate as auxiliary into SKlm and KED, respectively, obtained higher correlation coefficient and Nash‐Sutcliffe efficiency, and lower relative mean absolute error than other methods tested. The GIS program developed in this study can serve as an effective and efficient tool to implement advanced geostatistics methods that incorporate auxiliary information to improve spatial precipitation estimation for hydrologic models.     

Environmental mapping based on spatial variability

Environmental maps show the probable environmental states of different types of land use or development of landscape in a geographic context. Remotely sensed data are particularly efficient for environmental mapping in order to outline major environmental types. Multiple schemes of image classification used in environmental mapping are either traditionally statistical or heuristic. While the former methods do not take account of spatial variability in space and aerial data, the latter ones does not lend themselves to optimal solutions we present. Novel probabilistic models of piecewise-homogeneous images are used in environmental mapping to segment real images. The models consider both an image and a land cover map. Such a pair constitutes an example of a Markov random field specified by a joint Gibbs probability distribution of images and maps. Parameters of the model are estimated by using a stochastic approximation technique. Its convergence to the desired values is studied experimentally. Addition of spatial attributes appears to be necessary in most areas where the differences in spatial data between regions in the image occur. Experiments in generating the pairs of images and environmental maps and in segmenting the simulated as well as real images are discussed.     

MODELING HYDROLOGIC TIME SERIES FROM THE ARCTIC1

The general intervention model is applied to hydrologic and meteorologjc time series from the Canadian Arctic. The authors show how the model is able to account for environmental interventions, missing observations in the data, changes in data collection procedures, the effects of external inputs, as well as seasonality and autocorrelation. Methods for identifying transfer functions by making use of a physical understanding of the processes involved are demonstrated and sample applications of the general intervention model to Arctic data are shown.     

A COMPARATIVE ANALYSIS OF TECHNIQUES FOR SPATIAL INTERPOLATION OF PRECIPITATION1

One of the problems which often arises in engineering hydrology is to estimate data at a given site because either the data are missing or the site is ungaged. Such estimates can be made by spatial interpolation of data available at other sites. A number of spatial interpolation techniques are available today with varying degrees of complexity. It is the intent of this paper to compare the applicability of various proposed interpolation techniques for estimating annual precipitation at selected sites. The interpolation techniques analyzed include the commonly used Thiessen polygon, the classical polynomial interpolation by least-squares or Lagrange approach, the inverse distance technique, the multiquadric interpolation, the optimal interpolation and the Kriging technique. Thirty years of annual precipitation data at 29 stations located in the Region II of the North Central continental United States have been used for this study. The comparison is based on the error of estimates obtained at five selected sites. Results indicate that the Kriging and optimal interpolation techniques are superior to the other techniques. However, the multiquadric technique is almost as good as those two. The inverse distance interpolation and the Thiessen polygon gave fairly satisfactory results while the polynomial interpolation did not produce good results.     

Analyzing the Effects of Groundwater Pumping on an Urban Stream‐Aquifer System

When the cone of influence of a pumping well reaches a nearby river, the resulting hydraulic gradient can induce enhanced seepage of streamflow into the aquifer. The rate of seepage is often modeled using analytical solutions that are simple to apply but may not reproduce field data due to mathematical assumptions not being met in the field. Furthermore, the appropriateness of such models has not been investigated in detail due to difficulty in measuring streamflow loss in the field. In this study, a field experiment was conducted on a reach of the South Platte River near Denver, Colorado to estimate pumping‐induced streamflow loss. A network of stream gauges, monitoring wells, and in situ measurements was used to observe streamflow rates, groundwater levels, and temperature to assess if pumping wells have a significant impact on streamflow, and to compare observed streamflow depletion against analytical solutions. Data collected suggest that pumping wells have a noticeable impact on streamflow. The analytical solutions proved accurate if streamflow was low and constant but performed poorly if streamflow was high and variable. Therefore, for this reach, the use of analytical solutions to predict streamflow may only be appropriate under low‐flow, constant‐flow conditions. Methods and results can be used to guide other streamflow depletion studies and to inform cases of pumping‐induced streamflow depletion, particularly in regard to water rights.     

Exploring Persistence in Streamflow Forecasting

In this study, the authors explore three persistence approaches in streamflow forecasting motivated by the need for forecasting model skill evaluation. The authors use streamflow observations with 15 min resolution from the year 2008 to 2017 at 140 United States Geological Survey streamflow gauges monitoring the streams and rivers over the State of Iowa. The spatial scale of the basins ranges from about 7 to 37,000 km². The study explores three approaches: simple persistence, gradient persistence, and anomaly persistence. The study shows that persistence forecasts skill has strong dependence on basin scales and weaker but non‐negligible dependence on geometric properties of the river network for a given basin. Among the three approaches explored, anomaly persistence shows highest skill especially for small basins, under about 500 km². The anomaly persistence can serve as a benchmark for model evaluations considering the effect of basin scales and geometric properties of river network of the basin. This study further reiterates that persistence forecasts are hard‐to‐beat methods for larger basin scales at short to medium forecast range.     

Evaluation of the Global Precipitation Measurement (GPM) Satellite Rainfall Products over the Lower Colorado River Basin,Texas

Quality of precipitation products from the Integrated Multi‐satellitE Retrievals for Global Precipitation Measurement mission (IMERG) was evaluated over the Lower Colorado River Basin of Texas. Observations of several rainfall events of a wide range of magnitudes during May 2015 by a very dense network of 241 rain gauges over the basin were used as a reference. The impact of temporal and spatial downscaling of different satellite products (near/post‐real‐time) on their accuracy was studied. Generally, all IMERG products perform better when the temporal and spatial resolutions are downscaled. The Final product shows relatively better performance compared to the near‐real‐time products in terms of basic performance measures; however, regarding rainfall detection, all products show nearly similar performance. When considering rainfall detection, IMERG adequately captures the precipitation events; however, in terms of spatial patterns and accuracy, more improvements are needed. IMERG products analysis results may help developers gain insight into the regional performance of the product, improve the product algorithms, and provide information to end users on the products’ suitability for potential hydrometeorological applications. Overall, the IMERG products, even the uncalibrated product at its finest resolution, showed reasonable performance indicating their great potential for applications such as water resources management, prevention of natural disasters, and flood forecasting.     

Resampling methods for evaluating classification accuracy of wildlife habitat models

Predictive models of wildlife-habitat relationships often have been developed without being tested The apparent classification accuracy of such models can be optimistically biased and misleading. Data resampling methods exist that yield a more realistic estimate of model classification accuracy These methods are simple and require no new sample data. We illustrate these methods (cross-validation, jackknife resampling, and bootstrap resampling) with computer simulation to demonstrate the increase in precision of the estimate. The bootstrap method is then applied to field data as a technique for model comparison We recommend that biologists use some resampling procedure to evaluate wildlife habitat models prior to field evaluation.     

A geophysical perspective of value of information: examples of spatial decisions for groundwater sustainability

The value of information (VOI) can be used to determine what kind of spatial information maybe relevant and useful for groundwater sustainability decisions. In this paper, the unique challenges for applying VoI to spatial information from geophysical data are described. The uncertainty regarding the spatial structure or continuity of the subsurface properties can be described with geostatistical sample models. Using these models, one can quantify the prior value given our present state of uncertainty and a set of decision alternatives and outcomes. Because geophysical techniques are a type of remote-sensing data, assuming “perfect” information is not realistic since the techniques usually are indirectly sampling the aquifer properties. Therefore, the focus of this paper is describing how the data reliability (the measure of imperfectness) can be quantified. One of the foremost considerations is the non-unique relationship between geological parameters (which determine groundwater flow) and geophysical observables (what determines the response of the technique). Another is to have the information in a form such that it is useful for spatial decisions. This will often require inversion and interpretation of the geophysical data. Inversion reconstructs an image of the subsurface from the raw geophysical data. How closely the image reproduces the true subsurface structure or property of interest depends on the particular technique’s resolution, depth of investigation and sensor locations. Lastly, in some cases, interpretation of the geophysical data or inversion will be necessary to link the data to the variables that determine the outcome of the decision. Three examples are provided that illustrate different approaches and methods for addressing these challenges. In the examples, time-domain electromagnetic and electrical resistivity techniques are evaluated for their ability to assist in spatial decisions for aquifer management. The examples considered address these three situations: aquifer vulnerability to surface–borne contaminants, managed aquifer recharge and CO₂/brine leakage (related to CO₂ geologic sequestration activities). The methods presented here are transferable to other subsurface sciences and decisions that involve risk. Recent work has been applied to geothermal well-siting using electromagnetic techniques. These approaches can also be applied for oil and mining spatial decisions, and they offer advantages over previous VOI work done for oil applications: they explicitly include the geologic uncertainty modeling and simulate the physics of the considered geophysical technique.