Global surface net-radiation at 5 km from MODIS Terra

Reliable and fine resolution estimates of surface net-radiation are required for estimating latent and sensible heat fluxes between the land surface and the atmosphere. However, currently, fine resolution estimates of net-radiation are not available and consequently it is challenging to develop multi-year estimates of evapotranspiration at scales that can capture land surface heterogeneity and are relevant for policy and decision-making. We developed and evaluated a global net-radiation product at 5 km and 8-day resolution by combining mutually consistent atmosphere and land data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra. Comparison with net-radiation measurements from 154 globally distributed sites (414 site-years) from the FLUXNET and Surface Radiation budget network (SURFRAD) showed that the net-radiation product agreed well with measurements across seasons and climate types in the extratropics (Wilmott’s index ranged from 0.74 for boreal to 0.63 for Mediterranean sites). Mean absolute deviation between the MODIS and measured net-radiation ranged from 38.0 ± 1.8 W∙m−2 in boreal to 72.0 ± 4.1 W∙m−2 in the tropical climates. The mean bias was small and constituted only 11%, 0.7%, 8.4%, 4.2%, 13.3%, and 5.4% of the mean absolute error in daytime net-radiation in boreal, Mediterranean, temperate-continental, temperate, semi-arid, and tropical climate, respectively. To assess the accuracy of the broader spatiotemporal patterns, we upscaled error-quantified MODIS net-radiation and compared it with the net-radiation estimates from the coarse spatial (1° × 1°) but high temporal resolution gridded net-radiation product from the Clouds and Earth’s Radiant Energy System (CERES). Our estimates agreed closely with the net-radiation estimates from the CERES. Difference between the two was less than 10 W•m−2 in 94% of the total land area. MODIS net-radiation product will be a valuable resource for the science community studying turbulent fluxes and energy budget at the Earth’s surface.

Monitoring of active layer dynamics at two transects on Svalbard using multi-channel ground-penetrating radar

Multi-channel ground-penetrating radar is used to investigate the late-summer evolution of the thaw depth and the average soil water content of the thawed active layer at a high-arctic continuous permafrost site on Svalbard, Norway. Between mid of August and mid of September 2008, five surveys have been conducted over transect lengths of 130 and 175 m each. The maximum thaw depths range from 1.6 m to 2.0 m, so that they are among the deepest thaw depths recorded for Svalbard so far. The thaw depths increase by approximately 0.2 m between mid of August and beginning of September and subsequently remain constant until mid of September. The thaw rates are approximately constant over the entire length of the transects within the measurement accuracy of about 5 to 10 cm. The average volumetric soil water content of the thawed soil varies between 0.18 and 0.27 along the investigated transects. While the measurements do not show significant changes in soil water content over the first four weeks of the study, strong precipitation causes an increase in average soil water content of up to 0.04 during the last week. These values are in good agreement with evapotranspiration and precipitation rates measured in the vicinity of the the study site. While we cannot provide conclusive reasons for the detected spatial variability of the thaw depth at the study site, our measurements show that thaw depth and average soil water content are not directly correlated. The study demonstrates the potential of multi-channel ground-penetrating radar for mapping thaw depth in permafrost areas. The novel non-invasive technique is particularly useful when the thaw depth exceeds 1.5 m, so that it is hardly accessible by manual probing. In addition, multi-channel ground-penetrating radar holds potential for mapping the latent heat content of the active layer and for estimating weekly to monthly averages of the ground heat flux during the thaw period.

Land cover classification of Samoylov Island and Landsat subpixel water cover of Lena River Delta, Siberia, with links to ESRI grid files

Ignoring small-scale heterogeneities in Arctic land cover may bias estimates of water, heat and carbon fluxes in large-scale climate and ecosystem models. We investigated subpixel-scale heterogeneity in CHRIS/PROBA and Landsat-7 ETM+ satellite imagery over ice-wedge polygonal tundra in the Lena Delta of Siberia, and the associated implications for evapotranspiration (ET) estimation. Field measurements were combined with aerial and satellite data to link fine-scale (0.3 m resolution) with coarse-scale (upto 30 m resolution) land cover data. A large portion of the total wet tundra (80%) and water body area (30%) appeared in the form of patches less than 0.1 ha in size, which could not be resolved with satellite data. Wet tundra and small water bodies represented about half of the total ET in summer. Their contribution was reduced to 20% in fall, during which ET rates from dry tundra were highest instead. Inclusion of subpixel-scale water bodies increased the total water surface area of the Lena Delta from 13% to 20%. The actual land/water proportions within each composite satellite pixel was best captured with Landsat data using a statistical downscaling approach, which is recommended for reliable large-scale modelling of water, heat and carbon exchange from permafrost landscapes.

