A two‐step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis

The Messinian Salinity Crisis is well known to have resulted from a significant drop of the Mediterranean sea level. Considering both onshore and offshore observations, the subsequent reflooding is generally thought to have been very sudden. We present here offshore seismic evidence from the Gulf of Lions and re‐visited onshore data from Italy and Turkey that lead to a new concept of a two‐step reflooding of the Mediterranean Basin after the Messinian Salinity Crisis. The refilling was first moderate and relatively slow accompanied by transgressive ravinement, and later on very rapid, preserving the subaerial Messinian Erosional Surface. The amplitude of these two successive rises of sea level has been estimated at ≤500 m for the first rise and 600–900 m for the second rise. Evaporites from the central Mediterranean basins appear to have been deposited principally at the beginning of the first step of reflooding. After the second step, which preceeded the Zanclean Global Stratotype Section and Point, successive connections with the Paratethyan Dacic Basin, then the Adriatic foredeep, and finally the Euxinian Basin occurred, as a consequence of the continued global rise in sea level. A complex morphology with sills and sub‐basins led to diachronous events such as the so‐called ‘Lago Mare’.This study helps to distinguish events that were synchronous over the entire Mediterranean realm, such as the two‐step reflooding, from those that were more local and diachronous. In addition, the shoreline that marks the transition between these two steps of reflooding in the Provence Basin provides a remarkable palaeogeographical marker for subsidence studies.     

Field study on drainage densities and rescaled width functions in a high‐altitude alpine catchment

This paper addresses the effect of accurately mapping spatially heterogeneous drainage densities in high‐altitude alpine basins on Rescaled Width Functions (RWFs), used in some applications as a minimalist model of the hydrologic response. The channel network and 373 of its channel heads were mapped in the field in a high mountain catchment in the Swiss Alps. The mapped channel network is characterized by highly uneven drainage density, here described by the distribution of the length to the first channelized site computed along steepest descent from any unchannelled site. Various channel networks were extracted from a 1 m lidar‐derived digital terrain model and compared with the field‐mapped channel network using geomorphologic parameters, hillslope‐to‐channel distance and RWFs. Our results show that the channel network derived by statistical analysis of surface morphology is consistent with the field‐mapped network. Larger discrepancies were observed when the channel network was obtained with classical threshold‐based approaches relying on cumulative drainage area and local slope. The actual arrangement of the drainage densities has a significant impact on the RWFs. The discrepancy was largest between RWFs derived from classical extraction methods and RWFs derived with the field‐mapped network, indicating an inappropriate extraction of the channelled portion of the high‐altitude catchment that is a reflection of the variety of channel initiation processes. Our results suggest that spatial heterogeneity of the drainage density might play an important role in modelling streamflow generation. Copyright © 2016 John Wiley & Sons, Ltd.     

Volcanic and structural evolution of Pico do Fogo,Cape Verde

In recent months the media have drawn attention to the Cape Verde archipelago, with particular focus on the island of Fogo, the only island presently active and with an eruption that began on 23 November 2014, finally ceasing on 7 February 2015. The Monte Amarelo conical shield forms most of the 476 km² almost circular island of Fogo. After attaining a critical elevation of about 3500 m, the Monte Amarelo shield volcano was decapitated by a giant landslide that formed a caldera‐like depression (Cha das Caldeiras), which was subsequently partially filled by basaltic nested volcanism. This younger eruptive activity culminated in the construction of the 2829 m‐high Pico do Fogo stratocone, apparently entirely made of layers of basaltic lapilli. Continued growth of the Pico do Fogo summit eruptions was interrupted in 1750, most likely after the stratocone reached a critical height. Since then, at least eight eruptions have taken place inside the landslide depression at the periphery of the Pico do Fogo cone, including the 2014–2015 eruptive event. Strong geological similarities with the Canary Islands, 1400 km to the north, have been frequently noted, probably as a consequence of a common process of origin and evolution associated with a mantle hot‐spot. These similarities are particularly evident when comparing Fogo with the Teide Volcanic Complex on Tenerife, where a lateral collapse of the Las Cañadas stratovolcano also formed a large depression (the Caldera de Las Cañadas), now partially filled with the 3718 m‐high Teide stratocone. However, important geological differences also exist and probably relate to the contrasting evolutionary stages of both islands. The Las Cañadas volcano on Tenerife formed at a late post‐erosional stage, with predominantly evolved (trachyte and phonolite) magmas, while at Fogo basaltic volcanism is still dominant.     

Evaluation of the MODIS snow cover fraction product

