Spatial and temporal patterns of dry spell lengths in the Iberian Peninsula for the second half of the twentieth century

Summary Daily pluviometric records of 43 meteorological stations across the Iberian Peninsula have permitted a detailed analysis of dry spell patterns for the period 1951–2000 by distinguishing daily amount thresholds of 0.1, 1.0, 5.0 and 10.0 mm/day. The analyses are based on three annual series, namely the number of dry spells, N, the average dry spell length, L, and the extreme dry spell length, L _max. First, the statistical significance of local trends for the annual series of N, L and L _max has been investigated by means of the Mann-Kendall test and significant field trends have been established by means of Monte Carlo simulations. Clear signs of negative field trends are detected for N (1.0 and 10.0 mm/day) and L (0.1 mm/day). Second, the Weibull model fits well the empirical distributions of dry spell lengths for all the rain gauges, whatever the daily amount threshold, with a well ranged spatial distribution of their parameters u and k. On the basis of the Weibull distribution, return period maps for 2, 5, 10, 25 and 50 years have been obtained for dry spell lengths with respect to the four daily threshold levels. While for 0.1 and 1.0 mm/day the longest dry spells are expected at the south of the Iberian Peninsula, for 5.0 and 10.0 mm/day they are mostly detected at the southeast. Finally, the elapsed time between consecutive dry spells has been analysed by considering the same rain amount thresholds and different dry spell lengths at increasing intervals of 10 days. This analysis makes evident a significant negative field trend of the elapsed time between consecutive dry spells of lengths ranging from 10 to 20 days for daily amount thresholds of 1.0, 5.0 and 10.0 mm/day. Authors’ addresses: X. Lana, C. Serra, Departament de Física i Enginyeria Nuclear, ETSEIB, Universitat Politècnica de Catalunya, Av. Diagonal 647 planta 11, 08028 Barcelona, Spain; M. D. Marínez, Departament de Física Aplicada, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain; A. Burgue?o, Departament de Meteorologia i Astronomia, Universitat de Barcelona, 08028 Barcelona, Spain; J. Martín-Vide, L. Gómez, Grup de Climatologia, Universitat de Barcelona, 08028 Barcelona, Spain.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

Geometry and accuracy of reflecting points in bistatic satellite altimetry

The analysis of the time and space distribution of specular (reflecting) points in bistatic altimetry between GPS and CHAMP satellites or SAC-C (taken as examples) is extended from Wagner and Klokočník (2003 J. Geod 77: 128–138). We demonstrate a significantly higher number and density of reflecting points in bistatic altimetry in comparison with traditional monostatic altimetry. After an outline of our older accuracy assessment for the vertical position of the reflecting point, we add a new independent derivation and compare both approaches. We account for orbit errors of both the transmitters (GPS) and receiver (CHAMP) satellites, and the measurement (delay) error. We found that the accuracy of the vertical position of the reflecting point decreases only slowly with increasing off-nadir angle and that the orbit errors must be accounted for if decimeter and better accuracy is required. In this paper, we do not study errors such as state of the ocean, technical parameters of the receiving system, and atmospheric corrections.     

The Huelva (Spain) Tsunami-Ready Station and Its Interaction with Storm Emma (March 2018)

Abstract

The sea level station operating since 1996 at Mazagón (Huelva, Spain) has been progressively upgraded to fit tsunami warning requirements, due to its location in one of the main regions at risk. Its radar water level sensor was complemented in 2017, with the addition of a pressure sensor. The performance of both sea level sensors and their response to sea level oscillations, at different frequencies, is assessed. Particular emphasis is put on the effect of extreme events, such as Storm Emma, when alternative methods to obtain 1-min data are tested, in contrast to the one based on arithmetic means. The overall differences are small, for the whole period of study (centered-root-mean-square-error below 1?cm, for 5-min, and hourly data; similar tidal parameters and sea level oscillations with periods between 30?s and 5?min). However, during Storm Emma, the pressure sensor presents sensibly lower readings than the radar, with the centered-root-mean-square-error rising to 80?mm on the March 2^nd 2018. A new method to compute 1-min data, based on medians, reduced this value to 10?mm for the same day.     

Modeling rhamnolipids production by Pseudomonas aeruginosa from immiscible carbon source in a batch system

A mathematical model to predict the rhamnolipids production by Pseudomonas aeruginosa from oleic acid in a two phase liquid-liquid batch reaction system, was developed in this study. The model was based on two theoretical assumptions: 1) the convective oleic acid mass transfer is coupled to a bioreaction in the aqueous liquid bulk, and 2) the volume of the immiscible oleic acid drops and the saturation concentration at the interface are a function of rhamnolipids production. The model was able to accurately predict the experimental growth of the Pseudomonas aeruginosa strain, and the rhamnolipids production data with oleic acid as carbon source. This mathematical approach indicated a high correspondence between the saturation dimensionless profiles of oleic acid at the interface and the experimental profiles of surface tension difference. This modeling approach may constitute a useful tool in the design and scaling-up of bioreactors applied to the production of biosurfactants with immiscible carbon sources.     

Groundwater flow systems in turbidites of the Northern Apennines (Italy): natural discharge and high speed railway tunnel drainage

Turbidites crop out extensively in the Northern Apennine mountains (Italy). The huge amounts of groundwater drained by tunnels, built for the high speed railway connection between Bologna and Florence, demonstrate the aquifer-like behaviour of these units, up to now considered as aquitards. A conceptual model of groundwater flow systems (GFS) in fractured aquifers of turbidites is proposed, taking into account both system natural state and the perturbation induced by tunnel drainage. Analysis of hydrological data (springs, streams and tunnel discharge), collected over 10 years, was integrated with analysis of hydrochemical and isotopic data and a stream-tunnel tracer test. Hydrologic recession analysis of undisturbed conditions is a key tool in studying turbiditic aquifer hydrogeology, permitting the discrimination of GFS, the estimation of recharge relative to the upstream reach portion and the identification of springs most vulnerable to tunnel drainage impacts. The groundwater budgeting analysis provides evidence that the natural aquifer discharge was stream-focused through GFS, developed downslope or connected to main extensional tectonic lineaments intersecting stream beds; now tunnels drain mainly active recharge groundwater and so cause a relevant stream baseflow deplenishment (approximately two-thirds of the natural value), possibly resulting in adverse effects on local ecosystems.     

Testing optically stimulated luminescence dating on sand-sized quartz of deltaic deposits from the Sperchios delta plain,central Greece

This study reports on the first investigation into the potential of luminescence dating to establish a chronological framework for the depositional sequences of the Sperchios delta plain, central Greece. A series of three borehole cores(20 m deep) and two shallow cores(4 m deep), from across the delta plain, were extracted, and samples were collected for luminescence dating. The luminescence ages of sand-sized quartz grains were obtained from small aliquots of quartz, using the Single-Aliquot Regenerative-dose(SAR) protocol.The equivalent dose determination included a series of tests and the selection of the Minimum Age Model(MAM) as the most appropriate statistical model. This made it possible to confirm the applicability of quartz Optically Stimulated Luminescence(OSL) dating to establish absolute chronology for deltaic sediments from the Sperchios delta plain.Testing age results of the five cores showed that the deltaic sediments were deposited during the Holocene.A relatively rapid deposition is implied for the top ~14 m possibly as a result of the deceleration in the rate of the sea-level rise and the transition to terrestrial conditions, while on the deeper parts, the reduced sedimentation rate may indicate a lagoonal or coastal environment.