Eigen analysis of tree-ring records: part 2, posing the eigen problem

The technique of expanding meteorological fields on eigenvectors of the field covariation matrix is popular. In this paper, we propose for the first time to use a mathematically similar technique to solve the main problem of dendrochronology: classifying variations in tree-ring records as either age- and microenvironment-dependent or climate-induced. Applying this technique to a sample of very long-lived Qilian junipers (Sabina przewalskii Kom.) from the Dulan region in western China, we demonstrate that the ring-width variations projected on the first eigenvector are age-dependent, but those projected on several of the first subsequent vectors are mainly climate-induced. In particular, the second and third projections capture multi-centennial climatic variations, and the variations projected on the fourth through seventh eigenvectors show periodic variations that are probably induced by the 178-year solar cycle. The projections on the smallest eigenvectors seem to be negligible.     

Ecohydrological modelling in a deciduous boreal forest: Model evaluation for application in non-stationary climates

Soil moisture is an important driver of growth in boreal Alaska, but estimating soil hydraulic parameters can be challenging in this data-sparse region. Parameter estimation is further complicated in regions with rapidly warming climate, where there is a need to minimize model error dependence on interannual climate variations. To better identify soil hydraulic parameters and quantify energy and water balance and soil moisture dynamics, we applied the physically based, one-dimensional ecohydrological Simultaneous Heat and Water (SHAW) model, loosely coupled with the Geophysical Institute of Permafrost Laboratory (GIPL) model, to an upland deciduous forest stand in interior Alaska over a 13-year period. Using a Generalized Likelihood Uncertainty Estimation parameterisation, SHAW reproduced interannual and vertical spatial variability of soil moisture during a five-year validation period quite well, with root mean squared error (RMSE) of volumetric water content at 0.5 m as low as 0.020 cm³/cm³. Many parameter sets reproduced reasonable soil moisture dynamics, suggesting considerable equifinality. Model performance generally declined in the eight-year validation period, indicating some overfitting and demonstrating the importance of interannual variability in model evaluation. We compared the performance of parameter sets selected based on traditional performance measures such as the RMSE that minimize error in soil moisture simulation, with one that is designed to minimize the dependence of model error on interannual climate variability using a new diagnostic approach we call CSMP, which stands for Climate Sensitivity of Model Performance. Use of the CSMP approach moderately decreases traditional model performance but may be more suitable for climate change applications, for which it is important that model error is independent from climate variability. These findings illustrate (1) that the SHAW model, coupled with GIPL, can adequately simulate soil moisture dynamics in this boreal deciduous region, (2) the importance of interannual variability in model parameterisation, and (3) a novel objective function for parameter selection to improve applicability in non-stationary climates.     

Temporal variability in SO<Subscript>2</Subscript> exposures at Hawai’i Volcanoes National Park,USA

Sulfur dioxide measured at the Visitor’s Center of Hawai’i Volcanoes National Park, USA exhibits seasonal and diurnal patterns and is sensitive to rainfall, but the relationship between wind direction and SO₂ is too inconsistent to support deterministic predictions of hourly SO₂ based on hourly wind direction. Although SO₂ at the Visitor Center has usually been below regulatory levels, high SO₂ levels and adverse health effects remain a concern. This investigation identified patterns in hourly SO₂ based on wind direction, time of day, month, and rainfall occurrence using 4 years of hourly data. Empirical probabilities were investigated using a Bayesian approach. Winds from the volcanic vents were rare, but when they did occur SO₂ was elevated about half of the time. Conversely, half of the hours with elevated SO₂ occurred when volcanic vents were not directly upwind. Episodes of elevated SO₂ tended to occur during the months of November–March and between 8:00 a.m. and 5:00 p.m. Rainfall was associated with a marked reduction in SO₂ (29–81% depending on wind direction). Individuals that wish to avoid exposure to volcanic fumes can take these patterns into account.     

Numerical modeling of 3D fully nonlinear potential periodic waves

