Transport of Snow by the Wind: A Comparison Between Observations in Ad��lie Land, Antarctica, and Simulations Made with the Regional Climate Model MAR

For the first time a simulation of blowing snow events was validated in detail using one-month long observations (January 2010) made in Adélie Land, Antarctica. A regional climate model featuring a coupled atmosphere/blowing snow/snowpack model is forced laterally by meteorological re-analyses. The vertical grid spacing was 2 m from 2 to 20 m above the surface and the horizontal grid spacing was 5?km. The simulation was validated by comparing the occurrence of blowing snow events and other meteorological parameters at two automatic weather stations. The Nash test allowed us to compute efficiencies of the simulation. The regional climate model simulated the observed wind speed with a positive efficiency (0.69). Wind speeds higher than 12 m s ^?1 were underestimated. Positive efficiency of the simulated wind speed was a prerequisite for validating the blowing snow model. Temperatures were simulated with a slightly negative efficiency (?0.16) due to overestimation of the amplitude of the diurnal cycle during one week, probably because the cloud cover was underestimated at that location during the period concerned. Snowfall events were correctly simulated by our model, as confirmed by field reports. Because observations suggested that our instrument (an acoustic sounder) tends to overestimate the blowing snow flux, data were not sufficiently accurate to allow the complete validation of snow drift values. However, the simulation of blowing snow occurrence was in good agreement with the observations made during the first 20 days of January 2010, despite the fact that the blowing snow flux may be underestimated by the regional climate model during pure blowing snow events. We found that blowing snow occurs in Adélie Land only when the 30-min wind speed value at 2 m a.g.l. is >10 m s ^?1. The validation for the last 10 days of January 2010 was less satisfactory because of complications introduced by surface melting and refreezing.     

Identification and Explanation of a Change in the Groundwater Regime using Time Series Analysis

Time series analysis is applied to identify and analyze a transition in the groundwater regime in the aquifer below the sand ridge of Salland in the Netherlands, where groundwater regime refers to the range of head variations throughout the seasons. Standard time series analysis revealed a discrepancy between modeled and observed heads in several piezometers indicating a possible change in the groundwater regime. A new time series modeling approach is developed to simulate the transition from the initial regime to the altered regime. The transition is modeled as a weighted sum of two responses, one representing the initial state of the system, the other representing the altered state. The inferred timing and magnitude of the change provided strong evidence that the transition was the result of significant dredging works that increased the river bed conductance of the main river draining the aquifer. The plausibility of this explanation is corroborated by an analytical model. This case study and the developed approach to identify a change in the groundwater regime are meant to stimulate a more systematic application of time series analysis to detect and understand changes in groundwater systems which may easily go unnoticed in groundwater flow modeling.     

Temporary and definitive fixation of atmospheric lead in deep-sea sediments of the Western Mediterranean Sea

Most lead brought to the Mediterranean Sea has an anthropogenic origin and is mainly transported through the atmosphere. Atmospheric Pb was continuously collected at Cap Ferrat in 1986 and 1987. From this study, the estimation of the anthropogenic Pb flux on the whole Western Mediterranean was, averaged on 1986 and 1987 data, 4080 t. Assuming that the atmospheric anthropogenic Pb input varied in this course of time similarly to the consumption of Pb added to gasolines in France, the mean annual flux could be calculated: 3.95 kg km⁻² yr⁻¹, that is an annual input of 3360 t yr⁻¹. Reaching the sea, this metal seems to become rapidly bound to phytoplankton. Grazing by zooplankton leads to the production of faecal pellets which frequently contain rather high metal concentrations. The sinking rate of pellets of various zooplankton species is high; within a few days pellets may reach deep-sea sediments. After deposition, Pb is released from this organic-rich material during early diagenesis. In most cases, it, therefore, returns to the overlaying water body by ascending diffusion. But, in a deep-sea area of approximately 80 000 km² where Mn oxide precipitation occurs in surficial sediments, Pb seems to remain stored by coprecipitation processes. By considering the lead stored in ‘excess’ in the surficial sediment of the deep-sea area, we estimate that a mean annual anthropogenic Pb amount ranging from 800 up to 1080 t was stored every year from 1950. On the same area, taking into account the Pb loss at the straits, the ‘direct’ atmospheric input to the sea bottom is, on average, 184 t yr⁻¹. The remaining part, that is (800–1080)−184=(616–896) t yr⁻¹, corresponds to an additional ‘indirect’ Pb flux in water due to Pb released from sediments of the surrounding areas where it does not remain stored.     

Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods

