A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill
of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and
sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic,
hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It
also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust,
fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence
of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality
constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They
are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a
proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and
efficiency on two real basin test cases. 相似文献
FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation
tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov
scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go
phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into
FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley
and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and
published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian
Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses
revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data.
Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow
behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered
satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material.
The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation
strongly improves the simulation results and better represents steep fronts of final deposits. 相似文献
Analysis of strain in Jurassic argillites forming part of the folded and thrusted sedimentary succession of the Lagonegro basin (southern Italian Apennines) has been carried out using ellipsoid-shaped reduction spots as strain markers. Most of the determined finite strain ellipsoids are of oblate type and show a peculiar distribution of the maximum extension direction (X), with maxima either subparallel or subperpendicular to the local fold axes. Using the strain matrix method, two different deformation histories have been considered to assist the interpretation of the observed finite strain pattern. A first deformation history involved vertical compaction followed by horizontal shortening (occurring by a combination of true tectonic strain and volume loss), whereby all strain is coaxial and there is no change in the intermediate axis of the strain ellipsoid. By this type of deformation sequence, which produces a deformation path where total strain moves from the oblate to the prolate strain field and back to the oblate field, prolate strain ellipsoids can be generated and may be recorded where tectonic deformation has not been large enough to reverse pretectonic compaction. This type of deformation history may be of local importance within the study area (i.e. it may characterize some fold hinge regions) and, more generally, is probably of limited occurrence in deformed pelitic rocks. A second deformation sequence considered the superposition of pre-tectonic compaction and tectonic strain consisting of initial layer-parallel shortening followed by layer-parallel shear (related to flexural folding). Also in this instance, volume change during tectonic deformation and tectonic plane strain have been assumed. For geologically reasonable amounts of volume loss due to compaction and of initial layer-parallel shortening, this type of deformation history is capable of producing a deformation path entirely lying within the oblate strain field, but still characterized by a changeover, during deformation, of the maximum extension axis (X) from a position parallel to the fold axis to one perpendicular to it. This type of deformation sequence may explain the main strain features observed in the study area, where most of the measured finite strain ellipsoids, determined from the limb regions of flexural folds, display an oblate shape, irrespective of the orientation of their maximum extension direction (X) with respect to the local structural trends. More generally, this type of deformation history provides a mechanism to account for the predominance of oblate strains in deformed pelitic rocks. 相似文献
We have identified an extinct E–W spreading center in the northern Natal valley on the basis of magnetic anomalies which was active from chron M11 (133 Ma) to 125.3 Ma, just before chron M2 (124 Ma) in the Early Cretaceous. Seafloor spreading in the northern Natal valley accounts for approximately 170 km of north–south motion between the Mozambique Ridge and Africa. This extension resolves the predicted overlap of the continental (central and southern) Mozambique Ridge and Antarctica in the chron M2 to M11 reconstructions from Mesozoic finite rotation parameters for Africa and Antarctica. In addition, the magnetic data reveal that the Mozambique Ridge was an independent microplate from at least 133 to 125 Ma. The northern Natal valley extinct spreading center connects to the spreading center separating the Mozambique Basin and the Riiser-Larsen Sea to the east. It follows that the northern Mozambique Ridge was either formed after the emplacement of the surrounding oceanic crust or it is the product of a very robust spreading center. To the west the extinct spreading center connects to the spreading center separating the southern Natal valley and Georgia Basin via a transform fault. Prior to chron M11, there is still a problem with the overlap of Mozambique Ridge if it is assumed to be fixed with respect to either the African or Antarctic plates. Some of the overlap can be accounted for by Jurassic deformation of the Mozambique Ridge, Mozambique Basin, and Dronning Maud land. It appears though that the Mozambique Ridge was an independent microplate from the breakup of Gondwana, 160 Ma, until it became part of the African plate, 125 Ma. 相似文献
The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.
The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.
The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.
The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.
The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins. 相似文献