Moisture transport and seasonal variations in the stable isotopic composition of rainfall in Central American and Andean Páramo during El Niño conditions (2015–2016)

High‐elevation tropical grassland systems, called Páramo, provide essential ecosystem services such as water storage and supply for surrounding and lowland areas. Páramo systems are threatened by climate and land use changes. Rainfall generation processes and moisture transport pathways influencing precipitation in the Páramo are poorly understood but needed to estimate the impact of these changes, particularly during El Niño conditions, which largely affect hydrometeorological conditions in tropical regions. To fill this knowledge gap, we present a stable isotope analysis of rainfall samples collected on a daily to weekly basis between January 2015 and May 2016 during the strongest El Niño event recorded in history (2014–2016) in two Páramo regions of Central America (Chirripó, Costa Rica) and the northern Andes (Cajas, south Ecuador). Isotopic compositions were used to identify how rainfall generation processes (convective and orographic) change seasonally at each study site. Hybrid Single Particle Lagrangian Integrated Trajectory model (HYSPLIT) air mass back trajectory analysis was used to identify preferential moisture transport pathways to each Páramo site. Our results show the strong influence of north‐east trade winds to transport moisture from the Caribbean Sea to Chirripó and the South American low‐level jet to transport moisture from the Amazon forest to Cajas. These moisture contributions were also related to the formation of convective rainfall associated with the passage of the Intertropical Convergence Zone over Costa Rica and Ecuador during the wetter seasons and to orographic precipitation during the transition and drier seasons. Our findings provide essential baseline information for further research applications of water stable isotopes as tracers of rainfall generation processes and transport in the Páramo and other montane ecosystems in the tropics.     

Simulating damage evolution and fracture propagation in sandstone containing a preexisting 3-D surface flaw under uniaxial compression

Preexisting flaws and rock heterogeneity have important ramifications on the process of rock fracturing and on rock stability in many applications. Therefore, there is great interest in numerical modelling of rock fracture and the underlying mechanisms. We simulated damage evolution and fracture propagation in sandstone specimens containing a preexisting 3-D surface flaw under uniaxial compression. We applied the linear elastic damage model based on the unified strength theory following the rock failure process analysis code. However, in contrast to the rock failure process analysis code, we used the finite element method with tetrahedron elements on unstructured meshes. It provided higher geometrical flexibility and allowed for a more accurate representation of the disk-shaped flaw with various flaw depths, angles, and lengths through locally adapted meshes. The rock heterogeneity was modelled by sampling the initial local Young's modulus from a Weibull distribution over a cubic grid. The values were then interpolated to the computational finite element method mesh. This method introduced an additional length scale for the rock heterogeneity represented by the cell size in the sampling grid. The generation of three typical surface cracking patterns, called wing cracks, anti-wing cracks, and far-field cracks, were identified in the simulation results. These depend on the geometry of the preexisting surface flaw. The simulated fracture propagation, coalescence types, and failure modes for the specimens with preexisting surface flaw show good agreement with recent experimental studies.     

Effect of land cover and hydro‐meteorological controls on soil water DOC concentrations in a high‐elevation tropical environment

Páramo soils store high amounts of organic carbon. However, the effects of climate change and changes in land cover and use (LC/LU) in this high‐elevation tropical ecosystem may cause a decrease in their carbon storage capacity. Therefore, better understanding of the factors influencing the Páramo soils' carbon storage and export is urgently needed. To fill this knowledge gap, we investigated the differences in dissolved organic carbon (DOC) content in the soil water of four LC/LU types (tussock grass, natural forest, pine plantations, and pasture) and the factors controlling its variability in the Quinuas Ecohydrological Observatory in south Ecuador. Weekly measurements of soil water DOC concentrations, meteorological variables, soil water content, and temperature from various depths and slope positions were monitored within the soils' organic and mineral horizons between October 2014 and January 2017. These data were used to generate regression trees and random forest statistical models to identify the factors controlling soil water DOC concentrations. From high to low concentrations, natural forest depict the highest DOC concentrations followed by pasture, tussock grass, and pine forest. For all LC/LU types, DOC concentrations increase with decreasing soil moisture. Our results also show that LC/LU is the most important predictor of soil water DOC concentrations, followed by sampling depth and soil moisture. Interestingly, atmospheric variables and antecedent evapotranspiration and precipitation conditions show only little influence on DOC concentrations during the monitoring period. Our findings provide unique information that can help improve the management of soil and water resources in the Páramo and other peat dominated ecosystems elsewhere.     

Recent progress in the study of accessory minerals

Mineralogy and Petrology -     

Profile convexities in bedrock and alluvial streams

