Geomorphological constraining of tsunami (?) run-up in the Alcantarilha coastal lowland (Central Algarve, Portugal)

The Alcantarilha lowland, partly barred by a well developed barrier, including foredunes covering Pleistocene-Holocene beachrock and aeolianite, develops across the Alcantarilha infilled estuary, the beach-dune extending further SE until the Salgados lagoon. A topographic and coring survey revealed a peculiar feature at the leeward toe of the dune ridge close to the inlet area: a sandy fan with location, shape and morphology suggesting emplacement by single or multiple overwash of the barrier tip rather than tidal forcing. Its storm or tsunami origin and age are under investigation, and the only time-constrain available at present is that it should post-date ca. 6600 cal BP, the most recent in situ aeolianite (Moura et al., 2007) dated so far. METHODS, DATA SET AND RESULTS The fan boundaries are distinctive in aerial photos and satellite images: it is roughly ellipsoidal, ~200 m wide and ~300 m elongated paralleling the shoreline, rising ~ 0.9-1.2 m above the surrounding floodplain surface. Detailed topography shows that its short axis aligns with SW-NE elongated (though irregular) depressions in the dune crest, which link the beach with the fan. This could have favoured funnelling of, or erosion by, water overtopping the barrier but, in either case, the fan should correspond to extreme and abrupt event(s) of coastal flooding. 18 trenches and cores were performed in the exposed area of the fan and nearby flood plain to obtain samples and data on its sedimentology, lithostratigraphy and geometry. The fan consists of well sorted and rounded sand (Fig. 2). It thins away and wedges out landwards of the apex (located near Alc29T) where it is partly covered by dune sand. Its lower boundary is undulating and marked by textural contrast between sand (fan) and underlying mud (alluvial/lagoonal); an accumulation of marine-sourced perforated pebbles showing limited lateral continuity may pinpoint this boundary near the foredune (core Alc 25, ca. 80 m westward of profile in Fig.1); mud-balls were also observed immediately above this surface in cores and trenches. As the washover was probably emplaced in a barred lagoonal/estuarine floodplain setting, the fan’s northern outer belt is enclosed by low-energy sediments (not shown in Fig. 2).

Phase-field analysis of finite-strain plates and shells including element subdivision

With the theme of fracture of finite-strain plates and shells based on a phase-field model of crack regularization, we introduce a new staggered algorithm for elastic and elasto-plastic materials. To account for correct fracture behavior in bending, two independent phase-fields are used, corresponding to the lower and upper faces of the shell. This is shown to provide a realistic behavior in bending-dominated problems, here illustrated in classical beam and plate problems. Finite strain behavior for both elastic and elasto-plastic constitutive laws is made compatible with the phase-field model by use of a consistent updated-Lagrangian algorithm. To guarantee sufficient resolution in the definition of the crack paths, a local remeshing algorithm based on the phase- field values at the lower and upper shell faces is introduced. In this local remeshing algorithm, two stages are used: edge-based element subdivision and node repositioning. Five representative numerical examples are shown, consisting of a bi-clamped beam, two versions of a square plate, the Keesecker pressurized cylinder problem, the Hexcan problem and the Muscat-Fenech and Atkins plate. All problems were successfully solved and the proposed solution was found to be robust and efficient.

Central Algarve karst system tufa-related dynamics, Portugal.

