Application of multi-sensor advanced DInSAR analysis to severe land subsidence recognition: Alto Guadalentín Basin (Spain)

Multi-sensor advanced DInSAR analyses have been performed and compared with two GPS station measurements, in order to evaluate the land subsidence evolution in a 20-year period, in the Alto Guadalentín Basin where the highest rate of man-induced subsidence (> 10 cm yr−1) of Europe had been detected. The control mechanisms have been examined comparing the advanced DInSAR data with conditioning and triggering factors (i.e. isobaths of Plio-Quaternary deposits, soft soil thickness and piezometric level).

Twenty-year advanced DInSAR analysis of severe land subsidence: The Alto Guadalentín Basin (Spain) case study

A twenty-year period of severe land subsidence evolution in the Alto Guadalentín Basin (southeast Spain) is monitored using multi-sensor SAR images, processed by advanced differential interferometric synthetic aperture radar (DInSAR) techniques. The SAR images used in this study consist of four datasets acquired by ERS-1/2, ENVISAT, ALOS and COSMO-SkyMed satellites between 1992 and 2012. The integration of ground surface displacement maps retrieved for different time periods allows us to quantify up to 2.50 m of cumulated displacements that occurred between 1992 and 2012 in the Alto Guadalentín Basin. DInSAR results were locally compared with global positioning system (GPS) data available for two continuous stations located in the study area, demonstrating the high consistency of local vertical motion measurements between the two different surveying techniques. An average absolute error of 4.6 ± 4 mm for the ALOS data and of 4.8 ± 3.5 mm for the COSMO-SkyMed data confirmed the reliability of the analysis. The spatial analysis of DInSAR ground surface displacement reveals a direct correlation with the thickness of the compressible alluvial deposits. Detected ground subsidence in the past 20 years is most likely a consequence of a 100–200 m groundwater level drop that has persisted since the 1970s due to the overexploitation of the Alto Guadalentín aquifer system. The negative gradient of the pore pressure is responsible for the extremely slow consolidation of a very thick (> 100 m) layer of fine-grained silt and clay layers with low vertical hydraulic permeability (approximately 50 mm/h) wherein the maximum settlement has still not been reached.

Imaging Land Subsidence Induced by Groundwater Extraction in Beijing (China) Using Satellite Radar Interferometry

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.

Land use image classification through optimum-path forest clustering

Land use classification has been paramount in the last years, since we can identify illegal land use and also to monitor deforesting areas. Although one can find several research works in the literature that address this problem, we propose here the land use recognition by means of Optimum-Path Forest Clustering (OPF), which has never been applied to this context up to date. Experiments among Optimum-Path Forest, Mean Shift and K-Means demonstrated the robustness of OPF for automatic land use classification of images obtained by CBERS-2B and Ikonos-2 satellites. © 2011 IEEE.

Analysis of the long-term subsidence in the urban area of Bologna

Land subsidence in urban areas represents a widespread geological hazard and a pressing challenge for modern society. This research focuses on the subsidence process affecting the city of Bologna (Italy). Since the 1960s, Bologna has experienced ground deformation due to aquifers overexploitation that peaked during the 1970s with rates of 10 cm/year. Despite a general reduction in these rates over the subsequent decades, thanks to groundwater regulations policies, recent data underscore a substantial subsidence resurgence. To reconstruct the subsurface stratigraphic architecture of Bologna’s urban area and generate a 3D geological model, a multidisciplinary approach centred on a stratigraphic analysis relying on the lithofacies criterion was adopted. The convergence of the analyses within this framework resulted in partitioning the study area into three geological domains exhibiting unique morphological features and depositional stacking patterns. Subsequently, since long-term data are crucial for a comprehensive understanding of ongoing subsidence, a methodology was developed to generate cumulative ground displacement time series and maps by integrating ground-based and remotely-sensed measurements. While the reconstructed long-term subsidence field consistently aligns with the primary geological variations summarised in the 3D model, the generated cumulative displacement curves systematically match pluriannual trends observed in groundwater level and pumping monitoring data. Lastly, to evaluate the expression of the observed relationships from a geotechnical perspective, a series of one-dimensional subsidence calculations were conducted considering the mechanical properties of the investigated deposits and piezometric data. These analyses provided valuable insight into the overall mechanical behaviour of the existing soils, as well as the post-pumping groundwater level and pore pressure distributions, consistent with field data. The methodological approach employed enables a comprehensive analysis of land subsidence in urban areas, allowing the exploration of individual factors governing the deformation process and their interactions, even within complex stratigraphic and hydrogeological environments such as Bologna’s urban area.

