Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter

The ensemble Kalman filter (EnKF) is now widely used in diverse disciplines to estimate model parameters and update model states by integrating observed data. The EnKF is known to perform optimally only for multi-Gaussian distributed states and parameters. A new approach, the normal-score EnKF (NS-EnKF), has been recently proposed to handle complex aquifers with non-Gaussian distributed parameters. In this work, we aim at investigating the capacity of the NS-EnKF to identify patterns in the spatial distribution of the model parameters (hydraulic conductivities) by assimilating dynamic observations in the absence of direct measurements of the parameters themselves. In some situations, hydraulic conductivity measurements (hard data) may not be available, which requires the estimation of conductivities from indirect observations, such as piezometric heads. We show how the NS-EnKF is capable of retrieving the bimodal nature of a synthetic aquifer solely from piezometric head data. By comparison with a more standard implementation of the EnKF, the NS-EnKF gives better results with regard to histogram preservation, uncertainty assessment, and transport predictions. © 2011 International Association for Mathematical Geosciences.The authors gratefully acknowledge the financial support by the Spanish Ministry of Science and Innovation through project CGL2011-23295. The first author appreciates the financial aid from China Scholarship Council (CSC No. [2007]3020).Zhou, H.; Li, L.; Hendricks Franssen, H.; Gómez-Hernández, JJ. (2012). Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter. Mathematical Geosciences. 44(2):169-185. https://doi.org/10.1007/s11004-011-9372-3S169185442Arulampalam MS, Maskell S, Gordon N, Clapp T (2002) A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. IEEE Trans Signal Process 50(2):174–188Bertino L, Evensen G, Wackernagel H (2003) Sequential data assimilation techniques in oceanography. Int Stat Rev 71(2):223–241Burgers G, Jan van Leeuwen P, Evensen G (1998) Analysis scheme in the ensemble Kalman filter. Mon Weather Rev 126(6):1719–1724Carrera J, Neuman SP (1986b) Estimation of aquifer parameters under transient and steady state conditions: 2. Uniqueness, stability, and solution algorithms. Water Resour Res 22(2):211–227Chen Y, Zhang D (2006) Data assimilation for transient flow in geologic formations via ensemble Kalman filter. Adv Water Resour 29:1107–1122Delhomme JP (1979) Spatial variability and uncertainty in groundwater flow parameters: a geostatistical approach. Water Resour Res 15(2):269–280Evensen G (1994) Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics. J Geophys Res 99(C5):10143–10162Evensen G (2007) Data assimilation: the ensemble Kalman filter. Springer, Berlin, 279 ppFernàndez-Garcia D, Illangasekare T, Rajaram H (2005) Differences in the scale dependence of dispersivity and retardation factors estimated from forced-gradient and uniform flow tracer tests in three-dimensional physically and chemically heterogeneous porous media. Water Resour Res 41(3):W03012Gómez-Hernández JJ, Journel AG (1993) Joint sequential simulation of multi-Gaussian fields. In: Soares A (ed) Geostatistics Tróia ’92, vol 1. Kluwer Academic, Dordrecht, pp 85–94Gómez-Hernández JJ, Wen XH (1998) To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology. Adv Water Resour 21(1):47–61Gu Y, Oliver DS (2006) The ensemble Kalman filter for continuous updating of reservoir simulation models. J Energy Resour Technol 128:79–87Harbaugh AW, Banta ER, Hill MC, McDonald MG (2000) MODFLOW-2000, the U.S. geological survey modular ground-water model—user guide to modularization concepts and the ground-water flow process. Tech rep. Open-File Report 00-92, U.S. Department of the Interior, U.S. Geological Survey. Reston, Virginia, 121 ppHendricks Franssen HJ, Kinzelbach W (2008) Real-time groundwater flow modeling with the Ensemble Kalman Filter: joint estimation for states and parameters and the filter inbreeding problem. Water Resour Res 44:W09408Hendricks Franssen HJ, Kinzelbach W (2009) Ensemble Kalman filtering versus sequential self-calibration for inverse modelling of dynamic groundwater flow systems. J Hydrol 365(3–4):261–274Houtekamer PL, Mitchell HL (2001) A sequential ensemble Kalman filter for atmospheric data assimilation. Mon Weather Rev 129:123–137Journel AG, Deutsch CV (1993) Entropy and spatial disorder. Math Geol 25(3):329–355Li L, Zhou H, Gómez-Hernández JJ (2011a) A comparative study of three-dimensional hydraulic conductivity upscaling at the macrodispersion experiment (MADE) site, Columbus air force base, Mississippi (USA). J Hydrol 404(3–4):278–293Li L, Zhou H, Gómez-Hernández JJ (2011b) Transport upscaling using multi-rate mass transfer in three-dimensional highly heterogeneous porous media. Adv Water Resour 34(4):478–489Moradkhani H, Sorooshian S, Gupta HV, Houser PR (2005) Dual state-parameter estimation of hydrological models using ensemble Kalman filter. Adv Water Resour 28:135–147Naevdal G, Johnsen L, Aanonsen S, Vefring E (Mar. 2005) Reservoir monitoring and continuous model updating using ensemble Kalman filter. SPE J 10(1):66–74Pardo-Igúzquiza E, Dowd PA (2003) CONNEC3D: a computer program for connectivity analysis of 3D random set models. Comput Geosci 29:775–785Schöniger A, Nowak W, Hendricks Franssen HJ (2011) Parameter estimation by ensemble Kalman filters with transformed data: approach and application to hydraulic tomography. Water Resour Res (submitted)Simon E, Bertino L (2009) Application of the Gaussian anamorphosis to assimilation in a 3-D coupled physical-ecosystem model of the North Atlantic with the EnKF: a twin experiment. Ocean Sci 5:495–510Stauffer D, Aharony A (1994) Introduction to percolation theory. Taylor and Francis, London. 181 ppStrébelle S 2000. Sequential simulation drawing structures from training images. PhD thesis, Stanford University. 187 ppStrebelle S (2002) Conditional simulation of complex geological structures using multiple-point statistics. Math Geol 34(1):1–21Wen X, Chen W (2006) Real-time reservoir model updating using ensemble Kalman filter: the confirming approach. SPE J 11(4):431–442Wen X, Chen W (2007) Some practical issues on real time reservoir updating using ensemble Kalman filter. SPE J 12(2):156–166Zhou H, Gómez-Hernández JJ, Hendricks Franssen H-J, Li L (2011) An approach to handling non-gaussianity of parameters and state variables in ensemble Kalman filtering. Adv Water Resour 34(7):844–864Zinn B, Harvey C (2003) When good statistical models of aquifer heterogeneity go bad: a comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields. Water Resour Res 39(3):105

