Modeling Postfire Response and Recovery using the Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS)1

Abstract: This paper investigates application of the Army Corps of Engineers’ Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS) to a burned watershed in San Bernardino County, California. We evaluate the HEC‐HMS’ ability to simulate discharge in prefire and postfire conditions in a semi arid watershed and the necessary parameterizations for modeling hydrologic response during the immediate, and subsequent recovery, period after a wildfire. The model is applied to City Creek watershed, which was 90% burned during the Old Fire of October 2003. An optimal spatial resolution for the HEC‐HMS model was chosen based on an initial sensitivity analysis of subbasin configurations and related model performance. Five prefire storms were calibrated for the selected model resolution, defining a set of parameters that reasonably simulate prefire conditions. Six postfire storms, two from each of the following rainy (winter) seasons were then selected to simulate postfire response and evaluate relative changes in parameter values and model behavior. There were clear trends in the postfire parameters [initial abstractions (I_a), curve number (CN), and lag time] that reveal significant (and expected) changes in watershed behavior. CN returns to prefire (baseline) values by the end of Year 2, while I_a approaches baseline by the end of the third rainy season. However, lag time remains significantly lower than prefire values throughout the three‐year study period. Our results indicate that recovery of soil conditions and related runoff response is not entirely evidenced by the end of the study period (three rainy seasons postfire). Understanding the evolution of the land surface and related hydrologic properties during the highly dynamic postfire period, and accounting for these changes in model parameterizations, will allow for more accurate and reliable discharge simulations in both the immediate, and subsequent, rainy seasons following fire.     

Development of a Quantitative Pasture Phosphorus Management Tool Using the SWAT Model1

Abstract: Assessment tools to evaluate phosphorus loss from agricultural lands allow conservation planners to evaluate the impact of management decisions on water quality. Available tools to predict phosphorus loss from agricultural fields are either: (1) qualitative indices with limited applicability to address offsite water quality standards, or (2) models which are prohibitively complex for application by most conservation planners. The purpose of this research was to develop a simple interface for a comprehensive hydrologic/water quality model to allow its usage by farmers and conservation planners. The Pasture Phosphorus Management (PPM) Calculator was developed to predict average annual phosphorus (P) losses from pastures under a variety of field conditions and management options. PPM Calculator is a vastly simplified interface for the Soil and Water Assessment Tool (SWAT) model that requires no knowledge of SWAT by the user. PPM Calculator was validated using 33 months of data on four pasture fields in northwestern Arkansas. This tool has been extensively applied in the Lake Eucha/Spavinaw Basin in northeastern Oklahoma and northwestern Arkansas. PPM Calculator allows conservation planners to take advantage of the predictive capacity of a comprehensive hydrologic water quality model typically reserved for use by hydrologists and engineers. This research demonstrates the applicability of existing water quality models in the development of user friendly P management tools.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Watershed Regressions for Pesticides (warp) Models for Predicting Atrazine Concentrations in Corn Belt Streams1

Stone, Wesley W. and Robert J. Gilliom, 2012. Watershed Regressions for Pesticides (WARP) Models for Predicting Atrazine Concentrations in Corn Belt Streams. Journal of the American Water Resources Association (JAWRA) 48(5): 970‐986. DOI: 10.1111/j.1752‐1688.2012.00661.x Abstract: Watershed Regressions for Pesticides (WARP) models, previously developed for atrazine at the national scale, are improved for application to the United States (U.S.) Corn Belt region by developing region‐specific models that include watershed characteristics that are influential in predicting atrazine concentration statistics within the Corn Belt. WARP models for the Corn Belt (WARP‐CB) were developed for annual maximum moving‐average (14‐, 21‐, 30‐, 60‐, and 90‐day durations) and annual 95th‐percentile atrazine concentrations in streams of the Corn Belt region. The WARP‐CB models accounted for 53 to 62% of the variability in the various concentration statistics among the model‐development sites. Model predictions were within a factor of 5 of the observed concentration statistic for over 90% of the model‐development sites. The WARP‐CB residuals and uncertainty are lower than those of the National WARP model for the same sites. Although atrazine‐use intensity is the most important explanatory variable in the National WARP models, it is not a significant variable in the WARP‐CB models. The WARP‐CB models provide improved predictions for Corn Belt streams draining watersheds with atrazine‐use intensities of 17 kg/km² of watershed area or greater.     

Comparison of Stream Invertebrate Response Models for Bioassessment Metrics1

