Current state and trends in Canadian Arctic marine ecosystems: II. Heterotrophic food web, pelagic-benthic coupling, and biodiversity

As part of the Canadian contribution to the International Polar Year (IPY), several major international research programs have focused on offshore arctic marine ecosystems. The general goal of these projects was to improve our understanding of how the response of arctic marine ecosystems to climate warming will alter food web structure and ecosystem services provided to Northerners. At least four key findings from these projects relating to arctic heterotrophic food web, pelagic-benthic coupling and biodiversity have emerged: (1) Contrary to a long-standing paradigm of dormant ecosystems during the long arctic winter, major food web components showed relatively high level of winter activity, well before the spring release of ice algae and subsequent phytoplankton bloom. Such phenological plasticity among key secondary producers like zooplankton may thus narrow the risks of extreme mismatch between primary production and secondary production in an increasingly variable arctic environment. (2) Tight pelagic-benthic coupling and consequent recycling of nutrients at the seafloor characterize specific regions of the Canadian Arctic, such as the North Water polynya and Lancaster Sound. The latter constitute hot spots of benthic ecosystem functioning compared to regions where zooplankton-mediated processes weaken the pelagic-benthic coupling. (3) In contrast with another widely shared assumption of lower biodiversity, arctic marine biodiversity is comparable to that reported off Atlantic and Pacific coasts of Canada, albeit threatened by the potential colonization of subarctic species. (4) The rapid decrease of summer sea-ice cover allows increasing numbers of killer whales to use the Canadian High Arctic as a hunting ground. The stronger presence of this species, bound to become a new apex predator of arctic seas, will likely affect populations of endemic arctic marine mammals such as the narwhal, bowhead, and beluga whales.     

A procedure for automated quality control and homogenization of historical daily temperature and precipitation data (APACH): part 1: quality control and application to the Argentine weather service stations

The present paper describes the quality-control component of an automatic procedure (APACH: A Procedure for Automated Quality Control and Homogenization of Weather Station Data) developed to control quality and homogenize the historical daily temperature and precipitation data from meteorological stations. The quality-control method is based on a set of decision-tree algorithms analyzing separately precipitation and minimum and maximum temperature. All our tests are non-parametric and therefore are potentially useful in regions or countries presenting different climates as those observed in Argentina. The method is applied to the 1959–2005 historical daily database of the Argentine National Weather Service. Our results are coherent with the history of the Weather Service and more specifically with the history of implementation of systematized quality control processes. In temperature, our method detects a larger number of suspect values before 1967 (when there was no quality control) and after 1997 (when only real-time quality control had been applied). In precipitation, the detection of error in extreme precipitations is complex, but our method clearly detected a strong decrease in the number of potential outliers after 1976 when the National Weather Service was militarized, and the network was strongly reduced, focusing more on airport weather stations. Also in precipitation, we analyze in detail the long dry sequences and are able to identify potential long erroneous sequences. This is important for the use of the data for hydrological or agricultural impact studies. Finally, all the data are flagged with codes representing the path followed by the record in our decision-tree algorithms. While each code is associated to one of the categories (“Useful”, “Need-Check”, “Doubtful” or “Suspect”), the final user is free to redefine such category-assignment.     

Dynamics of sea-surface temperature anomalies in the Southern Ocean diagnosed from a 2D mixed-layer model

We analyze the processes responsible for the generation and evolution of sea-surface temperature anomalies observed in the Southern Ocean during a decade based on a 2D diagnostic mixed-layer model in which geostrophic advection is prescribed from altimetry. Anomalous air–sea heat flux is the dominant term of the heat budget over most of the domain, while anomalous Ekman heat fluxes account for 20–40% of the variance in the latitude band 40°?60°S. In the ACC pathway, lateral fluxes of heat associated with anomalous geostrophic currents are a major contributor, dominating downstream of several topographic features, reflecting the influence of eddies and frontal migrations. A significant fraction of the variability of large-scale SST anomalies is correlated with either ENSO or the SAM, each mode contributing roughly equally. The relation between the heat budget terms and these climate modes is investigated, showing in particular that anomalous Ekman and air–sea heat fluxes have a co-operating effect (with regional exceptions), hence the large SST response associated with each mode. It is further shown that ENSO- or SAM-locked anomalous geostrophic currents generate substantial heat fluxes in all three basins with magnitude comparable with that of atmospheric forcings for ENSO, and smaller for the SAM except for limited areas. ENSO-locked forcings generate SST anomalies along the ACC pathway, and advection by mean flows is found to be a non-negligible contribution to the heat budget, exhibiting a wavenumber two zonal structure, characteristic of the Antarctic Circumpolar Wave. By contrast SAM-related forcings are predominantly zonally uniform along the ACC, hence smaller zonal SST gradients and a lesser role of mean advection, except in the SouthWest Atlantic. While modeled SST anomalies are significantly correlated with observations over most of the Southern Ocean, the analysis of the data-model discrepancies suggests that vertical ocean physics may play a significant role in the nonseasonal heat budget, especially in some key regions for mode water formation.     

