A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.     

A simulation study of capillary transport,preferential retention and distribution of salts in historic sandstone buildings

Salt crystallisation is a major problem of deterioration in historic stone buildings, monuments and sculptures. The capillary rise of soil water is one of the primary sources of salts in stone structures, which evaporates leaving the salts behind. It has been noted that the spatial distribution profile of different species of salts crystallised in historic stone buildings is not homogeneous, i.e. different salts crystallise at different locations. The capillary transport and inhomogeneous spatial distribution of different salts in the porous building materials has been considered to be a result of solubility-dependent crystallisation; however, the factors responsible for this phenomenon are not clearly known. This paper aims to investigate the factors influencing the differential distribution of salts during capillary rise of soil water. In this study, the capillary transport of salts was simulated on two different sandstones—Locharbriggs, a Permo–Triassic, red sandstone and Stoke Hall, a Carboniferous, buff sandstone. The experiments were carried out under controlled environmental conditions to eliminate the possibility of evaporation-driven crystallisation of salts depending on their solubilities. The results indicate that fractionation or differential distribution of salts takes place even in the absence of evaporation and crystallisation. The sandstones exhibit properties like an ion exchange column, and ionic species present in the salt solution show differential distribution within the porous network of sandstone.     

Markov Chain Random Fields for Estimation of Categorical Variables

Multi-dimensional Markov chain conditional simulation (or interpolation) models have potential for predicting and simulating categorical variables more accurately from sample data because they can incorporate interclass relationships. This paper introduces a Markov chain random field (MCRF) theory for building one to multi-dimensional Markov chain models for conditional simulation (or interpolation). A MCRF is defined as a single spatial Markov chain that moves (or jumps) in a space, with its conditional probability distribution at each location entirely depending on its nearest known neighbors in different directions. A general solution for conditional probability distribution of a random variable in a MCRF is derived explicitly based on the Bayes’ theorem and conditional independence assumption. One to multi-dimensional Markov chain models for prediction and conditional simulation of categorical variables can be drawn from the general solution and MCRF-based multi-dimensional Markov chain models are nonlinear.     

一种改进的基于非高斯性最大化的预测反褶积算法

The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.     

REXEL: computer aided record selection for code-based seismic structural analysis

In code-based seismic design and assessment it is often allowed the use of real records as an input for nonlinear dynamic analysis. On the other hand, international seismic guidelines, concerning this issue, have been found hardly applicable by practitioners. This is related to both the difficulty in rationally relating the ground motions to the hazard at the site and the required selection criteria, which do not favor the use of real records, but rather various types of spectrum matching signals. To overcome some of these obstacles a software tool for code-based real records selection was developed. REXEL, freely available at the website of the Italian network of earthquake engineering university labs (), allows to search for suites of waveforms, currently from the European Strong-motion Database, compatible to reference spectra being either user-defined or automatically generated according to Eurocode 8 and the recently released new Italian seismic code. The selection reflects the provisions of the considered codes and others found to be important by recent research on the topic. In the paper, record selection criteria are briefly reviewed first, and then the algorithms implemented in the software are discussed. Finally, via some examples, it is shown how REXEL can effectively be a contribution to code-based real records selection for seismic structural analysis.     

Lava flow hazard at Nyiragongo Volcano,DRC

Mt. Nyiragongo is one of the most dangerous volcanoes in the world for the risk associated with the propagation of lava flows. In 2002 several vents opened along a huge system of fractures, pouring out lava which reached and destroyed a considerable part of Goma, a town of about 500,000 inhabitants on the shore of Lake Kivu. In a companion paper (Favalli et al. in Bull Volcanol, this issue, 2008) we employed numerical simulations of probable lava flow paths to evaluate the lava flow hazard on the flanks of the volcano, including the neighbouring towns of Goma (DRC) and Gisenyi (Rwanda). In this paper we use numerical simulations to investigate the possibility of significantly reducing the lava flow hazard in the city through the construction of protective barriers. These barriers are added to the DEM of the area as additional morphological elements, and their effect is evaluated by repeating numerical simulations with and without the presence of barriers. A parametric study on barrier location, size, shape and orientation led to the identification of barriers which maximize protection while minimizing their impact. This study shows that the highest hazard area corresponding to eastern Goma, which was largely destroyed by lava flows in 2002, cannot be effectively protected from future lava flows towards Lake Kivu and should be abandoned. On the contrary, the rest of the town can be sheltered from lava flows by means of two barriers that deviate or contain the lava within the East Goma sector. A proposal for the future development of the town is formulated, whereby “new” Goma is completely safe from the arrival of lava flows originating from vents outside its boundaries. The proposal minimizes the risk of further destruction in town due to future lava flows.     