Longitudinal profiles of bedrock streams in central Kentucky, and of coastal plain streams in southeast Texas, were analyzed to determine the extent to which they exhibit smoothly concave profiles and to relate profile convexities to environmental controls. None of the Kentucky streams have smoothly concave profiles. Because all observed knickpoints are associated with vertical joints, if they are migrating it either occurs rapidly between vertical joints, or migrating knickpoints become stalled at structural features. These streams have been adjusting to downcutting of the Kentucky River for at least 1.3 Ma, suggesting that the time required to produce a concave profile is long compared to the typical timescale of environmental change. A graded concave longitudinal profile is not a reasonable prediction or benchmark condition for these streams. The characteristic profile forms of the Kentucky River gorge area are contingent on a particular combination of lithology, structure, hydrologic regime, and geomorphic history, and therefore do not represent any general type of equilibrium state. Few stream profiles in SE Texas conform to the ideal of the smoothly, strongly concave profile. Major convexities are caused by inherited topography, geologic controls, recent and contemporary geomorphic processes, and anthropic effects. Both the legacy of Quaternary environmental change and ongoing changes make it unlikely that consistent boundary conditions will exist for long. Further, the few exceptions within the study area–i.e., strongly and smoothly concave longitudinal profiles–suggest that ample time has occurred for strongly concave profiles to develop and that such profiles do not necessarily represent any mutual adjustments between slope, transport capacity, and sediment supply. The simplest explanation of any tendency toward concavity is related to basic constraints on channel steepness associated with geomechanical stability and minimum slopes necessary to convey flow. This constrained gradient concept (CGC) can explain the general tendency toward concavity in channels of sufficient size, with minimal lithological constraints and with sufficient time for adjustment. Unlike grade- or equilibrium-based theories, the CGC results in interpretations of convex or low-concavity profiles or reaches in terms of local environmental constraints and geomorphic histories rather than as “disequilibrium” features.     

New impact‐melt rock from the Roter Kamm impact structure,Namibia: Further constraints on impact age,melt rock chemistry,and projectile composition

Abstract— A new locality of in situ massive impact‐melt rock was discovered on the south‐southwestern rim of the Roter Kamm impact structure. While the sub‐samples from this new locality are relatively homogeneous at the hand specimen scale, and despite being from a nearby location, they do not have the same composition of the only previously analyzed impact‐melt rock sample from Roter Kamm. Both Roter Kamm impact‐melt rock samples analyzed to date, as well as several suevite samples, exhibit a granitic‐granodioritic precursor composition. Micro‐chemical analyses of glassy matrix and Al‐rich orthopyroxene microphenocrysts demonstrate rapid cooling and chemical disequilibrium at small scales. Platinum‐group element abundances and ratios indicate an ordinary chondritic composition for the Roter Kamm impactor. Laser argon dating of two sub‐samples did not reproduce the previously obtained age of 3.7 ± 0.3 (1s?) for this impact event, based on ⁴⁰Ar/³⁹Ar dating of a single vesicular impact‐melt rock. Instead, we obtained ages between 3.9 and 6.3 Ma, with an inverse isochron age of 4.7 ± 0.3 Ma for one analyzed sub‐sample and 5.1 ± 0.4 Ma for the other. Clearly a post‐5 Ma impact at Roter Kamm remains indicated, but further analytical work is required to better constrain the currently best estimate of 4–5 Ma. Both impactor and age constraints are clearly obstructed by the inherent microscopic heterogeneity and disequilibrium melting and cooling processes demonstrated in the present study.     

Forest ecosystems, disturbance, and climatic change in Washington State, USA

Climatic change is likely to affect Pacific Northwest (PNW) forests in several important ways. In this paper, we address the role of climate in four forest ecosystem processes and project the effects of future climatic change on these processes across Washington State. First, we relate Douglas-fir growth to climatic limitation and suggest that where Douglas-fir is currently water-limited, growth is likely to decline due to increased summer water deficit. Second, we use existing analyses of climatic controls on tree species biogeography to demonstrate that by the mid twenty-first century, climate will be less suitable for key species in some areas of Washington. Third, we examine the relationships between climate and the area burned by fire and project climatically driven regional and sub-regional increases in area burned. Fourth, we suggest that climatic change influences mountain pine beetle (MPB) outbreaks by increasing host-tree vulnerability and by shifting the region of climate suitability upward in elevation. The increased rates of disturbance by fire and mountain pine beetle are likely to be more significant agents of changes in forests in the twenty-first century than species turnover or declines in productivity, suggesting that understanding future disturbance regimes is critical for successful adaptation to climate change.     

EXOTIME: searching for planets around pulsating subdwarf B stars

