The role of horizontal resolution in simulating drivers of the global hydrological cycle

The role of atmospheric general circulation model (AGCM) horizontal resolution in representing the global energy budget and hydrological cycle is assessed, with the aim of improving the understanding of model uncertainties in simulating the hydrological cycle. We use two AGCMs from the UK Met Office Hadley Centre: HadGEM1-A at resolutions ranging from 270 to 60 km, and HadGEM3-A ranging from 135 to 25 km. The models exhibit a stable hydrological cycle, although too intense compared to reanalyses and observations. This over-intensity is explained by excess surface shortwave radiation, a common error in general circulation models (GCMs). This result is insensitive to resolution. However, as resolution is increased, precipitation decreases over the ocean and increases over the land. This is associated with an increase in atmospheric moisture transport from ocean to land, which changes the partitioning of moisture fluxes that contribute to precipitation over land from less local to more non-local moisture sources. The results start to converge at 60-km resolution, which underlines the excessive reliance of the mean hydrological cycle on physical parametrization (local unresolved processes) versus model dynamics (large-scale resolved processes) in coarser HadGEM1 and HadGEM3 GCMs. This finding may be valid for other GCMs, showing the necessity to analyze other chains of GCMs that may become available in the future with such a range of horizontal resolutions. Our finding supports the hypothesis that heterogeneity in model parametrization is one of the underlying causes of model disagreement in the Coupled Model Intercomparison Project (CMIP) exercises.     

Long‐term demographic and reproductive trends in Ciona intestinalis sp. A

Ciona intestinalis sp. A (Tunicata, Ascidiacea) is a marine invertebrate with a widespread distribution and high invasive potential, yet little is known about its ecology. Here, the macrodynamics of demographic and reproductive trends were first determined using archive data covering 2002–2012. Species abundance in the central Tyrrhenian Sea (Campania, Italy ‐ NW Mediterranean) displayed Poisson distribution across 2002–2008, massive die‐off in 2009/2010 and population recovery in 2011/2012. We defined the seasonal ranges of sea surface temperature in which sampling size, animal size and reproductive status reached their highest values. In 2003/2004, self and non‐self fertilization appeared to follow reverse trends across the year. The long‐term recording run here on Ciona intestinalis sp. A provides novel insights into its population biology and ecology, and represents a valuable tool for sampling management and research planning.     

Improving Bayesian analysis for LISA Pathfinder using an efficient Markov Chain Monte Carlo method

We present a parameter estimation procedure based on a Bayesian framework by applying a Markov Chain Monte Carlo algorithm to the calibration of the dynamical parameters of the LISA Pathfinder satellite. The method is based on the Metropolis-Hastings algorithm and a two-stage annealing treatment in order to ensure an effective exploration of the parameter space at the beginning of the chain. We compare two versions of the algorithm with an application to a LISA Pathfinder data analysis problem. The two algorithms share the same heating strategy but with one moving in coordinate directions using proposals from a multivariate Gaussian distribution, while the other uses the natural logarithm of some parameters and proposes jumps in the eigen-space of the Fisher Information matrix. The algorithm proposing jumps in the eigen-space of the Fisher Information matrix demonstrates a higher acceptance rate and a slightly better convergence towards the equilibrium parameter distributions in the application to LISA Pathfinder data. For this experiment, we return parameter values that are all within ～1σ of the injected values. When we analyse the accuracy of our parameter estimation in terms of the effect they have on the force-per-unit of mass noise, we find that the induced errors are three orders of magnitude less than the expected experimental uncertainty in the power spectral density.     

New insights on diapirism in the Adriatic Sea: the Tremiti salt structure (Apulia offshore,southeastern Italy)

The reinterpretation of public seismic profiles in the Adriatic offshore of Gargano (Apulia, southern Italy) allowed the detection of a kilometre‐scale salt‐anticline, the Tremiti diapir, within the larger Tremiti Structure. This anticline was generated by diapirism of Upper Triassic anhydrites within a thick Mesozoic to Quaternary sedimentary succession. Both internal stratal patterns and shapes of Plio‐Quaternary units, and the occurrence of an angular unconformity between early Tortonian and Pliocene rocks on the Tremiti Islands, suggest that halokinesis began during the late Miocene and is still active today. An ancient extensional SE‐dipping fault, cutting an older Mesozoic low‐amplitude anhydritic ridge, played an important role during salt mobilization, which was promoted by NW‐SE shortening. The diapir grew in the footwall of this fault, causing its upward propagation. In some places, the ancient fault served as a preferential channel for the upward migration of the anhydrites.     

