Data assimilation methods for estimating a heterogeneous conductivity field by assimilating transient solute transport data via ensemble Kalman filter

An ensemble Kalman filter (EnKF) is developed to identify a hydraulic conductivity distribution in a heterogeneous medium by assimilating solute concentration measurements of solute transport in the field with a steady‐state flow. A synthetic case with the mixed Neumann/Dirichlet boundary conditions is designed to investigate the capacity of the data assimilation methods to identify a conductivity distribution. The developed method is demonstrated in 2‐D transient solute transport with two different initial instant solute injection areas. The influences of the observation error and model error on the updated results are considered in this study. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating solute concentration measurements. The larger area of the initial distribution and the more observed data obtained, the better the calculation results. When the standard deviation of the observation error varies from 1% to 30% of the solute concentration measurements, the simulated results by the data assimilation method do not change much, which indicates that assimilation results are not very sensitive to the standard deviation of the observation error in this study. When the inflation factor is more than 1.0 to enlarge the model error by increasing the forecast error covariance matrix, the updated results of the hydraulic conductivity by the data assimilation method are not good at all. Copyright © 2012 John Wiley & Sons, Ltd.     

Using data assimilation method to calibrate a heterogeneous conductivity field conditioning on transient flow test data

A data assimilation method is developed to calibrate a heterogeneous hydraulic conductivity field conditioning on transient pumping test data. The ensemble Kalman filter (EnKF) approach is used to update model parameters such as hydraulic conductivity and model variables such as hydraulic head using available data. A synthetical two-dimensional flow case is used to assess the capability of the EnKF method to calibrate a heterogeneous conductivity field by assimilating transient flow data from observation wells under different hydraulic boundary conditions. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating continuous hydraulic head measurements and the hydraulic boundary condition will significantly affect the simulation results. For our cases, after a few data assimilation steps, the assimilated conductivity field with four Neumann boundaries matches the real field well while the assimilated conductivity field with mixed Dirichlet and Neumann boundaries does not. We found in our cases that the ensemble size should be 300 or larger for the numerical simulation. The number and the locations of the observation wells will significantly affect the hydraulic conductivity field calibration.     

Managing Groundwater to Ensure Ecosystem Function

Groundwater is a critical resource not only for human communities but also for many terrestrial, riparian, and aquatic ecosystems and species. Yet groundwater planning and management decisions frequently ignore or inadequately address the needs of these natural systems. As a consequence, ecosystems dependent on groundwater have been threatened, degraded, or eliminated, especially in arid regions. There is growing acknowledgment that governmental protections for these ecological resources are necessary, but current legal, regulatory and voluntary provisions are often inadequate. Groundwater management premised on “safe yield,” which aims to balance human withdrawals with natural recharge rates, typically provides little to no consideration for water needed by ecosystems. Alternatively, the “sustainable yield” concept aims to integrate social, economic and environmental needs for groundwater, but the complexity of groundwater systems creates substantial uncertainty about the impact that current or future groundwater withdrawals will have on ecosystems. Regardless of the legal or regulatory framework, guidance is needed to help ensure environmental water needs will be met, especially in the face of pressure to increase human uses of groundwater resources. In this paper, we describe minimum provisions for planning, managing, and monitoring groundwater that collectively can lower the risk of harm to groundwater-dependent ecosystems and species, with a special emphasis on arid systems, where ecosystems and species may be especially reliant upon and sensitive to groundwater dynamics.     

AGU新出版地球系统模拟杂志

美国地球物理联合会（AGU）很高兴地宣布,AGU不久将成为《地球系统模拟进展杂志》（Journal of Advances in Modeling Earth Systems,JAMES）新的出版机构。JAMES是一份同行评议、开放获取、全电子化的杂志,通过提供高质量的学术研究论文推动地球系统模拟科学的发展。     

Family or Enterprise? What shapes the business structures of Australian farming?

