A Multiple Training Image Approach for Spatial Modeling of Geologic Domains

Characterization of complex geological features and patterns remains one of the most challenging tasks in geostatistics. Multiple point statistics (MPS) simulation offers an alternative to accomplish this aim by going beyond classical two-point statistics. Reproduction of features in the final realizations is achieved by borrowing high-order spatial statistics from a training image. Most MPS algorithms use one training image at a time chosen by the geomodeler. This paper proposes the use of multiple training images simultaneously for spatial modeling through a scheme of data integration for conditional probabilities known as a linear opinion pool. The training images (TIs) are based on the available information and not on conceptual geological models; one image comes from modeling the categories by a deterministic approach and another comes from the application of conventional sequential indicator simulation. The first is too continuous and the second too random. The mixing of TIs requires weights for each of them. A methodology for calibrating the weights based on the available drillholes is proposed. A measure of multipoint entropy along the drillholes is matched by the combination of the two TIs. The proposed methodology reproduces geologic features from both TIs with the correct amount of continuity and variability. There is no need for a conceptual training image from another modeling technique; the data-driven TIs permit a robust inference of spatial structure from reasonably spaced drillhole data.     

Projections of severe droughts in future climate in Southeast Brazil: a case study in Southern Minas Gerais State,Brazil

Theoretical and Applied Climatology - South of Minas Gerais state, in Southeast Brazil, is known for the coffee crop production (more than 30% of country’s production) and hydroelectricity...     

Analysis of aquifer vulnerability and water quality using SINTACS and geographic weighted regression

Aquifer vulnerability and water quality were assessed in the Central Valleys of Oaxaca (Mexico) using the SINTACS method, based on a geographic information system. SINTACS layers were prepared using data such as climate (rainfall and temperature), water table, hydraulic conductivity, geology, soil type and topographic model. Maps for water quality index (WQI), contamination index and pollution sources index (PSI) were also obtained by this work. Groundwater quality in the Central Valleys may be affected by two factors, those with an anthropogenic origin and those with natural origin. High vulnerability values are located in the valleys of the basin, where granular sediments are exposed. Low vulnerability values are distributed in the basin??s ranges, where metamorphic rocks are found. Given that many of the zones with the highest groundwater vulnerability values correspond to zones with the greatest PSI values, there is great risk of groundwater contamination for the area of study because external (indicated by PSI) and internal (indicated by SINTACS) factors that cause pollution can be frequently observed in the same place. Geographic weighted regression (GWR) is used to test the dependency between WQI as dependent variable and SINTACS, PSI, Urban localities, Agriculture, Pastures and Rivers as predictors. The results indicate the non-stationary behavior of the dependent variable with respect to the predictors. While the obtained GWR models used to model WQI cannot be used in practical situations to predict the behavior of said variable, they can be used to estimate the degree to which the predictors influence the variable of interest.     

Iron(II) oxide determination in rocks and minerals

The determination of FeO of geologic materials by modern instrumental methods (such as atomic absorption spectroscopy (AAS), inductively coupled plasma atomic emission spectroscopy (ICP-AES), X-ray fluorescence (XRF), etc.) cannot distinguish between different oxidation states of elements. In many cases, the oxidation state of Fe has to be known in order to perform several chemical calculations (norms, etc.) and discuss the reactions that occur during weathering, hydrothermal alteration and other processes. A modified Wilson method is proposed, giving reproducible results in a much shorter time than the classical method. Back-titration with potassium dichromate and an Fe(II) and ammonia sulphate solution is used, after dissolution of the sample powder in a heated HF/H₃PO₄ mixture and an ammonium vanadate solution. This modified method, tested with several international reference materials, gives reliable results, equivalent to the ones cited in the literature for the reference materials.     

Estimation of transit trip production using Factorial Kriging with External Drift: an aggregated data case study

Studies in transportation planning routinely use data in which location attributes are an important source of information. Thus, using spatial attributes in urban travel forecasting models seems reasonable. The main objective of this paper is to estimate transit trip production using Factorial Kriging with External Drift (FKED) through an aggregated data case study of Traffic Analysis Zones in São Paulo city, Brazil. The method consists of a sequential application of Principal Components Analysis (PCA) and Kriging with External Drift (KED). The traditional Linear Regression (LR) model was adopted with the aim of validating the proposed method. The results show that PCA summarizes and combines 23 socioeconomic variables using 4 components. The first component is introduced in KED, as secondary information, to estimate transit trip production by public transport in geographic coordinates where there is no prior knowledge of the values. Cross-validation for the FKED model presented high values of the correlation coefficient between estimated and observed values. Moreover, low error values were observed. The accuracy of the LR model was similar to FKED. However, the proposed method is able to map the transit trip production in several geographical coordinates of non-sampled values.     

Multiparametric Traveltime Inversion

In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method.     

Late- to post-orogenic basins of the Pan-African – Brasiliano collision orogen in southern Africa and southern Brazil

