Geomorphological and sedimentological features in Quaternary fluvial systems affected by solution-induced subsidence (Ebro Basin, NE-Spain)

The Quaternary evolution and the morpho-sedimentary features of some of the most important rivers in Spain (Ebro and Tagus rivers among others) have been controlled by subsidence due to alluvial karstification of the evaporitic bedrock. The subsidence mechanism may range from catastrophic collapse to slow sagging of the alluvium by passive bending. In the Ebro Basin, the mechanisms and processes involved in karstic subsidence were studied through the analysis of present-day closed depressions as well as through old subsidence depressions (palaeocollapses and solution-induced basins) and associated deformations recorded in the Quaternary alluvial sediments. The Gállego–Ebro river system is presented as a case study of channel adjustments and geomorphic and sedimentary evolution of fluvial systems in dissolution-induced subsidence areas. In this fluvial system, evaporite dissolution during particular Quaternary time intervals (namely early and middle Pleistocene) have lead to the development of a solution-induced basin, approximately 30 km-long by 8 km-wide, filled by Quaternary deposits with a total thickness in excess of 190 m. The main river response to balance the subsidence in the alluvial plain was aggradation in the central reach of the subsiding area, and degradation both in the upstream reach and in the valley sides where alluvial fans and covered pediments may prograde over the fluvial sediments. The main sinking areas are recognized in the sedimentary record by anomalous thickenings in the alluvial deposits and fine-grained sediments deposited in backswamp and ponded areas.     

Palaeoenvironmental interpretation of Late Pleistocene-Holocene morphosedimentary record in the Valsalada saline wetlands (Central Ebro Basin, NE Spain)

This work presents a palaeoenvironmental interpretation of the Upper Pleistocene-Holocene sedimentary sequence recorded in the Valsalada saline wetland system (Monegros, Central Ebro Basin). This morphosedimentary system developed on karstified, gypsiferous bedrock and was mainly fed by local saline groundwater. Based on geomorphological, sedimentological, palynological and radiocarbon data, three depositional units have been differentiated: 1) a lacustrine unit, which accumulated under cold/cool and humid environmental conditions at 41-40 kyr cal BP, 2) a fluvial unit that was deposited under arid climatic conditions with cold episodic periods from 14 to 3.5 kyr cal BP, and 3) an alluvial slope unit associated with arid climate with human influence during the Late Roman Period (1.5 kyr cal BP). The morphosedimentary evolution of the Valsalada system is strongly influenced by the geomorphic setting (small catchment with karstic depressions and low gradient slopes on gypsum), the hydrological availability (runoff and groundwater discharges) related to climate variability during the Late Pleistocene-Holocene, and the human activities developed in the area during historical times. The response of the Valsalada system displays a low sensitivity to short-term climatic changes but responds to major long-term climatic conditions, improving the regional paleoenvironmental picture of the Late Quaternary in NE Spain.     

Upper Oligocene-Lower Miocene shallowing-upwardlacustrine sequences controlled by periodic and non-periodic processes(Ebro Basin, northeastern Spain)

The lacustrine facies from two sections (Candasnos and Fraga) ofthe Oligocene-Miocene Torrente de Cinca lithostratigraphic Unit in thecentral part of the Ebro Basin (Spain) have been analysed to determine theinfluence of orbital parameters in lacustrine sedimentation. The unit ispredominantly composed of limestones and marls, and represents a shallowlacustrine freshwater system. The sedimentological features of the faciesstudied demonstrate that the lower part of the Candasnos section representsoffshore lacustrine subenvironments whereas the upper part, and the whole ofthe Fraga section, characterise marginal lacustrine areas. Series of stratalthickness variations of limestone, marl, and limestone/marl couplets fromboth sections have been analysed using spectral analysis. This shows thatinformation corresponding to periodic cycles only appears in the offshorefacies, that is to say, in the lower part of Candasnos section, and disappearsin the marginal facies where non-periodic cycles exist. Furthermore, thespectral analysis of the offshore facies highlights the existence of a peak inthe power spectrum with a period of around 7 (6.8 to 7.8) that can berecognised in the field as shallowing-upward lacustrine sequences.Magnetostratigraphic data from the Candasnos section allow us to establish atime span of 2,808 years for the limestone/marl couplet from the lower partof this section, and between 19,000 and 22,000 years for the periodic cycleidentified, thus representing the climatic precession cycle. Shallowingsequences from marginal areas do not correspond with any periodiccycle.     

Internal structure of the aeolian sand dunes of El Fangar spit, Ebro Delta (Tarragona, Spain)

