ERA-interim reanalysis debiased at FLUXNET sites

(preliminary) Exchanges of carbon, water and energy between the land surface and the atmosphere are monitored by eddy covariance technique at the ecosystem level. Currently, the FLUXNET database contains more than 500 sites registered and up to 250 of them sharing data (Free Fair Use dataset). Many modelling groups use the FLUXNET dataset for evaluating ecosystem model's performances but it requires uninterrupted time series for the meteorological variables used as input. Because original in-situ data often contain gaps, from very short (few hours) up to relatively long (some months), we develop a new and robust method for filling the gaps in meteorological data measured at site level. Our approach has the benefit of making use of continuous data available globally (ERA-interim) and high temporal resolution spanning from 1989 to today. These data are however not measured at site level and for this reason a method to downscale and correct the ERA-interim data is needed. We apply this method on the level 4 data (L4) from the LaThuile collection, freely available after registration under a Fair-Use policy. The performances of the developed method vary across sites and are also function of the meteorological variable. On average overall sites, the bias correction leads to cancel from 10% to 36% of the initial mismatch between in-situ and ERA-interim data, depending of the meteorological variable considered. In comparison to the internal variability of the in-situ data, the root mean square error (RMSE) between the in-situ data and the un-biased ERA-I data remains relatively large (on average overall sites, from 27% to 76% of the standard deviation of in-situ data, depending of the meteorological variable considered). The performance of the method remains low for the Wind Speed field, in particular regarding its capacity to conserve a standard deviation similar to the one measured at FLUXNET stations.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

Near-bottom zooplankton aggregations in Kongsfjorden with link to images from the optical zooplankton sensor MOKI

Near-bottom zooplankton communities have rarely been studied despite numerous reports of high zooplankton concentrations, probably due to methodological constraints. In Kongsfjorden, Svalbard, the near-bottom layer was studied for the first time by combining daytime deployments of a remotely operated vehicle (ROV), the optical zooplankton sensor moored on-sight key species investigation (MOKI), and Tucker trawl sampling. ROV data from the fjord entrance and the inner fjord showed high near-bottom abundances of euphausiids with a mean concentration of 17.3 ± 3.5 n/100 m**3. With the MOKI system, we observed varying numbers of euphausiids, amphipods, chaetognaths, and copepods on the seafloor at six stations. Light-induced zooplankton swarms reached densities in the order of 90,000 (euphausiids), 120,000 (amphipods), and 470,000 ind/m**3 (chaetognaths), whereas older copepodids of Calanus hyperboreus and C. glacialis did not respond to light. They were abundant at the seafloor and 5 m above and showed maximum abundance of 65,000 ind/m**3. Tucker trawl data provided an overview of the seasonal vertical distribution of euphausiids. The most abundant species Thysanoessa inermis reached near-bottom concentrations of 270 ind/m**3. Regional distribution was neither related to depth nor to location in the fjord. The taxa observed were all part of the pelagic community. Our observations suggest the presence of near-bottom macrozooplankton also in other regions and challenge the current view of bentho-pelagic coupling. Neglecting this community may cause severe underestimates of the stock of elagic zooplankton, especially predatory species, which link secondary production with higher trophic levels.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.

G-force measurements at Mount Nyiragongo

To evaluate the mechanical stress on the volcanic edifice that results from lava lake level variations, we deployed a self-recording, differential capacitance (MEMS Inertial Sensor STMicroelectronics LIS3LV02DQ), 3-axis X6-1A accelerometer (Gulf Coast Data Concepts, LLC) at a distance of ~100m from the center of the Nyiragongo lava lake on freshly erupted lava flows. The device range was used in high (12-bit) resolution mode, which corresponds to a sensitivity of about 1 mg. The device was set to high-sensitivity mode with four additional bits to improve resolution, yet with a much lower signal-noise ratio. Once in position, the accelerometer continuously recorded data for three-day periods in June 2010. The system was oriented so that the X- and Y-axes form a plain parallel to the lava lake. During data collection, we did not attempt to calibrate the precision of the angle because relative G-force measurements were required instead of absolute G-force measurements. To distinguish the tiny accelerations caused by temperature differentials of the atmosphere, from the forces caused by magma movements, the temperature of the X6-1A device was continuously recorded. Temperature variations were corrected for by applying a de-correlation method to the recorded signal. Data was collected at 20 Hz, regrouped into batches that cover 1 hour per observation and associated with one averaged temperature measurement. This method was reproducible because diurnal temperature variations were the main cause for heating and cooling.

Predicting changes in abundance and juvenile proportions of earthworms in the field using glmmADMB in R

The study aimed at investigating effects of three differently acting biocides; the insecticide esfenvalerate, the fungicide picoxystrobin and the bactericide triclosan, applied individually and as a mixture, on an earthworm community in the field. A concentration-response design was chosen and results were analyzed using univariate and multivariate approaches. Effects on juvenile proportions were less pronounced and more variable than effects on abundance, but effects in general were species- and chemical-specific, and temporal variations distinct. Esfenvalerate and picoxystrobin appeared to elicit stronger effects than triclosan at laboratory-based ECx values, which is in accordance with our previous laboratory study on Eisenia fetida. The mixture affected abundance and juvenile proportions, but the latter only at high mixture concentrations. Esfenvalerate and picoxystrobin appeared to be the main drivers for the mixture's toxicity. Species-specific toxicity patterns question the reliability of mixture toxicity predictions derived on E. fetida for field earthworms. Biocide concentrations equaling EC50s (reproduction) for E. fetida provoked effects on the field earthworms mainly exceeding 50%, indicating effect intensification from the laboratory to field as well as the influence of indirect effects produced by species interactions. The differing results of the present field study and the previous laboratory study imply that lower- and higher-tier studies may not be mutually exclusive, but to be used in complementary.

Occurrence and Range-through database of functional diversity of marine ecosystems after the Late Permian mass extinction event

The Late Permian mass extinction event about 252 million years ago was the most severe biotic crisis of the past 500 million years and occurred during an episode of global warming. The loss of around two-thirds of marine genera is thought to have had substantial ecological effects, but the overall impacts on the functioning of marine ecosystems and the pattern of marine recovery are uncertain. Here we analyse the fossil occurrences of all known benthic marine invertebrate genera from the Permian and Triassic periods, and assign each to a functional group based on their inferred lifestyle. We show that despite the selective extinction of 62-74% of these genera, all but one functional group persisted through the crisis, indicating that there was no significant loss of functional diversity at the global scale. In addition, only one new mode of life originated in the extinction aftermath. We suggest that Early Triassic marine ecosystems were not as ecologically depauperate as widely assumed. Functional diversity was, however, reduced in particular regions and habitats, such as tropical reefs; at these smaller scales, recovery varied spatially and temporally, probably driven by migration of surviving groups. We find that marine ecosystems did not return to their pre-extinction state, and by the Middle Triassic greater functional evenness is recorded, resulting from the radiation of previously subordinate groups such as motile, epifaunal grazers.