Time Evolution of Deformation Using Time Series of Differential Interferograms: Application to La Palma Island (Canary Islands)

Differential interferometry is a very powerful tool for detecting changes in the Earth’s crust where coherence conditions are good, but is difficult to employ in some volcanic areas due to dense vegetation. We apply two differential InSAR methods using the time series associated with the interferograms to perform a phase analysis on a data set for La Palma island (Canary Islands) from the ERS-1 and ERS-2 European Space Agency (ESA) satellites for the time period 1992 to 2000. Both methods involve choosing a master image from the database and creating a series of interferograms with respect to this image. The “Coherent Pixel Time Series” (CPTS) technique chooses pixels with good average coherence, aligns the unwrapped interferograms with a stable area and then performs an inversion to calculate the linear velocity to quantify the deformation. The Coherent Target Modeling (CTM) method calculates the temporal coherence of each pixel to identify stable targets and then determines the best velocity for each pixel by using a linear fit that maximizes the temporal coherence. Using these two methods we have been able to detect deformation on La Palma Island that has been previously undetectable by conventional InSAR methods. There is a roughly circular region on the Southern part of the island that is actively deforming at ~ −4 to −8 mm/yr. This region is located near the Teneguia valcano, the host of the last known eruption on La Palma in 1971. A thorough investigation of the possible sources for this deformation revealed that it was most likely created by a subsurface thermal source.     

On the Mathematical Analysis of an Elastic-gravitational Layered Earth Model for Magmatic Intrusion: The Stationary Case

In the early eighties Rundle (1980, 1981a,b, 1982) developed the techniques needed for calculations of displacements and gravity changes due to internal sources of strain in layered linear elastic-gravitational media. The approximation of the solution for the half space was obtained by using the propagator matrix technique. The Earth model considered is elastic-gravitational, composed of several homogeneous layers overlying a bottom half space. Two dislocation sources can be considered, representing magma intrusions and faults. In recent decades theoretical and computational extensions of that model have been developed by Rundle and co-workers (e.g., Fernández and Rundle, 1994a,b; Fernández et al., 1997, 2005a; Tiampo et al., 2004; Charco et al., 2006, 2007a,b). The source can be located at any depth in the media. In this work we prove that the perturbed equations representing the elastic-gravitational deformation problem, with the natural boundary and transmission conditions, leads to a well-posed problem even for varied domains and general data. We present constructive proof of the existence and we show the uniqueness and the continuous dependence with respect to the data of weak solutions of the coupled elastic-gravitational field equations.     

Evaluating methods to quantify sediment volumes trapped behind check dams,Salda(n)a badlands (Spain)

The determination of sediment yield in catchments based on the sediment trapped by check dams is becoming a subject of interest.In fact,several methods have been developed in recent years to estimate the sediment retained by check dams.The complexity,precision and accuracy of each method vary greatly.In this study,we evaluate the sediment trapped by check dams comparing the Sections method respect to the Prism,Pyramid,DTMs and Trapezoid methods.We analyzed a sample of 25 check dams(α-1=90%;ε=10%) in the Saldana badlands(Spain).The results showed that the Sections method offered a volume of retained sediment between the others,which gave an absolute variation from 22%to40%.The high variability of the check dam and sediment wedge sizes made necessary to compare methods in groups combining both characteristics.No significant differences in sediment volumes could be found between the methods for very small(height(h):1.8-2.3 m;trapped volume(V):6-102 m~3) or large check dams(h 2.3 m;V:165-387 m~3),while significant differences are found for small(h 1.5 m;V:1-229 m~3) or middle-sized check dams(h:2.2-3.2 m;V:65-235 m3).Nevertheless,volume differences between groups ranged up to + 25%.For these reason,the size of the check dam,the shape of the sediment wedge and the accuracy of the measure methods must considered when selecting an appropriate method to obtain the volume of retained sediment by check dams.A correct estimation of the sediment retention is needed to evaluate the role and efficiency of check dams in restoration projects or to estimate sediment yields.     

Factors determining subsidence in urbanized floodplains: evidence from MT‐InSAR in Seville (southern Spain)