Australian farmers navigate their contemporary circumstances through the use of different business and legal arrangements that are shaped by the commercial realities of farming and the aspirations of farm-owning households. In posing the question ‘Family or Enterprise?’, this paper examines the extent to which various household and farm business indicators are associated with different forms of farm ownership, namely sole proprietorships, partnerships, trusts and companies. Results from a postal survey of farm enterprises in Victoria, Australia suggest that both household and enterprise factors contribute to the business structure used, although the strongest determinants appear to be those factors that are less well understood in the rural geographical and sociological literature: household composition, farmer age and farm size. Greater scrutiny of the business instruments deployed by farmers to manage family and enterprise pressures should inform expectations of the fate of family farming in advanced financialised economies.     

SN 1994W: an interacting supernova or two interacting shells?

We present a multi-epoch quantitative spectroscopic analysis of the Type IIn supernova (Type IIn SN) 1994W, an event interpreted by Chugai et al. as stemming from the interaction between the ejecta of a SN and a  0.4 M_⊙  circumstellar shell ejected 1.5 yr before core collapse. During the brightening phase, our models suggest that the source of optical radiation is not unique, perhaps associated with an inner optically thick cold dense shell and outer optically thin shocked material. During the fading phase, our models support a single source of radiation, an hydrogen-rich optically thick layer with a near-constant temperature of ∼7000 K that recedes from a radius of  4.3 × 10¹⁵  at a peak to  2.3 × 10¹⁵ cm  40 d later. We reproduce the hybrid narrow-core broad-wing line profile shapes of SN 1994W at all times, invoking an optically thick photosphere exclusively (i.e. without any external optically thick shell). In SN 1994W, slow expansion makes scattering with thermal electrons a key escape mechanism for photons trapped in optically thick line cores, and allows the resulting broad incoherent electron-scattering wings to be seen around narrow-line cores. In SNe with larger expansion velocities, the thermal broadening due to incoherent scattering is masked by the broad profile and the dominant frequency redshift occasioned by bulk motions. Given the absence of broad lines at all times and the very low ⁵⁶Ni yields, we speculate whether SN 1994W could have resulted from an interaction between two ejected shells without core collapse. The high conversion efficiency of kinetic to thermal energy may not require a SN-like energy budget for SN1994W.     

Surface water resources for prehistoric peoples in western Papaguer?́a of the North American south-west

Water for personal use was a limiting factor for the prehistoric and historic indigenous peoples who inhabited the extremely arid North American region of south-western Arizona and north-western Sonora, known as Papaguer?́a. This paper lists and quantifies the perennial and intermittent water resources available to them. In the 16,824 km²study area of the Sonoran Desert, there were 128 waterholes, most of them in the mountains and, when full, these offered at least 3,099,141 l of water. At today's precipitation levels, virtually all points of the region were accessible to inhabitants within a 40 km radius of one or more water sources.     

A primer on upscaling tools for porous media

Over the last few decades a number of powerful approaches have been developed to intelligently reduce the number of degrees of freedom in very complex heterogeneous environs, e.g. mathematical homogenization, mixture and hybrid mixture theory, spatial averaging, moment methods, central limit or Martingale methods, stochastic-convective approaches, various other Eulerian and Lagrangian perturbation schemes, projection operators, renormalization group techniques, variational approaches, space transformational methods, continuous time random walks, and etc. In this article we briefly review many of these approaches as applied to specific examples in the hydrologic sciences.     

A New Magnetically Active Binary System Discovered in Yunnan-Hong Kong Wide Field Survey

We present a newly discovered magnetically active binary system detected by Yunnan-Hong Kong wide field survey, with an orbital period of 0.60286 days. Two color photometry for the system was performed using the 1 m Cassegrain telescope of Yunnan Observatories with its CCD (Charge-Coupled Device) camera. In the observed light curves, there are clearly different light maxima existed in the out-of-eclipse regions. We made spectroscopic observations for the binary system using the 2.4 m telescope and YFOSC (Yunnan Faint Object Spectrograph and Camera) of Lijiang station of Yunnan Observatories, China. The radial velocity curve was derived for primary star of the binary system. The primary star exhibited strong chromospheric activity, which confirms that the distortion of the light curves results from the starspot activity on the primary star. Through analyzing the light curves and RV (Radial Velocity) curve mentioned above by means of the W-D (Wilson-Devinney) code, orbital parameters and starspot configuration of the binary system are obtained. Finally, we have discussed the properties of the binary system, and given the prospects on the future work.