fsr: An R package for fuzzy spatial data handling

Spatial libraries are core components in many geographic information systems, spatial database systems, and spatial data science projects. These libraries provide the implementation of spatial type systems that include spatial data types and a large diversity of geometric operations. Their focus relies on handling crisp spatial objects, which are characterized by an exact location and a precisely defined extent, shape, and boundary in space. However, there is an increasing interest in analyzing spatial phenomena characterized by fuzzy spatial objects, which have inexact locations, vague boundaries, and/or blurred interiors. Unfortunately, available spatial libraries do not provide support for fuzzy spatial objects. In this article, we describe the R package named fsr, which is based on the Spatial Plateau Algebra and is publicly available at https://cran.r-project.org/package=fsr . Our tool provides methods for building fuzzy spatial objects as spatial plateau objects and conducting exploratory spatial data analysis by using fuzzy spatial operations.     

Barnacles as biomonitors of metal contamination in coastal waters

The use of barnacles as biomonitors of metal contamination in coastal waters worldwide is reviewed as a critique compilation of the reported studies and presents resume-tables of available data for future reference. The barnacle body reflects both short and long-term metal level environmental variations and the metal bioaccumulation occurs mainly in their granules (relatively inactive pools). The barnacle body is considered as good biomonitoring material and different barnacle species could bioaccumulate metal concentration ranges of 40–153,000 μg/g of Zn, 20–22,230 μg/g de Fe, 1.5–21,800 μg/g of Cu, 5.9–4742 μg/g of Mn, 0.1–1000 μg/g of Pb, 0.7–330 μg/g of Cd, 0.4–99 μg/g of Ni and 0.2–49 μg/g of Cr. However, as the plates (‘shells’) of barnacle exoskeletons can be affected by metal levels in coastal waters, mainly in their composition and morphology, they are not considered good biomonitoring material. Despite this, the use of a specific barnacle species or group of species in a specific region must firstly be carefully validated and the interpretation of the contaminant bioaccumulation levels should involve specific environmental variations of the region, physiological parameters of the barnacle species and the relationship between the potential toxicity of the contaminant for the environment and their significance for the barnacle species. Barnacles, particularly a widespread cosmopolitan species such as Amphibalanus amphitrite, have a great potential as biomonitors of anthropogenic contamination in coastal waters and have been used worldwide, including Europe (United Kingdom, Turkey, Poland, Croatia, Spain and Portugal), Asia (India and China), Oceania (Australia), North America (Florida, Massachusetts and Mexico) and South America (Brazil). The use of barnacle species as biomonitors of metal contamination in coastal waters is considered an important and valuable tool to evaluate and predict the ecological quality of an ecosystem.     

Improvement of the JRC Calculation Using Different Parameters Obtained Through a New Survey Method Applied to Rock Discontinuities

This paper presents a new, fast and economical method for digitizing a rock surface profile to determine its joint roughness coefficient (JRC) value (Barton and Choubey, Rock Mech Rock Eng 10:1–54, 1977) using a set of existing parameters, the appropriateness of which has been tested. The new procedure was applied to four rock road cuts in the Alpujarra area of Granada province (Andalusia, Spain), supported by 117 digitized profiles. The method avoids any personal subjectivity bias. The usual subjectivity in the visualization and assessment of the JRC value is estimated here by a survey test answered by 90 trained users. The survey test consisted of a set of JRC visual comparisons of 12 roughness profiles, selected from the 117 profiles surveyed during the field research in the Alpujarra region, with the Barton typical profiles. The results of the survey clearly show systematic inaccuracies in the traditional procedure that are reduced if the advice presented in the conclusions of this paper is taken in account when performing the visual estimation.     

The DRASTIC-Sg model: an extension to the DRASTIC approach for mapping groundwater vulnerability in aquifers subject to differential land subsidence,with application to Mexico City

Mexico City relies on groundwater for most of its domestic supply. Over the years, intensive pumping has led to significant drawdowns and, subsequently, to severe land subsidence. Tensile cracks have also developed or reactivated as a result. All such processes cause damage to urban infrastructure, increasing the risk of spills and favoring contaminant propagation into the aquifer. The effects of ground deformation are frequently ignored in groundwater vulnerability studies, but can be relevant in subsiding cities. This report presents an extension to the DRASTIC methodology, named DRASTIC-Sg, which focuses on evaluating groundwater vulnerability in urban aquifers affected by differential subsidence. A subsidence parameter is developed to represent the ground deformation gradient (Sg), and then used to depict areas where damage risk to urban infrastructure is higher due to fracture propagation. Space-geodetic SqueeSAR data and global positioning system (GPS) validation were used to evaluate subsidence rates and gradients, integrating hydrogeological and geomechanical variables into a GIS environment. Results show that classic DRASTIC approaches may underestimate groundwater vulnerability in settings such as the one at hand. Hence, it is concluded that the Sg parameter is a welcome contribution to develop reliable vulnerability assessments in subsiding basins.     

优视摄影测量方法及精度分析EI北大核心CSCD

无人机最优路径规划和精准数据采集是复杂城市场景三维重建的关键问题。优视摄影测量利用被摄对象的概略模型,结合密集采样的初始视点生成和采样点可重建性约束的视点优化技术,实现最小数据采集代价下的精确三维模型重建。本文首先研究了优视摄影测量的基本原理;然后利用优视摄影测量技术进行真实场景的无人机影像采集,并结合激光扫描点云进行空三精度分析和Mesh模型质量验证。试验结果表明,优视摄影测量显著提高了Mesh模型的重建质量;对于倾斜摄影难以观测的部分建筑物立面区域,精度提升达到3~5倍;与倾斜摄影相比,优视摄影测量能建立稳健影像连接,实现较高的绝对定位精度。     