Waite, Ian R., Jonathan G. Kennen, Jason T. May, Larry R. Brown, Thomas F. Cuffney, Kimberly A. Jones, and James L. Orlando, 2012. Comparison of Stream Invertebrate Response Models for Bioassessment Metrics. Journal of the American Water Resources Association (JAWRA) 48(3): 570-583. DOI: 10.1111/j.1752-1688.2011.00632.x Abstract: We aggregated invertebrate data from various sources to assemble data for modeling in two ecoregions in Oregon and one in California. Our goal was to compare the performance of models developed using multiple linear regression (MLR) techniques with models developed using three relatively new techniques: classification and regression trees (CART), random forest (RF), and boosted regression trees (BRT). We used tolerance of taxa based on richness (RICHTOL) and ratio of observed to expected taxa (O/E) as response variables and land use/land cover as explanatory variables. Responses were generally linear; therefore, there was little improvement to the MLR models when compared to models using CART and RF. In general, the four modeling techniques (MLR, CART, RF, and BRT) consistently selected the same primary explanatory variables for each region. However, results from the BRT models showed significant improvement over the MLR models for each region; increases in R² from 0.09 to 0.20. The O/E metric that was derived from models specifically calibrated for Oregon consistently had lower R² values than RICHTOL for the two regions tested. Modeled O/E R² values were between 0.06 and 0.10 lower for each of the four modeling methods applied in the Willamette Valley and were between 0.19 and 0.36 points lower for the Blue Mountains. As a result, BRT models may indeed represent a good alternative to MLR for modeling species distribution relative to environmental variables.     

Critical Steps for Continuing Advancement of Satellite Rainfall Applications for Surface Hydrology in the Nile River Basin1

Gebremichael, Mekonnen, Emmanouil N. Anagnostou, and Menberu M. Bitew, 2010. Critical Steps for Continuing Advancement of Satellite Rainfall Applications for Surface Hydrology in the Nile River Basin. Journal of the American Water Resources Association (JAWRA) 46(2):361-366. DOI: 10.1111/j.1752-1688.2010.00428.x. Abstract: Given the increasingly higher resolution and data accessibility, satellite precipitation products could be useful for hydrological application in the Nile River Basin, which is characterized by lack of reasonably dense hydrological in situ sensors and lack of access to the existing dataset. However, in the absence of both extreme caution and research results for the Nile basin, the satellite rainfall (SR) products may not be used, or may even be used erroneously. We identify two steps that are critical to enhance the value of SR products for hydrological applications in the Nile basin. The first step is to establish representative validation sites in the Nile basin. The validation site will help to quantify the errors in the different kinds of SR products, which will be used to select the best products for the Nile basin, include the errors in decision making, and design strategies to minimize the errors. Using rainfall measurements collected from the unprecedented high-density rain gauge network over a small region within the Nile basin, we indicate that SR estimates could be subject to significant errors, and quantification of estimation errors by way of establishing validation sites is critically important in order to use the SR products. The second step is to identify the degree of hydrologic model complexity required to obtain more accurate hydrologic simulation results for the Nile basin when using SR products as input. The level of model complexity may depend on basin size and SR algorithm, and further research is needed to spell out this dependence for the Nile basin.     

Inhalation dose assessment of indoor radon progeny using biokinetic and dosimetric modeling and its application to Jordanian population

High indoor radon concentrations in Jordan result in internal exposures of the residents due to the inhalation of radon and its short-lived progeny. It is therefore important to quantify the annual effective dose and further the radiation risk to the radon exposure. This study describes the methodology and the biokinetic and dosimetric models used for calculation of the inhalation doses exposed to radon progeny. The regional depositions of aerosol particles in the human respiratory tract were firstly calculated. For the attached progeny, the activity median aerodynamic diameters of 50 nm, 230 nm and 2500 nm were chosen to represent the nucleation, accumulation and coarse modes of the aerosol particles, respectively. For the unattached progeny, the activity median thermodynamic diameter of 1 nm was chosen to represent the free progeny nuclide in the room air. The biokinetic models developed by the International Commission on Radiological Protection (ICRP) were used to calculate the nuclear transformations of radon progeny in the human body, and then the dosimetric model was applied to estimate the organ equivalent doses and the effective doses with the specific effective energies derived from the mathematical anthropomorphic phantoms. The dose conversion coefficient estimated in this study was 15 mSv WLM⁻¹ which was in the range of the values of 6-20 mSv WLM⁻¹ reported by other investigators. Implementing the average indoor radon concentration in Jordan, the annual effective doses were calculated to be 4.1 mSv y⁻¹ and 0.08 mSv y⁻¹ due to the inhalation of radon progeny and radon gas, respectively. The total annual effective dose estimated for Jordanian population was 4.2 mSv y⁻¹. This high annual effective dose calculated by the dosimetric approach using ICRP biokinetic and dosimetric models resulted in an increase of a factor of two in comparison to the value by epidemiological study. This phenomenon was presented by the ICRP in its new published statement on radon.     