Changement identitaire et revendications régionalistes du Kansaï au Japon

Depuis quelque temps, au Japon, on utilise de plus en plus le terme « Kansaï » pour désigner la région du Kinki (littéralement « le voisinage de la capitale »). Cette thèse propose d’analyser l’émergence de cette entité régionale et de son discours culturel dans le but de pallier le manque de recherches antérieures sur la diversité socioculturelle et le régionalisme au Japon. Il y existe, d’une part, une volonté de considérer le Japon comme une entité homogène, et d’autre part, un contexte dans lequel le Japon lui-même prône l’homogénéité de son peuple. Historiquement, ces énoncés ont été réfutés à plusieurs reprises par différents chercheurs et organismes. Entre-temps, sur le plan régional, la diversité devient de moins en moins clairement observable dû à l’urbanisation, aux moyens de transport, à la migration interne et au développement des médias de masse. Cette situation à l'époque post-industrielle a engendré aujourd’hui le discours régionaliste du Kansaï. Dans ce contexte, cette étude porte spécifiquement sur le discours culturel concernant la région et la population du Kansaï, c’est-à-dire la région Kinki, où étaient situés les anciennes capitales et le berceau de l’État japonais du Yamato. On observe une modification et une transformation de cette région depuis l’époque Tokugawa. À partir de l’époque Meiji, l'intégration spatiale de l’archipel japonais est devenue indissociable de l’émergence de l’État soi-disant « moderne ». En outre, une distinction existe toujours entre le Japon de l’Ouest (Kansaï) et le Japon de l’Est (Kantō) qui repose sur des différences de coutumes et de mentalités, ainsi que sur des variations linguistiques : une dichotomie mieux représentée de nos jours par l’opposition entre les villes d’Osaka et de Tokyo. Aujourd’hui, le Japon permettre une centralisation continuelle à Tokyo et l’équilibre du pouvoir sur le plan économique s’en trouve fragilisé. Dans cette thèse, j’examine l’émergence de l’entité Kansaï dans ce contexte socio-économique, depuis l’arrivée du phénomène que les Kansaïens appellent « l’affaissement de terrain » du Kansaï, le jibanchinka, jusqu’aux revendications récentes pour l’introduction d’un système quasi-fédéraliste, le dōshū-sei, dans le contexte du développement régional déséquilibré du pays. En m'appuyant sur mon enquête effectuée sur terrain auprès des gens du Kansaï, je soutiens que leur discours régionaliste est bel et bien existant, mais ne repose pas sur l’homogénéité de la région. Il repose plutôt sur la position du Kansaï en tant qu’antithèse à la tendance centralisatrice perçue par les Kansaïens comme étant plutôt de nature tokyoïte. Leur discours met l’accent sur la diversité existant à l’intérieur même de la région tout en soulignant que celle-ci constitue l’entité kansaïenne. Mots-clés : Japon, Kinki, Kansai, Osaka, Nihonjinron, région, villes, discours culturel, État-nation, multiculturalisme, Oda Sakunosuke, Tanizaki Jun’ichiro.

Late Quaternary sedimentation in the Paraiba Basin, Northeastern Brazil: Landform, sea level and tectonics in Eastern South America passive margin

