Peat depth as a control on Sphagnum moisture stress during seasonal drought

Peatlands are globally important long-term sinks of carbon, however there is concern that enhanced peat decomposition and moss moisture stress due to climate change mediated drought will reduce moss productivity making these ecosystems vulnerable to carbon loss and associated long-term degradation. Peatlands are resilient to summer drought moss stress because of negative ecohydrological feedbacks that generally maintain a wet peat surface, but where feedbacks may be contingent on peat depth. We tested this ‘survival of the deepest’ hypothesis by examining water table (WT) position, near-surface moisture content, and soil water tension in peatlands that differ in size, peat depth, and catchment area during a summer drought. All shallow sites (<40 cm depth) lost their WT (i.e., the groundwater well was dry) for considerable time during the drought period. Near-surface soil water tension increased dramatically at shallow sites following WT loss, increasing ~5–7.5× greater at shallow sites compared to deep sites (≥40 cm depth). During a mid-summer drought intensive field survey, we found that 60–67% of plots at shallow sites exceeded a 100 mb tension threshold used to infer moss water stress. Unlike the shallow sites, tension typically did not exceed this 100 mb threshold at the deep sites. Using species dependent water content – chlorophyll fluorescence thresholds and relations between volumetric water content and WT depth, Monte Carlo simulations suggest that moss had nearly twice the likelihood of being stressed at shallow sites (0.38 ± 0.24) compared to deep sites (0.22 ± 0.18). This study provides evidence that mosses in shallow peatland may be particularly vulnerable to warmer and drier climates in the future, but where species composition may play an important role. We argue that a critical ‘threshold’ peat depth specific for different hydrogeological and hydroclimatic regions can be used to assess what peatlands are especially vulnerable to climate change mediated drought.     

Stratigraphy and Mesozoic–Cenozoic tectonic history of northern Sierra Los Ajos and adjacent areas,Sonora, Mexico

Geologic mapping in the northern Sierra Los Ajos reveals new stratigraphic and structural data relevant to deciphering the Mesozoic–Cenozoic tectonic evolution of the range. The northern Sierra Los Ajos is cored by Proterozoic, Cambrian, Devonian, Mississippian, and Pennsylvanian strata, equivalent respectively to the Pinal Schist, Bolsa Quartzite and Abrigo Limestone, Martin Formation, Escabrosa Limestone, and Horquilla Limestone. The Proterozoic–Paleozoic sequence is mantled by Upper Cretaceous rocks partly equivalent to the Fort Crittenden and Salero Formations in Arizona, and the Cabullona Group in Sonora, Mexico.Absence of the Upper Jurassic–Lower Cretaceous Bisbee Group below the Upper Cretaceous rocks and above the Proterozoic–Paleozoic rocks indicates that the Sierra Los Ajos was part of the Cananea high, a topographic highland during the Late Jurassic and Early Cretaceous. Deposition of Upper Cretaceous rocks directly on Paleozoic and Proterozoic rocks indicates that the Sierra Los Ajos area had subsided as part of the Laramide Cabullona basin during Late Cretaceous time. Basal beds of the Upper Cretaceous sequence are clast-supported conglomerate composed locally of basement (Paleozoic) clasts. The conglomerate represents erosion of Paleozoic basement in the Sierra Los Ajos area coincident with development of the Cabullona basin.The present-day Sierra Los Ajos reaches elevations of greater than 2600 m, and was uplifted during Tertiary basin-and-range extension. Upper Cretaceous rocks are exposed at higher elevations in the northern Sierra Los Ajos and represent an uplifted part of the inverted Cabullona basin. Tertiary uplift of the Sierra Los Ajos was largely accommodated by vertical movement along the north-to-northwest-striking Sierra Los Ajos fault zone flanking the west side of the range. This fault zone structurally controls the configuration of the headwaters of the San Pedro River basin, an important bi-national water resource in the US-Mexico border region.     

Nanophase magnetite in matrix of anomalous EL3 chondrite Northwest Africa (NWA) 8785

NWA 8785 is a remarkable EL3 chondrite with a high abundance (~34 vol%) of an Fe-rich matrix. This is the highest matrix abundance known among enstatite chondrites (ECs) and more similar to the matrix abundances in some carbonaceous and Rumuruti chondrites. X-ray diffraction and TEM data indicate that the fine-grained portion of the NWA 8785 matrix consists of nanoscale magnetite mixed with a noncrystalline silicate material and submicron-sized enstatite and plagioclase grains. This is the first report of magnetite nanoparticles in an EL3. The Si content of the metal (0.7 wt%), presence of ferroan alabandite, and its O isotopic composition indicate NWA 8785 is EL3-related. Having more abundant matrix than in other ECs, and that the matrix is rich in magnetite nanoparticles, which are not present in any other EC, suggest classification as an EL3 anomalous. Although we cannot completely exclude any of the mechanisms or environments for formation of the magnetite, we find a secondary origin to be the most compelling. We suggest that the magnetite formed due to hydrothermal activity in the meteorite parent body. Although ECs are relatively dry and likely formed within the nebular snow line, ices may have drifted inward from just beyond the snow line to the region where the EL chondrites were accreting, or more likely the snow line migrated inward during the early evolution of the solar system. This may have resulted in the condensation of ices and provided an ice-rich region for accretion of the EL3 parent body. Thus, the EL3 parent body may have had hydrothermal activity and if Earth formed near the EC accretion zone, similar bodies may have contributed to the Earth's water supply. NWA 8785 greatly extends the range of known characteristics of ECs and EC parent body processes.     

