Downscaling of surface moisture flux and precipitation in the Ebro Valley (Spain) using analogues and analogues followed by random forests and multiple linear regression

In this paper, reanalysis fields from the ECMWF have been statistically downscaled to predict from large-scale atmospheric fields, surface moisture flux and daily precipitation at two observatories (Zaragoza and Tortosa, Ebro Valley, Spain) during the 1961-2001 period. Three types of downscaling models have been built: (i) analogues, (ii) analogues followed by random forests and (iii) analogues followed by multiple linear regression. The inputs consist of data (predictor fields) taken from the ERA-40 reanalysis. The predicted fields are precipitation and surface moisture flux as measured at the two observatories. With the aim to reduce the dimensionality of the problem, the ERA-40 fields have been decomposed using empirical orthogonal functions. Available daily data has been divided into two parts: a training period used to find a group of about 300 analogues to build the downscaling model (1961-1996) and a test period (19972001), where models' performance has been assessed using independent data. In the case of surface moisture flux, the models based on analogues followed by random forests do not clearly outperform those built on analogues plus multiple linear regression, while simple averages calculated from the nearest analogues found in the training period, yielded only slightly worse results. In the case of precipitation, the three types of model performed equally. These results suggest that most of the models' downscaling capabilities can be attributed to the analogues-calculation stage.

Detecting changes in surface moisture and water table position with spectral changes in surface vegetation in northern peatlands

Due to warmer and drier conditions, wildland fire has been increasing in extent into peatland ecosystems during recent decades. As such, there is an increasing need for broadly applicable tools to detect surface peat moisture, in order to ascertain the susceptibility of peat burning, and the vulnerability of deep peat consumption in the event of a wildfire. In this thesis, a field portable spectroradiometer was used to measure surface reflectance of two Sphagnum moss dominated peatlands. Relationships were developed correlating spectral indices to surface moisture as well as water table position. Spectral convolutions were also applied to the high resolution spectra to represent spectral sensitivity of earth observing sensors. Band ratios previously used to monitor surface moisture with these sensors were assessed. Strong relationships to surface moisture and water table position are evident for both the narrowband indices as well as broadened indices. This study also found a dependence of certain spectral relationships on changes in vegetation cover by leveraging an experimental vegetation manipulation. Results indicate broadened indices employing the 1450-1650 nm region may be less stable under changing vegetation cover than those located in the 1200 nm region.

Simulation and diagnosis of large-scale atmospheric moisture and surface hydrology over North America

We describe a preliminary investigation into large-scale atmospheric and surface moisture variations over North America. We compare large-scale hydrologic budgets in the Los Alamos general circulation model (GCM) to observed precipitation and vertically integrated atmospheric moisture fluxes derived from the National Meteorological Center's operational analyses. THe GCM faithfully simulates the integrated flux divergence and P-E differences. However, the integrated moisture content is too low, and precipitation and evaporation are too high. The model produces summertime soil moisture dryness, which supports previous studies showing increased droughts under warmer conditions.

Pore, surface and fillet moisture in centrifuged beds of coarse coal

The moisture content of the coarse coking coal product from the centrifuges of preparation plants was investigated to evaluate the contribution of three types of water: that held internally in pores, that in fillets at points of contacts between the particles, and the moisture covering the surface. A standardised laboratory centrifuge test was used to measure the total non-centrifugable moisture (NCM) content and also the quantity held in internal pores, called NCMi. The fillet moisture NCMf was estimated by means of a formulation which relies on experimentally measured holdup volumes, supplemented by a physical model. The surface moisture NCMs could then be derived by difference. The NCMf, which depends on the body force, the particle size and the surface tension and contact angle of the liquid, ranges from effectively zero for large particles to 10% for fines. The surface moisture NCMs is of the order of 0.5% for high rank coals and increases to 4.5% for lower rank coals. (C) 2002 Elsevier Science Ltd. All rights reserved.

Evaporative Moisture Loss from Heterogeneous Stone: Material- Environment Interactions During Drying