Late Quaternary deposits in the northeastern Brazil have been scarcely investigated, despite their relevance to the discussion of the post-rift evolution of the South American passive margin within the context of landform, sea level and tectonic deformation. Sedimentological, stratigraphic and morphological characterization of these deposits, referred as Post-Barreiras Sediments, led to their distinction from underlying Early/Middle Miocene strata. Based on optically stimulated luminescence (OSL) dating, two sedimentary units (PB1 and PB2) were recognized and related to the time intervals between 74.8 +/- 9.3 and 30.8 +/- 6.9 ka, and 8.8 +/- 0.9 and 1.8 +/- 0.2 ka, respectively. Unit PB1 consists of indurated sandstones and breccias either with massive bedding or complex types of soft sediment deformation structures generated by contemporaneous seismic activity. Unit PB2 is composed of massive sands or sands related to structures developed by dissipation of dunes. The present work, focusing on the Post-Barreiras Sediments, discusses landform, sea level and tectonics of the eastern South American passive margin during the latest Quaternary. Non-deposition and sub-aerial exposure related to the Tortonian worldwide low sea level combined with tectonic quiescence followed the Miocene transgression. Tectonic deformation in the latest Pleistocene created space to accommodate unit PB1 in downthrown faulted blocks and, perhaps, also synclines produced by strike-slip deformation. Although deposition of this unit was simultaneous with the progressive fall in sea level that followed the Last Interglacial Maximum, punctuated rises combined with land subsidence led to marine deposition close to the modern coastline. Renewed subsidence in the Holocene gave rise to accommodation of the Post-Barreiras Sediments. Most of unit PB2 was deposited during the Holocene Transgression, but it is not composed of marine sediments, which suggests either an insignificant rise in relative sea level or aeolian reworking of thin transgressive sands. The data presented here lead to a review of the evolution of the South American passive margin based on assumptions of uniform sedimentation and undeformed planation surfaces over a wide coastal area of the northeastern Brazil. (C) 2011 Elsevier B.V. All rights reserved.

Wave height (regular and irregular waves) reduction in 2m water depth over a 40-m-long salt marsh test section in a large scale laboratory flume

Coastal communities around the world face increasing risk from flooding as a result of rising sea level, increasing storminess, and land subsidence. Salt marshes can act as natural buffer zones, providing protection from waves during storms. However, the effectiveness of marshes in protecting the coastline during extreme events when water levels and waves are highest is poorly understood. Here, we experimentally assess wave dissipation under storm surge conditions in a 300-m-long wave flume that contains a transplanted section of natural salt marsh. We find that the presence of marsh vegetation causes considerable wave attenuation, even when water levels and waves are high. From a comparison with experiments without vegetation, we estimate that up to 60% of observed wave reduction is attributed to vegetation. We also find that although waves progressively flatten and break vegetation stems and thereby reduce dissipation, the marsh substrate remained remarkably stable and resistant to surface erosion under all conditions.The effectiveness of storm wave dissipation and the resilience of tidal marshes even at extreme conditions suggest that salt marsh ecosystems can be a valuable component of coastal protection schemes.

New assessments of global sea level rise from tide gauges

Estimates of global sea-level change rates based on observations from Tide Gauges (TGs) show a long-term global mean sea-level rise (GMSLR) of 1÷2 mm/yr for the 20th century. The considerable scatter in these estimates is mainly attributable to the uneven distribution of the TG sites and to several physical phenomena that cause local sea level to deviate from the global mean, or to affect the TG record through land subsidence or uplift. The main cause of vertical ground motion on a regional space scale is the response of the Earth to past ice loads, called Glacial Isostatic Adjustment (GIA), which is often modelled and corrected for. In this work, a simple average approach was used to revisit two past estimates based on small sets of long, high-quality TG records in view of the longer record available, employing a newer GIA model (ICE-6G) from Peltier et al. [2015]. The value of GMSLR obtained from both sets is (1.5±0.4) mm/yr. In addition, a much larger set of TGs was used to estimate the contemporary (post 1993) GMSLR using satellite estimates from Cazenave et al. [2018] as a benchmark, in an attempt to understand how a simple average approach could perform for larger sets. The resulting estimate of (3.4÷3.5)±0.2 mm/yr (depending on the GIA correction applied) is comparable to the satellite result of (3.1±0.3) mm/yr.

Use of synthetic aperture radar for recognition of Coastal Geomorphological Features, land-use assessment and shoreline changes in Bragança coast, Pará, Northern Brazil