Implications of simultaneity in a physical damage function

A modeler must often rely on highly simplified representations of complex physical systems when analyzing associated economic issues. Herein, we consider a management problem in which a bioeconomic system exhibits simultaneity in processes governing productivity and damage. In this case, it may benefit the producer to sacrifice productivity to reduce the costs associated with increased damage. We specify empirically a structural damage relationship that explains the biological process by which an invasive species damages a host and estimate the structural model and its reduced form with an exceptional dataset on infestation of olives by the olive fruit fly. We contrast the results of these models with the approach typically taken in the economic literature, which expresses damage as a function of pest density. The population-based approach introduces significantly greater bias into the individual grower's choice of damage-control inputs than estimates based on the structural model.     

Characterizing sources of nitrate leaching from an irrigated dairy farm in Merced County, California

Dairy farms comprise a complex landscape of groundwater pollution sources. The objective of our work is to develop a method to quantify nitrate leaching to shallow groundwater from different management units at dairy farms. Total nitrate loads are determined by the sequential calibration of a sub-regional scale and a farm-scale three-dimensional groundwater flow and transport model using observations at different spatial scales. These observations include local measurements of groundwater heads and nitrate concentrations in an extensive monitoring well network, providing data at a scale of a few meters and measurements of discharge rates and nitrate concentrations in a tile-drain network, providing data integrated across multiple farms. The various measurement scales are different from the spatial scales of the calibration parameters, which are the recharge and nitrogen leaching rates from individual management units. The calibration procedure offers a conceptual framework for using field measurements at different spatial scales to estimate recharge N concentrations at the management unit scale. It provides a map of spatially varying dairy farming impact on groundwater nitrogen. The method is applied to a dairy farm located in a relatively vulnerable hydrogeologic region in California. Potential sources within the dairy farm are divided into three categories, representing different manure management units: animal exercise yards and feeding areas (corrals), liquid manure holding ponds, and manure irrigated forage fields. Estimated average nitrogen leaching is 872 kg/ha/year, 807 kg/ha/year and 486 kg/ha/year for corrals, ponds and fields respectively. Results are applied to evaluate the accuracy of nitrogen mass balances often used by regulatory agencies to assess groundwater impacts. Calibrated leaching rates compare favorably to field and farm scale nitrogen mass balances. These data and interpretations provide a basis for developing improved management strategies.     

Kinetics of carbon dioxide absorption into mixed aqueous solutions of MDEA and MEA using a laminar jet apparatus and a numerically solved 2D absorption rate/kinetics model

New comprehensive numerically solved 1D and 2D absorption rate/kinetics models have been developed, for the first time, to interpret the experimental kinetic data obtained with a laminar jet apparatus for the absorption of carbon dioxide (CO₂) in CO₂ loaded mixed solutions of mixed amine system of methyldiethanolamine (MDEA) and monoethanolamine (MEA). Three MDEA/MEA weight ratios ranging from 27/03 to 23/07, over a concentration range of 2.316–1.996 kmol/m³ for MDEA and of 0.490–1.147 kmol/m³ for MEA were studied. The models take into account the coupling between chemical equilibrium, mass transfer, and the chemical kinetics of all possible chemical reactions involved in the CO₂ reaction with MDEA/MEA solvent. The partial differential equations of the 1D model were solved by two numerical techniques; the finite difference method (FDM) based on in-house coded Barakat–Clark scheme and the finite element method (FEM) based on COMSOL software. The FEM comprehensive model was then used to solve the set of partial differential equations in the 2D cylindrical coordinate system setting. Both FDM and FEM produced very accurate results for both the 1D and 2D models, which were much better than our previously published simplified model. The reaction rate constant obtained for MEA blended into MDEA at 298–333 K was k_MEA = 5.127 × 10⁸ exp(−3373.8/T). In addition, the 2D model, for the first time, has provided the concentration profiles of all the species in both the radial and axial directions of the laminar jet, thereby enabling an understanding of the correct sequence in which the reaction steps involved in the reactive absorption of CO₂ in aqueous mixed amines occur.