The complexities of evaporation from structurally and mineralogically heterogeneous sandstone (Locharbriggs Sandstone) are investigated through a laboratory-based experiment in which a variety of environmental conditions are simulated. Data reported demonstrate the significance of material-environment interactions on the spatial and temporal variability of evaporative dynamics. Evaporation from porous stone is determined by the interplay between environmental, material and solution properties, which govern the rate and mode by which water is transmitted to, and subsequently removed from, an evaporating surface. Initially evaporation is marked by high rates of moisture loss controlled by external atmospheric conditions; then, when a critical level of surface moisture content is reached, hydraulic continuity between the stone surface and subsurface is disrupted and the drying front recedes
beneath the surface, evaporation rates decrease and are controlled by the ability of the material to transport water vapour to the surface. Pore size distribution and connectivity, as well as other material properties, control the timing of each stage of evaporation and the nature of the transition.

These experimental data highlight the complexity of evaporation, demonstrating that different regions of the same stone can exhibit varying moisture dynamics during drying and that the rate and nature of evaporative loss differs under different environmental conditions. The results identify the importance of material-environment interactions during drying and that stone micro-environmental conditions cannot be inferred from ambient data alone.
These data have significance for understanding the spatial distribution of stone surface weathering-related morphologies in both the natural and built environments where mineralogical and/or structural heterogeneity creates differences in moisture flux and hence variable drying rates. Such differences may provide a clearer explanation for the initiation and subsequent development of complex weathering responses where areas of significant deterioration can be found alongside areas that exhibit little or no evidence surface breakdown.

Stationary Eddies and Zonal Variations of the Global Hydrological Cycle in a Changing Climate

This thesis advances our physical understanding of the sensitivity of the hydrological cycle to global warming. Specifically, it focuses on changes in the longitudinal (zonal) variation of precipitation minus evaporation (P - E), which is predominantly controlled by planetary-scale stationary eddies. By studying idealized general circulation model (GCM) experiments with zonally varying boundary conditions, this thesis examines the mechanisms controlling the strength of stationary-eddy circulations and their role in the hydrological cycle. The overarching goal of this research is to understand the cause of changes in regional P - E with global warming. An understanding of such changes can be useful for impact studies focusing on water availability, ecosystem management, and flood risk.

Based on a moisture-budget analysis of ERA-Interim data, we establish an approximation for zonally anomalous P - E in terms of surface moisture content and stationary-eddy vertical motion in the lower troposphere. Part of the success of this approximation comes from our finding that transient-eddy moisture fluxes partially cancel the effect of stationary-eddy moisture advection, allowing divergent circulations to dominate the moisture budget. The lower-tropospheric vertical motion is related to horizontal motion in stationary eddies by Sverdrup and Ekman balance. These moisture- and vorticity-budget balances also hold in idealized and comprehensive GCM simulations across a range of climates.

By examining climate changes in the idealized and comprehensive GCM simulations, we are able to show the utility of the vertical motion P - E approximation for splitting changes in zonally anomalous P - E into thermodynamic and dynamic components. Shifts in divergent stationary-eddy circulations dominate changes in zonally anomalous P - E. This limits the local utility of the "wet gets wetter, dry gets drier” idea, where existing P - E patterns are amplified with warming by the increase in atmospheric moisture content, with atmospheric circulations held fixed. The increase in atmospheric moisture content manifests instead in an increase in the amplitude of the zonally anomalous hydrological cycle as measured by the zonal variance of P - E. However, dynamic changes, particularly the slowdown of divergent stationary-eddy circulations, limit the strengthening of the zonally anomalous hydrological cycle. In certain idealized cases, dynamic changes are even strong enough to reverse the tendency towards "wet gets wetter, dry gets drier” with warming.

Motivated by the importance of stationary-eddy vertical velocities in the moisture budget analysis, we examine controls on the amplitude of stationary eddies across a wide range of climates in an idealized GCM with simple topographic and ocean-heating zonal asymmetries. An analysis of the thermodynamic equation in the vicinity of topographic forcing reveals the importance of on-slope surface winds, the midlatitude isentropic slope, and latent heating in setting the amplitude of stationary waves. The response of stationary eddies to climate change is determined primarily by the strength of zonal surface winds hitting the mountain. The sensitivity of stationary-eddies to this surface forcing increases with climate change as the slope of midlatitude isentropes decreases. However, latent heating also plays an important role in damping the stationary-eddy response, and this damping becomes stronger with warming as the atmospheric moisture content increases. We find that the response of tropical overturning circulations forced by ocean heat-flux convergence is described by changes in the vertical structure of moist static energy and deep convection. This is used to derive simple scalings for the Walker circulation strength that capture the monotonic decrease with warming found in our idealized simulations.