Detrital provenance of the Grenvillian Oaxacan Complex,southern Mexico: a zircon perspective

The Oaxacan Complex is the largest exposure of Grenvillian-age rocks in Mexico, constituting the backbone of the Oaxaquia microcontinent. Whereas the main rock-forming events were previously established at 1,150–1,200 Ma (charnockite–syenite–gabbros), 1,020 Ma (AMCG suite), 990 Ma (granulite-facies metamorphism), and ca. 970 Ma post-tectonic pegmatites, no data are yet available to establish provenance links with other Grenville-age terranes. In this work, we studied detrital zircons belonging to 12 samples, all metamorphosed under granulite facies but variably affected by retrogression. Laser ablation inductively coupled plasma mass spectrometry U–Pb geochronology was employed on selected zircons to determine their crystallization age and geochemistry. The results of the analysis of about 100 crystals per sample show that the studied zircons range between ca. 940 and 1,400 Ma, with only three samples having zircons between 1,400 and 1,600 Ma, and only one showing older zircons up to ca. 1,775 Ma. Whereas some of the slightly discordant (1–5 %) zircons in several samples show ages younger than the granulite metamorphism (probably as a result of Pb loss), and thus a disturbed geochemical pattern (abnormal enrichment in LREE, decreasing HREE), a few metamorphic zircons show flat and depleted HREE patterns, contrasting with the igneous pattern of older zircons (positive Ce anomaly, negative Eu anomaly, enriched HREE pattern). The main distributions observed using the kernel density estimator diagrams fall in the range 975–995 Ma (six samples), 1,100 Ma (four samples) and 1,120–1,170 Ma (six samples). Only the southernmost sample shows a marked peak at ca. 1,400 Ma. The application of the Kolmogorov–Smirnov (K–S) statistical test to the studied samples and particularly the comparison of obtained P values yield interesting similarities. Overall, two sample groups show internal similarities, i.e., they may belong to the same source area, whereas only one sample is dissimilar, failing to pass the K–S test. Comparison of these data with the timing of comparable events in the Sveconorwegian orogens, the Sunsas and Rondonia-San Ignacio belts of Amazonia, and some of the Precambrian massifs cropping out in the Andes help to constrain possible Mesoproterozoic conjugate margins of Oaxaquia.     

Accurate 3D surface measurement of mountain slopes through a fully automated image-based technique

In this paper an automated procedure for surface reconstruction that can be used for surveying and monitoring rock and ground slopes is presented. This method has been developed for geological and engineering applications, where accuracy and completeness are factors of primary importance for the final 3D model, which must provide a detailed metric survey and not only a visual reconstruction of the scene. The proposed procedure integrates two image matching techniques. The first one is used to automatically extract a set of tie points that are needed for computing the exterior orientation parameters of all images through a photogrammetric bundle adjustment. Such tie oints are also exploited to obtain a preliminary seed model that is then enriched based on Multi-photo Least Squares Matching. During this second stage, the surface model is improved in terms of point density and accuracy. Different strategies were implemented to successfully combine both techniques, along with some new improvements. The presented procedure has been tested in two different applications: geometric modelling of rock cliffs and evaluation of weathering of a ground slope. In both cases the obtained results presented accuracy sufficient for geological investigation. Moreover, outcomes were comparable to the ones from laser scanning surveying and other photogrammetric implementations.     

Geometric and kinematic characterization of landslides affecting urban areas: the Lungro case study (Calabria,Southern Italy)

The geometric and kinematic characterization of landslides affecting urban areas is a challenging goal that is routinely pursued via geological/geomorphological method and monitoring of ground displacements achieved by geotechnical and, more recently, advanced differential interferometric synthetic aperture radar (A-DInSAR) data. Although the integration of all the above-mentioned methods should be planned a priori to be more effective, datasets resulting from the independent use of these different methods are commonly available, thus making crucial the need for their standardized a posteriori integration. In this regard, the present paper aims to provide a contribution by introducing a procedure that, taking into account the specific limits of geological/geomorphological analyses and deep/surface ground displacement monitoring via geotechnical and A-DInSAR data, allows the a posteriori integration of the results by exploiting their complementarity for landslide characterization. The approach was tested in the urban area of Lungro village (Calabria region, southern Italy), which is characterized by complex geological/geomorphological settings, widespread landslides and peculiar urban fabric. In spite of the different level of information preliminarily available for each landslide as result of the independent use of the three methods, the implementation of the proposed procedure allowed a better understanding and typifying of the geometry and kinematics of 50 landslides. This provided part of the essential background for geotechnical landslide models to be used for slope stability analysis within landslide risk mitigation strategies.     