In 2007, a companion with planetary mass was found around the pulsating subdwarf B star V391 Pegasi with the timing method, indicating that a previously undiscovered population of substellar companions to apparently single subdwarf B stars might exist. Following this serendipitous discovery, the EXOTIME (http://www.na.astro.it/~silvotti/exotime/) monitoring program has been set up to follow the pulsations of a number of selected rapidly pulsating subdwarf B stars on time scales of several years with two immediate observational goals:

1. determine $\dot{P}$ of the pulsational periods P
2. search for signatures of substellar companions in O–C residuals due to periodic light travel time variations, which would be tracking the central star’s companion-induced wobble around the centre of mass

These sets of data should therefore, at the same time, on the one hand be useful to provide extra constraints for classical asteroseismological exercises from the $\dot{P}$ (comparison with “local” evolutionary models), and on the other hand allow one to investigate the preceding evolution of a target in terms of possible “binary” evolution by extending the otherwise unsuccessful search for companions to potentially very low masses. While timing pulsations may be an observationally expensive method to search for companions, it samples a different range of orbital parameters, inaccessible through orbital photometric effects or the radial velocity method: the latter favours massive close-in companions, whereas the timing method becomes increasingly more sensitive toward wider separations. In this paper we report on the status of the on-going observations and coherence analysis for two of the currently five targets, revealing very well-behaved pulsational characteristics in HS?0444+0458, while showing HS?0702+6043 to be more complex than previously thought.     

城市环境气候图的发展及其应用现状

城市环境气候图最初于20世纪70年代由德国气候研究者开发制作,通过近40年的发展,目前在世界范围内已有20个国家开展了相关研究与应用项目,从而为改善城市气候环境与提高人居生活条件提供决策依据。该文首先对城市环境气候图的发展及其制作方法进行了一个总体性的回顾,然后选取德国、日本和中国香港地区的研究案例用以考察目前该领域研究和应用现状与存在问题,进一步详细探讨了未来发展的方向,最后指出了在我国城市化进程中为应对诸多气候环境问题所迫切需要开展的相关研究问题。     

Majoritic garnet grains within shock-induced melt veins in amphibolites from the Ries impact crater suggest ultrahigh crystallization pressures between 18 and 9 GPa

Shock-induced melt veins in amphibolites from the Nördlinger Ries often have chemical compositions that are similar to bulk rock (i.e., basaltic), but there are other veins that are confined to chlorite-rich cracks that formed before the impact and these are poor in Ca and Na. Majoritic garnets within the shock veins show a broad chemical variation between three endmembers: (1) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}{\text{Al}}_{2} ({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (normal garnet, Grt), (2) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}[{\text{M}}^{2 + } ({\text{Si,Ti}})]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\)  (majorite, Maj), and (3) \({}^{\text{VIII}}({{\text {Na} {\text M}^{2+}}_2}) {}^{\text{VI}}[ ({\text{Si,Ti}}){\text {Al}}]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (Na-majorite₅₀Grt₅₀), whereby M²⁺ = Mg²⁺, Fe²⁺, Mn²⁺, Ca²⁺. In particular, we observed a broad variation in ^VI(Si,Ti) which ranges from 0.12 to 0.58 cations per formula unit (cpfu). All these majoritic garnets crystallized during shock pressure release at different ultrahigh pressures. Those with high ^VI(Si,Ti) (0.36–0.58 cpfu) formed at high pressures and temperatures from amphibole-rich melts, while majoritic garnets with lower ^VI(Si,Ti) of 0.12–0.27 cpfu formed at lower pressures and temperatures from chlorite-rich melts. Furthermore, majoritic garnets with intermediate values of ^VI(Si,Ti) (0.24–0.39) crystallized from melts with intermediate contents of Ca and Na. To the best of our knowledge the ‘MORB-type’ Ca–Na-rich majoritic garnets with maximum contents of 2.99 wt% Na₂O and calculated crystallisation pressures of 16–18 GPa are the most extreme representatives ever found in terrestrial shocked materials. At the Ries, the duration of the initial contact and compression stage at the central location of impact is estimated to only ~ 0.1 s. We used a ~ 200-µm-thick shock-induced vein in a moderately shocked amphibolite to model its pressure–temperature–time (P–T–t) path. The graphic model manifests a peak temperature of ~ 2600 °C for the vein, continuum pressure lasting for ~ 0.02 s, a quench duration of ~ 0.02 s and a shock pulse of ~ 0.038 s. The small difference between the continuum pressure and the pressure of majoritic garnet crystallization underlines the usefulness of applying crystallisation pressures of majoritic garnets from metabasites for calculation of dynamic shock pressures of host rocks. Majoritic garnets of chlorite provenance, however, are not suitable for the determination of continuum pressure since they crystallized relatively late during shock release. An extraordinary glass- and majorite-bearing amphibole fragment in a shock-vein of one amphibolite documents the whole unloading path.