Condensation versus air density change as a major cause of airflow: experimental evidence, link to data files and videos

The dominant model of atmospheric circulation posits that hot air rises, creating horizontal winds. A second major driver has recently been proposed by Makarieva and Gorshkov in their biotic pump theory (BPT), which suggests that evapotranspiration from natural closed-canopy forests causes intense condensation, and hence winds from ocean to land. Critics of the BPT argue that air movement to fill the partial vacuum caused by condensation is always isotropic, and therefore causes no net air movement (Bunyard, 2015, hdl:11232/397). This paper explores the physics of water condensation under mild atmospheric conditions, within a purpose-designed square-section 4.8 m-tall closed-system structure. Two enclosed vertical columns are connected at top and bottom by two horizontal tunnels, around which 19.5 m**3 of atmospheric air can circulate freely, allowing rotary airflows in either direction. This air can be cooled and/or warmed by refrigeration pipes and a heating mat, and changes in airflow, temperature, humidity and barometric pressure measured in real time. The study investigates whether the "hot-air-rises" or an implosive condensation model can better explain the results of more than 100 experiments. The data show a highly significant correlation (R2 >0.96, p value <0.001) between observed airflows and partial pressure changes from condensation. While the kinetic energy of the refrigerated air falls short of that required in bringing about observed airflows by a factor of at least 30, less than a tenth of the potential kinetic energy from condensation is shown to be sufficient. The assumption that condensation of water vapour is always isotropic is therefore incorrect. Condensation can be anisotropic, and in the laboratory does cause sustained airflow.

Hydrogen isotope compositions of terrestrial plant-wax n-alkanes analysed on a transect of marine surface sediment samples from 1°N to 28°S off southwest Africa

A transect of marine surface sediment samples from 1° N to 28° S off southwest Africa was analysed to verify the application of hydrogen isotope compositions of terrestrial plant-wax n-alkanes preserved in ocean sediments as a proxy for continental hydrological conditions. Conditions on the adjacent continent range from humid evergreen forests to deciduous forests, wood- and shrub land and further to arid grasslands and deserts. The hydrogen isotope values for the dominant n-alkane homologues (C29, C31 and C33) vary from -123 per mil to -141 per mil VSMOW and correlate with the modelled hydrogen isotope composition of mean annual and growing season precipitation of postulated continental source areas (r up to 0.8, p < 0.01). The apparent hydrogen isotope fractionation between alkanes and mean annual precipitation is remarkably uniform (-109 per mil on average, Sigma <= 5 per mil, n = 27). Potentially, effects of aridity on the apparent hydrogen isotope fractionation are concealed by the contribution of different plants (C3 dicotyledons vs C4 grasses). Thus, isotope ratios of leaf wax n-alkanes preserved in ocean margin sediments in these and similar tropical regions may be directly converted to dD ratios of ancient precipitation by employing a constant hydrogen isotope fractionation.

Younger Dryas-Holocene sedimentary plant-wax hydrogen isotopes from the southern European Alps