Petrophysical and mechanical properties of Euganean trachyte and implications for dimension stone decay and durability performance

Euganean trachyte is a subvolcanic porphyritic rock extracted in northern Italy with an age-old tradition of use as dimension stone, historically linked, in particular, to the fervent building activity brought by the Roman Empire and, later on, the Republic of Venice Serenissima. The results of a comprehensive petrophysical and mechanical characterization of Euganean trachyte from the most representative quarries are discussed here, involving the following properties: density, porosity, water absorption, capillary water uptake, hygroscopic water adsorption, hydric/hygric dilatation, water vapor diffusion, thermal expansion, and resistance to salt attack and abrasion. The different trachyte varieties, although belonging to the same quarry basin, exhibit a relatively wide array of technical performances, which are strongly dependent on pore volume, size, size distribution, shape, and degree of interconnection, controlling modes and rates of water transport and retention. Therefore, indications are provided for evaluating durability performance of the stone, with stress on water-driven weathering. Complementary information is finally given on the possible criteria followed in the antiquity for properly selecting the trachyte quarries to be exploited, and by a comparison with the properties of the most important trachytes extracted in Europe.     

Geochemistry, petrofabrics and seismic properties of eclogites from the Chinese Continental Scientific Drilling boreholes in the Sulu UHP terrane, eastern China

We present an integrated study of geochemistry, petrofabrics and seismic properties of strongly sheared eclogites from the Chinese Continental Scientific Drilling (CCSD) project in the Sulu ultrahigh-pressure (UHP) metamorphic terrane, eastern China. First, geochemical data characterize diverse protoliths of the studied eclogites. The positive Eu- and Sr-anomalies, negative Nb anomaly and flat portion of heavy rare earth elements in coarse-grained rutile eclogites (samples B270 and B295) suggest a cumulate origin in the continental crust, whereas the negative Nb anomaly and enrichment of light rare earth elements in retrograde eclogites (samples B504, B15 and B19) imply an origin of continental basalts or island arc basalts. Second, P-wave velocities (V_p) of three typical eclogite samples were measured under confining pressures up to 500 MPa and temperatures to 700 °C. At 500 MPa and room temperature, the mean V_p reaches 8.50-8.53 km/s in samples B270 and B295 but drops to 7.86 km/s in sample B504, and the P-wave anisotropy changes from 1.7-2.7% to 5.5%, respectively. The pressure and temperature derivatives of V_p are larger in the retrograde eclogite than in fresh ones. Third, the electron backscatter diffraction (EBSD) measurements of the eclogites reveal random crystal preferred orientation (CPO) of garnet and pronounced CPO of omphacite, which is characterized by a strong concentration of [001]-axes sub-parallel to the lineation and of (010)-poles perpendicular to the foliation. The asymmetric CPO of omphacite in sample B270 recorded a top-to-the-south shear event during subduction of the Yangtze plate. The calculated fastest V_p is generally sub-parallel to the lineation, but a different deformation environment during exhumation could form second-order variations in omphacite CPO and affect the V_p distribution in eclogites (e.g., the fastest V_p is at ~ 35° from the foliation in sample B295). Comparison between measured and calculated seismic properties indicates that the CPO of omphacite controls the seismic anisotropy of eclogites at high pressure, and compositional layering and retrograde minerals will increase the anisotropy. Calculated P-wave velocities agree well with velocities measured at 500 MPa and room temperature for fresh eclogites, but much higher than those of retrograde eclogite. As a case study, the laboratory-derived V_p-P and V_p-T relationships were used to estimate P-wave velocities of eclogites and peridotites beneath the Western Superior Province, Canada. The results indicate that besides the fabric-induced anisotropy, the direction dependence of pressure and temperature derivatives of V_p can significantly increase seismic anisotropy of eclogites with depth, which results in eclogites being an important candidate for the seismic anisotropy in the upper mantle. Due to their very high density and velocity, garnet-rich eclogites within peridotite could be detected in seismic reflections in subduction zones.