Geomorphological mapping is a powerful instrument improving the geomorphological interpretation and understanding of the processes and forms used in landscape studies, with the ability of organizing different thematic layers in the same map. The presented map provide relevant information about the different geomorphological units of the central Algarve (i.e. the Carboniferous flysch mountains; the Barrocal, with marly and karstified subunits), where a karst system is prominent. Solution karst morphologies and large dry areas are common in the elevated areas of the Barrocal, suggesting deep circulation of groundwater. These recharge areas feed the perched aquifers of the area, where discharge is controlled by the impervious lithologies (clay-rich strata of the turbidites, marls and argilites) in the valley bottoms or other leaks in dammed aquifers. In springs related to the main aquifers tufa are actively being formed and, close coupled to spring location, different tufa depositional systems develop.

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA Pedro P. Cunha 1, António A. Martins 2, Jan-Pieter Buylaert 3,4, Andrew S. Murray 4, Luis Raposo 5, Paolo Mozzi 6, Martin Stokes 7 1 MARE - Marine and Environmental Sciences Centre, Department of Earth Sciences, University of Coimbra, Portugal: pcunha@dct.uc.pt 2 MARE - Marine and Environmental Sciences Centre, Dep. Geociências, University of Évora, Portugal; aam@uevora.pt 3 Centre for Nuclear Technologies, Technical University of Denmark, Risø Campus, Denmark; jabu@dtu.dk 4 Nordic Laboratory for Luminescence Dating, Aarhus University, Risø DTU, Denmark; anmu@dtu.dk 5 Museu Nacional de Arqueologia, Lisboa, Portugal; 3raposos@sapo.pt 6 Department of Geosciences, University of Padova, Italy; paolo.mozzi@unipd.it 7 School of Geography, Earth and Environmental Sciences, University of Plymouth, UK; m.stokes@plymouth.ac.uk The stratigraphic units that record the evolution of the Tagus River in Portugal (study area between Vila Velha de Ródão and Porto Alto villages; Fig. 1) have different sedimentary characteristics and lithic industries (Cunha et al., 2012): - a culminant sedimentary unit (the ancestral Tagus, before the drainage network entrenchment) – SLD13 (+142 to 262 m above river bed – a.r.b.; with probable age ca. 3,6 to 1,8 Ma), without artefacts; - T1 terrace (+84 to 180 m; ca. 1000? to 900 ka), without artefacts; - T2 terrace (+57 to 150 m; top deposits with a probable age ca. 600 ka), without artefacts; - T3 terrace (+43 to 113 m; ca. 460 to 360? ka), without artefacts; - T4 terrace (+26 to 55 m; ca. 335 a 155 ka), Lower Paleolithic (Acheulian) at basal and middle levels but early Middle Paleolithic at top levels; - T5 terrace (+5 to 34 m; 135 to 73 ka), Middle Paleolithic (Mousterian; Levallois technique); - T6 terrace (+3 to 14 m; 62 to 32 ka), late Middle Paleolithic (late Mousterian); - Carregueira Sands (aeolian sands) and colluvium (+3 a ca. 100 m; 32 to 12 ka), Upper Paleolithic to Epipaleolithic; - alluvial plain (+0 to 8 m; ca. 12 ka to present), Mesolithic and more recent industries. The differences in elevation (a.r.b.) of the several terrace staircases results from differential uplift due to active faults. Longitudinal correlation with the terrace levels indicates that a graded profile ca. 200 km long was achieved during terrace formation periods and a strong control by sea base level was determinant for terrace formation. The Neogene sedimentary units constituted the main source of sediments for the fluvial terraces (Fig. 2). Geomorphological mapping, coupled with lithostratigraphy, sedimentology and luminescence dating (quartz-OSL and K-feldspar post-IRIR290) were used in this study focused on the T4 terrace, which comprises a Lower Gravels (LG) unit and an Upper Sand (US) unit. The thick, coarse and dominantly massive gravels of the LG unit indicate deposition by a coarse bed-load braided river, with strong sediment supply, high gradient and fluvial competence, during conditions of rapidly rising sea level. Luminescence dating only provided minimum ages but it is probable that the LG unit corresponds to the earlier part of the MIS9 (ca. 335 to 325 ka), immediately postdating the incision promoted by the very low sea level (reaching ca. -140 m) during MIS10 (362 to 337 ka), a period of relatively cold climate conditions with weak vegetation cover on slopes and low sea level. Fig. 1. Main Portuguese reaches in which the Tagus River can be divided (Lower Tagus Basin): I – from the Spanish border to Arneiro (a general E–W trend, mainly consisting of polygonal segments); II – from Arneiro to Gavião (NE–SW); III – from Gavião to Arripiado (E–W); IV – from Arripiado to Vila Franca de Xira (NNE-SSW); V – from Vila Franca de Xira to the Atlantic shoreline. The faults considered to be the limit of the referred fluvial sectors are: F1 – Ponsul-Arneiro fault (WSW-ENE); F2 – Gavião fault (NW-SE); F3 – Ortiga fault (NW-SE); F4 – Vila Nova da Barquinha fault (W-E); F5 – Arripiado-Chamusca fault (NNE-SSW). 1 – estuary; 2 – terraces; 3 – faults; 4 – Tagus main channel. The main Iberian drainage basins are also represented (inset). The lower and middle parts of the US unit, comprising an alternation of clayish silts with paleosols and minor sands to the east (flood-plain deposits) and sand deposits to the west (channel belt), have a probable age of ca. 325 to 200 ka. This points to formation during MIS9 to MIS7, under conditions of high to medium sea levels and warm to mild conditions. The upper part of the US unit, dominated by sand facies and with OSL ages of ca. 200 to 154 ka, correlates with the early part of the MIS6. During this period, progradation resulted from climate deterioration and relative depletion of vegetation that promoted enhanced sediment production in the catchment, coupled with initiation of sea-level lowering that increased the longitudinal slope. The Vale do Forno and Vale da Atela archaeological sites (Alpiarça, central Portugal) document the earliest human occupation in the Lower Tagus River, well established in geomorphological and environmental terms, within the Middle Pleistocene. The Lower Palaeolithic sites were found on the T4 terrace (+26 m, a.r.b.). The oldest artefacts previously found in the LG unit, display crude bifacial forms that can be attributed to the Acheulian, with a probable age of ca. 335 to 325 ka. The T4 US unit has archaeological sites stratigraphically documenting successive phases of an evolved Acheulian, that probably date ca. 325 to 300 ka. Notably, these Lower Palaeolithic artisans were able to produce tools with different sophistication levels, simply by applying different strategies: more elaborated reduction sequences in case of bifaces and simple reduction sequences to obtain cleavers. Fig. 2. . Simplified geologic map of the Lower Tagus Cenozoic basin, adapted from the Carta Geológica de Portugal, 1/500000, 1992). The study area (comprising the Vale do Forno and Vale de Atela sites) is located on the more upstream sector of the Lower Tagus River reach IV, between Arripiado and Chamusca villages. 1 – alluvium (Holocene); 2 – terraces (Pleistocene); 3 – sands, silts and gravels (Paleogene to Pliocene); 4 – Sintra Massif (Cretaceous); 5 – limestones, marls, silts and sandstones (Mesozoic); 6 – quartzites (Ordovician); 7 – basement (Proterozoic to Palaeozoic); 8 – main fault. The main Portuguese reaches of the Tagus River are identified (I to V). The VF3 site (Milharós), containing a Final Acheulian industry, with fine and elaborated bifaces) found in a stratigraphic level located between the T4 terrace deposits and a colluvium associated with Late Pleistocene aeolian sands (32 to 12 ka), has an age younger than ca. 154 ka but much older than 32 ka. In the study area, the sedimentary units of the T4 terrace seem to record the river response to sea-level changes and climatically-driven fluctuations in sediment supply. REFERENCES Cunha P. P., Almeida N. A. C., Aubry T., Martins A. A., Murray A. S., Buylaert J.-P., Sohbati R., Raposo L., Rocha L., 2012, Records of human occupation from Pleistocene river terrace and aeolian sediments in the Arneiro depression (Lower Tejo River, central eastern Portugal). Geomorphology, vol. 165-166, pp. 78-90.

