Smoothing and high risk areas detection in space-time disease mapping: a comparison of P-splines,autoregressive, and moving average models

Recently, several models have been proposed for smoothing risks in disease mapping. These models consider different ways of introducing both spatial and temporal dependence as well as spatio-temporal interactions. In this work, a comparison among some autoregressive, moving average, and P-spline models is performed. Firstly, brain cancer mortality data are used to analyze the degree of smoothness introduced by these models. Secondly, two separate simulation studies (one model-based and the other model-free) are carried out to evaluate the model performance in terms of bias, variability, sensitivity, and specificity. We conclude that P-spline models seem to be a good alternative to autoregressive and moving average models when analyzing highly sparse disease mapping data.     

Seasonal copula models for the analysis of glacier discharge at King George Island,Antarctica

Modelling glacier discharge is an important issue in hydrology and climate research. Glaciers represent a fundamental water resource when melting of ice and snow contributes to runoff. Glaciers are also studied as natural global warming sensors. GLACKMA association has implemented one of their Pilot Experimental Catchment areas at the King George Island in the Antarctica which records values of the liquid discharge from Collins glacier. In this paper, we propose the use of time-varying copula models for analyzing the relationship between air temperature and glacier discharge, which is clearly non constant and non linear through time. A seasonal copula model is defined where both the marginal and copula parameters vary periodically along time following a seasonal dynamic. Full Bayesian inference is performed such that the marginal and copula parameters are estimated in a one single step, in contrast with the usual two-step approach. Bayesian prediction and model selection is also carried out for the proposed model such that Bayesian credible intervals can be obtained for the conditional glacier discharge given a value of the temperature at any given time point. The proposed methodology is illustrated using the GLACKMA real data where there is, in addition, a hydrological year of missing discharge data which were not possible to measure accurately due to problems in the sounding.     

A stochastic methodology for risk assessment of a large earthquake when a long time has elapsed

We propose a stochastic methodology for risk assessment of a large earthquake when a long time has elapsed from the last large seismic event. We state an approximate probability distribution for the occurrence time of the next large earthquake, by knowing that the last large seismic event occurred a long time ago. We prove that, under reasonable conditions, such a distribution is exponential with a rate depending on the asymptotic slope of the cumulative intensity function corresponding to a nonhomogeneous Poisson process. As it is not possible to obtain an empirical cumulative distribution function of the waiting time for the next large earthquake, an estimator of its cumulative distribution function based on existing data is derived. We conduct a simulation study for detecting scenario in which the proposed methodology would perform well. Finally, a real-world data analysis is carried out to illustrate its potential applications, including a homogeneity test for the times between earthquakes.     

Frontal dynamics boost primary production in the summer stratified Mediterranean sea

Bio-physical glider measurements from a unique process-oriented experiment in the Eastern Alboran Sea (AlborEx) allowed us to observe the distribution of the deep chlorophyll maximum (DCM) across an intense density front, with a resolution (～ 400 m) suitable for investigating sub-mesoscale dynamics. This front, at the interface between Atlantic and Mediterranean waters, had a sharp density gradient (Δρ ～ 1 kg/m³ in ～ 10 km) and showed imprints of (sub-)mesoscale phenomena on tracer distributions. Specifically, the chlorophyll-a concentration within the DCM showed a disrupted pattern along isopycnal surfaces, with patches bearing a relationship to the stratification (buoyancy frequency) at depths between 30 and 60 m. In order to estimate the primary production (PP) rate within the chlorophyll patches observed at the sub-surface, we applied the Morel and Andrè (J Geophys Res 96:685–698 1991) bio-optical model using the photosynthetic active radiation (PAR) from Argo profiles collected simultaneously with glider data. The highest production was located concurrently with domed isopycnals on the fresh side of the front, suggestive that (sub-)mesoscale upwelling is carrying phytoplankton patches from less to more illuminated levels, with a contemporaneous delivering of nutrients. Integrated estimations of PP (1.3 g C m^?2d^?1) along the glider path are two to four times larger than the estimations obtained from satellite-based algorithms, i.e., derived from the 8-day composite fields extracted over the glider trip path. Despite the differences in spatial and temporal sampling between instruments, the differences in PP estimations are mainly due to the inability of the satellite to measure DCM patches responsible for the high production. The deepest (depth > 60 m) chlorophyll patches are almost unproductive and probably transported passively (subducted) from upper productive layers. Finally, the relationship between primary production and oxygen is also investigated. The logarithm of the primary production in the DCM interior (chlorophyll (Chl) > 0.5 mg/m³) shows a linear negative relationship with the apparent oxygen utilization, confirming that high chlorophyll patches are productive. The slope of this relationship is different for Atlantic, mixed interface waters and Mediterranean waters, suggesting the presence of differences in planktonic communities (whether physiological, population, or community level should be object of further investigation) on the different sides of the front. In addition, the ratio of optical backscatter to Chl is high within the intermediate (mixed) waters, which is suggestive of large phytoplankton cells, and lower within the core of the Atlantic and Mediterranean waters. These observations highlight the relevance of fronts in triggering primary production at DCM level and shaping the characteristic patchiness of the pelagic domain. This gains further relevance considering the inadequacy of optical satellite sensors to observe DCM concentrations at such fine scales.     

