Advanced interpretation of interferometric SAR data to detect,monitor and model ground subsidence: outcomes from the ESA-GMES Terrafirma project

Land subsidence is a common natural hazard striking extensive areas worldwide, with severe impacts on economy and environment. Subsidence has been recognized as one of geohazards needing research efforts and knowledge transfer at international level, especially when urban fabrics and infrastructures are directly involved in the land settling. Policies and solutions for land subsidence management can be different. Despite this variability, where mitigation methods need to be adopted, mapping, monitoring and simulation of subsidence have to precede their design and implementation. In this framework, Earth Observation (EO) and remote sensing have a major role to play. Satellite Synthetic Aperture Radar Interferometry, thanks to its wide spatial coverage and its millimeter accuracy, provides a valuable contribution in the management of hazard posed by subsidence-related deformation. The ESA-GMES Terrafirma project (2003–2014) has worked for the promotion of the persistent scatterer interferometry, a family of techniques ideally suited for the assessment of magnitude of surface deformations associated with subsidence phenomena. Within the Terrafirma Project a series of products, based on the integration of EO technologies and in situ data, has been established and delivered to a wide community of end user. Three case studies, outcomes of the Terrafirma project, are presented: the wide area of Rome (Italy), the Anthemountas basin and the Kalochori village (Greece). These case studies have been selected with the purpose of showing the essential contribution of interferometric data during the main activities that must be covered when dealing with geohazard investigations (i.e., mapping, monitoring and modeling). These three case studies are meant to be representative of the suite of services delivered by the Terrafirma project to specific end users with the legal mandated to manage the geohazard.     

Basin filling evolution of the central basins of Mallorca since the Pliocene

A new compilation of data from 436 drill cores using decompaction and backstripping techniques was used to reconstruct the basin filling history from the Pliocene until the present day in the Palma, Inca and Sa Pobla Basins on the island of Mallorca (Spain). Calcareous rocks dominate the source area and provide a limited amount of clastic input to the basins that has resulted in an average accumulation rate of between 5 and 20 m/Ma during the last 5.3 Ma. Carbonate sediment production dominated the basin filling history during early‐mid Pliocene, but during the Quaternary, the sedimentation processes in the Palma Basin were probably enhanced by an evolution in the drainage network that increased the sediment supply and the accumulated thickness caused by stream capture. However, the maximum sedimentation rate filling the depocentres of the three basins has been decreasing since the Pliocene, showing that not only the catchment transport efficiency but also the relative sea level have been controlling the sediment accumulation in these carbonate basins. The isopach cross‐sections support the idea that a palaeorelief was generated during the Messinian sea level drop and that heterogeneities were filled in from the Pliocene to the Quaternary. We conclude that the central basins of Mallorca were filled heterogeneously due to tectonic and geomorphic processes that controlled sediment transport and production, resulting in different average sedimentation thicknesses that decreased since the Pliocene as the accommodation space became filled and the relative sea level dropped.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

Responsible fisheries: An overview of FAO policy developments (1945–1994)

The paper outlines the evolution of fisheries development policies (aquaculture excluded), from an international perspective as seen within the various UN and FAO fishery committees. It begins with the first FAO Technical Committee on Fisheries in 1945 and extends through the various periods corresponding to post-war reconstruction (1945–1958), geographical expansion (1959–1972), establishment of a new economic order of the oceans (1973–1982), and transition towards global concerns (1983–1992). It underlines the gaps between government resolutions at the international level and effective implementation at national policy level. This identifies the elements to be considered in responsible fishing policies. Finally, the FAO International Code of Conduct for Responsible fisheries is presented as the conceptual framework for achieving sustainable fisheries development and UNCED's objectives within the framework of the 1992 Convention.     

Volatile contents of phlogopite micas from South African kimberlite

