Soil production and hillslope transport in mid‐latitudes during the last glacial–interglacial cycle: a combined data and modelling approach in northern Ardennes

The relative efficiency of various hillslope processes through Quaternary glacial–interglacial cycles in the mid‐latitudes is not yet well constrained. Based on a unique set of topographic and soil thickness data in the Ardennes (Belgium), we combine the new CLICHE model of climate‐dependent hillslope evolution with an inversion algorithm in order to get deeper insight into the ways and timing of hillslope dynamics under one such climatic cycle. We simulate the evolution of a synthetic hill reproducing the slope, curvature, and contributing area distributions of the hillslopes of a ~ 2500 km² real area under a simple two‐stage 120‐kyr‐long climatic scenario with linear transitions between cold and warm stages. The inversion method samples a misfit function in the model parameter space, based on estimates of the fit of topographic derivative distributions in classes of soil thickness and of the relative frequencies of the predicted soil thickness classes. Though the inversion results show remarkable convergence patterns for most parameters, no unique solution emerges. We obtain five clusters of good fits, whose centroids are taken as acceptable model solutions. Based on the predicted time series of average denudation rate and soil thickness, plus snapshots of the soil distribution at characteristic times, we discuss these solutions and, comparing them with independent data not involved in the misfit function, we identify the most realistic scenario. Beyond providing first‐order estimates of several parameters that compare well with published data, our results show that denudation rates increase dramatically for a short time at both warm–cold and cold–warm transitions, when the mean annual temperature passes through the [0, ?5 °C] range. We also point to the overwhelming importance of solifluction in shaping hillslopes and transporting soil, and the role of depth‐dependent creep (including frost creep) throughout the climatic cycle, whereas the contributions of simple creep and overland flow are minor. Copyright © 2016 John Wiley & Sons, Ltd.     

Path planning in uncertain flow fields using ensemble method

An ensemble-based approach is developed to conduct optimal path planning in unsteady ocean currents under uncertainty. We focus our attention on two-dimensional steady and unsteady uncertain flows, and adopt a sampling methodology that is well suited to operational forecasts, where an ensemble of deterministic predictions is used to model and quantify uncertainty. In an operational setting, much about dynamics, topography, and forcing of the ocean environment is uncertain. To address this uncertainty, the flow field is parametrized using a finite number of independent canonical random variables with known densities, and the ensemble is generated by sampling these variables. For each of the resulting realizations of the uncertain current field, we predict the path that minimizes the travel time by solving a boundary value problem (BVP), based on the Pontryagin maximum principle. A family of backward-in-time trajectories starting at the end position is used to generate suitable initial values for the BVP solver. This allows us to examine and analyze the performance of the sampling strategy and to develop insight into extensions dealing with general circulation ocean models. In particular, the ensemble method enables us to perform a statistical analysis of travel times and consequently develop a path planning approach that accounts for these statistics. The proposed methodology is tested for a number of scenarios. We first validate our algorithms by reproducing simple canonical solutions, and then demonstrate our approach in more complex flow fields, including idealized, steady and unsteady double-gyre flows.     

Variability of nutrients and phytoplankton primary production in a shallow macrotidal coastal ecosystem (Arcachon Bay,France)

Seasonal and spatial variations of phytoplankton primary production were studied using a high frequency sampling strategy in the external (ENW) and internal (INW) part of Arcachon Bay, during 2002 and 2003. In order to better assess the availability of nutrients and their relationship with phytoplankton primary production, nutrient variability was studied in relation to environmental conditions and phytoplankton production. During winter, when primary production rates were the lowest, nutrient concentrations were maximal but did not show excessive levels compared to highly urbanised areas. Seasonal and spatial variations of nutrient concentrations (especially DIN-nitrate + nitrite + ammonium- and Si) were largely influenced by Leyre River loads coupled with high tidal exchange with the Atlantic Ocean creating a nutrient gradient between the INW and ENW. By February, diatom growth leads to an early severe nutrient depletion in the entire bay. Examination of nutrient ratios showed that the potential limiting nutrient during spring was P in 2003, and Si in 2002. During summer 2003, N and Si concentrations reached their lowest values, and nutrient ratios revealed a N-deficient environment, more pronounced in the INW. The high Si:N ratios during this period might be explained by (1) important N-uptake by all autotroph communities and (2) benthic-pelagic coupling with high Si regeneration. This study shows that nutrient levels in Arcachon Bay seem to play an important role in the control of phytoplankton primary production rates during the productive period and explain their spatial, seasonal and inter-annual variability. Our estimates of annual integrated phytoplankton primary production (103 g C m⁻² y⁻¹) place this bay within the low to moderate phytoplankton primary production systems.     

