Estimation of acceptable mercury intake from fish in Taiwan

The US Food and Drug Administration (USFDA) and Environmental Protection Agency (USEPA) are advising women and young children to avoid eating fish that contain high levels of mercury (Hg). However in Taiwan, the annual Black Fin Tuna Festival encourages the public to consume fish. The aim of this study was to assess fish intake in relation to the health risks of mercury exposure and calculate the acceptable and safe intake of fish in children and women of childbearing age. From the Monte Carlo simulation, based on USEPA's reference dose (0.1 microgkg(-1)d(-1)), we found that 21.6%-24.3% and 45.6%-57.4% of the daily mercury dose estimates exceeded the reference dose for typical and high-seafood consumers. The acceptable ingestion rates are <50 g d(-1) (children) and 90.8 +/- 15.7 g d(-1) (women of childbearing). Sensitivity analysis suggests that Hg concentration in fish may be a key parameter to aid governments as they offer guidance for risk management.     

A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions

This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing N(E) equilibrium reactions and a set of reactive transport equations of M-N(E) kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

Technology innovations and experience curves for nitrogen oxides control technologies

This paper reviews the regulatory history for nitrogen oxides (NOx) pollutant emissions from stationary sources, primarily in coal-fired power plants. Nitrogen dioxide (NO2) is one of the six criteria pollutants regulated by the 1970 Clean Air Act where National Ambient Air Quality Standards were established to protect public health and welfare. We use patent data to show that in the cases of Japan, Germany, and the United States, innovations in NOx control technologies did not occur until stringent government regulations were in place, thus "forcing" innovation. We also demonstrate that reductions in the capital and operation and maintenance (O&M) costs of new generations of high-efficiency NOx control technologies, selective catalytic reduction (SCR), are consistently associated with the increasing adoption of the control technology: the so-called learning-by-doing phenomena. The results show that as cumulative world coal-fired SCR capacity doubles, capital costs decline to approximately 86% and O&M costs to 58% of their original values. The observed changes in SCR technology reflect the impact of technological advance as well as other factors, such as market competition and economies of scale.     

Geostatistical analysis of soil hydrologic properties in a field plot

A 91-m transect was set up in an irrigated field near Las Cruces, New Mexico to obtain soil water tension and water content data to investigate their spatial variability. A total of 455 sampling points were monitored along a grid consisting of 91 stations placed 1 m apart by 5 depths per station (at 0.3, 0.6, 0.9, 1.2 and 1.5 m below the surface). Post-irrigation tension and wetness measurements were recorded over 45 days at 11 time periods. Soil water tension was measured with tensiometers using a hand-held pressure transducer. A neutron probe was used to obtain volumetric water content. Using the observed wetness and tension data, unsaturated hydraulic conductivity values were derived (using a cubic spline function to estimate the gradient), and an exponential model was used to fit the calculated conductivity-tension curves to obtain hydraulic conductivity parameter values. The spatial and temporal variability of wetness, tension, saturated hydraulic conductivity and pore-size distribution parameters, and texture at the 0.3-m depth were examined using geostatistical techniques.The exponential model was found to inadequately describe the hydraulic conductivity/tension relationship for the full range of tension, particularly in the tension range near saturation. The derived values of the saturated hydraulic conductivity parameter were much greater than expected and do not correspond to reasonable saturated hydraulic conductivity values.All of the soil parameters studied exhibited large spatial variability horizontally and vertically in the field. Ranges of dependence determined from semivariogram analysis over the 44-day drainage period are 3–32 m for wetness, 6–34 m for soil water tension, 5–35 m for natural log of saturated hydraulic conductivity parameters, 5–11 m for pore-size distribution parameter, and 8–24 m for percent sand, silt and clay at the 0.3 m depth. An alternate hole-effect model is suggested to describe the texture semivariograms.It was determined that the variance of volumetric water content generally increased at each depth over the measured time periods, which is consistent with certain past field studies and a stochastic analysis of unsaturated flow in heterogeneous soils. Future research is recommended relating soil texture to soil hydrologic parameters with the goal of predicting soil behavior with less extensive sampling schemes.     

MODEL DEVELOPMENT FOR CONJUNCTIVE USE STUDY OF THE SAN JACINTO BASIN,CALIFORNIA1

ABSTRACT: A two-layered confined-unconfined numerical model for flow and mass transport is developed for the San Jacinto Basin. The model structure is determined by the geological structure of the Basin and model parameters are calibrated using 20 years of historical records. The total number of historical head observations used for the flow model calibration is 1,117 and the total number of the estimated parameters is 91. The two-layered transport model is also calibrated using historical water quality records. Sensitivity analysis of the flow model shows that only 68 parameters (out of a total of 91) are relatively sensitive and reliable. However, the unreliable parameters (23 of them) are found to be insensitive and thus not significant to the prediction and management of conjunctive use of surface water and ground water. The developed flow model has been used to study the two proposed artificial recharge scenarios for the San Jacinto Basin. We have found that during a relatively dry condition, an artificial recharge rate of 80 acre-ft/day can be achieved during the recharge period October through January. However, for a relatively wet condition, only 80 percent of the proposed rate can be effectively stored in the Basin during these months.     

Karst Hydrological Processes and Grey System Model1

Hao, Yonghong, Jiaojuan Zhao, Huamin Li, Bibo Cao, Zhongtang Li, and Tian‐Chyi J. Yeh, 2012. Karst Hydrological Processes and Grey System Model. Journal of the American Water Resources Association (JAWRA) 48(4): 656‐666. DOI: 10.1111/j.1752‐1688.2012.00640.x Abstract: The karst hydrological processes are the response of karst groundwater system to precipitation. This study provided a concept model of karst hydrological processes. The hydraulic response time of spring discharge to precipitation includes the time that precipitation penetrates through the vadose zone, and the subsequent groundwater pressure wave propagates to a spring outlet. Due to heterogeneities in karst aquifers, the hydraulic response time is different in different areas. By using grey system theory, we proposed a karst hydrological model that offers a calculation of hydraulic response time, and a response model of spring discharge to precipitation. Then, we applied the models to the Liulin Springs Basin, China. In the south part of the Liulin Springs Basin, where large fields of carbonate rocks outcrop with intensive karstification, the hydraulic response time is one year. In the north, where the Ordovician karst aquifer is covered by Quaternary loess sediments, the response time is seven years. The grey system GM(1,3) response model of spring discharge to precipitation was applied in consideration of the hydraulic response time. The model calibration showed that the average error was 6.55%, and validation showed that the average error was 12.19%.