Abstract

This paper compares the performance of three geostatistical algorithms, which integrate elevation as an auxiliary variable: kriging with external drift (KED); kriging combined with regression, called regression kriging (RK) or kriging after detrending; and co-kriging (CK). These three methods differ by the way by in which the secondary information is introduced into the prediction procedure. They are applied to improve the prediction of the monthly average rainfall observations measured at 106 climatic stations in Tunisia over an area of 164 150 km² using the elevation as the auxiliary variable. The experimental sample semivariograms, residual semivariograms and cross-variograms are constructed and fitted to estimate the rainfall levels and the estimation variance at the nodes of a square grid of 20 km?×?20 km resolution and to develop corresponding contour maps. Contour diagrams for KED and RK were similar and exhibited a pattern corresponding more closely to local topographic features when (a) the network is sparse and (b) the rainfall–elevation correlation is poor, while CK showed a smooth zonal pattern. Smaller prediction variances are obtained for the RK algorithm. The cross-validation showed that the RMSE obtained for CK gave better results than for KED or RK.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Feki, H., Slimani, M., and Cudennec, C., 2012. Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods. Hydrological Sciences Journal, 57 (7), 1294–1314.     

Panta Rhei 2013–2015: global perspectives on hydrology,society and change

ABSTRACT

In 2013, the International Association of Hydrological Sciences (IAHS) launched the hydrological decade 2013–2022 with the theme “Panta Rhei: Change in Hydrology and Society”. The decade recognizes the urgency of hydrological research to understand and predict the interactions of society and water, to support sustainable water resource use under changing climatic and environmental conditions. This paper reports on the first Panta Rhei biennium 2013–2015, providing a comprehensive resource that describes the scope and direction of Panta Rhei. We bring together the knowledge of all the Panta Rhei working groups, to summarize the most pressing research questions and how the hydrological community is progressing towards those goals. We draw out interconnections between different strands of research, and reflect on the need to take a global view on hydrology in the current era of human impacts and environmental change. Finally, we look back to the six driving science questions identified at the outset of Panta Rhei, to quantify progress towards those aims.

Editor D. Koutsoyiannis; Associate editor not assigned     

In situ oxygen uptake rates by coastal sediments under the influence of the Rhône River (NW Mediterranean Sea)

The influence of riverine inputs on biogeochemical cycling and organic matter recycling in sediments on the continental shelf off the Rhône River mouth (NW Mediterranean Sea) was investigated by measuring sediment oxygen uptake rates using a combination of in situ and laboratory techniques. Four stations were investigated during two cruises in June 2001 and June 2002, with depths ranging from 9 to 192 m and over a distance to the Rhône River mouth ranging from 4 to 36 km. Diffusive oxygen uptake (DOU) rates were determined using an in situ sediment microprofiler and total oxygen uptake (TOU) rates were measured using sediment core incubations. There was good agreement between these two techniques which indicates that the non-diffusive fraction of the oxygen flux was minimal at the investigated stations. DOU rates ranged from 3.7±0.4 mmol O₂ m⁻² d⁻¹ at the continental shelf break to 19.3±0.5 mmol O₂ m⁻² d⁻¹ in front of the Rhône River mouth. Sediment oxygen uptake rates mostly decreased with increasing depth and with distance from the Rhône mouth. The highest oxygen uptake rate was observed at 63 m on the Rhône prodelta, corresponding to intense remineralization of organic matter. This oxygen uptake rate was much larger than expected for the increasing bathymetry, which indicates that biogeochemical cycles and benthic deposition are largely influenced by the Rhône River inputs. This functioning was also supported by the detailed spatial distribution of total organic carbon (TOC), total nitrogen (TN) and C/N atomic ratio in surficial sediments. Sediments of the Rhône prodelta are enriched in organic carbon (2–2.2%) relative to the continental shelf sediments (<1%) and showed C/N ratios exceeding Redfield stoichiometry for fresh marine organic matter. A positive exponential correlation was found between DOU and TOC contents (r²=0.98, n=4). South-westward of the Rhône River mouth, sediments contained highly degraded organic matter of both terrestrial and marine origin, due to direct inputs from the Rhône River, sedimentation of marine organic matter and organic material redeposition after resuspension events.     