Through the work of this thesis, the advances made in understanding the amplitude of stationary-waves in a changing climate can be directly applied to better understand and predict changes in the zonally anomalous hydrological cycle.

Testing solar forcing of pervasive Holocene climate cycles

The temporal and spatial extent of Holocene climate change is an area of considerable uncertainty, with solar forcing recently proposed to be the origin of cycles identified in the North Atlantic region. To address these issues we have developed an annually resolved record of changes in Irish bog tree populations over the last 7468 years which, together with radiocarbon-dated bog and lake-edge populations, extend the dataset back to 9000 yr ago. The Irish trees underpin the internationally accepted radiocarbon calibration curve, used to derive a proxy of solar activity, and allow us to test solar forcing of Holocene climate change. Tree populations and age structures provide unambiguous evidence of major shifts in Holocene surface moisture, with a dominant cyclicity of 800 yr, similar to marine cycles in the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. The cycles, however, are not coherent with changes in solar activity (both being on the same absolute timescale), indicating that Holocene North Atlantic climate variability at the millennial and centennial scale is not driven by a linear response to changes in solar activity.

Holocene climate change and past Irish societal response

The extent to which North Atlantic Holocene climatic perturbations influenced past human societies is an area of considerable uncertainty and fierce debate. Ireland is ideally placed to help resolve this issue, being occupied for over 9000 yr and located on the eastern Atlantic seaboard, a region dominated by westerly airflow. Irish bog and lake tree populations provide unambiguous evidence of major shifts in surface moisture through the Holocene similar to cycles recorded in the marine realm of the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. To test for human response to these cycles we summed the probabilities of 465 radiocarbon ages obtained from Irish archaeological contexts and observe enhanced archaeological visibility during periods of sustained wet conditions. These results suggest either increasing density of human populations in key, often defensive locations, and/or the development of subsistence strategies to overcome changing conditions, the latter recently proposed as a significant factor in avoiding societal collapse. Regardless, we demonstrate environmental change is a significantly more important factor in influencing human activity in the landscape than has hitherto been acknowledged.

Granular fertilizer drying and agglomeration in storage

In the manufacture of granular NPK fertilizer the product is cooled before packaging and storage in moisture-proof bags. It has been shown that the temperature of the fertilizer prior to packing is significant in that at high temperatures, drying of the granules takes place in the bag which causes an increase in the humidity of the air surrounding the granules and thus an increase in moisture content at the granule - granule interface. This surface moisture was shown to increase the likelihood of agglomeration in the fertilizer by a capillary adhesion/unconfined yield stress model. An iterative model was set up to establish conditions that would prevent drying occurring, which takes into account fertilizer drying rate, fertilizer cooling rate cooling rate and the effect of coating oils on the drying mechanism.

The use and meanings of 'time of wetness' in understanding stone decay

Moisture is a well documented, and crucial, control on the nature of stone decay. The term time of wetness has frequently been adopted to describe how long a stone block is wet, with a view to understanding the impact of this on decay processes. Although this term has proved conceptually useful, it has been used in different ways, by different groups to mean mean quite different things. For example, the time of wetness for a stone block surface (the traditional understanding) may be quite different from that of a block interior, controlled by the different dynamics of wetting and drying in those zones. Thus, surface wetting will occur regularly (sometimes swiftly followed by drying, depending on the time of year), with block interior wetting requiring the accumulation of surface moisture to penetrate to depth (more likely in autumn and winter months), and drying out much more slowly. This relatively new but important perspective, framed in the context of climate change, is crucial to understanding the length of time stone may remain damp at depth following a period of prolonged precipitation. The nature and speed of drying is also relevant in quantifying time of wetness of both surfaces and the interior of building stones.
These ideas related to time of wetness have implications for decay processes, specifically how a prolonged time of deep wetness may re-focus the emphasis of salt weathering in natural building stones toward chemical action. Literature on chemical change is discussed, suggesting that chemical change occurring during periods of prolonged wetness is likely to be significant in itself, with implications for weakening the stone (in terms of, for example, cement dissolution or grain boundary weakening) and exacerbating physical damage from salt crystallisation when blocks finally dry out.