Enhanced seamount location database for the western and central Pacific Ocean: Screening and cross-checking of 20 existing datasets

Seamounts are habitats of considerable interest in terms of conservation and biodiversity, and in terms of fisheries for bentho-pelagic and pelagic species. Twenty previously compiled datasets including seamount/underwater feature lists, bathymetric maps and emerged feature maps from different sources (ship-derived and satellite altimetry-derived) at different spatial scales (from individual cruise to worldwide satellite data) were gathered in order to compile an enhanced list of underwater features for parts of the western and central Pacific Ocean (WCPO). The KL04 dataset [Kitchingman, A., and Lai, S., 2004. Inferences on potential seamount locations from mid-resolution bathymetric data. Fisheries Centre Research Reports 12 (5), 7–12], listing seamount positions and depths as calculated from satellite altimetry-derived bathymetry, provided the baseline data for this study as it covered the entire region of interest and included summit depth information. All KL04 potential seamounts were cross-checked with other datasets to remove any atolls and islands that had been incorrectly classified as seamounts, to add seamounts undetected by KL04, to update the overall database (geolocation, depth, elevation, and name) and to compile a 12-class typology of the different types of underwater features. Of the 4626 potential seamounts identified in KL04, 719 were multiple identifications of the same large underwater features and 373 (10%) were actually emerged banks, atolls and islands, leaving 3534 actual underwater features. Conversely, 487 underwater features were documented in other datasets but not registered by KL04. The screening of all the potential WCPO seamounts produced a final list of 4021 underwater features with agreed upon position and information. This enhanced list should have many applications in oceanography, biodiversity conservation and studies of the influence of seamounts on pelagic ecosystems and fisheries.     

Downscaling ability of one-way nested regional climate models: the Big-Brother Experiment

A methodology is developed for testing the downscaling ability of nested regional climate models (RCMs). The proposed methodology, nick-named the Big-Brother Experiment (BBE), is based on a "perfect-prognosis" approach and hence does not suffer from model errors nor from limitations in observed climatologies. The BBE consists in first establishing a reference climate by performing a large-domain high-resolution RCM simulation: this simulation is called the Big Brother. This reference simulation is then degraded by filtering short scales that are unresolved in today's global objective analyses (OA) and/or global climate models (GCMs) when integrated for climate projections. This filtered reference is then used to drive the same nested RCM (called the Little Brother), integrated at the same high-resolution as the Big Brother, but over a smaller domain that is embedded in the Big-Brother domain. The climate statistics of the Little Brother are then compared with those of the Big Brother over the Little-Brother domain. Differences can thus be attributed unambiguously to errors associated with the nesting and downscaling technique, and not to model errors nor to observation limitations. The results of the BBE applied to a one-winter-month simulation over eastern North America at 45-km grid-spacing resolution show that the one-way nesting strategy has skill in downscaling large-scale information to the regional scales. The time mean and variability of fine-scale features in a number of fields, such as sea level pressure, 975-hPa temperature and precipitation are successfully reproduced, particularly over regions where small-scale surface forcings are strong. Over other regions such as the ocean and away from the surface, the small-scale reproducibility is more difficult to achieve.     

Using seasonal hindcasts to understand the origin of the equatorial cold tongue bias in CGCMs and its impact on ENSO

The cold equatorial SST bias in the tropical Pacific that is persistent in many coupled OAGCMs severely impacts the fidelity of the simulated climate and variability in this key region, such as the ENSO phenomenon. The classical bias analysis in these models usually concentrates on multi-decadal to centennial time series needed to obtain statistically robust features. Yet, this strategy cannot fully explain how the models errors were generated in the first place. Here, we use seasonal re-forecasts (hindcasts) to track back the origin of this cold bias. As such hindcasts are initialized close to observations, the transient drift leading to the cold bias can be analyzed to distinguish pre-existing errors from errors responding to initial ones. A time sequence of processes involved in the advent of the final mean state errors can then be proposed. We apply this strategy to the ENSEMBLES-FP6 project multi-model hindcasts of the last decades. Four of the five AOGCMs develop a persistent equatorial cold tongue bias within a few months. The associated systematic errors are first assessed separately for the warm and cold ENSO phases. We find that the models are able to reproduce either El Niño or La Niña close to observations, but not both. ENSO composites then show that the spurious equatorial cooling is maximum for El Niño years for the February and August start dates. For these events and at this time of the year, zonal wind errors in the equatorial Pacific are present from the beginning of the simulation and are hypothesized to be at the origin of the equatorial cold bias, generating too strong upwelling conditions. The systematic underestimation of the mixed layer depth in several models can also amplify the growth of the SST bias. The seminal role of these zonal wind errors is further demonstrated by carrying out ocean-only experiments forced by the AOCGCMs daily 10-meter wind. In a case study, we show that for several models, this forcing is sufficient to reproduce the main SST error patterns seen after 1 month in the AOCGCM hindcasts.     