Major rivers have traditionally been linked with important human settlements throughout history. The growth of cities over recent river deposits makes necessary the use of multidisciplinary approaches to characterize the evolution of drainage networks in urbanized areas. Since under‐consolidated fluvial sediments are especially sensitive to compaction, their spatial distribution, thickness, and mechanical behavior must be studied. Here, we report on subsidence in the city of Seville (Southern Spain) between 2003 and 2010, through the analysis of the results obtained with the Multi‐Temporal InSAR (MT‐InSAR) technique. In addition, the temporal evolution of the subsidence is correlated with the rainfall, the river water column and the piezometric level. Finally, we characterize the geotechnical parameters of the fluvial sediments and calculate the theoretical settlement in the most representative sectors. Deformation maps clearly indicate that the spatial extent of subsidence is controlled by the distribution of under‐consolidated fine‐grained fluvial sediments at heights comprised in the range of river level variation. This is clearly evident at the western margin of the river and the surroundings of its tributaries, and differs from rainfall results as consequence of the anthropic regulation of the river. On the other hand, this influence is not detected at the eastern margin due to the shallow presence of coarse‐grain consolidated sediments of different terrace levels. The derived results prove valuable for implementing urban planning strategies, and the InSAR technique can therefore be considered as a complementary tool to help unravel the subsidence tendency of cities located over under‐consolidated fluvial deposits. Copyright © 2017 John Wiley & Sons, Ltd.     

Geomorphological changes in an arid transgressive coastal dune field due to natural processes and human impacts

Geomorphological changes in recent decades in an arid transgressive coastal dune system (Maspalomas, Gran Canaria, Canary Islands) are analyzed. The methodology used is based on the generation of two geomorphological maps (1961 and 2003) by interpretation of digital orthophotos. The overlay of both maps in a geographic information system (GIS) enabled the spatial and surface changes of the landforms to be determined, and the processes that generated these changes. Twelve cultural and geomorphological processes were identified from highest to lowest importance, namely: anthropization by urban occupation (114 ha changed), stabilization (92.5 ha), barchanization (37 ha), salinization/halophytication (15 ha), anthropization (12.4 ha), deflation (11.8 ha), dune loss/beach gain (11.3 ha), dune formation (9.6 ha), progradation (8 ha), retrogradation (7.7 ha), destabilization (2.7 ha) and flooding (0.7 ha). Geomorphological changes are associated with a combination of five main factors, three of anthropogenic origin and two natural ones. The natural factors are: (1) the arid climate, which favors changes occurring at high speed; (2) the existence of a progressive sedimentary deficit. Anthropogenic factors are: (3) construction of tourist urbanizations, infrastructures and facilities; (4) installation of equipment or infrastructure on the beaches; (5) the activities carried out by users. These human factors have altered the aeolian dynamics and reduced the area occupied by vegetation in some areas, causing changes in aeolian sedimentary processes. The geomorphological processes identified can be used as indicators of environmental change, allowing us to synthesize the changes in landforms detected, and group all combinations derived from the analysis by GIS and analyze them spatially. Thus, the environmental changes in the transgressive coastal dune systems could be interpreted more effectively. Copyright © 2018 John Wiley & Sons, Ltd.     

Modelling feedbacks between geomorphological and riparian vegetation responses under climate change in a Mediterranean context

Climate change is expected to alter temperatures and precipitation patterns, affecting river flows and hence riparian corridors. In this context we have explored the potential evolution of riparian corridors under a dryness gradient of flow regimes associated with climate change in a Mediterranean river. We have applied an advanced bio‐hydromorphodynamic model incorporating interactions between hydro‐morphodynamics and vegetation. Five scenarios, representing drier conditions and more extreme events, and an additional reference scenario without climate change, have been designed and extended until the year 2100. The vegetation model assesses colonization, growth and mortality of Salicaceae species. We analysed the lower course of the Curueño River, a free flowing gravel bed river (NW Spain), as a representative case study of the Mediterranean region. Modelling results reveal that climate change will affect both channel morphology and riparian vegetation in terms of cover, age distribution and mortality. Reciprocal interactions between flow conditions and riparian species as bio‐engineers are predicted to promote channel narrowing, which becomes more pronounced as dryness increases. Reductions in seedling cover and increases in sapling and mature forest cover are predicted for all climate change scenarios compared with the reference scenario, and the suitable area for vegetation development declines and shifts towards lower floodplain elevations. Climate change also leads to younger vegetation becoming more subject to uprooting and flooding. The predicted reduction in suitable establishment areas and the narrowing of vegetated belts threatens the persistence of the current riparian community. This study highlights the usefulness of advanced bio‐hydromorphodynamic modelling for assessing climate change effects on fluvial landscapes. It also illustrates the need to consider climate change in river management to identify appropriate adaptation measures for riparian ecosystems. Copyright © 2018 John Wiley & Sons, Ltd.     

Linear inversion via the discrete wavelet transform pseudoinverse

