Histogram and variogram inference in the multigaussian model

Several iterative algorithms are proposed to improve the histogram and variogram inference in the framework of the multigaussian model. The starting point is the variogram obtained after a traditional normal score transform. The subsequent step consists in simulating many sets of gaussian values with this variogram at the data locations, so that the ranking of the original values is honored. The expected gaussian transformation and the expected variogram are computed by an averaging operation over the simulated datasets. The variogram model is then updated and the procedure is repeated until convergence. Such an iterative algorithm can adapt to the case of tied data and despike the histogram. Two additional issues are also examined, referred to the modeling of the empirical transformation function and to the optimal pair weighting when computing the sample variogram.     

An approach for quantifying geomorphological impacts for EIA of transportation infrastructures: a case study in northern Spain

A methodological proposal for the assessment of impacts due to linear infrastructures such as motorways, railways, etc. is presented. The approach proposed includes a series of specific issues to be addressed for each geomorphological feature analysed—both ‘static’ and ‘dynamic’—as well as a series of steps to be followed in the process.Geomorphic characteristics potentially affected were initially identified on the basis of a conceptual activities/impacts model that helps to single out geomorphic impacts related to environmental concerns for the area. The following issues were addressed for each individual impact: nature of potential effects; indicators that can be used to measure impacts; criteria of ‘geomorphologic performance’; procedure for measurement/prediction of changes; translation of geomorphologic impacts into significant terms from the viewpoint of human concerns; possible mitigation and/or compensation measures.The procedure has been applied to a case study corresponding to a new motorway in the Basque Country, northern Spain. Geomorphological impacts considered in this analysis included: (1) consumable resources; (2) sites of geomorphological interest; (3) land units with high potential for use, high productivity or value for conservation; (4) visual landscape; (5) slope instability processes. The procedure has been designed for implementation in a Geographic Information System (GIS) environment. Details are given on the application of the method to each individual impact analysed and results are presented in both numerical and map form.Impacts assessed were initially expressed by means of heterogeneous magnitudes, depending on the geomorphological feature considered. Those geomorphological impacts were then translated into significant terms and homogeneous magnitudes. Integration was carried out on the basis of impact values thus obtained. Final integrated results were also expressed in numerical and map form.The method proposed enables comparison of alternatives as well as ‘prediction’ and assessment of impacts in terms directly related to geomorphic characteristics. It also facilitates the expression of those impacts in terms that allow integration with other types of environmental impacts.     

Groundwater flow regime under natural conditions as inferred from past evidence and contemporary field observations in a semi-arid basin: Cuenca de la Independencia, Guanajuato, México

This study was carried out in the Cuenca de la Independencia, a semi-arid basin in Central Mexico. The objective is to describe the main features of a groundwater flow regime under natural conditions, based on groundwater discharge manifestations. Information obtained from paleoecological, paleontological, archaeological and historical data suggests that, prior to heavy development (starting in the 1950s), the hydrogeologic regime was characterized by a larger groundwater availability in a more humid and colder climate. Manifestations associated to groundwater discharges are springs, lagoons, wetlands, saline soils, chalcedony deposits, phreatophytes, thermalism, and artesianism. The different types of manifestations and their position in the basin indicate the influence of groundwater flow systems hierarchically nested, forming concentric zones at ground level. The groundwater flow regime corresponds to a classical gravity-induced flow system with generation of local, intermediate and regional patterns. Integrating several types of data to establish the flow geometry and its dynamics has proven a useful tool to increase understanding of the original groundwater regimes. This approach can also be applied in other over-exploited semi-arid basins.     

Curved structures and interference fold patterns associated with lateral ramps in the Eastern Cordillera, Central Andes of Argentina

