Wind as a cooling agent: substrate temperatures are responsible for variable lithobiont‐induced weathering patterns on west‐ and east‐facing limestone bedrock of the Negev

Spatial variability in lithobiont‐induced weathering patterns on desert rocks is aspect‐dependent. While differences between the northern and southern aspects have been extensively studied, little is known concerning the differences between east‐facing (EF) and west‐facing (WF) aspects in deserts, including the Negev Desert. Whereas cobbles on both slopes are inhabited by endolithic lichens, epilithic lichens, which render the bedrock a smooth appearance, and free‐living cyanobacteria, which give the bedrock a rugged microrelief, predominate on WF and EF bedrock, respectively. Following previous research that regarded dew as the principal factor that determines lithobiont distribution, measurements of radiation, temperature, wind and dew were carried out during 2008–2009 in the Negev Desert. The data indicated that albeit slightly higher midday surface temperatures that characterize WF surfaces (cobbles and bedrock), nocturnal temperatures on these surfaces were significantly lower, therefore facilitating higher dew condensation. High amounts of dew result from the relatively rapid drop in temperatures (14:00–20:00) due to the afternoon northwesterly sea‐breeze wind (with a cooling rate of the WF bedrock being 52.9% higher than on EF bedrock, 2.6 °C h^?1 in comparison to only 1.7 °C h^?1), and facilitate the growth of high‐chlorophyll dew‐fed (and rain‐fed) epilithic lichens, which may act as bio‐protectors on WF bedrock. Lack of condensation on EF bedrock results in turn in the growth of rain‐fed free‐living cyanobacteria, responsible for high rock dissolution and subsequently for a rugged microrelief. By affecting the nocturnal bedrock temperatures, wind acts as a cooling agent, impacting in turn the amount of dew, and subsequently lithobiont composition and weathering patterns in the Negev Desert. Copyright © 2016 John Wiley & Sons, Ltd. Copyright © 2016 John Wiley & Sons, Ltd.     

The Llullaillaco volcano, northwest Argentina: construction by Pleistocene volcanism and destruction by sector collapse

Llullaillaco is one of a chain of Quaternary stratovolcanoes that defines the present Andean Central Volcanic Zone (CVZ), and marks the border between Chile and Argentina/Bolivia. The current edifice is constructed from a series of thick dacitic lava flows, forming the second tallest active volcano in the world (6739 m). K–Ar and new biotite laser ⁴⁰Ar/³⁹Ar step-heating dates indicate that the volcano was constructed during the Pleistocene (≤1.5 Ma), with a youngest date of 0.048±0.012 Ma being recorded for a fresh dacite flow that descends the southern flank. Additional ⁴⁰Ar/³⁹Ar measurements for andesitic and dacitic lava flows from the surrounding volcanic terrain yield dates of between 11.94±0.13 Ma and 5.48±0.07 Ma, corresponding to an extended period of Miocene volcanism which defines much of the landscape in this region. Major- and trace-element compositions of lavas from Llullaillaco are typical of Miocene–Pleistocene volcanic rocks from the western margin of the CVZ, and are related to relatively shallow-dipping subduction of the Nazca plate beneath northern Chile and Argentina.Oversteepening of the edifice by stacking of thick, viscous, dacitic lava flows resulted in collapse of its southeastern flank to form a large volcanic debris avalanche. Biotite ⁴⁰Ar/³⁹Ar dating of lava blocks from the avalanche deposit indicate that collapse occurred at or after 0.15 Ma, and may have been triggered by extrusion of a dacitic flow similar to the one dated at 0.048±0.012 Ma. The avalanche deposits are exceptionally well preserved due to the arid climate, and prominent levées, longitudinal ridges, and megablocks up to 20-m diameter are observed.The avalanche descended 2.8 km vertically, and bifurcated around an older volcano, Cerro Rosado, before debouching onto the salt flats of Salina de Llullaillaco. The north and south limbs of the avalanche traveled 25 and 23 km, respectively, and together cover an area of approximately 165 km². Estimates of deposit volume are hampered by a lack of thickness information except at the edges, but it is likely to be between 1 and 2 km³. Equivalent coefficients of friction of 0.11 and 0.12, and excess travel distances of 20.5 and 18.5 km, are calculated for the north and south limbs, respectively. The avalanche ascended 400 m where it broke against the western flank of Cerro Rosado, and a minimum flow velocity of 90 m s⁻¹ can be calculated at this point; lower velocities of 45 m s⁻¹ are calculated where distal toes ascend 200 m slopes.It is suggested that the remaining precipitous edifice has a high probability for further avalanche collapse in the event of renewed volcanism.     

An integral approach to bedrock river profile analysis

Bedrock river profiles are often interpreted with the aid of slope–area analysis, but noisy topographic data make such interpretations challenging. We present an alternative approach based on an integration of the steady‐state form of the stream power equation. The main component of this approach is a transformation of the horizontal coordinate that converts a steady‐state river profile into a straight line with a slope that is simply related to the ratio of the uplift rate to the erodibility. The transformed profiles, called chi plots, have other useful properties, including co‐linearity of steady‐state tributaries with their main stem and the ease of identifying transient erosional signals. We illustrate these applications with analyses of river profiles extracted from digital topographic datasets. Copyright © 2012 John Wiley & Sons, Ltd.     

Small and large scale fluctuations in atmospheric wind speeds

