Ephemeral rollover points and clinothem evolution in the modern Po Delta based on repeated bathymetric surveys

Reconstructions of ancient delta systems rely typically on a two-dimensional (2D) view of prograding clinothems but may miss their three-dimensional (3D) stratigraphic complexity which can, instead, be best documented on modern delta systems by integrating high-resolution geophysical data, historical cartography, core data and geomorphological reconstructions offshore. We quantitatively compare three precisely positioned, high-resolution multi-beam bathymetry maps in the delta front and pro delta sectors (0.3 to 10 m water depth) of Po di Pila, the most active of the modern Po Delta five branches. By investigating the detailed morphology of the prograding modern Po Delta, we shed new light on the mechanisms that control the topset to foreset transition in clinothems and show the temporal and spatial complexity of a delta and its pro delta slope, under the impact of oceanographic processes. This study documents the ephemeral nature of the rollover point at the transition between sandy topset (fluvial, delta plain to mouth-bar) and muddy seaward-dipping foreset deposits advancing, in this case, in >20 m of water depth. Three multibeam surveys, acquired between 2013 and 2016, document the complexity in space and time of the topset and foreset regions and their related morphology, a diagnostic feature that could not be appreciated using solely 2D, even very high-resolution, seismic profiles. In addition, the comparison of bathymetric surveys gathered with one-year lapses shows the migration of subaqueous sand dunes on the clinothem topset, the formation of ephemeral cut-and-fill features at the rollover point (few m below mean sea level), the presence of collapse depressions derived by sagging of sediments and fluid expulsion (possibly induced by storm waves) on the foreset, and splays of sand likely reflecting gravity flows on the lower foreset. Though the modern Po Delta is anthropogenic in many respects, its subaqueous clinothem can be studied as a scale model for ancient clinothems that are less resolved geometrically and far less constrained chronologically.     

Climate control on stacked paleosols in the Pleistocene of the Po Basin (northern Italy)

Paleosols are recurrent features in alluvial successions and provide information about past sedimentary dynamics and climate change. Through sedimentological analysis on six sediment cores, the mud-dominated succession beneath the medieval ‘Two Towers’ of Bologna was investigated down to 100 m depth. A succession of weakly developed paleosols (Inceptisols) was identified. Four paleosols (P1, P2, P3 and PH) were radiocarbon-dated to 40–10 cal ka bp . Organic matter and CaCO₃ determinations indicate low groundwater levels during soil development, which spanned periods < 5 ka. The development and burial of soils, which occurred synchronously in the Bologna region and in other sectors of the Po Plain, are interpreted to reflect climatic and eustatic variations. Climatic oscillations, at the scale of the Bond cycles, controlled soil development and burial during Marine Isotope Stage (MIS) 3 (P1 and P2). Rapid sea-level oscillations probably induced soil development at the MIS 3/2 transition (P3) and favored burial of PH after 10 ka bp . Weakly developed paleosols in alluvial successions can provide clues to millennial-scale climatic and environmental variations. In particular, the paleosol-bearing succession of the Po Plain represents an unprecedent record of environmental changes across the Late Pleistocene (MIS 3 and 2) in the Mediterranean region.     

Late‐glacial to Holocene depositional architecture of the Ombrone palaeovalley system (Southern Tuscany,Italy): Sea‐level,climate and local control in valley‐fill variability

The Ombrone palaeovalley was incised during the last glacial sea‐level fall and was infilled during the subsequent Late‐glacial to Holocene transgression. A detailed sedimentological and stratigraphic study of two cores along the palaeovalley axis led to reconstruction of the post‐Last Glacial Maximum valley‐fill history. Stratigraphic correlations show remarkable similarity in the Late‐glacial to early‐Holocene succession, but discrepancy in the Holocene portion of the valley fill. Above the palaeovalley floor, about 60 m below sea‐level, Late‐glacial sedimentation is recorded by an unusually thick alluvial succession dated back to ca 18 cal kyr bp . The Holocene onset was followed by the retrogradational shift from alluvial to coastal facies. In seaward core OM1, the transition from inner to outer estuarine environments marks the maximum deepening of the system. By comparison, in landward core OM2, the emplacement of estuarine conditions was interrupted by renewed continental sedimentation. Swamp to lacustrine facies, stratigraphically equivalent to the fully estuarine facies of core OM1, represent the proximal expression of the maximum flooding zone. This succession reflects location in a confined segment of the valley, just landward of the confluence with a tributary valley. It is likely that sudden sediment input from the tributary produced a topographic threshold, damming the main valley course and isolating its landward segment from the sea. The seaward portion of the Ombrone palaeovalley presents the typical estuarine backfilling succession of allogenically controlled incised valleys. In contrast, in the landward portion of the system, local dynamics completely overwhelmed the sea‐level signal, following marine ingression. This study highlights the complexity of palaeovalley systems, where local morphologies, changes in catchment areas, drainage systems and tributary valleys may produce facies patterns significantly different from the general stratigraphic organization depicted by traditional sequence‐stratigraphic models.     