The conspicuous curved structures located at the eastern front of the Eastern Cordillera between 25° and 26° south latitude is coincident with the salient recognized as the El Crestón arc. Major oblique strike-slip faults associated with these strongly curved structures were interpreted as lateral ramps of an eastward displaced thrust sheet. The displacement along these oblique lateral ramps generated the local N–S stress components responsible for the complex hanging wall deformation. Accompanying each lateral ramp, there are two belts of strong oblique fault and folding: the upper Juramento River valley area and El Brete area.On both margins of the Juramento River upper valley, there is extensive map-scale evidence of complex deformation above an oblique ramp. The N–S striking folds originated during Pliocene Andean orogeny were subsequently or simultaneously folded by E–W oriented folds. The lateral ramps delimiting the thrust sheet coincident with the El Crestón arc salient are strike-slip faults emplaced in the abrupt transitions between thick strata forming the salient and thin strata outside of it. El Crestón arc is a salient related to the pre-deformational Cretaceous rift geometry, which developed over a portion of this basin (Metán depocenter) that was initially thicker. The displacement along the northern lateral ramp is sinistral, whereas it is dextral in the southern ramp. The southern end of the Eastern Cordillera of Argentina shows a particular structure reflecting a pronounced along strike variations related to the pre-deformational sedimentary thickness of the Cretaceous basin.     

Current states of stress in the northern Andes as indicated by focal mechanisms of earthquakes

Systematic inversion of double couple focal mechanisms of shallow earthquakes in the northern Andes reveals relatively homogeneous patterns of crustal stress in three main regions. The first region, presently under the influence of the Caribbean plate, includes the northern segment of the Eastern Cordillera of Colombia and the western flank of the Central Cordillera (north of 4°N). It is characterized by WNW–ESE compression of dominantly reverse type that deflects to NW–SE in the Merida Andes of Venezuela, where it becomes mainly strike–slip in type. A major bend of the Eastern thrust front of the Eastern Cordillera, near its junction with the Merida Andes, coincides with a local deflection of the stress regime (SW–NE compression), suggesting local accommodation of the thrust belt to a rigid indenter in this area. The second region includes the SW Pacific coast of Colombia and Ecuador, currently under the influence of the Nazca plate. In this area, approximately E–W compression is mainly reverse in type. It deflects to WSW–ENE in the northern Andes south of 4°N, where it is accommodated by right-lateral displacement of the Romeral fault complex and the Eastern front of the northern Andes. The third, and most complex, region is the area of the triple junction between the South American, Nazca and Caribbean plates. It reveals two major stress regimes, both mainly strike–slip in type. The first regime involves SW–NE compression related to the interaction between the Nazca and Caribbean plates and the Panama micro-plate, typically accommodated in an E–W left-lateral shear zone. The second regime involves NW–SE compression, mainly related to the interaction between the Caribbean plate and the North Andes block which induces left-lateral displacement on the Uramita and Romeral faults north of 4°N.Deep seismicity (about 150–170 km) concentrates in the Bucaramanga nest and Cauca Valley areas. The inversion reveals a rather homogeneous attitude of the minimum stress axis, which dips towards the E. This extension is consistent with the present plunge of the Nazca and Caribbean slabs, suggesting that a broken slab may be torn under gravitational stresses in the Bucaramanga nest. This model is compatible with current blocking of the subduction in the western northern Andes, inhibiting the eastward displacement of slabs, which are forced to break and sink in to the asthenosphere under their own weight.     

Flexural and gravity modelling of the Mérida Andes and Barinas–Apure Basin, Western Venezuela

The kinematic evolution of the Barinas–Apure Basin and the southern Mérida Andes from Lower Miocene to the Present is numerically modelled using flexural isostatic theory and geophysical and geological data. Two published regional transects are used to build up a reference section, which is then used to constrain important parameters (e.g. shortenings and sedimentary thicknesses) for the flexural modelling. To control the location of the main fault system in the flexural model earthquake information is also used. The estimated flexural elastic thickness of the South American lithosphere beneath the Barinas–Apure Basin and the Mérida Andes Range is 25 km. The value for the final total shortening is 60 km. The flexural isostatic model shows that the Andean uplift has caused the South American lithosphere subsidence and the development of the Barinas–Apure Basin.In addition, gravity modelling was used to understand deep crustal features that could not be predicted by flexural theory. Consequently, the best-fit flexural model is used to build a gravity model across the Mérida Andes and the Barinas–Apure Basin preserving the best-controlled structural features from the flexural modelling (e.g. basin wavelength and depth) and slightly changing the main bodies density values and deep crustal structures. The final gravity model is intended to be representative of the major features affecting the gravity field in the study area. The predicted morphology in the lower crustal level of the final gravity model favours the hypothesis of a present delamination or megathrust of the Maracaibo crust over the South American Shield. This process would use the Conrad discontinuity as a main detachment surface within an incipient NW dipping continental subduction.     