Atmospheric wind speeds and their fluctuations at different locations (onshore and offshore) are examined. One of the most striking features is the marked intermittency of probability density functions (PDF) of velocity differences, no matter what location is considered. The shape of these PDFs is found to be robust over a wide range of scales which seems to contradict the mathematical concept of stability where a Gaussian distribution should be the limiting one. Motivated by the non-stationarity of atmospheric winds it is shown that the intermittent distributions can be understood as a superposition of different subsets of isotropic turbulence. Thus we suggest a simple stochastic model to reproduce the measured statistics of wind speed fluctuations.     

Causes of the Regional Variability in Observed Sea Level,Sea Surface Temperature and Ocean Colour Over the Period 1993–2011

We analyse the regional variability in observed sea surface height (SSH), sea surface temperature (SST) and ocean colour (OC) from the ESA Climate Change Initiative datasets over the period 1993–2011. The analysis focuses on the signature of the ocean large-scale climate fluctuations driven by the atmospheric forcing and do not address the mesoscale variability. We use the ECCO version 4 ocean reanalysis to unravel the role of ocean transport and surface buoyancy fluxes in the observed SSH, SST and OC variability. We show that the SSH regional variability is dominated by the steric effect (except at high latitude) and is mainly shaped by ocean heat transport divergences with some contributions from the surface heat fluxes forcing that can be significant regionally (confirming earlier results). This is in contrast with the SST regional variability, which is the result of the compensation of surface heat fluxes by ocean heat transport in the mixed layer and arises from small departures around this background balance. Bringing together the results of SSH and SST analyses, we show that SSH and SST bear some common variability. This is because both SSH and SST variability show significant contributions from the surface heat fluxes forcing. It is evidenced by the high correlation between SST and buoyancy-forced SSH almost everywhere in the ocean except at high latitude. OC, which is determined by phytoplankton biomass, is governed by the availability of light and nutrients that essentially depend on climate fluctuations. For this reason, OC shows significant correlation with SST and SSH. We show that the correlation with SST displays the same pattern as the correlation with SSH with a negative correlation in the tropics and subtropics and a positive correlation at high latitude. We discuss the reasons for this pattern.     

Soil water dynamics under different forest vegetation cover: Implications for hillslope stability

Though it is well known that vegetation affects the water balance of soils through canopy interception and evapotranspiration, its hydrological contribution to soil hydrology and stability is yet to be fully quantified. To improve understanding of this hydrological process, soil water dynamics have been monitored at three adjacent hillslopes with different vegetation covers (deciduous tree cover, coniferous tree cover, and grass cover), for nine months from December 2014 to September 2015. The monitored soil moisture values were translated into soil matric suction (SMS) values to facilitate the analysis of hillslope stability. Our observations showed significant seasonal variations in SMS for each vegetation cover condition. However, a significant difference between different vegetation covers was only evident during the winter season where the mean SMS under coniferous tree cover (83.6 kPa) was significantly greater than that under grass cover (41 kPa). The hydrological reinforcing contribution due to matric suction was highest for the hillslope with coniferous tree cover, while the hillslope with deciduous tree cover was second and the hillslope with grass cover was third. The greatest contributions for all cover types were during the summer season. During the winter season, the wettest period of the monitoring study, the additional hydrological reinforcing contributions provided by the deciduous tree cover (1.5 to 6.5 kPa) or the grass cover (0.9 to 5.4 kPa) were insufficient to avoid potential slope failure conditions. However, the additional hydrological reinforcing contribution from the coniferous tree cover (5.8 to 10.4 kPa) was sufficient to provide potentially stable hillslope conditions during the winter season. Our study clearly suggests that during the winter season the hydrological effects from both deciduous tree and grass covers are insufficient to promote slope stability, while the hydrological reinforcing effects from the coniferous tree cover are sufficient even during the winter season. Copyright © 2018 John Wiley & Sons, Ltd.     

Low dispersive modeling of Rayleigh waves on partly staggered grids

In elastic media, finite-difference (FD) implementations of free-surface (FS) boundary conditions on partly staggered grid (PSG) use the highly dispersive vacuum formulation (VPSG). The FS boundary is embedded into a “vacuum” grid layer (null Lame’s constants and negligible density values) where the discretized equations of motion allow computing surface displacements. We place a new set of compound (stress-displacement) nodes along a planar FS and use unilateral mimetic FD discretization of the zero-traction conditions for displacement computation (MPSG). At interior nodes, MPSG reduces to standard VPSG methods and applies fourth-order centered FD along cell diagonals for staggered differentiation combined with nodal second-order FD in time. We perform a dispersion analysis of these methods on a Lamb’s problem and estimate dispersion curves from the phase difference of windowed numerical Rayleigh pulses at two FS receivers. For a given grid sampling criterion (e.g., six or ten nodes per reference S wavelength λ ^S), MPSG dispersion errors are only a quarter of the VPSG method. We also quantify root-mean-square (RMS) misfits of numerical time series relative to analytical waveforms. MPSG RMS misfits barely exceed 10 % when nine nodes sample the minimum S wavelength $\lambda _{\text {MIN}}^{\mathrm {S}}$ in transit (along distances $\sim $ 145 $\lambda _{\text {MIN}}^{\mathrm {S}}$ ). In same tests, VPSG RMS misfits exceed 70 %. We additionally compare MPSG to a consistently fourth-order mimetic method designed on a standard staggered grid. The latter equates the former’s dispersion errors on grids twice denser and shows higher RMS precision only on grids with six or less nodes per $\lambda _{\text {MIN}}^{\mathrm {S}}$ .