We present a Younger Dryas-Holocene record of the hydrogen isotopic composition of sedimentary plant waxes (dDwax) from the southern European Alps (Lake Ghirla, N-Italy) to investigate its sensitivity to climatic forcing variations in this mid-latitude region (45°N). A modern altitudinal transect of dD values of river water and leaf waxes in the Lake Ghirla catchment is used to test present-day climate sensitivity of dDwax. While we find that altitudinal effects on dDwax are minor at our study site, temperature, precipitation amount, and evapotranspiration all appear to influence dDwax to varying extents. In the lake-sediment record, dDwax values vary between -134 and -180 per mil over the past 13 kyr. The long-term Holocene pattern of dDwax parallels the trend of decreasing temperature and is thus likely forced by the decline of northern hemisphere summer insolation. Shorter-term fluctuations, in contrast, may reflect both temperature and moisture-source changes. During the cool Younger Dryas and Little Ice Age (LIA) periods we observe unexpectedly high dDwax values relative to those before and after. We suggest that a change towards a more D-enriched moisture source is required during these intervals. In fact, a shift from northern N-Atlantic to southern N-Atlantic/western Mediterranean Sea sources would be consistent with a southward migration of the Westerlies with climate cooling. Prominent dDwax fluctuations in the early and middle Holocene are negative and potentially associated with temperature declines. In the late Holocene (<4 kyr BP), excursions are partly positive (as for the LIA) suggesting a stronger influence of moisture-source changes on dDwax variation. In addition to isotopic fractionations of the hydrological cycle, changes in vegetation composition, in the length of the growing season, and in snowfall amount provide additional potential sources of variability, although we cannot yet quantitatively assess these in the paleo-record. We conclude that while our dDwax record from the Alps does contain climatic information, it is a complicated record that would require additional constraints to be robustly interpreted. This also has important implications for other water-isotope-based proxy records of precipitation and hydro-climate from this region, such as cave speleothems.

Modélisation de l'évolution hydroclimatique des flux et stocks d'eau verte et d'eau bleue du bassin versant de la Garonne

La gestion intégrée de la ressource en eau implique de distinguer les parcours de l’eau qui sont accessibles aux sociétés de ceux qui ne le sont pas. Les cheminements de l’eau sont nombreux et fortement variables d’un lieu à l’autre. Il est possible de simplifier cette question en s’attardant plutôt aux deux destinations de l’eau. L’eau bleue forme les réserves et les flux dans l’hydrosystème : cours d’eau, nappes et écoulements souterrains. L’eau verte est le flux invisible de vapeur d’eau qui rejoint l’atmosphère. Elle inclut l’eau consommée par les plantes et l’eau dans les sols. Or, un grand nombre d’études ne portent que sur un seul type d’eau bleue, en ne s’intéressant généralement qu’au devenir des débits ou, plus rarement, à la recharge des nappes. Le portrait global est alors manquant. Dans un même temps, les changements climatiques viennent impacter ce cheminement de l’eau en faisant varier de manière distincte les différents composants de cycle hydrologique. L’étude réalisée ici utilise l’outil de modélisation SWAT afin de réaliser le suivi de toutes les composantes du cycle hydrologique et de quantifier l’impact des changements climatiques sur l’hydrosystème du bassin versant de la Garonne. Une première partie du travail a permis d’affiner la mise en place du modèle pour répondre au mieux à la problématique posée. Un soin particulier a été apporté à l’utilisation de données météorologiques sur grille (SAFRAN) ainsi qu’à la prise en compte de la neige sur les reliefs. Le calage des paramètres du modèle a été testé dans un contexte differential split sampling, en calant puis validant sur des années contrastées en terme climatique afin d’appréhender la robustesse de la simulation dans un contexte de changements climatiques. Cette étape a permis une amélioration substantielle des performances sur la période de calage (2000-2010) ainsi que la mise en évidence de la stabilité du modèle face aux changements climatiques. Par suite, des simulations sur une période d’un siècle (1960-2050) ont été produites puis analysées en deux phases : i) La période passée (1960-2000), basée sur les observations climatiques, a servi de période de validation à long terme du modèle sur la simulation des débits, avec de très bonnes performances. L’analyse des différents composants hydrologiques met en évidence un impact fort sur les flux et stocks d’eau verte, avec une diminution de la teneur en eau des sols et une augmentation importante de l’évapotranspiration. Les composantes de l’eau bleue sont principalement perturbées au niveau du stock de neige et des débits qui présentent tous les deux une baisse substantielle. ii) Des projections hydrologiques ont été réalisées (2010-2050) en sélectionnant une gamme de scénarios et de modèles climatiques issus d’une mise à l’échelle dynamique. L’analyse de simulation vient en bonne part confirmer les conclusions tirées de la période passée : un impact important sur l’eau verte, avec toujours une baisse de la teneur en eau des sols et une augmentation de l’évapotranspiration potentielle. Les simulations montrent que la teneur en eau des sols pendant la période estivale est telle qu’elle en vient à réduire les flux d’évapotranspiration réelle, mettant en évidence le possible déficit futur des stocks d’eau verte. En outre, si l’analyse des composantes de l’eau bleue montre toujours une diminution significative du stock de neige, les débits semblent cette fois en hausse pendant l’automne et l’hiver. Ces résultats sont un signe de l’«accélération» des composantes d’eau bleue de surface, probablement en relation avec l’augmentation des évènements extrêmes de précipitation. Ce travail a permis de réaliser une analyse des variations de la plupart des composantes du cycle hydrologique à l’échelle d’un bassin versant, confirmant l’importance de prendre en compte toutes ces composantes pour évaluer l’impact des changements climatiques et plus largement des changements environnementaux sur la ressource en eau.