The Multiscale TROPIcal CatchmentS critical zone observatory M-TROPICS dataset II: Land use,hydrology and sediment production monitoring in Houay Pano,northern Lao PDR

Mountain regions of the humid tropics are characterized by steep slopes and heavy rains. These regions are thus prone to both high surface runoff and soil erosion. In Southeast Asia, uplands are also subject to rapid land-use change, predominantly as a result of increased population pressure and market forces. Since 1998, the Houay Pano site, located in northern Lao PDR (19.85°N 102.17°E) within the Mekong basin, aims at assessing the long-term impact of the conversion of traditional slash-and-burn cultivation systems to commercial perennial monocultures such as teak tree plantations, on the catchment hydrological response and sediment yield. The instrumented site monitors hydro-meteorological and soil loss parameters at both microplot (1 m²) and small catchment (0.6 km²) scales. The monitored catchment is part of the network of critical zone observatories named Multiscale TROPIcal CatchmentS (M-TROPICS). The data shared by M-TROPICS in Houay Pano are (1) rainfall, (2) air temperature, air relative humidity, wind speed, and global radiation, (3) catchment land use, (4) stream water level, suspended particulate matter, bed particulate matter and stones, (5) soil surface features, and (6) soil surface runoff and soil detachment. The dataset has already been used to interpret suspended particulate matter and bed particulate matter sources and dynamics, to assess the impact of land-use change on catchment hydrology, soil erosion, and sediment yields, to understand bacteria fate and weed seed transport across the catchment, and to build catchment-scale models focused on hydrology and water quality issues. The dataset may be further used to, for example, assess the role of headwater catchments in large tropical river basin hydrology, support the interpretation of new variables measured in the catchment (e.g., contaminants other than faecal bacteria), and assess the relative impacts of both climate and land-use change on the catchment.     

Ice flow physical processes derived from the ERS-1 high-resolution map of the Antarctica and Greenland ice sheets

The ERS-1 satellite, launched in 1991, has provided altimetric observations of the Greenland Ice Sheet and 80 per cent of the Antarctica Ice Sheet north of 82°S. It was placed in a geodetic (168-day repeat) orbit between April 1994 and March 1995, yielding a 1.5  km across-track spacing at latitude 70° with a higher along-track sampling of 350  m. We have analysed the waveform altimetric data from this period to compute maps with a 1/30° grid size. Data processing consists of correcting for environmental factors and editing and retracking the waveforms. A further step consists of reducing the radial orbit error through crossover analysis and correcting the slope error to second order. The high-resolution topography of both ice sheets reveals numerous details. A kilometre-scale surface roughness running at 45° from the flow direction is the dominant topographic characteristic of both continents. Antarctica also exhibits many scars due to local flow anomalies. Several physical processes can be identified: abrupt transitions from deformation to sliding and vice versa, and impressive strike-slip phenomena, inducing en echelon folds.     

A pilot study of genetic differentiation between two phenotypes of a Mediterranean population of the bivalve Cerastoderma glaucum and genetic discrimination with other Cerastoderma glaucum and Cerastoderma edule populations outside the Mediterranean

The internal transcribed spacer 1 (ITS1) region of 29 individuals of Cerastoderma glaucum (from the Mediterranean Sea, the Atlantic and the North Sea) and 18 individuals of Cerastoderma edule (from the Atlantic and the North Sea) and the cytochrome c oxidase subunit I (COI) region of 43 individuals of C. glaucum (from the Mediterranean Sea, the Atlantic and the North Sea) and 9 individuals of C. edule (from the Atlantic and the North Sea) were PCR amplified and sequenced. The ITS1 sequences led us to describe 27 genetic haplotypes, while 21 genetic haplotypes were determined via the COI sequences. The ITS1 and the COI sequences revealed an important genetic variability within the Mediterranean population of C. glaucum, which contained two different phenotypes: orange foot and yellow foot. A genetic differentiation between the two phenotypes was revealed by the nucleotide diversity index (π) and is strengthened by the AMOVA analysis. This result leads us to suspect the presence of two different groups in the same population of C. glaucum. Nevertheless to be confirmed this hypothesis requires further studies using more locations and a larger number of sample sizes.