Crustal xenoliths from Calbuco Volcano,Andean Southern Volcanic Zone: implications for crustal composition and magma-crust interaction

Crustal xenoliths in the 1961 andesite flow of Calbuco Volcano, in the southern Southern Volcanic Zone (SSVZ) of the Andes, consist predominantly of pyroxene granulites and hornblende gabbronorites. The granulites contain plagioclase+pyroxene+magnetite±amphibole, and have pronounced granoblastic textures. Small amounts of relict amphibole surrounded by pyroxene-plagioclase-magnetite-glass symplectites are found in some specimens. These and similar textures in the gabbronorites are interpreted as evidence of dehydration melting. Mineral and bulk rock geochemical data indicate that the granulites are derived from an incompatible trace element depleted basaltic protolith that underwent two stages of metamorphism: a moderate pressure, high temperature stage accompanied by melting and melt extraction from some samples, followed by thermal metamorphism after entrainment in the Calbuco andesite lavas. High _Nd ^T values (+4.0 to +8.6), Nd-isotope model ages of 1.7–2.0 Ga, and trace element characteristics like chondrite normalized La/Yb< and La/Nb1 indicate that the protoliths were oceanic basalts. Similar oceanic metabasalts of greenschist to amphibolite facies are found in the Paleozoic metamorphic belt that underlies the Chilean coastal ranges. Mineral and bulk rock compositions of the gabbronorite xenoliths indicate that they are cognate, crystallizing from the basaltic andesite magma at Calbuco. Crystallization pressures for the gabbros based on total Al contents in amphibole are 6–8 kbar. These pressures point to middle to lower crustal storage of the Calbuco magma. Neither granulite nor gabbro xenoliths have the appropriate geochemical characteristics to be contaminants of Calbuco andesites, although an ancient sedimentary contaminant is indicated by the lava compositions. The presence of oceanic metabasaltic xenoliths, together with the sedimentary isotopic imprint, suggests that the lower crust beneath the volcano is analogous to the coastal metamorphic belt, which is an accretionary complex of intercalated basalts and sediments that formed along the Paleozoic Gondwanan margin. If this is the case, the geochemical composition of the lower and middle crust beneath the SSVZ is significantly different from that of most recent SSVZ volcanic rocks.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Nest microhabitat selection in larks from a European semi-arid shrub-steppe: the role of sunlight and predation

We analysed nest site selection in two shrub-steppe larks, the lesser short-toed larkCalandrella rufescensand the Thekla larkGalerida theklae. Both species nest on the ground next to a bush preferably facing north or north-east. Nestling growth and fledging weight are greater when the nest faces these directions. Both species prefer to nest next to camephytes with little foliage, and the lesser short-toed lark chooses sites with high visibility. We conclude that sunlight and nest predation are fundamental selective pressures determining nest microhabitat selection in larks in semi-arid areas.     

Groundwater contamination with arsenic and other trace elements in an area of the pampa,province of Córdoba,Argentina

A geochemical study of groundwater of the pampa in the province of Córdoba, Argentina, was performed; the area covered approximately 10,000 km².Physical-chemical parameters, dissolved solids, and seven trace elements were determined in 60 selected water samples. Systematic and accurate measurements of arsenic, flourine, and vanadium were performed for the first time. Three trace element contaminants not reported earlier were found: an important one, selenium, and two others of less known effects, uranium and molybdenum.Eighty-four percent of the water analyzed showed arsenic contents over 0.05 mg/L, maximum contaminant level established by the U.S. Environmental Protection Agency (1982). The frequency distribution of trace elements was analyzed, and its fit to the lognormal distribution was proved by means of the Pearson and Kolmogorov-Smirnov test; the geographic distribution of the seven trace elements was mapped and its correlation with the anion-cation composition of the water was studied.The maximum arsenic, fluorine, vanadium, and uranium contents were found in the western part of the area under study, in waters containing dominant alkali metals in the cation composition. Maximum selenium and antimony contents were found in the eastern part of the area, while molybdenum distribution does not show any relationship with the other two groups. In addition, the geographic distribution of the trace elements seems to be related to the subsurface structure, which has been inferred using interactive digital analysis of Landsat imagery. The movements of the subsoil have disturbed surface and subsurface drainage influencing the water salinity and trace element contents.In order to investigate the origin of the contamination, 54 loess samples were collected in wells at depths ranging from the surface down to the water table. This loess, which has a high proportion of volcanic components, mainly rhyolitic glass, exhibits a chemical composition corresponding to that of a dacite.The loess and the volcanic glass show anomalous contents of all contaminant trace elements, mainly arsenic and selenium. For this reason loess is considered the most important contamination source in the groundwater under study.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Twenty years of progresses in oceanic rogue waves: the role played by weakly nonlinear models

Here we discuss some of the progresses that have been made in the last 20 years in the field of oceanic rogue waves, focusing on the role played by leading order equations such as the nonlinear Schrödinger and the Korteweg-De Vries equations. For such equations, it is possible, as shown in Onorato et al. (Origin of heavy tail statistics in equations of the nonlinear Schrödinger type: an exact result, 2016. arXiv:1601.04317), to derive a very simple relation in which the variation of the third (for the KdV) and fourth (for the NLS) moment of the probability density function of the wave field can be related to the variation of the spectral bandwidth. These relations give some new perspectives on the formation of rogue waves in a random sea state.     

Population genetic diversity and structure of a dominant tropical seagrass,Cymodocea rotundata,in the Western Pacific region

Cymodocea rotundata is an ecologically important tropical pioneer seagrass species distributed in the Indo‐Pacific region. The population genetic diversity and structure of this species were analysed at 46 sites spanning the Philippines, Ryukyu Islands (northern limit) and Hainan Island, by using microsatellite simple sequence repeat (SSR) markers. Analyses revealed the persistence of C. rotundata likely relies on local population dynamics and fitness influenced by environmental gradients, with sexual reproduction prevalent in the Philippines while the Ryukyu and Hainan populations were predominantly established by clonal spread. Analysis of molecular variance showed significant genetic differentiation (P < 0.001) among three geographic regions: the Philippines, Ryukyu and Hainan. Furthermore, the mean fixation index value was very high (F_ST = 0.36), indicating poor dispersal potential or limited gene flow. Allelic richness and heterozygosity of C. rotundata was comparable in the Philippines and Ryukyu Islands populations. More private alleles were found in the Philippines and excess heterozygotes in the Ryukyu Islands. STRUCTURE analysis revealed that the Ryukyu and Hainan populations were mosaics of admixed alleles of individuals from the Philippines. An assignment test suggested that recruitment occurs from the Northeast Philippines to Ryukyu Islands. These results suggest that the Philippine populations are tropical seagrass hotspots and perhaps the origin of the populations in the Ryukyu and Hainan Islands. This scenario is most likely driven by the Kuroshio Current and island integration events during plate tectonic activities. The Philippine Archipelago is considered of high importance for conservation objectives and management plan of seagrasses. Likewise, the Ryukyu Islands are also important in terms of conserving regional diversity because locally adapted genotypes have important evolutionary potential in the face of environmental change.     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.