Deformation patterns of upper Quaternary strata and their relation to active tectonics,Po Basin,Italy

Despite increased application of subsurface datasets below the limits of seismic resolution, reconstructing near‐surface deformation of shallow key stratigraphic markers beneath modern alluvial and coastal plains through sediment core analysis has received little attention. Highly resolved stratigraphy of Upper Pleistocene to Holocene (Marine Isotope Stage 5e to Marine Isotope Stage 1) alluvial, deltaic and coastal depositional systems across the southern Po Plain, down to 150 m depth, provides an unambiguous documentation on the deformation of previously flat‐lying strata that goes back in time beyond the limits of morphological, historical and palaeoseismic records. Five prominent key horizons, accurately selected on the basis of their sedimentological characteristics and typified for their fossil content, were used as highly effective stratigraphic markers (M1 to M5) that can be tracked for tens of kilometres across the basin. A facies‐controlled approach tied to a robust chronology (102 radiocarbon dates) reveals considerable deformation of laterally extensive nearshore (M1), continental (M2 and M3) and lagoon (M4 and M5) marker beds originally deposited in a horizontal position (M1, M4 and M5). The areas where antiformal geometries are best observed are remarkably coincident with the axes of buried ramp anticlines, across which new seismic images reveal substantially warped stratal geometries of Lower Pleistocene strata. The striking spatial coincidence of fold crests with the epicentres of historic and instrumental seismicity suggests that deformation of marker beds M1 to M5 might reflect, in part at least, syntectonically generated relief and, thus, active tectonism. Precise identification and lateral tracing of chronologically constrained stratigraphic markers in the ¹⁴C time window through combined sedimentological and palaeoecological data may delineate late Quaternary subsurface stratigraphic architecture at an unprecedented level of detail, outlining cryptic stratal geometries at the sub‐seismic scale. This approach is highly reproducible in tectonically active Quaternary depositional systems and can help to assess patterns of active deformation in the subsurface of modern alluvial and coastal plains worldwide.     

Elasto‐plastic analysis of block structures through a homogenization method

The paper describes the development and numerical implementation of a constitutive relationship for modeling the elasto‐plastic behavior of block structures with periodic texture, regarded at a macroscopic scale as homogenized anisotropic media. The macroscopic model is shown to retain memory of the mechanical characteristics of the joints and of the shape of the blocks. The overall mechanical properties display anisotropy and singularities in the yield surface, arising from the discrete nature of the block structure and the geometrical arrangement of the units. The model is formulated in the framework of multi‐surface plasticity. It is implemented in an finite element (FE) code by means of two different algorithms: an implicit return mapping scheme and a minimization algorithm directly derived from the Haar–Karman principle. The model is validated against analytical and experimental results: the comparison between the homogenized continuum and the original block assembly shows a good agreement in terms of ultimate inelastic behavior, when the size of the block is small as compared with that of the whole assembly. Copyright © 2009 John Wiley & Sons, Ltd.     

Biosedimentary record of postglacial coastal dynamics: high‐resolution sequence stratigraphy from the northern Tuscan coast (Italy)