This paper presents an analysis of the dune field dynamics of El Fangar Spit in the Ebro Delta (Spain), associating it with the internal structure of dunes carried out with ground-penetrating radar and supported by data from topographic DGPS. These analyses are of great importance to ascertain the state of the internal structure of dunes as an important element in their stability and, therefore in their evolution. The internal structure shows accretion and progradation sequences of dunes over beach deposits, which depend on dune morphology (height, crest orientation) and location, as well as the processes acting on them.     

Basin analysis palaeo-landscape modelling: Testing the critical controls using experimental design constrained by a real 3D geological model,Gippsland Basin,Australia

The main controlling variables for palaeo-landscape evolution are investigated to assess their relative importance using the Gippsland Basin geological history. Palaeo-landscape reconstruction is a complicated process controlled and affected by multiple variables, including tectonic, palaeo-environment, sea-level change, rainfall, sediment erosion, transportation, deposition, etc. The Basin and Landscape Dynamics software (Badlands) software was used with an efficient experimental design (ED) to guide the selected scenarios, process the results, and generate the multi-variate equations that define and identify the important controlling variables. The ED was used to test and identify the main uncertainties and their possible ranges, based on actual field data, while at the same time ensuring that the full multi-dimensional space for those variables was covered to enable the computation of multivariate equations from the minimum number of scenario runs. A full suite of 3D forward palaeo-landscape models of the Gippsland Basin was built to reconstruct the basin history from its formation to the present (Early Cretaceous to Holocene, 137-0 Ma). The models are compared to the corresponding full 3D realistic structural and stratigraphic model of the basin that has been built in Petrel (Schlumberger software). This constrains the sedimentary, stratigraphic, burial and thermal histories to the relative subsidence rates and basin-fill for each geological sequence by using the model isopachs input to the Badlands modelling. The ED required only 22 scenarios to fit 12 identified variables and test for possible interactions with each other. The most significant variables are those that control sediment supply including non-marine erodibility, rainfall, (Rainfall × Area) exponent m, Slope and critical slope while maximum % Marine Deposition and marine dispersal are also required to fill the marine accommodation space. Sea Level and subsidence only become significant when rapid enough to outpace sediment supply. The controlling factors change over time with basin development from rift to post-rift phases and interactions are highly significant.     

Three-dimensional numerical modelling of thermal transfer through the Vittel aquifer (Vosges, France), including approximated topographical effects

Summary. A geothermal and hydrochemical anomaly was observed in the Lower Triassic sandstone aquifer in the Vittel region (Vosges, France). This anomaly was attributed to a hydrothermal spring under sedimentary cover. In order to localize the hot spring and to quantify more precisely its flow rate and temperature, it is necessary to consider the 3-D thermal problem by taking into account the complex geometry of the domain and the flow rate in the aquifer. A 3-D numerical model of thermal conductive and forced convective transfer, developed for hydrological problems including approximate geometrical and topographical effects, is used which (i) can be directly applied to geologic strata (aquifers) with varying thickness and top and base slope; and (ii) allows calculation of heat flux anomalies associated with fluid flow in such geologically realistic aquifers.
The heat transfer equation is formulated in an orthogonal curvilinear coordinate system. As most geometries dealing with geothermal phenomena in sedimentary basins are nearly horizontal, this formulation can be simplified, leading readily to numerical solution with a finite difference method. The application of the 3-D model to the Vittel aquifer gives temperature results in agreement with measurements. These results provide evidence for the importance of associated forced convection and topographical effects for ground temperature distribution, and show clearly that heat flow in many basins is interpretable only if careful hydrological and thermal studies are made.     

Erosion surfaces and Neogene landscape evolution in the NE Duero Basin (north-central Spain)

We present a chronological model of erosion surface development in the Iberian and Cantabrian Ranges of north-central Spain. We map four erosion surfaces and interpret these to be related to Duero basin continental sediments and tectonic activity from Upper Oligocene to Plio-Pleistocene. The oldest erosion surface, SE1, formed across Upper Oligocene–Lower Miocene synorogenic deposits; while surface SE2 was contemporaneous with the Middle Miocene alluvial systems, ending with an uplift stage in the Astaracian. The two most recent erosion surfaces, SE3 and SE4, developed during extensional tectonic episodes and are associated with the deposition of Upper Páramo sedimentary units at the end of the Miocene (Upper Turolian) and alluvial fan deposits, known as rañas (Plio-Pleistocene). With the exception of SE1, which seems to be associated with a relatively wet climate, the surfaces formed during periods of marked aridity and generally warm temperatures. Through geostatistical reconstruction of the best preserved surface (SE2), applying ordinary kriging method to the topography (DEM) of the erosion surface and its correlating sedimentary plains, we identify the deformation processes which occurred on this surface after its formation.     

