Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

GRASS GIS-embedded Decision Support Framework for Flood Simulation and Forecasting

This study presents a spatiotemporal analysis tool, called Shyska. This tool allows the simulation and prediction of flash floods in semiarid basins. Shyska has been developed by Geographical Information System (GIS)‐embedded functions, allowing the integration of hydrometeorological information from modern technologies of data acquisition in real time. A Digital Elevation Model (DEM) is used in order to obtain the relevant parameters from the integrated rainfall‐runoff models. Some of its most relevant modules and methodology employed for its development are described. Case studies in basins of south‐east Spain illustrate the applicability of the proposed techniques.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.     

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.     

Temporal variability in the maternal environment and its effect on seed size and seed quality in Flourensia cernua DC. (Asteraceae).

The biological phenomenon commonly referred to as maternal effect occurs whenever the environment exerts an influence upon the mother that is later expressed in characteristics of her offsprings. For example, the environmental conditions experienced by a mother plant during flowering and fruiting can modify the quality of her descendants (seeds), and even interrupt their development. Between-year variation in minimum temperature and the classic, between-year variation in precipitation represent an unpredictable environment for winter blooming plants in arid zones. In this study, we investigated the effect of maternal environment (temperature and precipitation) on seed size and seed quality in F. cernua, over a 5-year period. Results indicated that the proportion of filled seeds, as well as their size, increased with higher absolute minimum temperatures during seed formation, and that, to a lesser extent, both the precipitation occurred during the same period and the annual precipitation also have a positive influence on seed size. In this way, environmental conditions prevailing during seed formation exerted a strong maternal effect on the size and quality of seeds produced and, probably also, in the developmental possibilities of F. cernua seedlings. This effect was subject to important, between-year variation.     

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.     

Age, distribution, tectonics, and eustatic controls of the Paranense and Caribbean marine transgressions in southern Bolivia and Argentina

Marine transgression onto the South American continent took place at least twice in the Miocene along distinct paleogeographic corridors. The first event occurred between 15 and 13 Ma and the second between 10 and 5? Ma. Each event has particular dominant variables (tectonism, eustacy, sediment accumulation rate) that permitted the preservation of the record and development of the sea on the continent. The 15–13 Ma transgression was tectonically and eustatically controlled, flooding older sedimentary accommodation zones on the South American plate during a global high sea level, whereas the 105? Ma event was predominantly tectonically controlled, generated by tectonic loading created in the Cordillera Oriental fold-and-thrust belt. A new 7.72±0.31 Ma ⁴⁰Ar/³⁹Ar date from the Río Parapetí in Bolivia suggests that the 15–13 Ma transgression registered in Argentina produced no continental connection to the Caribbean transgression, registered in Bolivia, because of temporal constraints.