Integrating analysis of the benthic palaeoecological record with multivariate ordination techniques represents a powerful synergy able to provide an improved characterization of coastal depositional facies in a sequence stratigraphical perspective. Through quantitative analysis of benthic foraminifer, ostracod and mollusc associations from the postglacial succession of Core M3 (Arno coastal plain, Tuscany, Italy), and application of detrended correspondence analysis (DCA) to the mollusc sub‐data set, we offer a refined picture of stratigraphical variations in faunal content from a paralic depositional setting, and reconstruct the palaeoenvironmental gradients that account for such variations. Despite distinct ecological behaviours, and taphonomic and sedimentological constraints, a strong ecological control on meio‐ and macrofaunal biofacies and taxa turnover is documented across the study succession. Amongst all possible mechanisms that may play a role in ‘shaping’ fossil distribution, the ecological signal driven by salinity represents the most prominent factor controlling the composition of fossil associations in the cored succession. Molluscs can even provide outstanding quantitative estimates of palaeosalinity along the sampled core. When plotted stratigraphically, the three fossil sub‐data sets show consistent patterns of vertical evolution that enable prompt identification of the key surfaces for sequence stratigraphical interpretation in otherwise lithologically indistinguishable deposits. The concomitant maximum richness of species with strong marine affinity, paralleled by the highest DCA salinity estimates, allows recognition of the maximum flooding zone, dated to ～7.7 cal. ka BP, within a homogeneous succession of outer lagoon clays. These clays are sandwiched between early transgressive, swamp to inner lagoon deposits and overlying prograding coastal?alluvial plain facies.     

Glaucony from the Cretaceous of the Sierra de Guadarrama (Central Spain) and its application in a sequence-stratigraphic context

Detailed sedimentological and geochemical analyses of Upper Cretaceous glaucony-bearing deposits from the middle portion of the Castro de Fuentidueña Formation, in Central Spain, enable identification of a multiphase history of glaucony accumulation. Despite its relatively high maturity (K₂O > 7%), glaucony from the transgressive subtidal sandstones has anomalously low concentration in the host rock (generally <10%), suggesting remobilization from a different source. The remarkable thickness of the glaucony-bearing horizons, concurrently with concentration of glaucony in laminae, small size, and high degree of roundness and sorting provide further evidence for an allochthonous origin of the green grains. In contrast, authigenic glaucony from the overlying offshore clays exhibits higher abundance (up to 45%) and lower maturity (K₂O < 7%) and is interpreted to have formed in situ. Allochthonous glaucony originally developed in more distal areas during a prolonged period of sediment starvation. The green grains were then swept away from their place of origin and concentrated into tidal bars and channels within the upper transgressive systems tract of the third-order depositional sequence. Maximum concentration of autochthonous glaucony is recorded at the transition from tidal-influenced to offshore deposits: this glaucony, which is relatively less evolved due to lower time available for maturation, is interpreted to reflect the turnaround from transgressive to highstand conditions, marking the condensed section of the depositional sequence. As postulated by previous sequence-stratigraphic models, this study documents that allochthonous glaucony can be widespread throughout the transgressive systems tract (TST), while the condensed section (CS) typically hosts autochthonous glaucony. Contrary to the existing literature, however, this study shows that glaucony from the TST may be even more mature than glaucony from the CS, if enough time for maturation during transgression is allowed before the ultimate concentration of the green grains. Identification of spatial and temporal characteristics of glaucony, thus, is critical to a reliable sequence-stratigraphic interpretation of the glaucony-bearing deposits.     

Numerical modelling of the transverse dynamic behaviour of circular tunnels in clayey soils

In this paper, different approaches aimed at investigating the dynamic behaviour of circular tunnels in the transverse direction are presented. The analysed cases refer to a shallow tunnel built in two different clayey deposits. The adopted approaches include 1D numerical analyses performed modelling the soil as a single-phase visco-elastic non-linear medium, the results of which are then used to evaluate the input data for selected analytical solutions proposed in the literature (uncoupled approach), and 2D fully coupled FE simulations adopting visco-elastic and visco-elasto-plastic effective stress models for the soil (coupled approach). The results are proposed in terms of seismic-induced loads in the transverse direction of the tunnel lining. The different constitutive hypotheses adopted in the coupled numerical approach prove to play a significant role on the results. In particular, the plasticity-based analyses indicate that a seismic event can produce a substantial modification of loads acting in the lining, leading to permanent increments of both hoop force and bending moment.