Reply to Molina-Garza et al. (2019) “Discussion of: Ortega-Flores et al. (2018) provenance analysis of Oligocene sandstone from the Cerro Pelón area,southern Gulf of Mexico”

ABSTRACT

The origin of the Oligocene turbidites from the Cerro Pelón area in south Gulf Mexico proposed by Ortega-Flores et al. (2018) is in disagree with the interpretations made by Molina-Garza et al. (2019), which main criticism is based on U-Pb ages of detrital zircons from the matrix of a conglomerate unit, which they refer to as ‘Nanchital Conglomerate’, as well as on the presence of limestone, gabbros, and mafic protolith-derived clasts. Molina-Garza et al. (2019) basically interpret the Nanchital Conglomerate as Miocene in age, which was sourced mainly from metamorphic complexes including their sedimentary covers located to the west and south of the Cerro Pelón area. For some reason, Molina-Garza et al. (2019) suppose that the Nanchital Conglomerate should have the same provenance sources that the Oligocene turbidites from Cerro Pelón area, reported by Ortega-Flores et al. (2018). Based on the foregoing, we strongly disagree with Molina-Garza et al. (2019) considering that, from the beginning, they intend to compare two units of different age. Additionally, the scarce data reported from both the matrix and the clasts of the Nanchital Conglomerate are not determinant for interpreting the provenance of this conglomeratic unit and subsequently, to consider the same rock sources from the Oligocene through Miocene time.     

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.