Estimating marine biogeochemical rates of the carbonate pH system—A Kalman filter tested

Oxygen (O₂), nitrate (NO₃), dissolved inorganic carbon (DIC) or pCO₂, and pH or total alkalinity (TA), are useful indices of marine chemical, physical and biological processes operating on varying time-scales. Although these properties are increasingly being monitored at high frequency, they have not been extensively used for studying ecosystem dynamics. We test whether we can estimate time-evolving biogeochemical rates (e.g. primary production, respiration, calcification and carbonate dissolution, and nitrification) from synthetic high frequency time-series of O₂, NO₃, DIC, pCO₂, TA or pH. More specifically, a Kalman filter has been implemented in a very simplified biogeochemical model describing the dynamics of O₂, NO₃, DIC and TA and linking the concentration data to biogeochemical fluxes. Different sets of concentration data are assimilated and biogeochemical rates are estimated. The frequency of assimilation required to get acceptable results is investigated and is compared with the frequency of sampling in the field or in controlled experimental settings.Smoothing of the data to remove data noise before assimilation improves the estimation of the biogeochemical rates. The best estimated rates are obtained when assimilating O₂, NO₃ and TA although the assimilation of DIC instead of TA also gives satisfactory results. In case pH or pCO₂ is assimilated rather than DIC or TA, the linearization of the (now nonlinear) observation equation introduces perturbations and the Kalman filter behaves suboptimal. We conclude that, given the resolution of data required, the tool has potential to estimate biogeochemical rates of the carbonate system under controlled settings.     

CLIMATE FACTOR FOR SMALL-BASIN FLOOD FREQUENCY1

ABSTRACT: A climate factor, C_T, (T = 2–, 25-, and 100-year recurrence intervals) that delineates regional trends in small-basin flood frequency was derived using data from 71 long-term rainfall record sites. Values of C_T at these sites were developed by a regression analysis that related rainfall-runoff model estimates of T-year floods to a sample set of 50 model calibrations. C_T was regionalized via kriging to develop maps depicting its geographic variation for a large part of the United States east of the 105th meridian. Kriged estimates of C_T and basin-runoff characteristics were used to compute regionalized T-year floods for 200 small drainage basins. Observed T-year flood estimates also were developed for these sites. Regionalized floods are shown to account for a large percentage of the variability in observed flood estimates with coefficients of determination ranging from 0.89 for 2-year floods to 0.82 for 100-year floods. The relative importance of the factors comprising regionalized flood estimates is evaluated in terms of scale (size of drainage area), basin-runoff characteristics (rainfall. runoff model parameters), and climate (C_T).     

INSTANTANEOUS PEAK FLOW ESTIMATION PROCEDURES FOR NEWFOUNDLAND STREAMS1

ABSTRACT: A procedure of estimating instantaneous flood flows for various return periods on the Island of Newfoundland is presented. The procedure is based on annual maximum instantaneous flows rather than annual maximum daily-mean flows, as the latter requires the conversion of estimated daily-mean flows into instantaneous flows. Regression equations were developed for each of three homogeneous regions for the desired return periods. The flood flow estimation capability of the presented procedure is demonstrated to be better than any other currently available procedure on the Island.     

POLLUTANT LOADS IN URBAN STORMWATER: REVIEW OF METHODS FOR PLANNING-LEVEL ESTIMATES1

ABSTRACT: A review of methods for planning-level estimates of pollutant loads in urban stormwater focuses on transfer of charac. teristic runoff quality data to unmonitored sites, runoff monitoring, and simulation models. Load estimation by transfer of runoff quality data is the least expensive, but the accuracy of estimates is unknown. Runoff monitoring methods provide best estimates of existing loads, but cannot be used to predict load changes resulting from runoff controls, or other changes of the urban system. Simulation models require extensive calibration for reliable application. Models with optional formulations of pollutant build up, washoff, and transport can be better calibrated and the selection of options should be based on a statistical analysis of calibration data. Calibrated simulation models can be used for evaluation of control alternatives.     

PARAMETRIC ESTIMATION OF COVARIANCES OF REGIONALIZED VARIABLES1

ABSTRACT: The usefulness of stochastic models in describing the spatial variability of hydrogeologic quantities, such as permeability, storativity, piezometric head, seepage velocity, and solute concentrations is now widely recognized. In practice, these quantities are represented as the sum of a well-structured component, or drift, and a more erratic fluctuation component which is described statistically through its covariance function. This paper reviews some of the most recent and most promising methods for the estimation of parameters of these covariances from existing data. They are maximum likelihood, restricted maximum likelihood, minimum-variance unbiased quadratic estimation, and minimum-norm (weighted least squares) estimation. The applicability of such methods to conditional and unconditional probability problems is discussed.     

RUNOFF PREDICTION USING REMOTE SENSING IMAGERY1