Imagens de radar de abertura sintética (SAR) vem sendo bem mais utilizadas do que antes nas aplicações de geociências em regiões tropicais úmidas. Nesta investigação, uma imagem RADARSAT-1, na banda C, polarização HH adquirida em 1998 foi usada para o mapeamento costeiro e avaliação da cobertura da terra na área de Bragança, norte do Brasil. Imagem do radar aerotransportado GEMS-1000, na banda X, polarização HH, adquirida em 1972 durante o projeto RADAM foi também utilizada para avaliar as variações costeiras ocorridas nas últimas três décadas. A pesquisa tem confirmado a utilidade da imagem RADARSAT-1 para o mapeamento geomorfológico e avaliação da cobertura da terra, particularmente em costas de manguezal de macromaré. Além disso, um novo método para estimar as variações da linha de costa baseado na superposição de vetores extraídos de diferentes imagens SAR, com alta acurácia geométrica, tem mostrado que a planície costeira de Bragança tem estado sujeita a severa erosão responsável pelo recuo de aproximadamente 32 km² e acreção de 20 km², resultando em uma perda de área de manguezal de aproximadamente 12 km². Como perspectiva de aplicação, dados SAR orbitais e aerotransportados provaram ser uma importante fonte de informação tanto para o mapeamento geomorfológico, quando para o monitoramento de modificações costeiras em ambientes tropicais úmidos.

Facies architecture of the CRP-3 drillhole, Victoria Land Basin, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diametites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ("Motif A"). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as "Motif B. Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ("o", "p" and "q") correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

Laser-derived particle size data from CRP-3, Victoria Land Basin, Antarctica: Implications for sequence and seismic stratigraphy

Seven hundred and nineteen samples from throughout the Cainozoic section in CRP-3 were analysed by a Malvern Mastersizer laser particle analyser, in order to derive a stratigraphic distribution of grain-size parameters downhole. Entropy analysis of these data (using the method of Woolfe and Michibayashi, 1995) allowed recognition of four groups of samples, each group characterised by a distinctive grain-size distribution. Group 1, which shows a multi-modal distribution, corresponds to mudrocks, interbedded mudrock/sandstone facies, muddy sandstones and diamictites. Group 2, with a sand-grade mode but showing wide dispersion of particle size, corresponds to muddy sandstones, a few cleaner sandstones and some conglomerates. Group 3 and Group 4 are also sand-dominated, with better grain-size sorting, and correspond to clean, well-washed sandstones of varying mean grain-size (medium and fine modes, respectively). The downhole disappearance of Group 1, and dominance of Groups 3 and 4 reflect a concomitant change from mudrock- and diamictite-rich lithology to a section dominated by clean, well-washed sandstones with minor conglomerates. Progressive downhole increases in percentage sand and principal mode also reflect these changes. Significant shifts in grain-size parameters and entropy group membership were noted across sequence boundaries and seismic reflectors, as recognised in others studies.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Idealized model for changes in equilibrium temperature, mixed layer depth, and boundary layer cloud over land in a doubled CO2 climate

An idealized equilibrium model for the undisturbed partly cloudy boundary layer (BL) is used as a framework to explore the coupling of the energy, water, and carbon cycles over land in midlatitudes and show the sensitivity to the clear‐sky shortwave flux, the midtropospheric temperature, moisture, CO2, and subsidence. The changes in the surface fluxes, the BL equilibrium, and cloud cover are shown for a warmer, doubled CO2 climate. Reduced stomatal conductance in a simple vegetation model amplifies the background 2 K ocean temperature rise to an (unrealistically large) 6 K increase in near‐surface temperature over land, with a corresponding drop of near‐surface relative humidity of about 19%, and a rise of cloud base of about 70 hPa. Cloud changes depend strongly on changes of mean subsidence; but evaporative fraction (EF) decreases. EF is almost uniquely related to mixed layer (ML) depth, independent of background forcing climate. This suggests that it might be possible to infer EF for heterogeneous landscapes from ML depth. The asymmetry of increased evaporation over the oceans and reduced transpiration over land increases in a warmer doubled CO2 climate.