Développement d'un indice de vulnérabilité à l'érosion éolienne à partir d'images satellitales, dans le Bassin arachidier du Sénégal : cas de la région de Thiès.

L’érosion éolienne est un problème environnemental parmi les plus sévères dans les régions arides, semi-arides et les régions sèches sub-humides de la planète. L’érosion des sols accélérée par le vent provoque des dommages à la fois localement et régionalement. Sur le plan local, elle cause la baisse des nutriments par la mobilisation des particules les plus fines et de la matière organique. Cette mobilisation est une des causes de perte de fertilité des sols avec comme conséquence, une chute de la productivité agricole et une réduction de la profondeur de la partie arable. Sur le plan régional, les tempêtes de poussières soulevées par le vent ont un impact non négligeable sur la santé des populations, et la déposition des particules affecte les équipements hydrauliques tels que les canaux à ciel ouvert ainsi que les infrastructures notamment de transport. Dans les régions où les sols sont fréquemment soumis à l’érosion éolienne, les besoins pour des études qui visent à caractériser spatialement les sols selon leur degré de vulnérabilité sont grands. On n’a qu’à penser aux autorités administratives qui doivent décider des mesures à prendre pour préserver et conserver les potentialités agropédologiques des sols, souvent avec des ressources financières modestes mises à leur disposition. Or, dans certaines de ces régions, comme notre territoire d’étude, la région de Thiès au Sénégal, ces études font défaut. En effet, les quelques études effectuées dans cette région ou dans des contextes géographiques similaires ont un caractère plutôt local et les approches suivies (modèles de pertes des sols) nécessitent un nombre substantiel de données pour saisir la variabilité spatiale de la dynamique des facteurs qui interviennent dans le processus de l’érosion éolienne. La disponibilité de ces données est particulièrement problématique dans les pays en voie de développement, à cause de la pauvreté en infrastructures et des problèmes de ressources pour le monitoring continu des variables environnementales. L’approche mise de l’avant dans cette recherche vise à combler cette lacune en recourant principalement à l’imagerie satellitale, et plus particulièrement celle provenant des satellites Landsat-5 et Landsat-7. Les images Landsat couvrent la presque totalité de la zone optique du spectre exploitable par télédétection (visible, proche infrarouge, infrarouge moyen et thermique) à des résolutions relativement fines (quelques dizaines de mètres). Elles permettant ainsi d’étudier la distribution spatiale des niveaux de vulnérabilité des sols avec un niveau de détails beaucoup plus fin que celui obtenu avec des images souvent utilisées dans des études environnementales telles que AVHRR de la série de satellites NOAA (résolution kilométrique). De plus, l’archive complet des images Landsat-5 et Landsat-7 couvrant une période de plus de 20 ans est aujourd’hui facilement accessible. Parmi les paramètres utilisés dans les modèles d’érosion éolienne, nous avons identifiés ceux qui sont estimables par l’imagerie satellitale soit directement (exemple, fraction du couvert végétal) soit indirectement (exemple, caractérisation des sols par leur niveau d’érodabilité). En exploitant aussi le peu de données disponibles dans la région (données climatiques, carte morphopédologique) nous avons élaboré une base de données décrivant l’état des lieux dans la période de 1988 à 2002 et ce, selon les deux saisons caractéristiques de la région : la saison des pluies et la saison sèche. Ces données par date d’acquisition des images Landsat utilisées ont été considérées comme des intrants (critères) dans un modèle empirique que nous avons élaboré en modulant l’impact de chacun des critères (poids et scores). À l’aide de ce modèle, nous avons créé des cartes montrant les degrés de vulnérabilité dans la région à l’étude, et ce par date d’acquisition des images Landsat. Suite à une série de tests pour valider la cohérence interne du modèle, nous avons analysé nos cartes afin de conclure sur la dynamique du processus pendant la période d’étude. Nos principales conclusions sont les suivantes : 1) le modèle élaboré montre une bonne cohérence interne et est sensible aux variations spatiotemporelles des facteurs pris en considération 2); tel qu’attendu, parmi les facteurs utilisés pour expliquer la vulnérabilité des sols, la végétation vivante et l’érodabilité sont les plus importants ; 3) ces deux facteurs présentent une variation importante intra et inter-saisonnière de sorte qu’il est difficile de dégager des tendances à long terme bien que certaines parties du territoire (Nord et Est) aient des indices de vulnérabilité forts, peu importe la saison ; 4) l’analyse diachronique des cartes des indices de vulnérabilité confirme le caractère saisonnier des niveaux de vulnérabilité dans la mesure où les superficies occupées par les faibles niveaux de vulnérabilité augmentent en saison des pluies, donc lorsque l’humidité surfacique et la végétation active notamment sont importantes, et décroissent en saison sèche ; 5) la susceptibilité, c’est-à-dire l’impact du vent sur la vulnérabilité est d’autant plus forte que la vitesse du vent est élevée et que la vulnérabilité est forte. Sur les zones où la vulnérabilité est faible, les vitesses de vent élevées ont moins d’impact. Dans notre étude, nous avons aussi inclus une analyse comparative entre les facteurs extraits des images Landsat et celles des images hyperspectrales du satellite expérimental HYPERION. Bien que la résolution spatiale de ces images soit similaire à celle de Landsat, les résultats obtenus à partir des images HYPERION révèlent un niveau de détail supérieur grâce à la résolution spectrale de ce capteur permettant de mieux choisir les bandes spectrales qui réagissent le plus avec le facteur étudié. Cette étude comparative démontre que dans un futur rapproché, l’amélioration de l’accessibilité à ce type d’images permettra de raffiner davantage le calcul de l’indice de vulnérabilité par notre modèle. En attendant cette possibilité, on peut de contenter de l’imagerie Landsat qui offre un support d’informations permettant tout de même d’évaluer le niveau de fragilisation des sols par l’action du vent et par la dynamique des caractéristiques des facteurs telles que la couverture végétale aussi bien vivante que sénescente.