Eleven years of daily 500 m gridded Terra Moderate Resolution Imaging Spectroradiometer (MODIS) (MOD10A1) snow cover fraction (SCF) data are evaluated in terms of snow presence detection in Colorado and Washington states. The SCF detection validation study is performed using in‐situ measurements and expressed in terms of snow and land detection and misclassification frequencies. A major aspect addressed in this study involves the shifting of pixel values in time due to sensor viewing angles and gridding artifacts of MODIS sensor products. To account for this error, 500 m gridded pixels are grouped and aggregated to different‐sized areas to incorporate neighboring pixel information. With pixel aggregation, both the probability of detection (POD) and the false alarm ratios increase for almost all cases. Of the false negative (FN) and false positive values (referred to as the total error when combined), FN estimates dominate most of the total error and are greatly reduced with aggregation. The greatest POD increases and total error reductions occur with going from a single 500 m pixel to 3×3‐pixel averaged areas. Since the MODIS SCF algorithm was developed under ideal conditions, SCF detection is also evaluated for varying conditions of vegetation, elevation, cloud cover and air temperature. Finally, using a direct insertion data assimilation approach, pixel averaged MODIS SCF observations are shown to improve modeled snowpack conditions over the single pixel observations due to the smoothing of more error‐prone observations and more accurately snow‐classified pixels. Copyright © 2012 John Wiley & Sons, Ltd.     

Distribution and chronology of Pleistocene permafrost features in France: Database and first results

Numerous periglacial features (polygons, nets, soil stripes, ice‐wedge pseudomorphs and sand‐wedge casts, involutions) have been recorded in France by examining bibliographical sources and aerial photographs. These data show that a large part of France was affected by permafrost during the Pleistocene and only the southern Aquitaine Basin and Languedoc seem to have been beyond its maximum extent. The first OSL ages obtained from the aeolian infill of wedge structures indicate that at least two phases of thermal contraction cracking occurred in southwestern France between ～25 and 36 ka. Chronostratigraphical data from loess in northern France indicate that these episodes correspond to the formation of ice‐wedge networks associated with tundra gleys. In the latter region, two additional permafrost episodes probably occurred during the Last Glacial, the older one corresponding to the end of Marine Isotope Stage (MIS) 4 around 60 ka and the more recent one to MIS 2 around 19–16 ka. Although stratigraphical data indicate that these episodes were relatively short (about one millennium), relict permafrost may have existed for longer periods in northern France.     

DISECA – A Matlab code for dispersive waveform calculations

The propagation velocity of certain elastic waves (so-called surface waves) in vertically inhomogeneous materials is dispersive. The dispersion determination is however frequently ambiguous. To verify the signal analysis reliability and to estimate the dispersive character of the waves an algorithm for dispersive waveform calculation was designed. It is based on the summation of frequency components, with shifts corresponding to the velocity dispersion and distance. The knowledge of physical parameters of the medium is not required. The resulting waveform only contains an individual dispersive wave of the selected mode, thus being particularly suitable for testing of methodologies for dispersive wave analysis.     

TP or not TP,that is the question

We give here a comparative study on the mathematical analysis of two (classes of) discretization schemes for the computation of approximate solutions to incompressible two-phase flow problems in homogeneous porous media. The first scheme is the well-known finite volume scheme with a two-point flux approximation, classically used in industry. The second class contains the so-called approximate gradient schemes, which include finite elements with mass lumping, mixed finite elements, and mimetic finite differences. Both (classes of) schemes are nonconforming and can be expressed using discrete function and gradient reconstructions within a variational formulation. Each class has its specific advantages and drawbacks: monotony properties are natural with the two-point finite volume scheme, but meshes are restricted due to consistency issues; on the contrary, gradient schemes can be used on general meshes, but monotony properties are difficult to obtain.     

Low Tropospheric Layers Over Reunion Island in Lidar-Derived Observations and a High-Resolution Model

In November and December 2008, ground-based mobile lidar (GBML) measurements were carried out on Reunion Island (Indian Ocean, $21^{\circ }07^{\prime }\hbox {S}, 55^{\circ }32^{\prime }\hbox {E}$ , 700 km east of Madagascar) with an ultraviolet (355 nm) aerosol-backscatter lidar. Complex substructures were identified within the planetary boundary layer (PBL). A 500-m-resolution non-hydrostatic model was used to simulate the dynamics of the lower troposphere for two observation periods characteristic of the two main weather regimes in this season: the “trade-wind” regime and the “breeze” regime. The model captured the observed structures with a high degree of realism compared to the GBML. A complete diurnal cycle of the PBL along the south coast of the island during a “trade-wind” day was observed and simulated. The PBL depth was found to be anti-correlated with the wind speed. The model showed that the PBL along the coast behaved as a shallow-water flow in hydraulic theory. As the flow accelerated in response to lateral constriction, conversion of potential into kinetic energy forced the PBL top downwards. This favoured rapid transport of concentrated surface emissions within the contracted surface layer, with a possible impact on air quality. GBML observations were also conducted during the early morning of a “breeze” day on the western slope of the Maïdo mountain (2,200 m), at the top of which a new atmospheric observatory has been in operation since 2012. Both model and GBML revealed two superposed layers. The upper layer, higher than approximately 1,600 m above mean sea level, corresponded to free tropospheric air driven by the trade winds. Below, westerly counterflow advection of humid marine air occurred as a result of wake vortices in the lee of the island. The model suggests that free-tropospheric conditions prevail at the observatory from the second half of the night to mid-morning.