Dental Occlusion in a Split Amazon Indigenous Population: Genetics Prevails over Environment

Background: Studies examining human and nonhuman primates have supported the hypothesis that the recent increase in the occurrence of misalignment of teeth and/or incorrect relation of dental arches, named dental malocclusion, is mainly attributed to the availability of a more processed diet and the reduced need for powerful masticatory action. For the first time on live human populations, genetic and tooth wear influences on occlusal variation were examined in a split indigenous population. The Arara-Iriri people are descendants of a single couple expelled from a larger village. In the resultant village, expansion occurred through the mating of close relatives, resulting in marked genetic cohesion with substantial genetic differences. Methodology/Principal Findings: Dental malocclusion, tooth wear and inbreeding coefficient were evaluated. The sample examined was composed of 176 individuals from both villages. Prevalence Ratio and descriptive differences in the outcomes frequency for each developmental stage of the dentition were considered. Statistical differences between the villages were examined using the chi-square test or Fisher’s exact statistic. Tooth wear and the inbreeding coefficient (F) between the villages was tested with Mann-Whitney statistics. All the statistics were performed using two-tailed distribution at p#0.05. The coefficient inbreeding (F) confirmed the frequent incestuous unions among the Arara-Iriri indigenous group. Despite the tooth wear similarities, we found a striking difference in occlusal patterns between the two Arara villages. In the original village, dental malocclusion was present in about one third of the population; whilst in the resultant village, the occurrence was almost doubled. Furthermore, the morphological characteristics of malocclusion were strongly different between the groups. Conclusions/Significance: Our findings downplay the widespread influence of tooth wear, a direct evidence of what an individual ate in the past, on occlusal variation of living human populations. They also suggest that genetics plays the most important role on dental malocclusion etiology