A simple and exact numerical scheme for long-term simulations of 3D potential fully nonlinear periodic gravity waves is suggested. The scheme is based on the surface-following nonorthogonal curvilinear coordinate system. Velocity potential is represented as a sum of analytical and nonlinear components. The Poisson equation for the nonlinear component of velocity potential is solved iteratively. Fourier transform method, the second-order accuracy approximation of vertical derivatives on a stretched vertical grid and the fourth-order Runge–Kutta time stepping are used. The scheme is validated by simulation of steep Stokes waves. A one-processor version of the model for PC allows us to simulate evolution of a wave field with thousands degrees of freedom for hundreds of wave periods. The scheme is designed for investigation of nonlinear 2D surface waves, generation of extreme waves, and direct calculations of nonlinear interactions.     

Dense 3D residual moveout analysis as a tool for HTI parameter estimation

Three‐dimensional residual moveout analysis is the basic step in velocity model refinement. The analysis is generally carried out using horizontal and/or vertical semblances defined on a sparse set of in‐lines or cross‐lines with densely sampled source–receiver offsets. An alternative approach, which we call dense residual moveout analysis (DRMA), is to use all the bins of a three‐dimensional survey but sparsely sampled offsets. The proposed technique is very fast and provides unbiased and statistically efficient estimates of the residual moveout. Indeed, for the sparsest possible offset distribution, when only near‐ and far‐angle stacks are used, the variance of the residual moveout estimate is only 1.4 times larger than the variance of the least‐squares estimate obtained using all offsets. The high performance of DRMA makes it a useful tool for many applications, of which azimuthal velocity analysis is considered here. For a horizontal transverse isotropy (HTI) model, a deterministic procedure is proposed to define, at every point of residual moveout estimation, the azimuthal angle of the HTI axis of symmetry, the Thomsen anisotropy coefficients, and the interval (or root‐mean‐square) velocities in both the HTI isotropy and symmetry planes. The procedure is not restricted by DRMA assumptions; for example, it is also applicable to semblance‐based residual moveout estimates. The high resolution of the technique is illustrated by azimuthal velocity analysis over an oilfield in West Siberia.     

Models of the wave boundary layer

A general approach to model the structure of the wave boundary layer, based on the nonlinear Reynolds equations in a curvilinear system of coordinates, is described. Both spectral and numerical grid models are used. The energetic interactions between wind and wave in terms of Miles' parameter are studied. For waves outrunning or running against the wind, the range of the inverse flux of energy is found. For waves running slower than the wind, quadratic growth of is established. Vertical profiles of the wave momentum flux for different fetches are given. Following P. Janssen, a one-dimensional analytical model of the wave boundary layer is suggested. The effect of waves on the drag coefficient is analyzed.     

Differing responses of zircon,chevkinite-(Ce), monazite-(Ce) and fergusonite-(Y) to hydrothermal alteration: Evidence from the Keivy alkaline province,Kola Peninsula,Russia

Mineralogy and Petrology - A quartzolite from the Rova occurrence, Keivy alkali granite province, Kola Peninsula, Russia, is used to examine the differing responses of certain rare-metal minerals...     

Proterozoic mafic magmatism in Siberian craton: An overview and implications for paleocontinental reconstruction

We present a summary of late Paleoproterozoic to Neoproterozoic mafic magmatism in the Siberian craton, including recently published U–Pb and ⁴⁰Ar–³⁹Ar dates. These new precise ages suggest that at least some of the previously published K–Ar ages of Siberian mafic bodies should be ignored. The time–space geochronological chart, or the ‘barcode’ of mafic magmatic events shows significant differences between northern and southern Siberia. Both are characterized by ∼1900–1700 Ma magmatic events, but then there was an almost 1 Ga mafic magmatic ‘pause’ in south Siberia until ∼800 Ma. Meanwhile there are indications of multiple mafic magmatic events in North Siberia (Anabar shield and Olenek uplift) between ∼1600 and 1000 Ma. A series of magmatic events probably related to the breakup of Rodinia occurred in southern Siberia after ∼800 Ma. So far, there are no indications of late Neoproterozoic mafic magmatism in North Siberia. Ca. 1000–950 Ma mafic sills were reported from Meso- to Neo-Proterozoic sedimentary successions in the Sette-Daban area on the east side of the Siberian craton, but their tectonic setting is debated. Recent Ar–Ar dates of ∼1750 Ma for NW-trending dykes in the Aldan and Anabar shields, together with similar-age NNE-trending Baikal uplift dykes in south-eastern Siberia suggest the existence of a giant radial dyke swarm possibly related to a mantle plume centred in the Vilyui River area.