Twentieth century ENSO characteristics in the IPCC database

In this paper, we assess and compare to observations the spatial characteristics of the twentieth Century ENSO SST variability simulated by 23 models of the IPCC-AR4/CMIP3 database. The analysis is confined to the SST anomalies along the equatorial Pacific and is based on the use of a non-linear neural classification algorithm, the Self-Organizing Maps. Systematic biases include a larger than observed proportion for modelled ENSO maximum variability occurring in the Western Pacific. No clear relationship is found between this bias and the characteristics of the modelled mean state bias in the equatorial Pacific. This bias is mainly related to a misrepresentation of both El Niño and La Niña termination phases for most of the models. In contrast, the onset phase is quite well simulated. Modelled El Niño and La Niña peak phases display an asymmetric bias. Whereas the main bias of the modelled El Niño peak is to exhibit a maximum in the western Pacific, the simulated La Niña bias mainly occurs in the central Pacific. In addition, some models are able to capture the observed El Niño peak characteristics while none of them realistically simulate La Niña peaks. It also arises that the models closest to the observations score unevenly in reproducing the different phases, preventing an accurate classification of the models quality to reproduce the overall ENSO-like variability.     

Nuclear modification factor for quarkonium production in high energy proton-nucleus collisions within the parton saturation approach

In this work, we analyze nuclear effects in the inclusive production of quarkonium in proton-nucleus collisions at high energy regime. A theoretical framework that includes initial state effects like nuclear shadowing and gluon density saturation is considered. Numerical results for nuclear modification factor, $R_{pA}(y)$, as a function of meson rapidity in proton-nucleus collisions are presented. The parameter-free predictions are compared to the available data from the Large Hadron Collider. Discussion on the main theoretical uncertainties is made, with emphasis on the phenomenological models for the nuclear saturation scale.     

A parametric study of stability and resonances around L 4 in the elliptic restricted three-body problem

The size of the stable region around the Lagrangian point L ₄ in the elliptic restricted three-body problem is determined by numerical integration as a function of the mass parameter and eccentricity of the primaries. The size distribution of the stable regions in the mass parameter-eccentricity plane shows minima at certain places that are identified with resonances between the librational frequencies of motions around L ₄. These are computed from an approximate analytical equation of Rabe relating the frequency, mass parameter and eccentricity. Solutions of this equation are determined numerically and the global behaviour of the frequencies depending on the mass parameter and eccentricity is shown and discussed. The minimum sizes of the stable regions around L ₄ change along the resonances and the relative strength of the resonances is analysed. Applications to possible Trojan exoplanets are indicated. Escape from L ₄ is also investigated.     

Modelling of sediment supply from torrent catchments in the Western Alps using the sediment contributing area (SCA) approach

The production of coarse sediment in mountain landscapes depends mainly on the type and activity of geomorphic processes and topographic and natural conditions (e.g. vegetation cover) of these catchments. The supply of sediment from these slopes to mountain streams and its subsequent transport lead to sediment connectivity, which describes the integrated coupled state of these systems. Studies from the Northern Calcareous Alps show that the size of the sediment contributing area (SCA), a subset of the drainage area that effectively delivers sediment to the channel network, can be used as a predictor of sediment delivery to mountain streams. The SCA concept is delineated on a digital elevation model (DEM) using a set of rules related to the steepness and length of slopes directly adjacent to the channel network, the gradient of the latter and the vegetation cover. The present study investigates the applicability of this concept to the Western Alps to identify geomorphologically active areas and to estimate mean annual sediment yield (SY) in mainly debris-flow-prone catchments. We use a statistical approach that shows a parameter optimisation and a linear regression of SY on SCA extent. We use a dataset of ~25 years of assessed coarse sediment accumulation in 35 sediment retention basins. In the investigated catchments, sediment transport is governed by several factors, mainly by the extent of vegetation-free areas with a minimum slope of 23° that is coupled to the channel network with a very low gradient of the latter. With our improved framework, we can show that the SCA approach can be applied to catchments that are widely distributed, in a large spatial scale (hectare area) and very heterogeneous in their properties. In general, the investigated catchments show high connectivity, resulting in significant correlations between long-term average yield and the size of the SCA.     