Rainfall infiltration and soil hydrological characteristics below ancient forest, planted forest and grassland in a temperate northern climate

How rainfall infiltration rate and soil hydrological characteristics develop over time under forests of different ages in temperate regions is poorly understood. In this study, infiltration rate and soil hydrological characteristics were investigated under forests of different ages and under grassland. Soil hydraulic characteristics were measured at different scales under a 250-year-old grazed grassland (GL), 6-year-old (6yr) and 48-year-old (48yr) Scots pine (Pinus sylvestris) plantations, remnant 300-year-old individual Scots pine (OT) and a 4000-year-old Caledonian Forest (AF). In situ field-saturated hydraulic conductivity (Kfs) was measured, and visible root:soil area was estimated from soil pits. Macroporosity, pore structure and macropore connectivity were estimated from X-ray tomography of soil cores, and from water-release characteristics. At all scales, the median values for Kfs, root fraction, macroporosity and connectivity values tended to AF>OT>48yr>GL>6yr, indicating that infiltration rates and water storage increased with forest age. The remnant Caledonian Forest had a huge range of Kfs (12 to >4922mmh-1), with maximum Kfs values 7 to 15 times larger than those of 48-year-old Scots pine plantation, suggesting that undisturbed old forests, with high rainfall and minimal evapotranspiration in winter, may act as important areas for water storage and sinks for storm rainfall to infiltrate and transport to deeper soil layers via preferential flow. The importance of the development of soil hydrological characteristics under different aged forests is discussed.

Relations entre apports terrigènes et conchyliculture dans les Pertuis Charentais

This work presents interactions between quantitative and qualitative river freshwater inputs and the shellfish farming (oyster and mussel) in the Pertuis Charentais. The quantity of freshwater (i.e. salinity) seems to have a weak influence on the shellfish farming contrarily to its quality determined by particulate and dissolved matters contained in the water. In autumn and winter, large precipitations have a "globally positive" effect amending the coastal ecosystem. Associated dissolved nutriments and the organic matter largely determine the quality of the coming spring growth for bred shellfish, itself controlling in turn the annual yield efficiencies. However, in winter their effects are postponed because of strong mineral load, low luminosity and temperature, then limiting the primary production. The spring contributions, directly linked to territorial practices, agriculture and tourism are more variable in quantity and quality from one year to another. They often correspond to high-risk inflows since numerous substances from anthropogenic watersheds can be found diluted in the coastal zone as in the Pertuis Charentais. Their impacts on in situ estuarine ecosystems are still poorly known since these substances are mainly studied and estimated in laboratory in controlled conditions. Several studies showed anthropogenic contaminations (i.e. cadmium, pesticides) could have significant direct or indirect effects on shellfish farming. For instance, the "summer" mortalities between 1990 and 2000 in the South of the Marennes-Oléron bay (MOB), that induced environmental and physiological oyster disorders, could be linked to pesticide effects, measured during consecutive years on the oyster bed of Ronce Perquis in the South of the MOB. The weak results from the spring larval rearing of the IFREMER experimental hatchery in the South of the bay, and chromosomal abnormalities measured on the stocks of wild oysters of the Pertuis could confirm a high-risk spring environment for the shellfish farming. In summer terrestrial inputs are reduced by low precipitations, anthropogenic water removals (drinking water, irrigation) and by plant evapotranspiration. Consequently certain years, a significant salinity increase in water masses of the Pertuis Charentais is observed. However, based on long-term observations, the significant interannual variability noticed in freshwater contributions constitutes one of the most important facts of these last years. When contributions are weak (i.e. 1991 and 2011), the mean annual salinity is 34.5 in the MOB. To the contrary, other years (i.e. 1977, 1981, 1983 and 1988), the mean salinity reduced to 30.5 shows the significant freshwater contributions to the bay. Elsewhere, particularly in the mediterranean region, oyster breeding water conditions characterized by high salinity values show the freshwater does not seem to be necessary for biological functions of the Pacific oyster Crassostrea gigas. Indeed, the oyster embryonic life in particular is well adapted to high salinity values as long as trophic resources are substantial and temperatures remain high. These two factors firstly condition the embryonic survival before the water salinity. Besides, in the Pertuis Charentais, wind conditions and the geographical bloodstock position rather determine the success of the larvae capture than seawater physic-chemical conditions. Finally, a misunderstanding still remains on summer freshwater contributions to the oyster larvae food supply.

