spatialCopula: A Matlab toolbox for copula-based spatial analysis

The spatialCopula toolbox contains a set of Matlab functions that provides utilities for copula-based analysis of spatially referenced data, a topic which has re cently attracted much attention in spatial statistics. These tools have been developed to support the work flow in parameter estimation, spatial interpolation and visualization. They offer flexible and user-friendly software for dealing with non-Gaussian and extreme value data that possibly contain a spatial trend or geometric anisotropy. The objective of this paper is to give an introduction to the concept behind the software and to outline the functionality of the toolbox. We illustrate its usefulness by analyzing a data set here referred to as the Gomel data set, which includes moderately skewed radioactivity measurements in the region of Gomel, Belarus. The source codes are freely available in Matlab language on the author’s website (fam.tuwien.ac.at/~hakazian/software.html).     

One step back for a leap forward: toward operational measurements of elements at risk

This study presents the evaluation of 1 year of operational lightning forecasts provided for Europe, using the Weather Research and Forecasting model coupled with a cloud-top height-based lightning parameterization scheme. Three different convective parameterization schemes were employed for parameterizing sub-grid cloud-top heights and consequently driving the lightning scheme. Triggering of the lightning scheme was controlled by means of a model-resolved microphysics-based masking filter, while the formulation for deriving lightning flash rates was also modified, assuming a single “marine” equation instead of the original equations discriminating between continental and marine lightning. Gridded lightning observations were used for evaluating model performance on a dichotomous decision basis. Analysis showed that the lightning scheme is sensitive to the parameterization of convection. In particular, the Kain–Fritsch convective scheme was found to outperform the Grell–Devenyi and Grell–Freitas schemes, showing a statistically significant better performance with respect to lightning prediction. This was most evident during the warm season, while smaller differences among the schemes were recorded during the cold season. Further, for all examined convective schemes, it was found that the application of the masking filter is desirable for improving model performance in terms of lightning forecasting. Last, the reported results revealed that the refinement of the formulation of the lightning parameterization scheme, adhering to a “global” marine equation instead of distinguishing between land and sea lightning, may be necessary in order to obtain reliable lightning forecasts.     

A contribution to the crystal chemistry of the voltaite group: solid solutions, Mössbauer and infrared spectra, and anomalous anisotropy

Voltaite is a mineral of fumaroles, solfatares, coal-fire gas vents, and acid-mine drainage systems. The nominal composition is K₂Fe₅ ²⁺Fe₃ ³⁺Al(SO₄)₁₂·18H₂O and the nominal symmetry is cubic, $Fd\overline{3}c$ . The tetragonal (I4₁/acd) superstructure of voltaite is known as the mineral pertlikite. In this study, we investigated 22 synthetic voltaite samples in which Fe²⁺ was partially or completely replaced by Mg, Zn, Mn, or Cd, by single-crystal and powder X-ray diffraction (both in-house and synchrotron). Two samples contained NH₄ ⁺ instead of K⁺. The structure of voltaite is based on a framework defined by kröhnkite-like heteropolyhedral chains which host both M³⁺ and M²⁺ in octahedral coordination. Unit cell dimensions of the end-members scale almost linearly with the size of M²⁺. In the Fe²⁺-Mg-Zn solid solutions, the Fe²⁺-Mg and Fe²⁺-Zn solutions are linear (ideal) in terms of their lattice-parameter variations. The Mg-Zn solid solution, however, is strongly non-ideal. A detailed analysis of the topology of the chains showed that this behavior originates in expansion and contraction of individual M²⁺-O bonds within the chains. In the Mg-Zn solid solution, some of the M²⁺-O bonds expand while none contract. In the other solid solutions, expansion of some M²⁺-O bonds is always compensated by contraction of the other ones. Parts of the nominally cubic crystals are optically anisotropic and their symmetry is found to be tetragonal by single crystal X-ray diffraction measurements. The coexistence of cubic and tetragonal sectors within a single crystal without any detectable difference in their chemical composition is difficult to explain in terms of growth of such composite crystals. Mössbauer and infrared spectra collected on our synthetic crystals conform with previously published data.     

Does the complexity in temporal precipitation disaggregation matter for a lumped hydrological model?

ABSTRACT