Slab pull effects from a flexural analysis of the Tonga and Kermadec trenches (Pacific Plate)

Thin-plate flexure models have been frequently used to explain the mechanical behaviour of the lithosphere at oceanic trenches, but little attention has been paid to using them as a way to check the relative importance of different plate-driving mechanisms. A 2-D numerical algorithm accounting for the flexural deflection of the lithosphere controlled by multilayered elastic–plastic rheology (brittle–elastic–ductile) has been applied to the seaward side of the Tonga and Kermadec trenches. This approach gives a better fit to the bathymetry on both trenches than assuming classical homogeneous plate models, and allows the interplate coupling forces and the lithospheric strength profile to be constrained. Our results show that, in order to fit the observed deflection of the lithosphere, a regional tensile horizontal force must act in both regions. This tensile force and its flexural effects are discussed in terms of slab pull as a main plate-driving mechanism. The predicted stress and yielding distributions partially match the outer-rise earthquake hypocentres within the subducting plate, and thus do not invalidate the model.     

Three-dimensional geometrical ray theory and modelling of transmitted seismic energy of data from the Nevada Test Site

This paper presents a geometrically based algorithm for computing synthetic seismograms for energy transmitted through a 3-D velocity distribution. 3-D ray tracing is performed to compute the traveltimes and geometrical spreading (amplitude). The formulations of both kinematic and dynamic ray-tracing systems are presented. The two-point ray-tracing problem is solved by systematically updating the initial conditions and adjusting the ray direction until the ray intersects the specified endpoint. The amount of adjustment required depends on the derivatives of the position with respect to the given starting angles between consecutive rays. The algorithm uses derivatives to define the steepest-descent direction and to update the initial directions. The convergence rate depends on the complexity of the model.
Test seismograms compare favourably with those from a 2-D asymptotic ray theory algorithm and a 3-D Gaussian-beam algorithm. The algorithm is flexible in modelling arbitrary source and recorder geometries for various smoothly varying 3-D velocity distributions. The algorithm is further tested by simulating surface-to-tunnel vibroseis field data. Shear waves as well as compressional waves may be approximately included. Application of the algorithm to a data set from the Rainier Mesa of the Nevada Test Site produced a good fit to the transmitted (first arrival) traveltimes and amplitudes, with approximately 15 per cent variation in the local 3-D velocity.     

The response of drainage basins to the late Quaternary tectonics in the Sicilian side of the Messina Strait (NE Sicily)

Despite more than 40 yr of research attributing temporal changes in streambank erosion rates to subaerial processes, little quantitative information is available on the relationships between streambank erodibility (k_d) and critical shear stress (τ_c) and the environmental conditions and processes that enhance streambank erosion potential. The study goal was to evaluate temporal changes in k_d and τ_c from soil desiccation and freeze–thaw cycling. Soil erodibility and τ_c were measured monthly in situ using a multiangle, submerged jet test device. Soil moisture, temperature, and bulk density as well as precipitation, air temperature, and stream stage were measured continuously to determine changes in soil moisture content and state. Pairwise Mann–Whitney tests indicted k_d was 2.9 and 2.1 times higher (p < 0.0065) during the winter (December–March) than in the spring/fall (April–May, October–November) and the summer (June–September), respectively. Regression analysis showed 80% of the variability in k_d was explained by freeze–thaw cycling alone. Study results also indicated soil bulk density was highly influenced by winter weather conditions (r² = 0.86): bulk density was inversely related to both soil water content and freeze–thaw cycling. Results showed that significant changes in the resistance of streambank soils to fluvial erosion can be attributed to subaerial processes. Water resource professionals should consider the implications of increased soil erodibility during the winter in the development of channel erosion models and stream restoration designs.     

Frequency response modelling of seismic waves using finite difference time domain with phase sensitive detection (TD–PSD)

This paper describes an efficient approach for computing the frequency response of seismic waves propagating in 2- and 3-D earth models within which the magnitude and phase are required at many locations. The approach consists of running an explicit finite difference time domain (TD) code with a time harmonic source out to steady-state. The magnitudes and phases at locations in the model are computed using phase sensitive detection (PSD). PSD does not require storage of time-series (unlike a fast Fourier transform), reducing its memory requirements. Additionally, the response from multiple sources can be obtained from a single finite difference run by encoding each source with a different frequency. For 2-D models with many sources, this time domain phase sensitive detection (TD–PSD) approach has a higher arithmetic complexity than direct solution of the finite difference frequency domain (FD) equations using nested dissection re-ordering (FD–ND). The storage requirements for 2-D finite difference TD–PSD are lower than FD–ND. For 3-D finite difference models, TD–PSD has significantly lower arithmetic complexity and storage requirements than FD–ND, and therefore, may prove useful for computing the frequency response of large 3-D earth models.     