Evaluation of frozen ground conditions along a coastal topographic gradient at Byers Peninsula (Livingston Island, Antarctica) by geophysical and geoecological methods

Geophysical surveying and geoelectricalmethods are effective to study permafrost distribution and conditions in polar environments. Geoelectrical methods are particularly suited to study the spatial distribution of permafrost because of its high electrical resistivity in comparison with that of soil or rock above 0 °C. In the South Shetland Islands permafrost is considered to be discontinuous up to elevations of 20–40ma.s.l., changing to continuous at higher altitudes. There are no specific data about the distribution of permafrost in Byers Peninsula, in Livingston Island, which is the largest ice-free area in the South Shetland Islands. With the purpose of better understanding the occurrence of permanent frozen conditions in this area, a geophysical survey using an electrical resistivity tomography (ERT)methodologywas conducted during the January 2015 field season, combined with geomorphological and ecological studies. Three overlapping electrical resistivity tomographies of 78meach were done along the same profile which ran from the coast to the highest raised beaches. The three electrical resistivity tomographies are combined in an electrical resistivitymodel which represents the distribution of the electrical resistivity of the ground to depths of about 13malong 158m. Several patches of high electrical resistivity were found, and interpreted as patches of sporadic permafrost. The lower limits of sporadic to discontinuous permafrost in the area are confirmed by the presence of permafrost-related landforms nearby. There is a close correspondence between moss patches and permafrost patches along the geoelectrical transect.