Mesozoic climates: General circulation models and the rock record

General circulation models (GCMs) use the laws of physics and an understanding of past geography to simulate climatic responses. They are objective in character. However, they tend to require powerful computers to handle vast numbers of calculations. Nevertheless, it is now possible to compare results from different GCMs for a range of times and over a wide range of parameterisations for the past, present and future (e.g. in terms of predictions of surface air temperature, surface moisture, precipitation, etc.). GCMs are currently producing simulated climate predictions for the Mesozoic, which compare favourably with the distributions of climatically sensitive facies (e.g. coals, evaporites and palaeosols). They can be used effectively in the prediction of oceanic upwelling sites and the distribution of petroleum source rocks and phosphorites. Models also produce evaluations of other parameters that do not leave a geological record (e.g. cloud cover, snow cover) and equivocal phenomena such as storminess. Parameterisation of sub-grid scale processes is the main weakness in GCMs (e.g. land surfaces, convection, cloud behaviour) and model output for continental interiors is still too cold in winter by comparison with palaeontological data. The sedimentary and palaeontological record provides an important way that GCMs may themselves be evaluated and this is important because the same GCMs are being used currently to predict possible changes in future climate. The Mesozoic Earth was, by comparison with the present, an alien world, as we illustrate here by reference to late Triassic, late Jurassic and late Cretaceous simulations. Dense forests grew close to both poles but experienced months-long daylight in warm summers and months-long darkness in cold snowy winters. Ocean depths were warm (8 degrees C or more to the ocean floor) and reefs, with corals, grew 10 degrees of latitude further north and south than at the present time. The whole Earth was warmer than now by 6 degrees C or more, giving more atmospheric humidity and a greatly enhanced hydrological cycle. Much of the rainfall was predominantly convective in character, often focused over the oceans and leaving major desert expanses on the continental areas. Polar ice sheets are unlikely to have been present because of the high summer temperatures achieved. The model indicates extensive sea ice in the nearly enclosed Arctic seaway through a large portion of the year during the late Cretaceous, and the possibility of sea ice in adjacent parts of the Midwest Seaway over North America. The Triassic world was a predominantly warm world, the model output for evaporation and precipitation conforming well with the known distributions of evaporites, calcretes and other climatically sensitive facies for that time. The message from the geological record is clear. Through the Phanerozoic, Earth's climate has changed significantly, both on a variety of time scales and over a range of climatic states, usually baldly referred to as "greenhouse" and "icehouse", although these terms disguise more subtle states between these extremes. Any notion that the climate can remain constant for the convenience of one species of anthropoid is a delusion (although the recent rate of climatic change is exceptional). (c) 2006 Elsevier B.V. All rights reserved.