Classical least‐squares techniques (Moore–Penrose pseudoinverse) are covariance based and are therefore unsuitable for the solution of very large‐scale linear systems in geophysical inversion due to the need of diagonalisation. In this paper, we present a methodology to perform the geophysical inversion of large‐scale linear systems via the discrete wavelet transform. The methodology consists of compressing the linear system matrix using the interesting properties of covariance‐free orthogonal transformations, to design an approximation of the Moore–Penrose pseudoinverse. We show the application of the discrete wavelet transform pseudoinverse to well‐conditioned and ill‐conditioned linear systems. We applied the methodology to a general‐purpose linear problem where the system matrix has been generated using geostatistical simulation techniques and also to a synthetic 2D gravimetric problem with two different geological set‐ups, in the noise‐free and noisy cases. In both cases, the discrete wavelet transform pseudoinverse can be applied to the original linear system and also to the linear systems of normal equations and minimum norm. The results are compared with those obtained via the Moore–Penrose and the discrete cosine transform pseudoinverses. The discrete wavelet transform and the discrete cosine transform pseudoinverses provide similar results and outperform the Moore–Penrose pseudoinverse, mainly in the presence of noise. In the case of well‐conditioned linear systems, this methodology is more efficient when applied to the least‐squares system and minimum norm system due to their higher condition number that allows for a more efficient compression of the system matrix. Also, in the case of ill‐conditioned systems with very high underdetermined character, the application of the discrete cosine transform to the minimum norm solution provides very good results. Both solutions might differ on their regularity, depending on the wavelet family that is adopted. These methods have a general character and can be applied to solve any linear inverse problem arising in technology, particularly in geophysics, and also to non‐linear inversion by linearisation of the forward operator.     

Monoclinal flexure of an orogenic plateau margin during subduction,south Turkey

Geologic evidence across orogenic plateau margins enables the discrimination of the relative contributions of orogenic, epeirogenic and/or climatic processes that lead to growth and maintenance of those plateaus and their margins. Here, we discuss the mode of formation of the southern margin of the Central Anatolian Plateau (SCAP) and evaluate its time of formation using fieldwork in the onshore and seismic reflection data in the offshore. In the onshore, uplifted Miocene rocks in a dip‐slope topography show monocline flexure over >100 km, km‐scale asymmetric folds verging south, and outcrop‐scale syn‐sedimentary reverse faults. On the Turkish shelf, vertical faults transect the basal latest Messinian of a 10 km fold where on‐structure syntectonic wedges and synsedimentary unconformities indicate pre‐Pliocene uplift and erosion, followed by Pliocene and younger deformation. Collectively, Miocene rocks delineate a flexural monocline at plateau margin scale that is expressed along our on‐offshore sections as a kink‐band fold with a steep flank 20–25 km long. In these reconstructed sections, we estimate a relative vertical displacement of 3.8 km at rates of ca. 0.5 mm/y, and horizontal shortening values <1 %. We use this evidence together with our observations of shortening at outcrop, basin, plateau‐margin and forearc‐system scales to infer that the SCAP forms as a monoclinal flexure to accommodate deep‐seated thickening and shortening since >5 Ma, and to contextualize the plateau margin as the forearc high of the Cyprus subduction system.     

The storm tracks and the energy cycle of the Southern Hemisphere: sensitivity to sea-ice boundary conditions

The effect of sea-ice on various aspects of the Southern Hemisphere (SH) extratropical climate is examined. Two simulations using the LMD GCM are performed: a control run with the observed sea-ice distribution and an anomaly run in which all SH sea-ice is replaced by open ocean. When sea-ice is removed, the mean sea level pressure displays anomalies predominantly negatives near the Antarctic coast. In general, the meridional temperature gradient is reduced over most of the Southern Ocean, the polar jet is weaker and the sea level pressure rises equatorward of the control ice edge. The high frequency filtered standard deviation of both the sea level pressure and the 300-hPa geopotential height decreases over the southern Pacific and southwestern Atlantic oceans, especially to the north of the ice edge (as prescribed in the control). In contrast, over the Indian Ocean the perturbed simulation exhibits less variability equatorward of about 50°S and increased variability to the south. The zonal averages of the zonal and eddy potential and kinetic energies were evaluated. The effect of removing sea-ice is to diminish the available potential energy of the mean zonal flow, the available potential energy of the perturbations, the kinetic energy of the growing disturbances and the kinetic energy of the mean zonal flow over most of the Southern Ocean. The zonally averaged intensity of the subpolar trough and the rate of the baroclinic energy conversions are also weaker.