A novel two-stage discrete crack method based on the screened Poisson equation and local mesh refinement

We propose an alternative crack propagation algo- rithm which effectively circumvents the variable transfer procedure adopted with classical mesh adaptation algo- rithms. The present alternative consists of two stages: a mesh-creation stage where a local damage model is employed with the objective of defining a crack-conforming mesh and a subsequent analysis stage with a localization limiter in the form of a modified screened Poisson equation which is exempt of crack path calculations. In the second stage, the crack naturally occurs within the refined region. A staggered scheme for standard equilibrium and screened Poisson equa- tions is used in this second stage. Element subdivision is based on edge split operations using a constitutive quantity (damage). To assess the robustness and accuracy of this algo- rithm, we use five quasi-brittle benchmarks, all successfully solved.

Damage and fracture algorithm using the screened Poisson equation and local remeshing

We propose a crack propagation algorithm which is independent of particular constitutive laws and specific element technology. It consists of a localization limiter in the form of the screened Poisson equation with local mesh refinement. This combination allows the cap- turing of strain localization with good resolution, even in the absence of a sufficiently fine initial mesh. In addition, crack paths are implicitly defined from the localized region, cir- cumventing the need for a specific direction criterion. Observed phenomena such as mul- tiple crack growth and shielding emerge naturally from the algorithm. In contrast with alternative regularization algorithms, curved cracks are correctly represented. A staggered scheme for standard equilibrium and screened equations is used. Element subdivision is based on edge split operations using a given constitutive quantity (either damage or void fraction). To assess the robustness and accuracy of this algorithm, we use both quasi-brittle benchmarks and ductile tests.

How did the AD 1755 tsunami impact on sand barriers across the southern coast of Portugal?

Tsunamis are highly energetic events that may destructively impact the coast. Resolving the degree of coastal resilience to tsunamis is extremely difficult and sometimes impossible. In part, our understanding is constrained by the limited number of contemporaneous examples and by the high dynamism of coastal systems. In fact, longterm changes of coastal systems can mask the evidence of past tsunamis, leaving us a short or incomplete sedimentary archive. Here, we present a multidisciplinary approach involving sedimentological, geomorphological and geophysical analyses and numerical modelling of the AD 1755 tsunami flood on a coastal segment located within the southern coast of Portugal. In particular, the work focuses on deciphering the impact of the tsunami waves over a coastal sand barrier enclosing two lowlands largely inundated by the tsunami flood. Erosional features documented by geophysical data were assigned to the AD 1755 eventwith support of sedimentological and age estimation results. Furthermore, these features allowed the calibration of the simulation settings to reconstruct the local conditions and establish the run-up range of the AD 1755 tsunami when it hit this coast (6– 8 m above mean sea level). Our work highlights the usefulness of erosional imprints preserved in the sediment record to interpret the impact of the extreme events on sand barriers