Influence of vegetation on the local climate and hydrology in the tropics: sensitivity to soil parameters

Land use change with accompanying major modifications to the vegetation cover is widespread in the tropics, due to increasing demands for agricultural land, and may have significant impacts on the climate. This study investigates (1) the influence of vegetation on the local climate in the tropics; (2) how that influence varies from region to region; and (3) how the sensitivity of the local climate to vegetation, and hence land use change, depends on the hydraulic characteristics of the soil. A series of idealised experiments with the Hadley Centre atmospheric model, HadAM3, are described in which the influence of vegetation in the tropics is assessed by comparing the results of integrations with and without tropical vegetation. The sensitivity of the results to the soil characteristics is then explored by repeating the experiments with a differing, but equally valid, description of soil hydraulic parameters. The results have shown that vegetation has a significant moderating effect on the climate throughout the tropics by cooling the surface through enhanced latent heat fluxes. The influence of vegetation is, however, seasonally dependent, with much greater impacts during the dry season when the availability of surface moisture is limited. Furthermore, there are significant regional variations both in terms of the magnitude of the cooling and in the response of the precipitation. Not all regions show a feedback of vegetation on the local precipitation; this result has been related both to vegetation type and to the prevailing meteorological conditions. An important finding has been the sensitivity of the results to the specification of the soil hydraulic parameters. The introduction of more freely draining soils has changed the soil-moisture contents of the control, vegetated system and has reduced, significantly, the climate sensitivity to vegetation and by implication, land use change. Changes to the soil parameters have also had an impact on the soil hydrology and its interaction with vegetation, by altering the partitioning between fast and slow runoff processes. These results raise important questions about the representation of highly heterogeneous soil characteristics in climate models, as well as the potential influence of land use change on the soil characteristics themselves.

Numerical modelling of the energy and water cycle of the Baltic Sea

This paper will introduce the Baltex research programme and summarize associated numerical modelling work which has been undertaken during the last five years. The research has broadly managed to clarify the main mechanisms determining the water and energy cycle in the Baltic region, such as the strong dependence upon the large scale atmospheric circulation. It has further been shown that the Baltic Sea has a positive water balance, albeit with large interannual variations. The focus on the modelling studies has been the use of limited area models at ultra-high resolution driven by boundary conditions from global models or from reanalysis data sets. The programme has further initiated a comprehensive integration of atmospheric, land surface and hydrological modelling incorporating snow, sea ice and special lake models. Other aspects of the programme include process studies such as the role of deep convection, air sea interaction and the handling of land surface moisture. Studies have also been undertaken to investigate synoptic and sub-synoptic events over the Baltic region, thus exploring the role of transient weather systems for the hydrological cycle. A special aspect has been the strong interests and commitments of the meteorological and hydrological services because of the potentially large societal interests of operational applications of the research. As a result of this interests special attention has been put on data-assimilation aspects and the use of new types of data such as SSM/I, GPS-measurements and digital radar. A series of high resolution data sets are being produced. One of those, a 1/6 degree daily precipitation climatology for the years 1996–1999, is such a unique contribution. The specific research achievements to be presented in this volume of Meteorology and Atmospheric Physics is the result of a cooperative venture between 11 European research groups supported under the EU-Framework programmes.