GPS-measurements of recent crustal deformation in the junction zone of the rift segments in the central Baikal rift system

Based on multiyear measurements of present-day motions in the central area of the Baikal rift system, new data on the kinematics of horizontal motions, relative horizontal deformation rates, and rotation velocities in the area of junction of the South Baikal, North Baikal, and Barguzin rift basins have been obtained. This area is an intricate structure with two transfer zones: Ol’khon–Svyatoi Nos and Ust’-Barguzin.It is shown that crustal blocks are moving southeastward, normally to the structures of transfer zones and at an acute angle to the Baikal Rift strike, which corresponds to the right-lateral strike-slip extensional faulting along the major structure. The average horizontal velocities increase from 3.0 mm yr–1 in the northern South Baikal basin to 6.5 mm yr–1 in the Barguzin basin. The elongation axes prevailing in the study region are mainly of NW–SE direction. The areas of intense deformations are confined to structures with high seismic activity in the South Baikal and, partly, Barguzin basins. This confirms the existence of a present-day zone of the Earth’s crust destruction in the Baikal rift system, which is the most likely source of strong earthquakes in the future. Two zones with rotations in opposite directions are recognized in the rotation velocity field. Clockwise rotation is typical of structures of N–NE strike (Maloe More basin, southern North Baikal basin, Barguzin Ridge rise). Counterclockwise rotation is determined for NE-striking structures (northern South Baikal basin, southern Barguzin basin). In general, the obtained data show an intricate pattern of present-day horizontal dislocations and deformations in the area of junction of NE- and N–NE-striking rift structures. This suggests left- and right-lateral strike-slip faults, respectively, within them.     

Geostatistical Data Integration Model for Contamination Assessment

Soil contamination assessments can be improved with new methods aimed at the accurate estimation of the volume and extension of contaminated soil to be remediated. Geostatistical models that use secondary information to characterize soil contamination are incorporated into a new integration model to provide accurate three-dimensional maps. The proposed integration model is based on a stochastic inversion approach and uses sequential indicator simulation. A two-dimensional reference image representing the areal extension of the contamination is combined with local measurements of contamination in the vertical direction, to render a three-dimensional contamination map. To demonstrate how well the integration model performs, the case study presented focuses on geophysical data and how it can be integrated with soil contamination measurements to improve the characterization of a contaminated site. The results show that the model reproduces successfully the reference image thus providing an accurate three-dimensional contamination map.     

Impacts of wildfire and post-fire land management on hydrological and sediment processes in a humid Mediterranean headwater catchment

The extensive afforestation of the Mediterranean rim of Europe in recent decades has increased the number of wildfire disturbances on hydrological and sediment processes, but the impacts on headwater catchments is still poorly understood, especially when compared with the previous agricultural landscape. This work monitored an agroforestry catchment in the north-western Iberian Peninsula, with plantation forests mixed with traditional agriculture using soil conservation practices, for one year before the fire and for three years afterwards, during which period the burnt area was ploughed and reforested. During this period, continuous data was collected for meteorology, streamflow and sediment concentration at the outlet, erosion features were mapped and measured after major rainfall events, and channel sediment dynamics were monitored downstream from the agricultural and the burnt forest area. Data from 202 rainfall events with over 10 mm was analysed in detail. Results show that the fire led to a notable impact on sediment processes during the first two post-fire years, but not on streamflow processes; this despite the small size of the burnt area (10% of the catchment) and the occurrence of a severe drought in the first year after the fire. During this period, soil loss at the burnt forest slopes was much larger than that at most traditionally managed fields, and, ultimately, led to sediment exhaustion. At the catchment scale, storm characteristics were the dominant factor behind streamflow and sediment yield both before and after the fire. However, the data indicated a shift from detachment-limited sediment yield before the fire, to transport-limited sediment yield afterwards, with important increases in streamflow sediment concentration. This indicates that even small fires can temporarily change sediment processes in agroforestry catchments, with potential negative consequences for downstream water quality.     

Modelling the first probable two-plane lens system B2114+022: reproducing two compact radio cores A and D

We test possible lensing scenarios of the JVAS system B2114+022, in which two galaxies at different redshifts ('G1' at z ₁=0.3157 and 'G2' at z ₂=0.5883) are found within 2 arcsec of quadruple radio sources. For our investigation, we use possible lensing constraints derived from a wealth of data on the radio sources obtained with VLA, MERLIN, VLBA and EVN as well as HST imaging data on the two galaxies, which were presented recently in Augusto et al. In the present study, we focus on reproducing the widest separated, observationally similar radio components A and D as lensed images. We first treat G2 (which is the more distant one from the geometric centre) as a shear term, and then consider two-plane lensing explicitly including G2's potential at the z ₂ plane as the first case of two-plane lens modelling. Our modelling results not only support the hypothesis that the system includes gravitationally lensed images of a higher-redshift extragalactic object, but they also show that the explicit inclusion of G2's potential at the second lens plane is necessary in order to fit the data with astrophysically plausible galaxy parameters. Finally, we illustrate a natural consequence of a two-plane lens system, namely the prediction of distortion as well as shift and stretching of G2's isophotes by G1's potential, which can in principle be measured by subtracting out G1's light distribution in an image of high signal-to-noise ratio and good angular resolution, especially a multicolour one.