Scalar-Flux Similarity in the Layer Near the Surface Over Mountainous Terrain

The scaled standard deviations of temperature and humidity are investigated in complex terrain. The study area is a steep Alpine valley, with six measurement sites of different slope, orientation and roughness (i-Box experimental site, Inn Valley, Austria). Examined here are several assumptions forming the basis of Monin–Obukhov similarity theory (MOST), including constant turbulence fluxes with height and the degree of self-correlation between the involved turbulence variables. Since the basic assumptions for the applicability of the MOST approach—horizontally homogeneous and flat conditions—are violated, the analysis is performed based on a local similarity hypothesis. The scaled standard deviations as a function of local stability are compared with previous studies from horizontally homogeneous and flat terrain, horizontally inhomogeneous and flat terrain, weakly inhomogeneous and flat terrain, as well as complex terrain. As a reference, similarity relations for unstable and stable conditions are evaluated using turbulence data from the weakly inhomogeneous and flat terrain of the Cabauw experimental site in the Netherlands, and assessed with the same post-processing method as the i-Box data. Significant differences from the reference curve and also among the i-Box sites are noted, especially for data derived from the i-Box sites with steep slopes. These differences concern the slope and the magnitude of the best-fit curves, illustrating the site dependence of any similarity theory.     

Comment on Pisarenko et al., “Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory”

In this short note, I comment on the research of Pisarenko et al. (Pure Appl. Geophys 171:1599–1624, 2014) regarding the extreme value theory and statistics in the case of earthquake magnitudes. The link between the generalized extreme value distribution (GEVD) as an asymptotic model for the block maxima of a random variable and the generalized Pareto distribution (GPD) as a model for the peaks over threshold (POT) of the same random variable is presented more clearly. Inappropriately, Pisarenkoet al. (Pure Appl. Geophys 171:1599–1624, 2014) have neglected to note that the approximations by GEVD and GPD work only asymptotically in most cases. This is particularly the case with truncated exponential distribution (TED), a popular distribution model for earthquake magnitudes. I explain why the classical models and methods of the extreme value theory and statistics do not work well for truncated exponential distributions. Consequently, these classical methods should be used for the estimation of the upper bound magnitude and corresponding parameters. Furthermore, I comment on various issues of statistical inference in Pisarenkoet al. and propose alternatives. I argue why GPD and GEVD would work for various types of stochastic earthquake processes in time, and not only for the homogeneous (stationary) Poisson process as assumed by Pisarenko et al. (Pure Appl. Geophys 171:1599–1624, 2014). The crucial point of earthquake magnitudes is the poor convergence of their tail distribution to the GPD, and not the earthquake process over time.     

Can GPS-Derived Surface Loading Bridge a GRACE Mission Gap?

We investigated two ‘gap-filler’ methods based on GPS-derived low-degree surface loading variations (GPS-I and GPS-C) and a more simple method (REF-S) which extends a seasonal harmonic variation into the expected Gravity Recovery and Climate Experiment (GRACE) mission gap. We simulated two mission gaps in a reference solution (REF), which is derived from a joint inversion of GRACE (RL05) data, GPS-derived surface loading and simulated ocean bottom pressure. The GPS-I and GPS-C methods both have a new type of constraint applied to mitigate the lack of GPS station network coverage over the ocean. To obtain the GPS-C solution, the GPS-I method is adjusted such that it fits the reference solution better in a 1.5 year overlapping period outside of the gap. As can be expected, the GPS-I and GPS-C solutions contain larger errors compared to the reference solution, which is heavily constrained by GRACE. Within the simulated gaps, the GPS-C solution generally fits the reference solution better compared to the GPS-I method, both in terms of spherical harmonic loading coefficients and in terms of selected basin-averaged hydrological mass variations. Depending on the basin, the RMS-error of the water storage variations (scaled for leakage effects) ranges between 1.6 cm (Yukon) and 15.3 cm (Orinoco). In terms of noise level, the seasonal gap-filler method (REF-S) even outperforms the GPS-I and GPS-C methods, which are still affected by spatial aliasing problems. However, it must be noted that the REF-S method cannot be used beyond the study of simple harmonic seasonal variations.     

Transient Divergent Flow and Transport in an Infinite Anisotropic Porous Formation

When seeking to predict plume geometry resulting from fluid injection through partially penetrating wells, it is common to assume a steady-state spherically diverging flow field. In reality, the flow field is transient. The steady-flow assumption is likely to cause overestimation of injection plume radius since the accommodation of fluid by increases in porosity and fluid density is ignored. In this paper, a transient solution is developed, resulting in a nonlinear ordinary differential equation expressing plume radius as a function of time. It is shown that the problem can be fully described by one type curve. A critical time, t_c, is identified at which the percentage error of the steady-state flow solution compared to the fully dynamic problem is less than 1%. Only for large injection rates and low permeabilities, does t_c become greater than 1 h. Nevertheless, an improved approximate solution is obtained by a simple linearization procedure. The critical time, t_c for the new approximate solution is 0.3% of that required for the steady-state flow solution.     