Soil liquefaction during the Arequipa Mw 8.4, June 23, 2001 earthquake, southern coastal Peru

The Arequipa June 23, 2001, earthquake with a moment magnitude of Mw 8.4 struck southern Peru, northern Chile and western Bolivia. This shallow (29 km deep) interplate event, occurring in the coupled zone of the Nazca subduction next to the southeast of the subducting Nazca ridge, triggered very localized but widely outspread soil liquefaction. Although sand blows and lateral spreading of river banks and road bridge abutments were observed 390 km away from the epicenter in the southeast direction (nearing the town of Tacna, close to the Chile border), liquefaction features were only observed in major river valleys and delta and coastal plains in the meizoseismal area. This was strongly controlled by the aridity along the coastal strip of Southern Peru. From the sand blow distribution along the coastal area, a first relationship of isolated sand blow diameter versus epicentral distance for a single event is ever proposed. The most significant outcome from this liquefaction field reconnaissance is that energy propagation during the main June 23, 2001, event is further supported by the distribution and size of the isolated sand blows in the meizoseismal area. The sand blows are larger to the southeast of the epicenter than its northwestern equivalents. This can be stated in other words as well. The area affected by liquefaction to the northwest is less spread out than to the southeast. Implications of these results in future paleoliquefaction investigations for earthquake magnitude and epicentral determinations are extremely important. In cases of highly asymmetrical distribution of liquefaction features such as this one, where rupture propagation tends to be mono-directional, it can be reliably determined an epicentral distance (between earthquake and liquefaction evidence) and an earthquake magnitude only if the largest sand blow is found. Therefore, magnitude estimation using this uneven liquefaction occurrence will surely lead to underrating if only the shortest side of the meizoseismal area is unluckily studied, which can eventually be the only part exhibiting liquefaction evidence, depending on the earthquake location and the distribution of liquefaction-prone environments.     

Geochemical characteristics of mineral air particles emitted inside a clay atomisation plant which introduces waste recycling processes

The objective of this work is to assess the concentrations of three factions of air particles (settable particles, TSP and PM10) and the levels of several toxic elements in a clay atomisation industry through aerosol sampling at several points inside an industrial plant. Mechanical activities, which produce diffuse emissions, are the main process of discharge of particles in both indoor and outdoor workplace environments in the atomisation plant. The levels of As, Cd, Pb, Zn, Ba and Ni increase in the zones with higher concentrations of particles and lower ventilation. The concentrations of As and F are not influenced by the recycling processes. The levels of Cd and Pb do not show great enrichment in air particles collected inside the atomisation plant although the content of both elements is associated with ceramic muck recycling. Finally, the content of B in waste water is mainly transferred in gaseous phase to the atmosphere during the process of drying by atomisation.     

Detection and 2d modelling of cavities using pole–dipole array

The resistivity method is often used in cave prospecting. In this paper the pole–dipole array ability to detect cavities at different depths and with different water contents is investigated. The research was performed using analogical and numerical modelling. According to the results, empty caves at a depth less than four times its diameter can be easily detected. The ability of the pole–dipole array to detect water-filled caves reaches a minimum at filling percentages of 30–50%. Overburden effects research shows that low resistivity overburden reduces the resolution capability of the array. This study shows that equivalent results can be obtained by modelling the empty caves as infinite resistivity bodies or, alternatively, as very high resistivity ones. The analysis of field data acquired in the Maciço Calcário Estremenho (Portugal) shows the practical importance of the pole–dipole array in cavities prospecting.