Impact of historical tsunamis on a sandy coastal barrier: an example from the northern Gargano coast, southern Italy

International audienceThe Lesina coastal barrier is characterized by the presence of three wide washover fans. They were formed by three distinct tsunamis which struck the northern coast of the Gargano Promontory (Apulia, Italy) during historical times. A model for their formation is presented. It takes into account the geomorphological data collected and some reports about the effect of recent tsunamis on coastal barriers and beaches. Washover fans were produced by tsunami waves which ran through coseismic cracks developed on dune ridges shaping a narrow, straight and relatively deep trench which constitutes the fan throat. Moreover, each tsunami event most likely caused severe erosion of the coastal barrier, shaping erosive grooves across the dune ridges, causing beach cliffs and causing the nourishment of submarine offshore bars. After the tsunami, a phase of coastal barrier recovery began, forming new dune ridges and closing washover fan throats. Morphological, archeological and radiometric data indicate a pre-Roman age for the oldest event, which was dated at 2430 years BP. The second tsunami struck the Lesina coastal barrier with similar magnitude 1550 years BP; it was caused by the strong earthquake that occurred at Gargano Promontory in the year 493 AD as reported by a medieval sacred legend. The smallest and more recent fan formed following the tsunami that hit the northern coast of Gargano on 30 July 1627

GTE. A new software for gravitational terrain effect computation: theory and performances

The computation of the vertical attraction due to the topographic masses, the so-called Terrain Correction, is a fundamental step in geodetic and geophysical applications: it is required in high-precision geoid estimation by means of the remove–restore technique and it is used to isolate the gravitational effect of anomalous masses in geophysical exploration. The increasing resolution of recently developed digital terrain models, the increasing number of observation points due to extensive use of airborne gravimetry in geophysical exploration and the increasing accuracy of gravity data represents nowadays major issues for the terrain correction computation. Classical methods such as prism or point masses approximations are indeed too slow while Fourier based techniques are usually too approximate for the required accuracy. In this work a new software, called Gravity Terrain Effects (GTE), developed to guarantee high accuracy and fast computation of terrain corrections is presented. GTE has been thought expressly for geophysical applications allowing the computation not only of the effect of topographic and bathymetric masses but also those due to sedimentary layers or to the Earth crust-mantle discontinuity (the so-called Moho). In the present contribution, after recalling the main classical algorithms for the computation of the terrain correction we summarize the basic theory of the software and its practical implementation. Some tests to prove its performances are also described showing GTE capability to compute high accurate terrain corrections in a very short time: results obtained for a real airborne survey with GTE ranges between few hours and few minutes, according to the GTE profile used, with differences with respect to both planar and spherical computations (performed by prism and tesseroid respectively) of the order of 0.02 mGal even when using fastest profiles

Quantitative analysis of extensional joints in the southern Adriatic foreland (Italy), and the active tectonics of the Apulia region

The Adriatic foreland of the Apennines comes ashore only in Apulia (easternmost Italy). Its southern part, our study area, lacks any structural analysis devoted to define its recent-to-active tectonics. Throughout the Quaternary, this region was affected by mild brittle deformation with rare faults, characterized by small displacement, and widespread extension joints, frequently organized in sets. Therefore, we conducted a quantitative and systematic analysis of the joint sets affecting Quaternary deposits, by applying an inversion technique ad hoc to infer the orientation and ratio of the principal stress axes, R = (σ2 - σ3)/(σ1 - σ3). Within a general extensional regime, we recognized three deformational events of regional significance. The oldest event, constrained to the early and middle part of the Middle Pleistocene, is characterized by variable direction of extension and R between 0.64-0.99. The penultimate event, dated late Middle Pleistocene, is characterized by an almost uniaxial tension, with a horizontal σ3 striking ~N43°E; R is high, between 0.85-0.99. The most recent event is characterized by the lowermost R values, that never exceed 0.47 and are frequently <0.30, indicating a sort of horizontal „radial‟ extension. This event is not older than the Late Pleistocene and possibly reflects the active stress field still dominating the entire study area

Quantitative analysis of extensional joints in the southern Adriatic foreland (Italy), and the active tectonics of the Apulia region

The Adriatic foreland of the Apennines comes ashore only in Apulia (easternmost Italy). Its southern part, our study area, lacks any structural analysis devoted to define its recent-to-active tectonics. Throughout the Quaternary, this region was affected by mild brittle deformation with rare faults, characterized by small displacement, and widespread extension joints, frequently organized in sets. Therefore, we conducted a quantitative and systematic analysis of the joint sets affecting Quaternary deposits, by applying an inversion technique ad hoc to infer the orientation and ratio of the principal stress axes, R = (σ2 − σ3)/(σ1 − σ3). Within a general extensional regime, we recognized three deformational events of regional significance. The oldest event, constrained to the early and middle part of the Middle Pleistocene, is characterized by variable direction of extension and R between 0.64 and 0.99. The penultimate event, dated late Middle Pleistocene, is characterized by an almost uniaxial tension, with a horizontal σ3 striking ∼N43°E; R is high, between 0.85 and 0.99. The most recent event is characterized by the lowermost R values, that never exceed 0.47 and are frequently <0.30, indicating a sort of horizontal ‘radial’ extension. This event is not older than the Late Pleistocene and possibly reflects the active stress field still dominating the entire study area