Intensification of decadal and multi-decadal sea level variability in the western tropical Pacific during recent decades

Previous studies have linked the rapid sea level rise (SLR) in the western tropical Pacific (WTP) since the early 1990s to the Pacific decadal climate modes, notably the Pacific Decadal Oscillation in the north Pacific or Interdecadal Pacific Oscillation (IPO) considering its basin wide signature. Here, the authors investigate the changing patterns of decadal (10–20 years) and multidecadal (>20 years) sea level variability (global mean SLR removed) in the Pacific associated with the IPO, by analyzing satellite and in situ observations, together with reconstructed and reanalysis products, and performing ocean and atmosphere model experiments. Robust intensification is detected for both decadal and multidecadal sea level variability in the WTP since the early 1990s. The IPO intensity, however, did not increase and thus cannot explain the faster SLR. The observed, accelerated WTP SLR results from the combined effects of Indian Ocean and WTP warming and central-eastern tropical Pacific cooling associated with the IPO cold transition. The warm Indian Ocean acts in concert with the warm WTP and cold central-eastern tropical Pacific to drive intensified easterlies and negative Ekman pumping velocity in western-central tropical Pacific, thereby enhancing the western tropical Pacific SLR. On decadal timescales, the intensified sea level variability since the late 1980s or early 1990s results from the “out of phase” relationship of sea surface temperature anomalies between the Indian and central-eastern tropical Pacific since 1985, which produces “in phase” effects on the WTP sea level variability.     

Simulating Dispersion in the Evening-Transition Boundary Layer

We investigate dispersion in the evening-transition boundary layer using large-eddy simulation (LES). In the LES, a particle model traces pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent subgrid-scale motions. The LES is forced with both a sudden switch-off of the surface heat flux and also a more gradual observed evolution. The LES shows ‘lofting’ of plumes from near-surface releases in the pre-transition convective boundary layer; it also shows the subsequent ‘trapping’ of releases in the post-transition near-surface stable boundary layer and residual layer above. Given the paucity of observations for pollution dispersion in evening transitions, the LES proves a useful reference. We then use the LES to test and improve a one-dimensional Lagrangian Stochastic Model (LSM) such as is often used in practical dispersion studies. The LSM used here includes both time-varying and skewed turbulence statistics. It is forced with the vertical velocity variance, skewness and dissipation from the LES for particle releases at various heights and times in the evening transition. The LSM plume spreads are significantly larger than those from the LES in the post-transition stable boundary-layer trapping regime. The forcing from the LES was thus insufficient to constrain the plume evolution, and inclusion of the significant stratification effects was required. In the so-called modified LSM, a correction to the vertical velocity variance was included to represent the effect of stable stratification and the consequent presence of wave-like motions. The modified LSM shows improved trapping of particles in the post-transition stable boundary layer.     

An empirical model of tropical ocean dynamics

To extend the linear stochastically forced paradigm of tropical sea surface temperature (SST) variability to the subsurface ocean, a linear inverse model (LIM) is constructed from the simultaneous and 3-month lag covariances of observed 3-month running mean anomalies of SST, thermocline depth, and zonal wind stress. This LIM is then used to identify the empirically-determined linear dynamics with physical processes to gauge their relative importance to ENSO evolution. Optimal growth of SST anomalies over several months is triggered by both an initial SST anomaly and a central equatorial Pacific thermocline anomaly that propagates slowly eastward while leading the amplifying SST anomaly. The initial SST and thermocline anomalies each produce roughly half the SST amplification. If interactions between the sea surface and the thermocline are removed in the linear dynamical operator, the SST anomaly undergoes less optimal growth but is also more persistent, and its location shifts from the eastern to central Pacific. Optimal growth is also found to be essentially the result of two stable eigenmodes with similar structure but differing 2- and 4-year periods evolving from initial destructive to constructive interference. Variations among ENSO events could then be a consequence not of changing stability characteristics but of random excitation of these two eigenmodes, which represent different balances between surface and subsurface coupled dynamics. As found in previous studies, the impact of the additional variables on LIM SST forecasts is relatively small for short time scales. Over time intervals greater than about 9?months, however, the additional variables both significantly enhance forecast skill and predict lag covariances and associated power spectra whose closer agreement with observations enhances the validation of the linear model. Moreover, a secondary type of optimal growth exists that is not present in a LIM constructed from SST alone, in which initial SST anomalies in the southwest tropical Pacific and Indian ocean play a larger role than on shorter time scales, apparently driving sustained off-equatorial wind stress anomalies in the eastern Pacific that result in a more persistent equatorial thermocline anomaly and a more protracted (and predictable) ENSO event.     

The use of Lagrangian trajectories for the identification of the environmentally safe fairways

We propose and test a method for the optimisation of marine fairways to minimise the risk to high-value areas, based on statistical analysis of Lagrangian trajectories of current-driven pollution transport. The offshore areas are quantified according to the probability of pollution released in these areas to reach vulnerable regions. The method contains an eddy-resolving circulation model, a scheme for tracking of Lagrangian trajectories, a technique for the calculation of quantities characterising the potential of different sea areas to supply adverse impacts, and routines to construct the optimum fairway. The gain is expressed in terms of the probability of pollution transport to the nearshore and the associated time (particle age). The use of the optimum fairway would decrease the probability of coastal pollution by 40% or increase the average time of reaching the pollution to the coast from 5.3 to about 9 days in the Gulf of Finland, the Baltic Sea.