Tides and mixing in the northwestern East China Sea Part I: Rotating and reversing tidal flows

Bottom-mounted ADV and ADCP instruments in combination with CTD profiling measurements taken along the Chinese coast of the East China Sea were used to study the vertical structure of temperature, salinity, and velocity in reversing tidal currents on a shallow inner shelf and in rotating tidal flows over a deeper sloping bottom of the outer shelf. These two regimes of barotropic tide affect small-scale dynamics in the lower part of the water column differently. The reversing flow was superimposed by seiches of ∼2.3 h period generated in semienclosed Jiaozhou Bay located nearby. As the tidal vector rotates over the sloping bottom, the height of the near-bottom logarithmic layer is subjected to tidal-induced variations. A maximum of horizontal velocity U_max appears at the upper boundary of the log layer during the first half of the current vector rotation from the minor to the major axis of tidal ellipse. In rotating tidal flow, vertical shear generated at the seafloor, propagated slowly to the water interior up to the height of U_max, with a phase speed of ∼5 m/h. The time-shifted shear inside the water column, relative to the shear at the bottom, was associated with periodically changing increases and decreases of the tidal velocity above the log layer toward the sea surface. In reversing flows, the shear generated near the bottom and the shear at the upper levels were almost in phase.     

Late Carboniferous porphyry copper mineralization at La Voluntad,Neuquén,Argentina: Constraints from Re–Os molybdenite dating

The La Voluntad porphyry Cu–Mo deposit in Neuquén, Argentina, is one of several poorly known porphyry-type deposits of Paleozoic to Early Jurassic age in the central and southern Andes. Mineralization at La Voluntad is related to a tonalite porphyry from the Chachil Plutonic Complex that intruded metasedimentary units of the Piedra Santa Complex. Five new Re–Os molybdenite ages from four samples representing three different vein types (i.e., quartz–molybdenite, quartz–sericite–molybdenite and quartz–sericite–molybdenite ± chalcopyrite–pyrite) are identical within error and were formed between ~312 to ~316 Ma. Rhenium and Os concentrations range between 34 to 183 ppm and 112 to 599 ppb, respectively. The new Re–Os ages indicate that the main mineralization event at La Voluntad, associated to sericitic alteration, was emplaced during a time span of 1.7 ± 3.2 Ma and that the deposit is Carboniferous in age, not Permian as previously thought. La Voluntad is the oldest porphyry copper deposit so far recognized in the Andes and indicates the presence of an active magmatic arc, with associated porphyry style mineralization, at the proto-Pacific margin of Gondwana during the Early Pennsylvanian.     

Thermal expansion and high-temperature P21/c–C2/c phase transition in clinopyroxene-type LiFeGe2O6 and comparison to NaFe(Si,Ge)2O6