Monitoring the regional and temporal variability of winter snowfall in the arid Andes using digital NOAA/AVHRR data

Abstract

This paper deals with the spatial distribution and the temporal variability of snowfall in the most arid part of the Andes (18°‐ 28°S) during southern hemisphere winter (May‐September). As the official precipitation data is of poor quality, analyses were carried out by means of digital image processing techniques, using NOAA/AVHRR satellite‐data. Through analysis of 24 different snowfall events from six winters, a previously unknown spatial and temporal precipitation pattern in this remote and unexplored area was revealed. Snowfall is most abundant in the southernmost part of the research area and on the western side of the Andes, indicating the Pacific origin of the snowfall.

Nevertheless, the typical snowfall pattern is modified during different periods of the winter. Three typical time periods could be defined and distinguished from one another. Each of these three periods is characterized by typical weather conditions (cold fronts and “cut‐offs “) leading to a distinct snowfall pattern.

As this study is part of a broader paleoclimatic project, the results will serve as a basis for paleoclimatic reconstruction of past climate. Only by knowing the modern circulation and precipitation patterns is it possible to interpret paleoclimatic signals and archives found in the study area (e.g. paleosol, moraines) correctly.     

Empirical behaviour of tests for the beta distribution and their application in environmental research

The beta distribution is used in different models of environmental research. The power of the test for beta distribution of Raschke [Biased transformation and its application in goodness-of-fit tests for the beta and gamma distribution. Commun. Statist. B–Computa. Simula. 38 (2009): 1870–1890] is researched here. The power of the Kolmogorov–Smirnov, Kuiper, Cramér-von Mises, Watson and Anderson–Darling tests are researched for different sample sizes, levels of significance and parameters of the beta distribution. The limitation to these tests is discussed including the differences between previous publications. The empirical behaviour is investigated by extensive Monte Carlo simulations. The most powerful test for the beta distribution is the Anderson–Darling test for the considered constellations of alternative distribution, contamination or scaling. The second best test is the Cramér-von Mises test, followed by the Watson test. The analysis of relative humidity data of meteorology and of runoff coefficients of the hydrology demonstrates the advantages of the new tests and the necessity to test an assumption of beta distribution.     

On the breakdown of zircon upon “dry” thermal annealing

Zircon samples without and with secondary chemical alteration from diverse sources were subjected to heat treatment at 1400 °C for 96 h. Resulting new phases and textures suggest that decomposition of zircon into component oxides occurred in all experiments to various degrees. The crucible material was found to have a strong influence on the extent of breakdown, especially in the case of altered starting materials. In this study the progressive stages of the breakdown of zircon grains are described. The factors that may govern the decomposition are discussed, including radiation damage, secondary alteration and external reaction conditions (sample container, atmosphere). Alumina crucibles should generally be avoided in dry annealing of zircon, to minimise uncontrolled breakdown into oxides.     

S?o Francisco-Congo Craton break-up delimited by U-Pb-Hf isotopes and trace-elements of zircon from metasediments of the Ara?uaí Belt

Detrital zircon U-Pb geochronology combined with Hf isotopic and trace element data from metasedimentary rocks of the Aracuai Belt in southeastern Brazil provide evidence for break-up of the Congo-Sao Francisco Craton. The U-Pb age spectra of detrital zircons from metasediments of the Rio Doce Group(RDG) range from 900-650 Ma and define a maximum depositional age of ca. 650 Ma. Zircon trace element and whole rock data constrain an oceanic island arc as source for the deposition setting of the protoliths to the metasediments. Zircon ε_(Hf)(t) values from these rocks are positive between +1 and +15, supporting previous evidence of a Neoproterozoic extensional phase and oceanic crust formation in a precursor basin to the Aracuai Belt. Recrystallization of detrital zircon at ca. 630 Ma is compatible with a regional metamorphic event associated with terrane accretion to the Paleoproterozoic basement after transition from an extensional to a convergent regime. The juvenile nature, age spectra and trace element composition recorded in detrital zircons of metasediments from the Aracuai Belt correspond with zircons from metasedimentary rocks and oceanic crust remnants of other orogenic belts to its south. This suggests that rifting and oceanic crust formation of the entire orogenic system, the so-called Mantiqueira Province, was contemporaneous, most likely related to the opening of a large ocean. It further indicates that the cratonic blocks involved in the orogenic evolution of the Mantiqueira Province were spatially connected as early as 900 Ma.