Trend analysis of nutrient loadings in a large prairie catchment

Nutrient loadings in many river catchments continue to increase due to rapid expansion of agriculture, urban and industrial development, and population growth. Nutrient enrichment of water bodies has intensified eutrophication which degrades water quality and ecosystem health. In this study, we carried out a trend analysis of total phosphorus and total nitrogen loads in the South Saskatchewan River (SSR) catchment using a novel approach to analyse nutrient time series. Seasonal analysis of trends at each of the water quality stations was performed to determine the relationships between annual flow regimes and nutrient loads in the catchment, in particular, the influence of the high spring runoff on nutrient export. Decadal analysis was also performed to determine the long-term relationships of nutrients with anthropogenic changes in the catchment. Although it was found that seasonal and historical variability of nutrient load trends is mainly determined by streamflow regime changes, there is evidence that increases in nitrogen concentration can also be attributed to anthropogenic changes.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Atlas of susceptibility to pollution in marinas. Application to the Spanish coast

An atlas of susceptibility to pollution of 320 Spanish marinas is provided. Susceptibility is assessed through a simple, fast and low cost empirical method estimating the flushing capacity of marinas. The Complexity Tidal Range Index (CTRI) was selected among eleven empirical methods. The CTRI method was selected by means of statistical analyses because: it contributes to explain the system's variance; it is highly correlated to numerical model results; and, it is sensitive to marinas' location and typology. The process of implementation to the Spanish coast confirmed its usefulness, versatility and adaptability as a tool for the environmental management of marinas worldwide. The atlas of susceptibility, assessed through CTRI values, is an appropriate instrument to prioritize environmental and planning strategies at a regional scale.     

Geophysical evidence of the Eastern Marginal Fault of the Cheb Basin (Czech Republic)

The western part of the Bohemian Massif hosts an intersection of two regional fault zones, the SW-NE trending Oh?e/Eger Graben and the NNW-SSE trending Mariánské Lázně Fault, which has been reactivated several times in the geological history and controlled the formation of the Tertiary Cheb Basin. The broader area of the Cheb Basin is also related to permanent seismic activity of M_L 3+ earthquake swarms. The Eastern Marginal Fault of the Cheb Basin (northern segment of the Mariánské Lázně Fault) separates the basin sediments and underlying granites in the SW from the Kru?né Hory/Erzgebirge Mts. crystalline unit in the NE. We describe a detailed geophysical survey targeted to locating the Eastern Marginal Fault and determining its geometry in the depth. The survey was conducted at the Kopanina site near the Nový Kostel focal zone, which shows the strongest seismic activity of the whole Western Bohemia earthquake swarm region. Complex geophysical survey included gravimetry, electrical resistivity tomography, audiomagnetotellurics and seismic refraction. We found that the rocks within the Eastern Marginal Fault show low resistivity, low seismic velocity and density, which indicates their deep fracturing, weathering and higher water content. The dip of the fault in shallow depths is about 60° towards SW. At greater depths, the slope turns to subvertical with dip angle of about 80°. Results of geoelectrical methods show blocky fabric of the Cheb Basin and deep weathering of the granite bedrock, which is consistent with geologic models based on borehole surveys.     

Deterioration processes affecting prehistoric rock art engravings in granite in NW Spain

Prehistoric rock art sites are endangered despite conservation efforts. The lack of scientific documentation regarding weathering agents affecting rock art and the absence of specific diagnostic protocols hinder the development of conservation strategies. The aim of this research was to investigate active deterioration processes in a granite petroglyph site located in Mougás (Galicia, NW Spain) by characterizing the granite, conducting a geotechnical study of the outcrop and describing and analysing the main weathering processes. Two main deterioration factors were identified. First, water favours block disjunction at the massif scale and causes pitting and surface erosion at the millimetre scale that affects the readability of the engravings. Second, high temperatures associated with wildfires cause mineral transformations that increase the susceptibility of the rock to weathering. Identifying deterioration factors is a first step in developing appropriate preventive conservation measures, which should aim to reduce rock contact time with water (technically affordable in the short term) and to reduce the probability of wildfire occurrence (technically more complex and possibly with longer‐term results). Copyright © 2018 John Wiley & Sons, Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.