The Geology and Petroleum Resources of Kansas: A Review From Alpha to Omega or From the Pleistocene to the Precambrian

The geological story of Kansas is told through the rocks that are present. It is a simple story in generalities but complex in detail. Knowing the story, gives insight into understanding the occurrence and location of possible economic valuable minerals, such as petroleum. This is a brief review of Kansas geology with respect to the known occurrence of oil and gas. Kansas is part of the Midcontinent oil province with oil having been discovered 150 years ago and commercial production commencing in 1873. Although many prospects remain in Kansas, the state has gone from the number 1 producer in the U.S. in 1916 to 8th today. Exploration for new oil and gas production therefore is going to have to be more imaginative and utilize new approaches and techniques to find the elusive petroleum. There are possibilities however for the prospector who can search diligently. Although the big fields probably have been discovered, the prospects today are deeper, in more undetectable traps, and in essentially untested places.     

Unique isoprene oxidation products demonstrate chlorine atom chemistry occurs in the Houston,Texas urban area

As part of the 2000 Texas Air Quality Study (TexAQS), we studied the isoprene oxidation process under ambient conditions to discern the presence of chlorine atom (Cl) chemistry in the Houston, Texas urban area. By measuring chloromethylbutenone (CMBO) and an isomer of chloromethylbutenal (CMBA), we clearly observed sixteen episodes of active Cl chemistry during the 24-day experiment. Estimated median Cl concentration during each of these episodes was between the detection limit of ~10² atoms cm⁻³ and 50 _{- 30} ^{+ 70} ×10⁴ 50_{ - 30}^{ + 70} \times {10^4} atoms cm⁻³. Cl concentration during all the episodes averaged 7.6 _{- 2.0} ^{+ 4.7} ×10⁴ 7.6_{ - 2.0}^{ + 4.7} \times {10^4} atoms cm⁻³ and thus amounted to less than 3% of the OH concentration during the same periods. During the episodes, the fraction of oxidation chemistry initiated by Cl ranged from 3–43% and was strongly dependent on the quantity and type of hydrocarbons present in the atmosphere. Because of its intermittent presence and low concentration, Cl is not a broadly influential oxidant in the Houston, Texas urban area.     

Preface

Climate change is one of the most challenging issues faced by human beings. Although many remarkable achievements have been made in climate change research,there is still a tremendous amount of uncertainty in     

An Analysis of the Published Mineral Resource Estimates of the Haji-Gak Iron Deposit,Afghanistan

The Haji-Gak iron deposit of eastern Bamyan Province, eastern Afghanistan, was studied extensively and resource calculations were made in the 1960s by Afghan and Russian geologists. Recalculation of the resource estimates verifies the original estimates for categories A (in-place resources known in detail), B (in-place resources known in moderate detail), and C₁ (in-place resources estimated on sparse data), totaling 110.8 Mt, or about 6% of the resources as being supportable for the methods used in the 1960s. C₂ (based on a loose exploration grid with little data) resources are based on one ore grade from one drill hole, and P₂ (prognosis) resources are based on field observations, field measurements, and an ore grade derived from averaging grades from three better sampled ore bodies. C₂ and P₂ resources are 1,659.1 Mt or about 94% of the total resources in the deposit. The vast P₂ resources have not been drilled or sampled to confirm their extent or quality. The purpose of this article is to independently evaluate the resources of the Haji-Gak iron deposit by using the available geologic and mineral resource information including geologic maps and cross sections, sampling data, and the analog-estimating techniques of the 1960s to determine the size and tenor of the deposit.