Evaluating and comparing methods of sinkhole susceptibility mapping in the Ebro Valley evaporite karst (NE Spain)

Multiple sinkhole susceptibility models have been generated in three study areas of the Ebro Valley evaporite karst (NE Spain) applying different methods (nearest neighbour distance, sinkhole density, heuristic scoring system and probabilistic analysis) for each sinkhole type separately (cover collapse sinkholes, cover and bedrock collapse sinkholes and cover and bedrock sagging sinkholes). The quantitative and independent evaluation of the predictive capability of the models reveals that: (1) The most reliable susceptibility models are those derived from the nearest neighbour distance and sinkhole density. These models can be generated in a simple and rapid way from detailed geomorphological maps. (2) The reliability of the nearest neighbour distance and density models is conditioned by the degree of clustering of the sinkholes. Consequently, the karst areas in which sinkholes show a higher clustering are a priori more favourable for predicting new occurrences. (3) The predictive capability of the best models obtained in this research is significantly higher (12.5–82.5%) than that of the heuristic sinkhole susceptibility model incorporated into the General Urban Plan for the municipality of Zaragoza. Although the probabilistic approach provides lower quality results than the methods based on sinkhole proximity and density, it helps to identify the most significant factors and select the most effective mitigation strategies and may be applied to model susceptibility in different future scenarios.     

The influence of a stratified rheology on the flexural response of the lithosphere to (un)loading by extensional faulting

We present a two-layered finite difference model for the flexural response of the lithosphere to extensional faulting. The model allows for three modes of flexure: (1) fully coupled, with the upper crust and mantle welded together by the lower crust; (2) fully decoupled, with the upper crust and mantle behaving as independent layers; and (3) partly decoupled, signifying that the response of the upper crust to small-wavelength loads is superimposed on the response of the entire lithosphere to long-wavelength loads. Which of these modes of flexure is to be expected depends on the rheology and especially the thermal state of the lithosphere. Coupled behaviour is related to a cold and strong lithosphere. The Baikal Rift Zone provides a typical example for this mode of flexure. A fully decoupled lithosphere is an exceptional case, related to anomalous high temperatures in the lower crust, and is observed in the Basin and Range province. The most common case is a partly decoupled lithosphere, with the degree of decoupling depending on the thickness and viscosity of the lower crust. This is inferred, for example, for the Bay of Biscay margin.     

Tectonic and climatic controls on the sequential arrangement of an alluvial fan/fan‐delta complex (Montserrat,Eocene, Ebro Basin,NE Spain)

A magnetostratigraphy‐based chronological framework has been constructed in the Eocene sediments of the Montserrat alluvial fan/fan‐delta complex (southeast Ebro Basin), in order to unravel forcing controls on their sequential arrangement and to revise the tectonosedimentary history of the region. The palaeomagnetic study is based on 403 sites distributed along an 1880‐m‐thick composite section, and provides improved temporal constraints based on an independent correlation to the geomagnetic polarity time scale. The new chronological framework together with sequence stratigraphy and geohistory analysis allow us to investigate the interplay between factors controlling the sequential arrangement of the Montserrat complex at the different temporal scales and to test for orbitally driven climate forcing. The results suggest that the internal stacking pattern in transgressive and regressive sequences sets within the more than 1000‐m‐thick Milany Composite Megasequence can be explained as the result of subsidence‐driven accommodation changes under a general increase of sediment supply. Composite sequences (tens to hundreds of metres thick) likely reflect orbitally forced cyclicity related to the 400‐kyr eccentricity cycle, possibly controlled by climatically induced sea‐level fluctuations. This study also provides new insights on the deformational history of the area, and shows a correlation between (tectonic) subsidence and forelimb rotation measured on basin‐margin deformed strata. Integration of subsidence curves from different sectors of the eastern Ebro Basin allows us to estimate the variable contribution of tectonic loads from the two active basin margins: the Catalan Coastal Ranges and the Pyrenees. The results support the presence of a double flexure from Late Lutetian to Late Bartonian, associated with the two tectonically active margins. From Late Bartonian to Early Priabonian the homogenization of subsidence values is interpreted as the result of the coupling of the two sources of tectonic load.