ABSTRACT: Percent imperviousness is an important parameter in modeling the urban rainfall-runoff process and is usually determined using manual methods such as random sampling or conventional accounting methods. In this study two computerized methods are used for estimating the percent imperviousness of urban watersheds using high altitude remote sensing imagery. These methods include the Laser Image Processing Scanner and the Video-Tape Camera system. Imperviousness is directly estimated in the former method while in the latter it is estimated as a function of the statistics of the responses on emulsions of the imagery. The percent imperviousness computed by utilizing remote sensing imagery was used with the conceptual models of rainfall-runoff models. The models were applied to four urban watersheds and the runoff prediction results indicate that imperviousness determined by using remote sensing imagery was as accurate as that obtained by the manual methods, and that the use of remote sensing imagery requires significantly less time and money.     

STATISTICAL METHODS AND SAMPLING DESIGN FOR ESTIMATING STEP TRENDS IN SURFACE-WATER QUALITY1

ABSTRACT: This paper addresses two components of the problem of estimating the magnitude of step trends in surface water quality. The first is finding a robust estimator appropriate to the data characteristics expected in water-quality time series. The Hodges-Lehmann class of estimators is found to be robust in comparison to other nonparametric and moment-based estimators. A seasonal Hodges-Lehmann estimator is developed and shown to have desirable properties. Second, the effectiveness of various sampling strategies are examined using Monte Carlo simulation coupled with application of this estimator. The simulation is based on a large set of total phosphorus data from the Potomac River. To assure that the simulated records have realistic properties, the data are modeled in a multiplicative fashion incorporating flow, hysteresis, seasonal, and noise components. The results demonstrate the importance of balancing the length of the two sampling periods and balancing the number of data values between the two periods. The inefficiency of sampling at frequencies much in excess of 12 samples per year is demonstrated. Rotational sampling designs are discussed, and efficient designs, at least for this river and constituent, are shown to involve more than one year of active sampling at frequencies of about 12 per year.     

VELOCITY EQUATION FOR WATER QUALITY MODELING IN GEORGIA1

ABSTRACT: Tabletop water quality modeling still plays an important role in the water pollution control activities of the Georgia Environmental Protection Division. Tabletop models are those developed with out the aid of extensive field data. One important component of GEORGIA DOSAG, our basic water quality model, is the equation used to predict flow through velocity. However, Georgia is characterized by wide physiographic diversity which reduces the effectiveness of uncalibrated velocity equations. Using 15 years of accumulated time-of-travel studies, a series of empirical velocity equations were developed and calibrated to various physiographic conditions in Georgia. Equations are available for each major soil province and for three stream flow ranges within each province - Q<100 cfs, 100<Q<1000 cfs, and Q>1000 cfs. Now, in the absence of extensive field data, we have data based velocity equations which can be tailored to each site under study.     

KRIGED ESTIMATES OF TRANSMISSIVITY IN THE MESILLA BOLSON,NEW MEXICO1

ABSTRACT: Kriging utilizes a statistically based procedure of spatial interpolation that incorporates the spatial correlation structure of the phenomenon, and provides an error estimate. Kriging was applied to a total of 141 transmissivity values in an attempt to quantify the transmissivity distribution of the Santa Fe aquifer in Mesilla Bolson. New Mexico. The analysis produced contour maps of estimated transmissivity values and associated estimation variances. Through variogram analysis and fitting of an exponential variogrsm to 141 natural log of transmissivity (InT) values, the range was determined to be 3 miles, the average variance 2.74 (σ_InT= 1.65) with a mean of 8.65. Kriged estimates were generally lower when compared to estimates based on available transmissivity maps.     

MODELING FECAL COLIFORM CONTAMINATION IN THE RIO GRANDE1

This study examines sources of fecal coliform in Segment 2302 of the Rio Grande, located south of the International Falcon Reservoir in southern Texas. The watershed is unique because the contributing drainage areas lie in Texas and Mexico. Additionally, the watershed is mostly rural, with populated communities known as “colonias.” The colonias lack sewered systems and discharge sanitary water directly to the ground surface, thus posing an increased health hazard from coliform bacteria. Monitoring data confirm that Segment 2302 is not safe for contact recreation due to elevated fecal coliform levels. The goal of the study was to simulate the observed exceedences in Segment 2302 and evaluate potential strategies for their elimination. Fecal coliform contributions from ranching and colonia discharges were modeled using the Hydrologic Simulation Program‐Fortran (HSPF). Model results indicated that the regulatory 30‐day geometric mean fecal coliform concentration of 200 colony forming units (cfu) per 100 milliliters is exceeded approximately three times per year for a total of 30 days. Ongoing initiatives to improve wastewater facilities will reduce this to approximately once per year for 14 days. Best management practices will be necessary to reduce cattle access to streams and eliminate all exceedences. The developed model was limited by the relatively sparse flow and fecal coliform data.