Tectonic significance of the present relief of the Betic Cordillera

Tortonian calcarenites of the Betic Cordillera were deposited in coastal or very shallow marine environments and represent an ideal marker for estimating vertical movements from the late Miocene to the Present. A map showing the heights at which these Tortonian marine rocks are situated has a clear correlation with the present relief, indicating that today's relief has been formed since the Tortonian. There is also a good correlation between present relief and the Bouguer anomaly distribution in the Betic Cordillera, as well as with crustal thickness. Likewise, the present relief is directly related to the geodynamic setting of a horizontal N–S to NNW–SSE compression and an almost perpendicular extension, along with isostatic readjustment, existing in the Betic Cordillera from the Tortonian. As a result of these regional stresses, faults and folds have produced notable vertical movements. The highest rates of uplift of the Betic Cordillera coincide with large antiforms, in particular those of the Sierra Nevada and the Sierra Filabres. Several subsiding sectors also exist (for example, the Granada Basin or the Guadalquivir Basin). The foreland Guadalquivir Basin has a complex history because the uplift in its eastern sector and subsidence in the western sector coexisted during the late Tortonian. Today the whole Betic Cordillera is characterized by differential regional uplift, even in the aforementioned subsiding sectors.     

Effects of spatially structured vegetation patterns on hillslope erosion in a semiarid Mediterranean environment: a simulation study

A general trend of decreasing soil loss rates with increasing vegetation cover fraction is widely accepted. Field observations and experimental work, however, show that the form of the cover‐erosion function can vary considerably, in particular for low cover conditions that prevail on arid and semiarid hillslopes. In this paper the structured spatial distribution of the vegetation cover and associated soil attributes is proposed as one of the possible causes of variation in cover–erosion relationships, in particular in dryland environments where patchy vegetation covers are common. A simulation approach was used to test the hypothesis that hillslope discharge and soil loss could be affected by variation in the spatial correlation structure of coupled vegetation cover and soil patterns alone. The Limburg Soil Erosion Model (LISEM) was parameterized and verified for a small catchment with discontinuous vegetation cover at Rambla Honda, SE Spain. Using the same parameter sets LISEM was subsequently used to simulate water and sediment fluxes on 1 ha hypothetical hillslopes with simulated spatial distributions of vegetation and soil parameters. Storms of constant rainfall intensity in the range of 30–70 mm h^?1 and 10–30 min duration were applied. To quantify the effect of the spatial correlation structure of the vegetation and soil patterns, predicted discharge and soil loss rates from hillslopes with spatially structured distributions of vegetation and soil parameters were compared with those from hillslopes with spatially uniform distributions. The results showed that the spatial organization of bare and vegetated surfaces alone can have a substantial impact on predicted storm discharge and erosion. In general, water and sediment yields from hillslopes with spatially structured distributions of vegetation and soil parameters were greater than from identical hillslopes with spatially uniform distributions. Within a storm the effect of spatially structured vegetation and soil patterns was observed to be highly dynamic, and to depend on rainfall intensity and slope gradient. Copyright © 2005 John Wiley & Sons, Ltd.     

The role of vegetation patterns in structuring runoff and sediment fluxes in drylands

The dynamics of vegetation‐driven spatial heterogeneity (VDSH) and its function in structuring runoff and sediment fluxes have received increased attention from both geomorphological and ecological perspectives, particularly in arid regions with sparse vegetation cover. This paper reviews the recent findings in this area obtained from field evidence and numerical simulation experiments, and outlines their implications for soil erosion assessment. VDSH is often observed at two scales, individual plant clumps and stands of clumps. At the patch scale, the local outcomes of vegetated patches on soil erodibility and hydraulic soil properties are well established. They involve greater water storage capacity as well as increased organic carbon and nutrient inputs. These effects operate together with an enhanced capacity for the interception of water and windborne resources, and an increased biological activity that accelerates breakdown of plant litter and nutrient turnover rates. This suite of relationships, which often involve positive feedback mechanisms, creates vegetated patches that are increasingly different from nearby bare ground areas. By this way a mosaic builds up with bare ground and vegetated patches coupled together, respectively, as sources and sinks of water, sediments and nutrients. At the stand scale within‐storm temporal variability of rainfall intensity controls reinfiltration of overland flow and its decay with slope length. At moderate rainfall intensity, this factor interacts with the spatial structure of VDSH and the mechanism of overland flow generation. Reinfiltration is greater in small‐grained VDSH and topsoil saturation excess overland flow. Available information shows that VDSH structures of sources and sinks of water and sediments evolve dynamically with hillslope fluxes and tune their spatial configurations to them. Rainfall simulation experiments in large plots show that coarsening VDSH leads to significantly greater erosion rates even under heavy rainfall intensity because of the flow concentration and its velocity increase. Copyright © 2005 John Wiley & Sons, Ltd.     

A comparative study of concentrated plasticity models in dynamic analysis of building structures

Concentrated plasticity (CP) models are frequently used in static and dynamic building analysis and have been implemented in available commercial software. This investigation deals with three different CP‐models, a simplified macroelement model (SEM) for a complete building story, a frame element with elasto‐plastic interaction hinges (PH), and a frame element with fiber hinges (FB). The objectives of this work are to evaluate the quality of the earthquake responses predicted by these models and to identify important aspects of their implementation and limitations for their use in dynamic analysis. The three elements are tested in a single‐story asymmetric plan building and in a three‐story steel building. Results show that base shear and global response values are usually computed with better accuracy than interstory deformations and local responses. Besides, the main limitation of elasto‐plastic CP models is to control the displacement offsets that result from perfect elasto‐plastic behavior. On the other hand, calibration of the SEM‐model shows that global responses in steel structures may be computed within 20% error in the mean at a computational cost two orders of magnitude smaller than that of the other CP elements considered. However, the three element models considered lead to increasing levels of accuracy in the dynamic response and their use depends on the refinement of the analysis performed. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

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.     

Kamafugitic diatremes: their textures and field relationships with examples from the Goiás alkaline province, Brazil

Kamafugitic rocks intruded the Precambrian basement and Phanerozoic sediments at the northeast border of the Paraná basin as part of the Late Cretaceous Goiás alkaline province (GAP). Plutonic complexes dominate the north of the province, whereas lavas and pyroclastic rocks prevail in the south. The central GAP is characterized by kamafugitic diatremes, which may crop out continuously for up to 850 m and consist of a central breccia body, surrounded and overlain by lava flows and crosscut by dykes. The breccias contain some special spheroidal juvenile fragments—namely, accretionary and armored lapilli, frozen droplets, spinning droplets, and wrapped fragments—whose textural and mineralogical aspects are described in detail. Irregularly shaped tuff pockets that occur within the breccias contain textures and structures similar to those of subaerial surge deposits and formed in confined, high gas to solid+liquid ratio domains in the conduit. Diatreme emplacement affected the country rock through thermal metamorphism, development of columnar jointing, and formation of peperite-like mixtures. There is no evidence of phreatomagmatic activity in the diatremes, and CO₂, rather than H₂O, seems to have been the major volatile component of the kamafugitic magmas. This finding implies that features such as accretionary lapilli and peperites are not exclusively associated with H₂O-dominated processes.     

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.