Design of passive methane oxidation biosystems considering their response to the presence of capillary barrier effect

La construction des biosystèmes d’oxydation passive du méthane (BOPM) est une option économique et durable pour réduire les émissions de méthane des sites d’enfouissement de déchets et des effets subséquents du réchauffement climatique. Les BOPM sont constitués de deux couches principales: la couche d'oxydation du méthane (MOL) et la couche de distribution du gaz (GDL). L'oxydation du méthane se produit dans la MOL par les réactions biochimiques des bactéries méthanotrophes, et la GDL est construite sous la MOL pour intercepter et distribuer les émissions fugitives de biogaz à la base de la MOL. Fondamentalement, l'efficacité d'un BOPM est définie en fonction de l'efficacité d'oxydation du méthane dans la MOL. Par conséquent, il est indispensable de fournir des conditions adéquates pour les activités bactériennes des méthanotrophes. En plus des paramètres environnementaux, l'intensité et la distribution du biogaz influencent l'efficacité des BOPM, et ils peuvent rendre le matériau de la MOL - avec une grande capacité d'accueillir les activités bactériennes - inutilisables en termes d'oxydation du méthane sur place. L'effet de barrière capillaire le long de l'interface entre la GDL et la MOL peut provoquer des émissions localisées de méthane, due à la restriction ou la distribution non uniforme de l’écoulement ascendant du biogaz à la base de la MOL. L'objectif principal de cette étude est d'incorporer le comportement hydraulique non saturé des BOPM dans la conception des BOPM, afin d’assurer la facilité et la distribution adéquates de l'écoulement du biogaz à la base de la MOL. Les fonctions de perméabilité à l'air des matériaux utilisés pour construire la MOL des BOPM expérimentaux au site d’enfouissement des déchets de St Nicéphore (Québec, Canada), ainsi que celles d'autres de la littérature technique, ont été étudiés pour évaluer le comportement d'écoulement non saturé du gaz dans les matériaux et pour identifier le seuil de migration sans restriction du gaz. Ce dernier seuil a été introduit en tant que un paramètre de conception avec lequel le critère de conception recommandé ici, c’est à dire la longueur de la migration sans restriction de gaz (LMSG), a été défini. La LMSG est considérée comme la longueur le long de l'interface entre la GDL et la MOL où le biogaz peut migrer à travers la MOL sans restriction. En réalisant des simulations numériques avec SEEP/W, les effets de la pente de l'interface, des paramètres définissant la courbe de rétention d'eau, de la fonction de la conductivité hydraulique du matériau de la MOL sur la valeur de la LMSG (représentant la facilité d'écoulement du biogaz à l'interface) et de la distribution de l'humidité (et par conséquent celle du biogaz) ont été évalués. Selon les résultats des simulations, la conductivité hydraulique saturée et la distribution des tailles de pores du matériau de la MOL sont les paramètres les plus importants sur la distribution de l'humidité le long de l'interface. Ce dernier paramètre influe également sur la valeur du degré de saturation et donc la facilité du biogaz à la base de la MOL. La densité sèche du matériau de MOL est un autre paramètre qui contrôle la facilité d'écoulement ascendant du biogaz. Les limitations principales de la présente étude sont associées au nombre de matériaux de MOL testés et à l'incapacité de SEEP/W de considérer l'évapotranspiration. Toutefois, compte tenu des hypothèses raisonnables dans les simulations et en utilisant les données de la littérature, on a essayé de réduire ces limitations. En utilisant les résultats des expériences et des simulations numériques, des étapes et des considérations de conception pour la sélection du matériau de MOL et de la pente d'interface ont été proposées. En effet,le comportement hydraulique non saturé des matériaux serait intégré dans les nécessités de conception pour un BOPM efficace, de sorte que la capacité maximale possible d'oxydation du méthane du matériau de la MOL soit exploitée.