A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) Å, β = 108.844(6)°, V = 442.75(16) ų, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the $(\bar{1}\,0\,1) A synthetic clinopyroxene with composition LiFe³⁺Ge₂O₆, monoclinic s.g. P2₁/c, a = 9.8792(7), b = 8.8095(5), c = 5.3754(3) ?, β = 108.844(6)°, V = 442.75(16) ?³, has been studied by in situ low- and high-temperature single-crystal X-ray diffraction. The variation of lattice parameters and the intensity of the b-type reflections (h + k = 2n + 1, only present in the P-symmetry) with increasing temperature showed a displacive phase transition from space group P2₁/c to C2/c at a transition temperature T _tr = 789 K, first order in character, with a sudden volume increase of 1.6% and a decrease of β by 1° at the transition. This spontaneous dilatation is reversible, shows a limited hysteresis of ±10°C, and corresponds to the vanishing of the b-type reflections, thus indicating a symmetry increase to space group C2/c. Below T _tr an expansion is observed for all the cell parameters, while the β angle remained almost constant; at T > T _tr the thermal volume expansion is due to dilatation of the structure in the ([`1] 0 1)(\bar{1}\,0\,1) plane, mostly along [0 1 0], and pure shear in the (0 1 0) plane due to the decrease of β. From comparison with silicate analogues, the germanate clinopyroxenes are more expansible, while the P2₁/c expands more than the C2/c phase. The evolution of Q ² (calculated as the normalized intensity of b-type reflections) with T in the framework of the Landau theory has been done using a standard expression for a first order phase transition. We observe a jump of Q ₀² = 0.538(2) at T _tr, with T _c of 481(7) K, b/a = −2,290 K, and c/a = 3,192 K, and thus far from being tri-critical point. A closely related composition (LiFe³⁺Si₂O₆) shows an equivalent phase transition at 228 K, which is very close to the tri-critical point and 561 K cooler. This result indicates that a change in the composition of tetrahedral sites can have dramatic effects on the P2₁/c ↔ C2/c displacive phase transition in clinopyroxenes. The major changes observed in the evolution of the crystal structure with T are observed in the M2 polyhedron, with a volume decrease by ca. 13.3%, compared to ca. 1.3% observed in the M1 polyhedron. The tetrahedra behave as rigid units with neither a significant change of volume at T > T _tr (<1‰), nor a change of tilting of the basal plane. No change in coordination is observed at T > T _tr in the M2 polyhedron, which remains sixfold coordinated although a strong deformation of this polyhedron is observed. This deformation is related to a strong change by 51.4° at T _tr of the kinking angle (O3–O3–O3 angle) of the B-chain of tetrahedra, which switches from O-rotated to S-rotated [from 143.3(5)° to 194.7(6)°]. The A-chain is S-rotated at T < T _tr [206.8(5)° at 703 K] and extends by 12° at the transition.     

Distribuzione geografica dell'attivita' magnetica terrestre

Riassunto Richiamate le caratteristiche secondo cui le perturbazioni regolari ed irregolari del campo terrestre concorrono a definire l'attività geomagnetica, se ne analizza la distribuzione geografica mediante le grandezze:HR _H ,ZR _Z , e (10^–4.²). Di queste sono illustrati e discussi gli andamenti in funzione della latitudine geomagnetica, dedotti dai valori medi calcolati per quasi tutte le stazioni funzionanti durante l'Anno Polare 1932–33 ed il triennio successivo 1933–36. Si rileva particolarmente il massimo assoluto di attività della zone anulare delle aurore, e quello relativo delle regioni equatoriali. Segue un'indagine, a carattere piuttosto indicativo, sulla ripartizione geografica del grado di variazione dell'attività nel tempo, introducendo a tale scopo opportuni indici ed esponendone il comportamento. Conclude la Nota un cenno sugli andamenti stagionali dell'attività e della variabilità geomagnetica.     

An audio-magnetotelluric investigation in Terceira Island (Azores)

Ten audio-magnetotelluric soundings have been carried out along a profile crossing the Serra do Cume caldera in the eastern part of the Terceira Island (Azores). The main objectives of this investigation were to detect geoelectrical features related with tectonic structures and to characterize regional hydrological and hydrothermal aspects mainly those related to geothermal fluid dynamics.Three-dimensional numerical investigation showed that the data acquired at periods shorter than 1 s are not significantly affected by ocean effect. The data was analysed using the Smith's decomposition method in order to investigate possible distortions caused by superficial structures and to estimate a global regional strike. The results suggest that in general the soundings were not distorted. A regional N55°W strike was chosen for the two-dimensional data inversion.The low-resistivity zones (10–30 ohm-m) displayed in the central part of the 2-D geoelectrical model have been interpreted as caused by hydrothermal circulation. The low-resistivity anomalies at the ends of the profile might be attributed to alteration zones with interaction of seawater intrusion. High-resistivity (> 300 ohm-m) values have been related with less permeable zones in the SW of Cinco Picos and Guilherme Moniz caldera walls.     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

A New Combined Assessment of Mixed Uncertainty in Spatial Models: Conceptualization and Implementation

Uncertainty quantification is not often performed in spatial modeling applications, especially when there is a mixture of probabilistic and non‐probabilistic uncertainties. Furthermore, the effect of positional uncertainty is often not assessed, despite its relevance to geographical applications. Although there has been much work in investigating the aforementioned types of uncertainty in isolation, combined approaches have not been much researched. This has resulted in a lack of tools for conducting mixed uncertainty analyses that include positional uncertainty. This research addresses the issue by first presenting a new, flexible, simulation‐oriented conceptualization of positional uncertainty in geographic objects called F‐Objects. F‐Objects accommodates various representations of uncertainty, while remaining conceptually simple. Second, a new Python‐based framework is introduced, termed Wiggly and capable of conducting mixed uncertainty propagation using fuzzy Monte Carlo simulation (FMCS). FMCS combines both traditional Monte Carlo with fuzzy analysis in a so‐called hybrid approach. F‐Objects is implemented within the Wiggly framework, resulting in a tool capable of considering any combination of: (1) probabilistic variables; (2) fuzzy variables; and (3) positional uncertainty of objects (probabilistic/fuzzy). Finally, a realistic GIS‐based groundwater contamination problem demonstrates how F‐Objects and Wiggly can be used to assess the effect of positional uncertainty.     

Biomass expansion factor and root-to-shoot ratio for Pinus in Brazil

The Biomass Expansion Factor (BEF) and the Root-to-Shoot Ratio (R) are variables used to quantify carbon stock in forests. They are often considered as constant or species/area specific values in most studies. This study aimed at showing tree size and age dependence upon BEF and R and proposed equations to improve forest biomass and carbon stock. Data from 70 sample Pinus spp. grown in southern Brazil trees in different diameter classes and ages were used to demonstrate the correlation between BEF and R, and forest inventory data, such as DBH, tree height and age. Total dry biomass, carbon stock and CO₂ equivalent were simulated using the IPCC default values of BEF and R, corresponding average calculated from data used in this study, as well as the values estimated by regression equations. The mean values of BEF and R calculated in this study were 1.47 and 0.17, respectively. The relationship between BEF and R and the tree measurement variables were inversely related with negative exponential behavior. Simulations indicated that use of fixed values of BEF and R, either IPCC default or current average data, may lead to unreliable estimates of carbon stock inventories and CDM projects. It was concluded that accounting for the variations in BEF and R and using regression equations to relate them to DBH, tree height and age, is fundamental in obtaining reliable estimates of forest tree biomass, carbon sink and CO₂ equivalent.     

Real-time cycle slip detection in triple-frequency GNSS

The modernization of the global positioning system and the advent of the European project Galileo will lead to a multifrequency global navigation satellite system (GNSS). The presence of new frequencies introduces more degrees of freedom in the GNSS data combination. We define linear combinations of GNSS observations with the aim to detect and correct cycle slips in real time. In particular, the detection is based on five geometry-free linear combinations used in three cascading steps. Most of the jumps are detected in the first step using three minimum-noise combinations of phase and code observations. The remaining jumps with very small amplitude are detected in the other two steps by means of two-tailored linear combinations of phase observations. Once the epoch of the slip has been detected, its amplitude is estimated using other linear combinations of phase observations. These combinations are defined with the aim of discriminating between the possible combinations of jump amplitudes in the three carriers. The method has been tested on simulated data and 1-second triple-frequency undifferenced GPS data coming from a friendly multipath environment. Results show that the proposed method is able to detect and repair all combinations of cycle slips in the three carriers.     

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.