Flood peaks and volumes are essential design variables and can be simulated by precipitation–runoff (P–R) modelling. The high-resolution precipitation time series that are often required for this purpose can be generated by various temporal disaggregation methods. Here, we compare a simple method (M1, one parameter), focusing on the effective precipitation duration for flood simulations, with a multiplicative cascade model (M2, 32/36 parameters). While M2 aims at generating realistic characteristics of precipitation time series, M1 aims only at accurately reproducing flood variables by P–R modelling. Both disaggregation methods were tested on precipitation time series of nine Swiss mesoscale catchments. The generated high-resolution time series served as input for P–R modelling using a lumped HBV model. The results indicate that differences identified in precipitation characteristics of disaggregated time series vanish when introduced into the lumped hydrological model. Moreover, flood peaks were more sensitive than flood volumes to the choice of disaggregation method.     

Seasonal Water Storage Variations as Impacted by Water Abstractions: Comparing the Output of a Global Hydrological Model with GRACE and GPS Observations

Better quantification of continental water storage variations is expected to improve our understanding of water flows, including evapotranspiration, runoff and river discharge as well as human water abstractions. For the first time, total water storage (TWS) on the land area of the globe as computed by the global water model WaterGAP (Water Global Assessment and Prognosis) was compared to both gravity recovery and climate experiment (GRACE) and global positioning system (GPS) observations. The GRACE satellites sense the effect of TWS on the dynamic gravity field of the Earth. GPS reference points are displaced due to crustal deformation caused by time-varying TWS. Unfortunately, the worldwide coverage of the GPS tracking network is irregular, while GRACE provides global coverage albeit with low spatial resolution. Detrended TWS time series were analyzed by determining scaling factors for mean annual amplitude (f _GRACE) and time series of monthly TWS (f _GPS). Both GRACE and GPS indicate that WaterGAP underestimates seasonal variations of TWS on most of the land area of the globe. In addition, seasonal maximum TWS occurs 1 month earlier according to WaterGAP than according to GRACE on most land areas. While WaterGAP TWS is sensitive to the applied climate input data, none of the two data sets result in a clearly better fit to the observations. Due to the low number of GPS sites, GPS observations are less useful for validating global hydrological models than GRACE observations, but they serve to support the validity of GRACE TWS as observational target for hydrological modeling. For unknown reasons, WaterGAP appears to fit better to GPS than to GRACE. Both GPS and GRACE data, however, are rather uncertain due to a number of reasons, in particular in dry regions. It is not possible to benefit from either GPS or GRACE observations to monitor and quantify human water abstractions if only detrended (seasonal) TWS variations are considered. Regarding GRACE, this is mainly caused by the attenuation of the TWS differences between water abstraction variants due to the filtering required for GRACE TWS. Regarding GPS, station density is too low. Only if water abstractions lead to long-term changes in TWS by depletion or restoration of water storage in groundwater or large surface water bodies, GRACE may be used to support the quantification of human water abstractions.     

Rb-Sr provenance ages from weathered and stream transported quartz grains from the Harney Peak Granite,Black Hills,South Dakota

Weathered quartz grus and stream transported quartz of the Harney Peak Granite, Black Hills, South Dakota, contain low concentrations of Rb (generally 0.3–6.8 ppm) and Sr (0.2–2.0 ppm) and variable Sr isotopic ratios (0.759–1.070).Six of seven single grains of large composite quartz grus which recently entered the weathering environment define an apparent isochron age (about 1800 Myr) and initial ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio (0.7066) that approximate the whole-rock isochron age (1707 Myr) and initial ratio (0.7143) of the Harney Peak Granite. Apparently the Rb-Sr systematics of these grains were not significantly altered during initial weathering. Leached fluid inclusion material from a ca. 2 g aggregate of composite quartz grains contains very little Rb or Sr (0.019 and 0.17 μg, respectively) and has a very low ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio (0.739). The Rb and Sr content of the quartz grains appears to be concentrated in minute, heterogeneously-distributed mineral inclusions.Five aggregates of more completely weathered, small non-composite quartz grains produce a widely scattered pattern on an isochron diagram with all samples plotting below the 1707 Myr isochron. Examination by SEM of these grains shows solution and precipitation features on their relatively large effective surface areas. The differential precipitation of Rb is believed to have been the major perturbating chemical process during weathering.Three aggregates of stream quartz grains define an apparent isochron age of 1777 Myr and an initial ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio of 0.720 that suggest the initial ‘igneous’ Rb-Sr characteristics of the stream quartz were re-attained during their transportation, probably as a result of removal of the outer weathered surface by abrasion. The apparent resistance to chemical weathering of stream quartz and quartz which has just entered the weathering environment suggests that this mineral may be extremely useful in studies of provenance and the geochronology of strongly weathered terranes.     