Late- to post-orogenic basins formed on both sides of the Pan-African – Brasiliano orogen when the Congo and Kalahari Cratons collided with the Rio de la Plata Craton during the formation of western Gondwana. Trace fossil evidence and radiometric age dating indicate that deposits on both sides are coeval and span the Cambrian–Precambrian boundary. A peripheral foreland basin, the Nama Basin, developed on the subducting southern African plate. Lower, craton-derived fluviomarine clastics are overlain by marine platform carbonates and deltaic flysch derived in part from the rising subduction complex along the northern (Damara Belt) and western (Gariep Belt) orogenic margins. Rare, thin volcanic ash layers (tuffs and cherts) are present. Upper sediments consist of unconformable red molasse related to collisional orogenesis. Orogenic loading from the north and west led to crustal flexure and the formation of a remnant ocean that drained to the south and closed progressively from north to south. During final collision SE-, E- and NE-verging nappes overrode the active basin margins. Although younger than most of the post-orogenic magmatism, its setting on the cratonic edge of the subducting plate precluded marked volcanism or granitic intrusion, the only exception being the youngest intrusions of the Kuboos-Bremen Suite dated at 521±6 Ma to 491±8 Ma. Two foreland-type basins, perhaps faulted remnants of a much larger NE–SW elongated retroarc foreland basin, are found west of the Dom Feliciano Belt on the edge of the Rio de la Plata Craton in southern Brazil. In the southern Camaqua Basin, basal fluvial deposits are followed by cyclical marine and coarsening-up deltaic deposits with a notable volcanic and volcaniclastic component. This lower deformed succession, comprising mainly red beds, contain stratabound Cu and Pb–Zn deposits and is overlain unconformably by a fluviodeltaic to aeolian succession of sandstones and conglomerates (with minor andesitic volcanics), derived primarily from an eastern orogenic source and showing southerly longitudinal transport. In the northern Itajaí Basin, sediments range from basal fluvial and platform sediments to fining-up submarine fan and turbidite deposits with intercalated acid tuffs. The Brazilian basins had faulted margins off which alluvial fans were shed. They also overlie parts of the Ribeira Belt. Thrust deformation along the orogenic margin bordering the Dom Feliciano Belt was directed westward in the Camaqua and Itajaí basins, but reactivated strike-slip and normal faults are also present. Late- to post-orogenic granitoids and volcanics of the Dom Feliciano Belt, ranging in age from 568±6 Ma to 529±4 Ma, occur in the foreland basins and are geochemically related to some of the synsedimentary volcanics.     

Numerical homogenization-based seismic assessment of an English-bond masonry prototype: Structural level application

A multiscale strategy is evaluated at a structural level for the analysis of unreinforced masonry structures. The mechanical characterization of the masonry is deduced from homogenization-based micro-scale finite element (FE) models. The derived data are here employed at a structural level via a discrete FE model. The discrete FE model is composed of quadrilateral rigid plates interconnected through vertical and horizontal interfaces. On the interfaces, between adjoining discrete elements, a model that accounts for the in- and out-of-plane behavior of masonry, with damage and plasticity, is adopted. Such interfaces represent the material pre- and post-peak regimes, its orthotropy, and, depending on the micro-model assumed, account by three-dimensional shear effects that are especially important for multi-leaf walls and complex regular textures. The discrete model has been implemented in an advanced structural analysis software where powerful built-in features as the arc-length method, line-search algorithm, and implicit or explicit solver schemes are available. The multi-scale model is applied for the dynamic study of a small English-bond masonry house prototype subjected to a series of consecutive earthquake records. Detailed comparisons between the experimental and numerical data are presented, including the results obtained through a continuous total strain rotating crack model. Quasi-static and dynamic analyses are conducted. Results demonstrate that when enough experimental information is available on the masonry components under tension, shear, and compression regimes, the approach predicts well the seismic structural response in terms of time-history displacements, seismic capacity, and damage patterns. The required computational cost (CPU time) is very attractive.     

Water Availability of Sub-Basins with Restrictions on Use under Climate Change Scenarios in Brazil

Water Resources - The Queima Pé Stream Basin (QPSB) and Ararão River Basin (ARB) in central Brazil have had recurrent problems with water scarcity. This may get worse under anthropogenic...     

Camacho inlet southern Brazil - Instability causes and user conflicts

Inlet stabilization is a common practice in Brazil, as it is worldwide. In this paper we present results of the analysis of the stabilization practices of the Camacho Inlet, located in Santa Catarina State, southern Brazil. This inlet is maintained intermittently, in accordance with a regime designed to benefit both local fishermen and rice farmers. Fishermen require this inlet to be open as a navigation channel during the fishing season; rice farmers needs it closed during the irrigation period to allow the pumping of fresh water from the Congonhas River which drains into the Camacho’s backbarrier lagoon. An additional complication arises from the tendency of this inlet to close naturally due to infilling by dune migration. Aerial images spanning a 65 year time period (1938-2003) were analyzed, demonstrating the variability in the inlet’s width between 0 (closed) and 360 m. The analysis of inlet stability and user needs indicate that an appropriate solution is the implementation of an annual inlet dredging program in the austral autumn. This would allow for the inlet to remain open from May to September and then close naturally for the rest of the year, providing for a backbarrier fresh water reservoir during the irrigation period.