Late Quaternary deformation of the southern Adriatic foreland (southern Apulia) from mesostructural data: preliminary results

Southern Apulia (Adriatic foreland, Italy), has long been considered a “stable area” lying in between two active orogens, but in fact its tectonic framework is poorly known. To learn more about this topic, we carried out an original structural analysis on Pleistocene deposits. The results indicate that southern Apulia has been affected by mild but discernible brittle deformation throughout the Middle and Late Pleistocene. Joints prevail, whereas faults are rare and all characterized by small displacement. Horizontal extension dominates throughout the entire study area; the SW-NE to SSW-NNE direction is the most widespread. WNW-ESE extension prevails in the Adriatic side portion of the study area, but the dispersion of the measured plane directions is high, suggesting that the local strain field is not characterized by a strongly predominant trend. A Middle and Late Pleistocene, SW-NE to SSW-NNE– oriented maximum extension is not surprising for the study area, as it is compatible with most of the available geodynamic models, whereas the different state of deformation affecting the Adriatic side of the study area requires further investigations. We tentatively interpreted this anomaly as reflecting some regional variation of the general geodynamic frame, for instance as the farthest evidence of ongoing compressional deformation across the W-verging Albanide-Hellenide foldand- thrust belt

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the ob-served data. However, the assumption that the spatial dependence is constant through-out the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage fordealing with non-stationarity in geodetic data. We then compared stationary and non-stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC

Error sources and data limitations for the prediction ofsurface gravity: a case study using benchmarks

Gravity-based heights require gravity values at levelled benchmarks (BMs), whichsometimes have to be predicted from surrounding observations. We use EGM2008 andthe Australian National Gravity Database (ANGD) as examples of model and terrestrialobserved data respectively to predict gravity at Australian national levelling network(ANLN) BMs. The aim is to quantify errors that may propagate into the predicted BMgravity values and then into gravimetric height corrections (HCs). Our results indicatethat an approximate ±1 arc-minute horizontal position error of the BMs causesmaximum errors in EGM2008 BM gravity of ~ 22 mGal (~55 mm in the HC at ~2200 melevation) and ~18 mGal for ANGD BM gravity because the values are not computed atthe true location of the BM. We use RTM (residual terrain modelling) techniques toshow that ~50% of EGM2008 BM gravity error in a moderately mountainous regioncan be accounted for by signal omission. Non-representative sampling of ANGDgravity in this region may cause errors of up to 50 mGals (~120 mm for the Helmertorthometric correction at ~2200 m elevation). For modelled gravity at BMs to beviable, levelling networks need horizontal BM positions accurate to a few metres, whileRTM techniques can be used to reduce signal omission error. Unrepresentative gravitysampling in mountains can be remedied by denser and more representative re-surveys,and/or gravity can be forward modelled into regions of sparser gravity

Closed-Form transformation between geodetic and ellipsoidal coordinates

We present formulas for direct closed-form transformation between geodetic coordinates(Φ, λ, h) and ellipsoidal coordinates (β, λ, u) for any oblate ellipsoid of revolution.These will be useful for those dealing with ellipsoidal representations of the Earth's gravityfield or other oblate ellipsoidal figures. The numerical stability of the transformations for nearpolarand near-equatorial regions is also considered

The role of large earthquakes and tsunami in the Late Holocene evolution of Fortore River coastal plain (Apulia, Italy): a synthesis

Morphological analysis of the Fortore River coastal plain and the Lesina Lake coastal barrier integrated with radiocarbon age data indicates that the evolution of the coastal landscape has been strongly affected by a number of strong earthquakes and related tsunamis which occurred during the last 3000 years. The first seismic event struck this coastal area in the V century BC. It produced strong erosion of the Fortore River coastal plain and significant emersion of Punta delle Pietre Nere, as well as the large tsunami responsible for the development of the Sant'Andrea washover fan. The second event occurred in 493 AD; it induced severe erosion of the Fortore River coastal plain and triggered the large tsunami that hit the Lesina Lake coastal barrier, producing the Foce Cauto washover fan. Then later in 1627, an earthquake was responsible for the further coseismic uplift of Punta delle Pietre Nere, the subsidence of Lesina village area and the development of a tsunami which produced two washover fans. Morphological analysis points out that seismic events strong enough to control the morphological evolution of local coastal landscapes show a statistical return period of about 1000 years. These major events produced important coseismic vertical movements and large tsunamis. However, the correct identification of the tectonic structure responsible for the generation of these strong earthquakes is still an unsolved problem