ESR/U-series chronology of the Lower Palaeolithic palaeoanthropological site of Visogliano, Trieste, Italy

The Visogliano shelter, in north-eastern Italy, is an important Middle Pleistocene occupation site where human remains were found together with an archaic lithic industry, including choppers, chopping tools and a few protobifaces. It is of utmost importance to try to document this period, when a second wave of settlement colonised Western Europe, carrying new flaking techniques and tools.Combined ESR/U-series analyses, integrated with biostratigraphical and environmental data, define a chronological frame for the layers from which the artefacts were unearthed. The lower levels, including human remains, can be dated to the 350–500 kyr time span, in agreement with micromammal and stratigraphical studies.These data make Visogliano one of the oldest palaeoanthropological sites in Italy, where human remains are directly associated with protobifaces, choppers and chopping tools. In Western Europe, Visogliano is contemporaneous to the G soil of the Arago Cave, France, with which it shares several similarities in faunal assemblages and radiometric data, and which contains human remains also. These data make Visogliano as one of the oldest sites in Europe where the Acheulian culture is observed.     

Diagenesis and the reconstruction of paleoenvironments: A method to restore original δO values of carbonate and phosphate from fossil tooth enamel

Intra-tooth δ¹⁸O variations within the carbonate (δ¹⁸Oc) and phosphate (δ¹⁸Op) components of tooth apatite were measured for Miocene and Pliocene hypsodont mammals from Afghanistan, Greece and Chad in order to evaluate the resistance of enamel to diagenetic alteration. Application of water-apatite interaction models suggest that the different kinetic behaviours of the phosphate-water and carbonate-water systems can be used to detect subtle oxygen isotope disequilibria in fossil enamel when intra-individual variations are considered. Selective alteration of the oxygen isotope composition from the carbonate component of Afghan and Greek enamels suggests inorganic isotopic exchange processes. Microbially-induced isotopic exchange for phosphate is demonstrated for the first time in enamel samples from Chad, in association with extensive recrystallization. In Chad, δ¹⁸Op values were derived from partial isotopic exchange with fossil groundwater during early diagenesis. Mass balance calculations using average carbonate content in enamel as a proxy for recrystallization, and the lowest δ¹⁸Op value of dentine as a proxy for the isotopic composition of the diagenetic fluid, indicate that diagenesis can alter δ¹⁸Op by as much as 3‰ in some enamel samples. This diagenetic alteration is also responsible for a decrease in intra-individual variations of up to 1‰ in affected specimens. The effects of diagenesis on δ¹⁸Op values of fossil enamel are not systematic, however, and can only be estimated if sequential δ¹⁸Op and δ¹⁸Oc analyses are performed on fossil enamel and dentine. Reconstruction of large temporal- or spatial-scale paleoclimates based on δ¹⁸Op analyses from mammalian teeth cannot be considered valid if enamel has been affected by bacterial activity or if the data cannot be corrected for diagenetic effects.     

Projected changes in the tropical Pacific Ocean of importance to tuna fisheries

Future physical and chemical changes to the ocean are likely to significantly affect the distribution and productivity of many marine species. Tuna are of particular importance in the tropical Pacific, as they contribute significantly to the livelihoods, food and economic security of island states. Changes in water properties and circulation will impact on tuna larval dispersal, preferred habitat distributions and the trophic systems that support tuna populations throughout the region. Using recent observations and ocean projections from the CMIP3 and preliminary results from CMIP5 climate models, we document the projected changes to ocean temperature, salinity, stratification and circulation most relevant to distributions of tuna. Under a business-as-usual emission scenario, projections indicate a surface intensified warming in the upper 400 m and a large expansion of the western Pacific Warm Pool, with most surface waters of the central and western equatorial Pacific reaching temperatures warmer than 29 °C by 2100. These changes are likely to alter the preferred habitat of tuna, based on present-day thermal tolerances, and in turn the distribution of spawning and foraging grounds. Large-scale shoaling of the mixed layer and increases in stratification are expected to impact nutrient provision to the biologically active layer, with flow-on trophic effects on the micronekton. Several oceanic currents are projected to change, including a strengthened upper equatorial undercurrent, which could modify the supply of bioavailable iron to the eastern Pacific.