Abordagem rítmica probabilística aplicada em análises de fluxo e estabilidade de taludes

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Programa de Pós-Graduação em Geotecnia, 2015.

Transfer function-noise modeling and spatial interpolation to evaluate the risk of extreme (shallow) water-table levels in the Brazilian Cerrados

Water regimes in the Brazilian Cerrados are sensitive to climatological disturbances and human intervention. The risk that critical water-table levels are exceeded over long periods of time can be estimated by applying stochastic methods in modeling the dynamic relationship between water levels and driving forces such as precipitation and evapotranspiration. In this study, a transfer function-noise model, the so called PIRFICT-model, is applied to estimate the dynamic relationship between water-table depth and precipitation surplus/deficit in a watershed with a groundwater monitoring scheme in the Brazilian Cerrados. Critical limits were defined for a period in the Cerrados agricultural calendar, the end of the rainy season, when extremely shallow levels (< 0.5-m depth) can pose a risk to plant health and machinery before harvesting. By simulating time-series models, the risk of exceeding critical thresholds during a continuous period of time (e.g. 10 days) is described by probability levels. These simulated probabilities were interpolated spatially using universal kriging, incorporating information related to the drainage basin from a digital elevation model. The resulting map reduced model uncertainty. Three areas were defined as presenting potential risk at the end of the rainy season. These areas deserve attention with respect to water-management and land-use planning.

Cálculo da evapotranspiração de referência e da necessidade de irrigação ao longo do Projeto de Integração do Rio São Francisco

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Agronomia e Medicina Veterinária, Programa de Pós-Graduação em Agronomia, 2016.

Comparison of Hydrologic and Hydraulic Characteristics of the Anacostia River to Non-Urban Coastal Streams

Streams in urban areas often utilize channelization and other bank erosion control measures to improve flood conveyance, reduce channel migration, and overbank flooding. This leads to reductions in evapotranspiration and sediment storage on floodplains. The purpose of this study is to quantify the evapotranspiration and sediment transport capacity in the Anacostia Watershed, a large Coastal Plain urban watershed, and to compare these processes to a similar sized non-urban watershed. Times series data of hydrologic and hydraulic changes in the Anacostia, as urbanization progressed between 1939-2014, were also analyzed. The data indicates lower values of warm season runoff in the non-urban stream, suggesting a shift from evapotranspiration to runoff in urban streams. Channelization in the Anacostia also increased flow velocities and decreased high flow width. The high velocities associated with channelization and the removal of floodplain storage sites allows for the continued downstream transport of sediment despite stream bank stabilization.

The Westerly Index as complementary indicator of the North Atlantic oscillation in explaining drought variability across Europe

This paper analyses the influence of different atmospheric circulation indices on the multi-scalar drought variability across Europe by using the Standardized Precipitation Evapotranspiration Index (SPEI). The monthly circulation indices used in this study include the North Atlantic oscillation (NAO), the East Atlantic (EA), the Scandinavian (SCAN) and the East Atlantic-Western Russia (EA-WR) patterns, as well as the recently published Westerly Index (WI), defined as the persistence of westerly winds over the eastern north Atlantic region. The results indicate that European drought variability is better explained by the station-based NAO index and the WI than by any other combination of circulation indices. In northern and central Europe the variability of drought severity for different seasons and time-scales is strongly associated with the WI. On the contrary, the influence of the NAO on southern Europe droughts is stronger than that exerted by the WI. The correlation patterns of the NAO and WI with the SPEI show a spatial complementarity in shaping drought variability across Europe. Lagged correlations of the NAO and WI with the SPEI also indicate enough skill of both indices to anticipate drought severity several months in advance. As long as instrumental series of the NAO and WI are available, their combined use would allow inferring European drought variability for the last two centuries and improve the calibration and interpretation of paleoclimatic proxies associated with drought.