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.     

Interpretation procedures for the determination of sediment parameters from time-dependent flux inputs

A formulation for the effects of a time-dependent input flux to a sedimentation system consisting of a diffusive layer underlain by a non-diffusive medium has been developed. The intent is to provide interpretation procedures for the determination of sedimentation rate, diffusion coefficient and mixed layer thickness from observed tracer concentrations. They have been applied, for illustrative purposes, to published¹³⁷Cs and^239,240Pu measurements.     

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.     

Late Paleozoic orogeny in Alaska's Farewell terrane

Evidence is presented for a previously unrecognized late Paleozoic orogeny in two parts of Alaska's Farewell terrane, an event that has not entered into published scenarios for the assembly of Alaska. The Farewell terrane was long regarded as a piece of the early Paleozoic passive margin of western Canada, but is now thought, instead, to have lain between the Siberian and Laurentian (North American) cratons during the early Paleozoic. Evidence for a late Paleozoic orogeny comes from two belts located 100–200 km apart. In the northern belt, metamorphic rocks dated at 284–285 Ma (three ⁴⁰Ar/³⁹Ar white-mica plateau ages) provide the main evidence for orogeny. The metamorphic rocks are interpreted as part of the hinterland of a late Paleozoic mountain belt, which we name the Browns Fork orogen. In the southern belt, thick accumulations of Pennsylvanian-Permian conglomerate and sandstone provide the main evidence for orogeny. These strata are interpreted as the eroded and deformed remnants of a late Paleozoic foreland basin, which we name the Dall Basin. We suggest that the Browns Fork orogen and Dall Basin comprise a matched pair formed during collision between the Farewell terrane and rocks to the west. The colliding object is largely buried beneath Late Cretaceous flysch to the west of the Farewell terrane, but may have included parts of the so-called Innoko terrane. The late Paleozoic convergent plate boundary represented by the Browns Fork orogen likely connected with other zones of plate convergence now located in Russia, elsewhere in Alaska, and in western Canada.     

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 numerical assessment of simple airblast models of impact airbursts

Asteroids and comets 10–100 m in size that collide with Earth disrupt dramatically in the atmosphere with an explosive transfer of energy, caused by extreme air drag. Such airbursts produce a strong blastwave that radiates from the meteoroid's trajectory and can cause damage on the surface. An established technique for predicting airburst blastwave damage is to treat the airburst as a static source of energy and to extrapolate empirical results of nuclear explosion tests using an energy‐based scaling approach. Here we compare this approach to two more complex models using the iSALE shock physics code. We consider a moving‐source airburst model where the meteoroid's energy is partitioned as two‐thirds internal energy and one‐third kinetic energy at the burst altitude, and a model in which energy is deposited into the atmosphere along the meteoroid's trajectory based on the pancake model of meteoroid disruption. To justify use of the pancake model, we show that it provides a good fit to the inferred energy release of the 2013 Chelyabinsk fireball. Predicted overpressures from all three models are broadly consistent at radial distances from ground zero that exceed three times the burst height. At smaller radial distances, the moving‐source model predicts overpressures two times greater than the static‐source model, whereas the cylindrical line‐source model based on the pancake model predicts overpressures two times lower than the static‐source model. Given other uncertainties associated with airblast damage predictions, the static‐source approach provides an adequate approximation of the azimuthally averaged airblast for probabilistic hazard assessment.     

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.