Environmental controls, oceanography and population dynamics of pathogens and harmful algal blooms: connecting sources to human exposure

© 2008 Author et al. This is an open access article distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Environmental Health 7 (2008): S5, doi:10.1186/1476-069X-7-S2-S5.Coupled physical-biological models are capable of linking the complex interactions between environmental factors and physical hydrodynamics to simulate the growth, toxicity and transport of infectious pathogens and harmful algal blooms (HABs). Such simulations can be used to assess and predict the impact of pathogens and HABs on human health. Given the widespread and increasing reliance of coastal communities on aquatic systems for drinking water, seafood and recreation, such predictions are critical for making informed resource management decisions. Here we identify three challenges to making this connection between pathogens/HABs and human health: predicting concentrations and toxicity; identifying the spatial and temporal scales of population and ecosystem interactions; and applying the understanding of population dynamics of pathogens/HABs to management strategies. We elaborate on the need to meet each of these challenges, describe how modeling approaches can be used and discuss strategies for moving forward in addressing these challenges.The authors acknowledge the financial support for the NSF/NIEHS and NOAA Centers for Oceans and Human Healt

The impact of co-infections on fish: a review

International audienceAbstractCo-infections are very common in nature and occur when hosts are infected by two or more different pathogens either by simultaneous or secondary infections so that two or more infectious agents are active together in the same host. Co-infections have a fundamental effect and can alter the course and the severity of different fish diseases. However, co-infection effect has still received limited scrutiny in aquatic animals like fish and available data on this subject is still scarce. The susceptibility of fish to different pathogens could be changed during mixed infections causing the appearance of sudden fish outbreaks. In this review, we focus on the synergistic and antagonistic interactions occurring during co-infections by homologous or heterologous pathogens. We present a concise summary about the present knowledge regarding co-infections in fish. More research is needed to better understand the immune response of fish during mixed infections as these could have an important impact on the development of new strategies for disease control programs and vaccination in fish

Pregnancy-related pelvic girdle pain: an update

A large number of scientists from a wide range of medical and surgical disciplines have reported on the existence and characteristics of the clinical syndrome of pelvic girdle pain during or after pregnancy. This syndrome refers to a musculoskeletal type of persistent pain localised at the anterior and/or posterior aspect of the pelvic ring. The pain may radiate across the hip joint and the thigh bones. The symptoms may begin either during the first trimester of pregnancy, at labour or even during the postpartum period. The physiological processes characterising this clinical entity remain obscure. In this review, the definition and epidemiology, as well as a proposed diagnostic algorithm and treatment options, are presented. Ongoing research is desirable to establish clear management strategies that are based on the pathophysiologic mechanisms responsible for the escalation of the syndrome's symptoms to a fraction of the population of pregnant women

Biophysical interactions in tropical agroforestry systems

sequential systems, simultaneous systems Abstract. The rate and extent to which biophysical resources are captured and utilized by the components of an agroforestry system are determined by the nature and intensity of interac-tions between the components. The net effect of these interactions is often determined by the influence of the tree component on the other component(s) and/or on the overall system, and is expressed in terms of such quantifiable responses as soil fertility changes, microclimate modification, resource (water, nutrients, and light) availability and utilization, pest and disease incidence, and allelopathy. The paper reviews such manifestations of biophysical interactions in major simultaneous (e.g., hedgerow intercropping and trees on croplands) and sequential (e.g., planted tree fallows) agroforestry systems. In hedgerow intercropping (HI), the hedge/crop interactions are dominated by soil fertility improvement and competition for growth resources. Higher crop yields in HI than in sole cropping are noted mostly in inherently fertile soils in humid and subhumid tropics, and are caused by large fertility improvement relative to the effects of competition. But, yield increases are rare in semiarid tropics and infertile acid soils because fertility improvement does not offse