‘Box‐Profile’ Ion Implants as Geochemical Reference Materials for Electron Probe Microanalysis and Secondary Ion Mass Spectrometry

Although electron probe microanalysis and secondary ion mass spectrometry are widely used analytical techniques for geochemical and mineralogical applications, metrologically rigorous quantification remains a major challenge for these methods. Secondary ion mass spectrometry (SIMS) in particular is a matrix‐sensitive method, and the use of matrix‐matched reference materials (RMs) is essential to avoid significant analytical bias. A major problem is that the number of available RMs for SIMS is extremely small compared with the needs of analysts. One approach for the production of matrix‐specific RMs is the use of high‐energy ion implantation that introduces a known amount of a selected isotope into a material. We chose the more elaborate way of implanting a so‐called ‘box‐profile’ to generate a quasi‐homogeneous concentration of the implanted isotope in three dimensions, which allows RMs not only to be used for ion beam analysis but also makes them suitable for EPMA. For proof of concept, we used the thoroughly studied mineralogically and chemically ‘simple’ SiO₂ system. We implanted either ⁴⁷Ti or ⁴⁸Ti into synthetic, ultra‐high‐purity silica glass. Several ‘box‐profiles’ with mass fractions between 10 and 1000 μg g^?1 Ti and maximum depths of homogeneous Ti distribution between 200 nm and 3 μm were produced at the Institute of Ion Beam Physics and Materials Research of Helmholtz‐Zentrum Dresden‐Rossendorf. Multiple implantation steps using varying ion energies and ion doses were simulated with Stopping and Range of Ions in Matter (SRIM) software, optimising for the target concentrations, implantation depths and technical limits of the implanter. We characterised several implant test samples having different concentrations and maximum implantation depths by means of SIMS and other analytical techniques. The results show that the implant samples are suitable for use as reference materials for SIMS measurements. The multi‐energy ion implantation technique also appears to be a promising procedure for the production of EPMA‐suitable reference materials.     

Comprehensive study of the drying behavior of Boom clay: Experimental investigation and numerical modeling

This paper presents a thermo‐hydro‐mechanical framework to model the drying behavior of Boom clay. First, the experimental campaign conducted Noémie Prime is briefly presented because it is used to validate the model. The data acquisition and processing is emphasized because of the use of X‐ray microtomography to be able to more accurately compare experimental and numerical strain fields. The different submodels are introduced. Numerical simulations are performed to illustrate the capability of the proposed model to reproduce the observed behavior. Finally, a comprehensive sensitivity study on several key model parameters associated with the water retention curve, and the permeability of the medium, is performed to get a better understanding of the physics behind the coupled model.     

Solid–fluid transition modelling in geomaterials and application to a mudflow interacting with an obstacle

Given the contrasting behaviour observed for geomaterials, for example, during landslides of the flow type, this contribution proposes an original constitutive model, which associates both an elasto‐plastic relation and a Bingham viscous law linked by a mechanical transition criterion. This last is defined as the second‐order work sign for each material point, which is a general criterion for divergence instabilities. Finite element method with Lagrangian integration points is chosen as a framework for implementing the new model because of its well‐known ability to deal with both solid and fluid behaviours in large deformation processes. A first boundary model considering a sample of initially stable soil, a slope and an obstacle is performed. The results show the power of the constitutive model because the consistent evolution of initiation, propagation and arrest of the mudflow is described. A parametric study is led on various plastic and viscous parameters to determine their influence on the flow development and arrest. Finally, forces against the obstacle are compared with good agreement with those of other authors for the same geometry and a pure viscous behaviour. Copyright © 2014 John Wiley & Sons, Ltd.