Phlogopite micas from nodules in South African kimberlites were analyzed for major elements with the electron microprobe and for volatile contents by high temperature mass spectrometry. The micas are from primary- (deformed) and secondary- (undeformed) textured grains in perodotite xenoliths, glimmerites, MARID (mica-amphibole-rutile-ilmenite-diopside) suite nodules and a mica megacryst. The major element and volatile contents of micas exhibiting these modes of occurrence overlap to a greater extent than indicated in previous studies. Concentrations of volatile species occupying structurally defined crystallographic sites (H₂O, F, Cl) are greater for many of the micas than predicted on the basis of the mica formula, particularly for the glimmerite and MARID suite samples. A correlation exists between micas with tetrahedral and octahedral cation deficiencies and those with excess H₂O, F and Cl. Substitution of H⁺ for tetrahedral and possibly octahedral cations may be responsible for the excess H₂O in these micas. Except for one sample, the major element and volatile data for the peridotite, glimmerite and MARID suite micas indicate that they crystallized at oxygen fugacities below the quartz-fayalite-magnetite buffer. F and K₂O are in the correct proportion in the micas to provide the source for these elements in alkali basalts, but not in mid-ocean ridge basalts. Kaersutite amphibole is a more likely source of potassium and fluorine in mid-ocean ridge basalts.     

Analysis of daily rainfall processes in lower extremadura (Spain) and homogenization of the data

Summary In this paper we analyze, from the point of view of stochastic processes, daily rainfall data recorded at the Badajoz Observatory (Southwestern Spain) since the beginning of the century. We attempt to identify any periodicities or trends in the daily rainfall occurrences and their dependence structure, and attempt to select the appropriate point stochastic model for the daily rainfall series.Standard regression analysis, graphical methods and the Cramer statistic show a rise in the number of cases of light rain (between 0.1 and 5 mm/d) and a decline in the number of cases of moderate to heavy rain (> 5 mm/d) in the daily rainfall at least at the 5% significance level.That the homogenization process was satisfactory is shown by the mean interarrival time of the homogenized series and the test of the rate of homogenized daily rainfall occurrences. Our analysis also shows that the behavior of the spectra of the homogenized daily rainfall counts is completely different from that of a Poisson process, so that the hypothesis of a non-homogeneous Poisson process is rejected.With 7 Figures     

Geochemical modeling of magma mixing and magma reservoir volumes during early episodes of Kīlauea Volcano’s Pu‘u ‘Ō‘ō eruption

Geochemical modeling of magma mixing allows for evaluation of volumes of magma storage reservoirs and magma plumbing configurations. A new analytical expression is derived for a simple two-component box-mixing model describing the proportions of mixing components in erupted lavas as a function of time. Four versions of this model are applied to a mixing trend spanning episodes 3–31 of Kilauea Volcano’s Puu Oo eruption, each testing different constraints on magma reservoir input and output fluxes. Unknown parameters (e.g., magma reservoir influx rate, initial reservoir volume) are optimized for each model using a non-linear least squares technique to fit model trends to geochemical time-series data. The modeled mixing trend closely reproduces the observed compositional trend. The two models that match measured lava effusion rates have constant magma input and output fluxes and suggest a large pre-mixing magma reservoir (46±2 and 49±1 million m³), with little or no volume change over time. This volume is much larger than a previous estimate for the shallow, dike-shaped magma reservoir under the Puu Oo vent, which grew from ∼3 to ∼10–12 million m³. These volumetric differences are interpreted as indicating that mixing occurred first in a larger, deeper reservoir before the magma was injected into the overlying smaller reservoir. Electronic Supplementary Material Supplementary material is available at and is accessible for authorized users.     

The Role of the Magnetic Field in the Ascension of the 22nd and 23rd Cycles of Activity

The first small active regions of the last two new cycles seem to be formed, in its very first, preparatory stage, close to the equator with uncertain magnetic orientation, and as secondary products at the periphery of equatorial `magnetic bubbles' of the new magnetic flux still belonging to the old cycle. At the real beginning of the new cycle, at the periphery of its higher-latitude `magnetic bubbles', the opposite applies to the secondary regions of the ending cycle. It seems that two modes of magnetic flux supply the visible layers of the photosphere during the cycle transition phase: the emergence of the magnetic flux from the depth of the convective zone, and the mutual interaction of the earlier emerging fluxes, or the action of the `local dynamo', giving rise to the new local magnetic field concentrations.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.