New spatial dimensions of global cityscapes: From reviewing existing concepts to a conceptual spatial approach

Current global urbanisation processes are leading to new forms of massive urban constellations. The conceptualisations and classifications of these, however, are often ambiguous, overlap or lag behind in scientific literature. This article examines whether there is a common denominator to define and delimitate–and ultimately map–these new dimensions of cityscapes. In an extensive literature review we analysed and juxtaposed some of the most common concepts such as megacity, megaregion or megalopolis. We observed that many concepts are abstract or unspecific, and for those concepts for which physical parameters exist, the parameters are neither properly defined nor used in standardised ways. While understandably concepts originate from various disciplines, the authors identify a need for more precise definition and use of parameters. We conclude that often, spatial patterns of large urban areas resemble each other considerably but the definitions vary so widely that these differences may surpass any inconsistencies in the spatial delimitation process. In other words, today we have tools such as earth observation data and Geographic Information Systems to parameterise if clear definitions are provided. This appears not to be the case. The limiting factor when delineating large urban areas seems to be a commonly agreed ontology.     

The German Bight (North Sea) is a nursery area for both locally and externally produced sprat juveniles

To better understand the role of the German Bight (GB) as a nursery area for juvenile North Sea sprat Sprattus sprattus we sought to determine whether the area may receive only locally or also externally produced offspring. We sampled juveniles during 3 trawl surveys in the GB in August, September, and October 2004 and applied otolith microstructure analysis in order to reconstruct their distributions of the day-of-first-increment-formation (dif). These were contrasted with spatial and seasonal patterns of sprat egg abundance in the GB and its adjacent areas, observed during 6 monthly plankton surveys. It was found that the majority of juveniles originated mainly from April/May 2004, coinciding with high spawning activity west of the GB, whereas spawning and larval production inside the GB peaked notably later, in May/June. This indicated that a large proportion of juveniles was produced outside the GB and transported subsequently into it through passive and/or active migration. Shifts to later mean difs from one survey to the next and length distributions indicative of the simultaneous presence of multiple cohorts, supported the notion that the GB is a complex retention and nursery area for sprat offspring from different North Sea spawning grounds and times. Later born juveniles had significantly faster initial growth rates than earlier born conspecifics, which was likely temperature-mediated, given the strong correlation between back-calculated growth histories and sea surface temperature as a proxy for thermal histories of juveniles (r² = 0.52).     

Scales of columnar jointing in igneous rocks: field measurements and controlling factors

Columnar jointing is a common feature of solidified lavas, sills and dikes, but the factors controlling the characteristic stoutness of columns remain debated, and quantitative field observations are few in number. In this paper, we provide quantitative measurements on sizing of columnar joint sets and our assessment of the principal factors controlling it. We focus on (1) chemistry, as it is the major determinant of the physical (mechanical and thermal) properties of the lava, and (2) geology, as it influences the style of emplacement and lava geometry, setting boundary conditions for the cooling process and the rate of heat loss. In our analysis, we cover lavas with a broad range of chemical compositions (from basanite to phonolite, for six of which we provide new geochemical analyses) and of geological settings. Our field measurements cover 50 columnar jointing sites in three countries. We provide reliable, manually digitized data on the size of individual columns and focus the mathematical analysis on their geometry (23,889 data on side length, of which 17,312 are from full column sections and 3,033 data on cross-sectional area and order of polygonality). The geometrical observations show that the variation in characteristic size of columns between different sites exceeds one order of magnitude (side length ranging from 8 to 338 cm) and that the column-bounding polygons’ average order is less than 6. The network of fractures is found to be longer than required by a minimum-energy hexagonal configuration, indicating a non-equilibrium, geologically quick process. In terms of the development and characteristic sizing of columnar joint sets, our observations suggest that columns are the result of an interplay between the geological setting of emplacement and magma chemistry. When the geological setting constrains the geometry of the emplaced body, it exerts a stronger control on characteristic column stoutness. At unconstrained geometries (e.g. unconfined lava flows), chemistry plays the major role, resulting in stouter columns in felsic lavas and slenderer columns in mafic lavas.