A new stepwise carbon cycle data assimilation system using multiple data streams to constrain the simulated land surface carbon cycle

Acknowledgements. This work was mainly funded by the EU FP7 CARBONES project (contracts FP7-SPACE-2009-1-242316), with also a small contribution from GEOCARBON project (ENV.2011.4.1.1-1-283080). This work used eddy covariance data acquired by the FLUXNET community and in particular by the following networks: AmeriFlux (U.S. Department of Energy, Biological and Environmental Research, Terrestrial Carbon Program; DE-FG02-04ER63917 and DE-FG02-04ER63911), AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada (supported by CFCAS, NSERC, BIOCAP, Environment Canada, and NRCan), GreenGrass, KoFlux, LBA, NECC, OzFlux, TCOS-Siberia, USCCC. We acknowledge the financial support to the eddy covariance data harmonization provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Université Laval and Environment Canada and US Department of Energy and the database development and technical support from Berkeley Water Center, Lawrence Berkeley National Laboratory, Microsoft Research eScience, Oak Ridge National Laboratory, University of California-Berkeley, University of Virginia. Philippe Ciais acknowledges support from the European Research Council through Synergy grant ERC-2013-SyG-610028 “IMBALANCE-P”. The authors wish to thank M. Jung for providing access to the GPP MTE data, which were downloaded from the GEOCARBON data portal (https://www.bgc-jena.mpg.de/geodb/projects/Data.php). The authors are also grateful to computing support and resources provided at LSCE and to the overall ORCHIDEE project that coordinate the development of the code (http://labex.ipsl.fr/orchidee/index.php/about-the-team).

Quantifying wind and pressure effects on trace gas fluxes across the soil–atmosphere interface

Acknowledgements. We would like to acknowledge the manufacturers of the inner toroid: Mark Bentley and Steve Howarth from the University of York, Dept. of Biology, mechanical and electronics workshops respectively. Furthermore, we would like to acknowledge the Forestry Commission for access and aid at Wheldrake Forest, Mike Bailey and Natural Resources Wales for access and assistance at Cors Fochno, and Norrie Russell and the Royal Society for the Protection of Birds for access and aid at Forsinard. We would also like to thank Graham Hambley, James Robinson, and Elizabeth Donkin for equipment preparation and sampling. Phil Ineson is thanked for the loan of essential equipment, site suggestions, and accessible power supply. Funding was provided by the University of York, Dept. of Biology, and by a grant to YAT by the UK Natural Environment Research Council (NE/H01182X/1).

FCE LTER Mid-term Review field trip to TS/Ph-1

http://digitalcommons.fiu.edu/fce_lter_photos/1346/thumbnail.jpg

Firn line elevation of the Dickson Land peninsula

The Dickson Land peninsula is located in central West-Spitsbergen between the NNE branches of Isfjorden. The climatic firn line lying at 500 m causes plateau glaciers with outlet tongues which are characteristic of S-Dickson Land. The distribution of valley glaciers and the variations of the orographic firn line depend on wind direction. In comparing the firn lines established by the methods of LICHTENECKER (1938) and VISSER (1938), to the values calculated by the method of v. HÖFER (1879), differences of up to l07 m are found. These differences may depend on the inclination and distance relationships of the glaciers above and below the real firn lines. During the latest glacial advance, Dickson Land was located on the peripheries of two local glaciation centers. At that time an inland glaciation of West-Spitsbergen did not exist . The formation of a subglacial channel system dates back to the maximum extent of the late glacial phase before 17500 B.P, (+2000/-1375 years). A correlation of postglacial stadia and 14C dated marine terraces (FEYLING-HANSSEN & OLSSON, 1960; FEYLING-HANSSEN, 1965) is possible. Considering isostatic movement and the difference between calculated and real firn lines, a postglacial stadium at about 10400 B. P. can be reconstructed with a firn line lying 265 m above former sea level. On average, the absolute depression below the recent firn line amounted to 246 m. Stagnation at 9650 B.P. coincided with a firn line at 315 m above former sea level and a depression of 173 m. Around 1890 A.D., glacial fluctuations corresponded to a firn line at 415 m (depression: 64 m). To some extent the morphology of the main valleys appears to depend on structure and petrography. Therefore their value as indicators of former glaciations is questionable. The periglacial forms are shown on a large-scale map. At the time of the "Holocene warm interval", between 7000 and 2000 B.P. (FEYLING-HANSSEN, 1955a, 1965), an increase of periglacial activity seems likely. This can be explained by a simultaneous increase in the depth of the active layer in both soil and bedrock.

Calculation of electric field strengths at a sharp edge

Sharp edges were first used for field ionisation mass spectrometry by Beckey. Although Cross and Robertson found that etched metal foils were more effective than razor blades for field ionisation, blades are very convenient for determination of field ionisation mass spectra, as reported by Robertson and Viney. The electric field at the vertex of a sharp edge can be calculated by the method of the conformal transformation. Here we give some equations for the field deduced with the assumption that the edge surface can be approximated by a hyperbola. We also compare two hyperbolae with radii of curvature at the vertex of 500 Angstrom and 1000 Angstrom with the profile of a commercial carbon-steel razor blade.

Influence of climate variability and stocking strategies on greenhouse gas emissions (GHGE), production and profit of a northern Queensland beef cattle herd

Previous studies of greenhouse gas emissions (GHGE) from beef production systems in northern Australia have been based on models of ‘steady-state’ herd structures that do not take into account the considerable inter-annual variation in liveweight gain, reproduction and mortality rates that occurs due to seasonal conditions. Nor do they consider the implications of flexible stocking strategies designed to adapt these production systems to the highly variable climate. The aim of the present study was to quantify the variation in total GHGE (t CO2e) and GHGE intensity (t CO2e/t liveweight sold) for the beef industry in northern Australia when variability in these factors was considered. A combined GRASP–Enterprise modelling platform was used to simulate a breeding–finishing beef cattle property in the Burdekin River region of northern Queensland, using historical climate data from 1982–2011. GHGE was calculated using the method of Australian National Greenhouse Gas Inventory. Five different stocking-rate strategies were simulated with fixed stocking strategies at moderate and high rates, and three flexible stocking strategies where the stocking rate was adjusted annually by up to 5%, 10% or 20%, according to pasture available at the end of the growing season. Variation in total annual GHGE was lowest in the ‘fixed moderate’ (~9.5 ha/adult equivalent (AE)) stocking strategy, ranging from 3799 to 4471 t CO2e, and highest in the ‘fixed high’ strategy (~5.9 ha/AE), which ranged from 3771 to 7636 t CO2e. The ‘fixed moderate’ strategy had the least variation in GHGE intensity (15.7–19.4 t CO2e/t liveweight sold), while the ‘flexible 20’ strategy (up to 20% annual change in AE) had the largest range (10.5–40.8 t CO2e/t liveweight sold). Across the five stocking strategies, the ‘fixed moderate’ stocking-rate strategy had the highest simulated perennial grass percentage and pasture growth, highest average rate of liveweight gain (121 kg/steer), highest average branding percentage (74%) and lowest average breeding-cow mortality rate (3.9%), resulting in the lowest average GHGE intensity (16.9 t CO2e/t liveweight sold). The ‘fixed high’ stocking rate strategy (~5.9 ha/AE) performed the poorest in each of these measures, while the three flexible stocking strategies were intermediate. The ‘fixed moderate’ stocking strategy also yielded the highest average gross margin per AE carried and per hectare. These results highlight the importance of considering the influence of climate variability on stocking-rate management strategies and herd performance when estimating GHGE. The results also support a body of previous work that has recommended the adoption of moderate stocking strategies to enhance the profitability and ecological stability of beef production systems in northern Australia.

Improved production of chlorogenic acid from cell suspension cultures of Lonicera macranthoids

Purpose: To evaluate the potential of Lonicera macranthoids Hand. -Mazz. Yulei1 suspension culture system for enhanced production of the main secondary metabolite, chlorogenic acid. Methods: The callus of L. macranthoides Hand.-Mazz. “Yulei1” was suspension cultured in B5 liquid medium supplemented with different plant growth regulators. Biomass accumulation was calculated by weight method and chlorogenic acid production was measured using high performance liquid chromatography (HPLC). HPLC was carried out on C18 analytical column at 35 °C and the detection wavelength was set at 324 nm. Results: The results showed that maximum accumulation of biomass and chlorogenic acid were achieved 15 days after culture growth. The optimized conditions for biomass accumulation and chlorogenic acid production were 50 g/L of inoculum on fresh weight basis, B5 medium supplemented with plant growth regulators, 30 - 40 g/L sucrose and initial medium pH of 5.5. Maximum accumulation of chlorogenic acid and biomass were observed when the culture medium was supplemented with 2.0 mg/L6-BA. Optimal accumulation of chlorogenic acid was observed using combination of hormones 2.0 mg/L 6-Benzyladenine (BA) + 0.5 mg/L2, 4-Dichlorophenoxyacetic acid (2,4-D), while optimal accumulation of biomass was observed with 2.0 mg/L 6-BA + 2.0 mg/L2, 4-D. In addition, phenylalanine also contributed to the synthesis of chlorogenic acid at a concentration > 50 mg/L. Conclusion: Cell suspension cultures of L. macranthoides Hand.-Mazz. “Yulei1” have successfully been established. The findings provide a potential basis for large scale production of chlorogenic acid using cell suspension cultures of L. macranthoides.

Degradacao termica e catalitica da borra oleosa de Petroleo com materiais nanoestruturados al-mcm-41 e AL-SBA-15

Aiming to reduce and reuse waste oil from oily sludge generated in large volumes by the oil industry, types of nanostructured materials Al-MCM-41 and Al-SBA-15, with ratios of Si / Al = 50, were synthesized , and calcined solids used as catalysts in the degradation of oily sludge thermocatalytic oil from oilfield Canto do Amaro, in the state of Rio Grande do Norte. Samples of nanostructured materials were characterized by thermogravimetric analysis (TG / DTG), X-ray diffraction (XRD), scanning electron microscopy (SEM), absorption spectroscopy in the infrared Fourier transform (FT-IR) and adsorption nitrogen (BET). The characterization showed that the synthesized materials resulted in a catalyst nanostructure, and ordered pore diameter and surface area according to existing literature. The oily sludge sample was characterized by determining the API gravity and sulfur content and SARA analysis (saturates, aromatics, resins and asphaltenes). The results showed a material equivalent to the average oil with API gravity of 26.1, a low sulfur content and considerable amount of resins and asphaltenes, presented above in the literature. The thermal and catalytic degradation of the oily sludge oil was performed from room temperature to 870 ° C in the ratios of heating of 5, 10 and 20 ° C min-1. The curves generated by TG / DTG showed a more accelerated degradation of oily sludge when it introduced the nanostructured materials. These results were confirmed by activation energy calculated by the method of Flynn-Wall, in the presence of catalysts reduced energy, in particular in the range of cracking, showing the process efficiency, mainly for extraction of lightweight materials of composition of oily sludge, such as diesel and gasoline

Échanges d’énergie et d’eau des écosystèmes nordiques dans un contexte de changement climatique

Le réchauffement climatique affecte fortement les régions nordiques du Canada où le dégel du pergélisol discontinu à sa limite sud est accompagné du mouvement de la limite des arbres vers le nord en zone de pergélisol continu. Ces altérations faites aux paysages de la Taïga des Plaines sont le point de départ de plusieurs rétroactions puisque les changements apportés aux caractéristiques de la surface (au niveau de l’albédo, l’humidité du sol et la rugosité de la surface) vont à leur tour entraîner des modifications biophysiques et éventuellement influencer l’augmentation ou la diminution subséquente des températures et de l’humidité de l’air. Seulement, il y a un nombre important de facteurs d’influence qu’il est difficile de projeter toutes les boucles rétroactives qui surviendront avec les présents changements climatiques en régions nordiques. Dans le but de caractériser les échanges d’eau et d’énergie entre la surface et l’atmosphère de trois sites des Territoires du Nord-Ouest subissant les conséquences de l’augmentation des températures de l’air, la méthode micro-météorologique de covariance des turbulences fut utilisée en 2013 aux sites de Scotty Creek (forêt boréale et tourbière nordique en zone de pergélisol sporadique-discontinu), de Havikpak Creek (forêt boréale nordique en zone de pergélisol continu) et de Trail Valley Creek (toundra arctique en zone de pergélisol continu). En identifiant les procédés biotiques et abiotiques (ex. intensité lumineuse, disponibilité en eau, etc.) d’évapotranspiration aux trois sites, les contrôles par l’eau et l’énergie furent caractérisés et permirent ainsi de projeter une augmentation de la limitation en eau, mais surtout en énergie du site de Trail Valley Creek. La répartition de l’énergie projetée est semblable à celle de Havikpak Creek, avec une augmentation de la proportion du flux de chaleur sensible au détriment de celui latent suite aux modifications des caractéristiques de la surface (albédo, rugosité et humidité du sol). L’augmentation relative du flux d’énergie sensible laisse présager une boucle rétroactive positive de l’augmentation des températures de l’air à ce site. Ensuite, en comparant des données modelées de la hauteur de la couche limite planétaire et des données provenant de profils atmosphériques d’Environnement Canada entre les trois sites, les changements de hauteur de cette couche atmosphérique furent aussi projetés. Trail Valley Creek pourrait connaître une hausse de la hauteur de sa couche limite planétaire avec le temps alors que Scotty Creek connaîtrait une diminution de celle-ci. Ces changements au niveau des couches atmosphériques liés à la répartition des flux d’énergie dans les écosystèmes se répercuteraient alors sur le climat régional de façon difficile à déterminer pour l’instant. Les changements apportés désignent une boucle rétroactive positive des températures de l’air à Trail Valley Creek et l’inverse à Scotty Creek. Les deux axes d’analyse arrivent donc aux mêmes conclusions et soulignent aussi l’importance de l’influence mutuelle entre le climat et les caractéristiques spécifiques des écosystèmes à la surface.

Échanges d’énergie et d’eau des écosystèmes nordiques dans un contexte de changement climatique

Le réchauffement climatique affecte fortement les régions nordiques du Canada où le dégel du pergélisol discontinu à sa limite sud est accompagné du mouvement de la limite des arbres vers le nord en zone de pergélisol continu. Ces altérations faites aux paysages de la Taïga des Plaines sont le point de départ de plusieurs rétroactions puisque les changements apportés aux caractéristiques de la surface (au niveau de l’albédo, l’humidité du sol et la rugosité de la surface) vont à leur tour entraîner des modifications biophysiques et éventuellement influencer l’augmentation ou la diminution subséquente des températures et de l’humidité de l’air. Seulement, il y a un nombre important de facteurs d’influence qu’il est difficile de projeter toutes les boucles rétroactives qui surviendront avec les présents changements climatiques en régions nordiques. Dans le but de caractériser les échanges d’eau et d’énergie entre la surface et l’atmosphère de trois sites des Territoires du Nord-Ouest subissant les conséquences de l’augmentation des températures de l’air, la méthode micro-météorologique de covariance des turbulences fut utilisée en 2013 aux sites de Scotty Creek (forêt boréale et tourbière nordique en zone de pergélisol sporadique-discontinu), de Havikpak Creek (forêt boréale nordique en zone de pergélisol continu) et de Trail Valley Creek (toundra arctique en zone de pergélisol continu). En identifiant les procédés biotiques et abiotiques (ex. intensité lumineuse, disponibilité en eau, etc.) d’évapotranspiration aux trois sites, les contrôles par l’eau et l’énergie furent caractérisés et permirent ainsi de projeter une augmentation de la limitation en eau, mais surtout en énergie du site de Trail Valley Creek. La répartition de l’énergie projetée est semblable à celle de Havikpak Creek, avec une augmentation de la proportion du flux de chaleur sensible au détriment de celui latent suite aux modifications des caractéristiques de la surface (albédo, rugosité et humidité du sol). L’augmentation relative du flux d’énergie sensible laisse présager une boucle rétroactive positive de l’augmentation des températures de l’air à ce site. Ensuite, en comparant des données modelées de la hauteur de la couche limite planétaire et des données provenant de profils atmosphériques d’Environnement Canada entre les trois sites, les changements de hauteur de cette couche atmosphérique furent aussi projetés. Trail Valley Creek pourrait connaître une hausse de la hauteur de sa couche limite planétaire avec le temps alors que Scotty Creek connaîtrait une diminution de celle-ci. Ces changements au niveau des couches atmosphériques liés à la répartition des flux d’énergie dans les écosystèmes se répercuteraient alors sur le climat régional de façon difficile à déterminer pour l’instant. Les changements apportés désignent une boucle rétroactive positive des températures de l’air à Trail Valley Creek et l’inverse à Scotty Creek. Les deux axes d’analyse arrivent donc aux mêmes conclusions et soulignent aussi l’importance de l’influence mutuelle entre le climat et les caractéristiques spécifiques des écosystèmes à la surface.

TECTONIC CONTROLS ON FLUVIAL LANDSCAPE DEVELOPMENT IN CENTRAL-EASTERN PORTUGAL: INSIGHTS FROM LONG PROFILE TRIBUTARY STREAM ANALYSES

This study examines the long profiles of tributaries of the Tejo (Tagus) and Zêzere rivers in central eastern Portugal (West Iberia) in order to provide new insights into the patterns, timing and controls on drainage development during the Pleistocene to Holocene incision stage. The long profiles were extracted from lower order tributary streams associated with the trunk drainage of the Tejo River and one main tributary, the Zêzere River (Fig. 1). These streams flow through a landscape strongly influenced by variations in bedrock lithology (mainly granites and metasediments), fault structures delimiting crustal blocks with distinct uplift rates, and a base-level lowering history (tectonic uplift / eustatic). The long profiles of the tributaries of the Tejo and Zêzere rivers record a series of transient and permanent knickpoints. The permanent knickpoints have direct correlation with the bedrock strength, corresponding to the outcropping of very hard quartzites or to the transition from softer (slates/metagreywaques) to harder (granite) basement. The analyzed streams/rivers record also an older transient knickpoint/knickzone separating: a) an upstream relict graded profile, with lower steepness and higher concavity, that reflects a long period of quasi-equilibrium conditions reached after the beginning of the incision stage; and b) a downstream reach displaying a rejuvenated long profile, with steeper gradient and lower concavity, particularly for the final segment, which is often convex (Fig. 2). The rejuvenated reaches testify the upstream propagation of several incision waves that are the response of each stream to continuous or increasing crustal uplift and dominant periods of base-level lowering by the trunk drainages, coeval of low sea level conditions. The long profiles and their morphological configurations enabled spatial and relative temporal patterns of incision to be quantified for each individual tributary stream. The incision values of streams flowing in uplifted blocks of the Portuguese Central Range (PCR) (ca.380-280 m) indicate differential uplift and are higher than the incision values of streams flowing on the adjacent South Portugal planation surface – the Meseta (ca. 200 m). The normalized steepness index, calculated using the method of Wobus et al. (2006), proved to be sensitive to active tectonics, as lower ksn values were found in relict graded profiles of streams located in less uplifted blocks, (e.g. Sertã stream in the PCR), or in those flowing through tectonic depressions. Fig. 1 – Geological map of the study area. 1 – fluvial terraces (Pleistocene); 2 – sedimentary cover (Paleogene and Neogene); 3 – slates and metasandstones (Devonian); 4 – slates and quartzites (Silurian); 5 – quartzites (Ordovician); 6 – slates and metagreywackes (Precambrian to Cambrian); 7 – slates, metagreywackes and limestones (Precambrian); 8 – granites and ortogneisses; 9 – diorites and gabros; 10 - fault. SFf – Sobreira Formosa fault; Sf – Sertã fault; Pf – Ponsul fault; Gf – Grade fault. The differential uplift indicated by the distribution of the ksn values and by the fluvial incision was likely accumulated on a few major faults, as the Sobreira Formosa fault (SFf), thus corroborating the tectonic activity of these faults. Due to the fact that the relict graded profiles can be correlated with other geomorphic references documented in the study area, namely the T1 terrace of the Tagus River (with an age of ca. 1 Myr), the following incision rates can be estimated: a) for the studied streams located in uplifted blocks of the PCR, 0.38 m/kyr to 0.28 m/kyr; b) for the streams flowing on the South Portugal planation surface, 0.20 m/kyr. The differential uplift inferred between crustal blocks in the study area corroborates the neotectonic activity of the bordering faults, which has been proposed in previous studies based upon less robust data. Fig. 2 – Longitudinal profile of the Nisa stream a tributary of the Tejo River. Note the equilibrium relict profile upstream the older transient knickpoint (hatched line) and the downstream rejuvenated profile (continuous line). Legend: tKP – transient knickpoint; rKp – resistant knickpoint; Mt – schist and phyllite; Gr – granite; Hf – hornfels; Og – orthogneisse. In the inset Distance – Slope plots, fill circles correspond to the relict graded profile, crosses correspond to the rejuvenated profile located downstream the older transient knickpoint (tKP).

Lake–atmosphere interactions at Alqueva reservoir: a case study in the summer of 2014

The study of lake–atmosphere interactions was the main purpose of a 2014 summer experiment at Alqueva reservoir in Portugal. Near-surface fluxes of momentum, heat and mass [water vapour (H2O) and carbon dioxide (CO2)] were obtained with the new Campbell Scientific’s IRGASON Integrated Open-Path CO2/H2O Gas Analyser and 3D Sonic Anemometer between 2 June and 2 October. On average, the reservoir was releasing energy in the form of sensible and latent heat flux during the study period. At the end of the 75 d, the total evaporation was estimated as 490.26 mm. A high correlation was found between the latent heat flux and the wind speed (R = 0.97). The temperature gradient between air and water was positive between 12 and 21 UTC, causing a negative sensible heat flux, and negative during the rest of the day, triggering a positive sensible heat flux. The reservoir acted as a sink of atmospheric CO2 with an average rate of −0.026 mg m−2 s−1. However, at a daily scale we found an unexpected uptake between 0 and 9 UTC and almost null flux between 13 and 19 UTC. Potential reasons for this result are further discussed. The net radiation was recorded for the same period and water column heat storage was estimated using water temperature profiles. The energy balance closure for the analysed period was 81%. In-water solar spectral downwelling irradiance profiles were measured with a new device allowing measurements independent of the solar zenith angle, which enabled the computation of the attenuation coefficient of light in the water column. The average attenuation coefficient for the photosynthetically active radiation spectral region varied from 0.849 ± 0.025 m−1 on 30 July to 1.459 ± 0.007 m−1 on 25 September.

Canopy-scale biophysical controls of transpiration and evaporation in the Amazon Basin.

Canopy and aerodynamic conductances (gC and gA) are two of the key land surface biophysical variables that control the land surface response of land surface schemes in climate models. Their representation is crucial for predicting transpiration (λET) and evaporation (λEE) flux components of the terrestrial latent heat flux (λE), which has important implications for global climate change and water resource management. By physical integration of radiometric surface temperature (TR) into an integrated framework of the Penman?Monteith and Shuttleworth?Wallace models, we present a novel approach to directly quantify the canopy-scale biophysical controls on λET and λEE over multiple plant functional types (PFTs) in the Amazon Basin. Combining data from six LBA (Large-scale Biosphere-Atmosphere Experiment in Amazonia) eddy covariance tower sites and a TR-driven physically based modeling approach, we identified the canopy-scale feedback-response mechanism between gC, λET, and atmospheric vapor pressure deficit (DA), without using any leaf-scale empirical parameterizations for the modeling. The TR-based model shows minor biophysical control on λET during the wet (rainy) seasons where λET becomes predominantly radiation driven and net radiation (RN) determines 75 to 80 % of the variances of λET. However, biophysical control on λET is dramatically increased during the dry seasons, and particularly the 2005 drought year, explaining 50 to 65 % of the variances of λET, and indicates λET to be substantially soil moisture driven during the rainfall deficit phase. Despite substantial differences in gA between forests and pastures, very similar canopy?atmosphere "coupling" was found in these two biomes due to soil moisture-induced decrease in gC in the pasture. This revealed the pragmatic aspect of the TR-driven model behavior that exhibits a high sensitivity of gC to per unit change in wetness as opposed to gA that is marginally sensitive to surface wetness variability. Our results reveal the occurrence of a significant hysteresis between λET and gC during the dry season for the pasture sites, which is attributed to relatively low soil water availability as compared to the rainforests, likely due to differences in rooting depth between the two systems. Evaporation was significantly influenced by gA for all the PFTs and across all wetness conditions. Our analytical framework logically captures the responses of gC and gA to changes in atmospheric radiation, DA, and surface radiometric temperature, and thus appears to be promising for the improvement of existing land?surface?atmosphere exchange parameterizations across a range of spatial scales.

Canopy-scale biophysical controls of transpiration and evaporation in the Amazon Basin.

Canopy and aerodynamic conductances (gC and gA) are two of the key land surface biophysical variables that control the land surface response of land surface schemes in climate models. Their representation is crucial for predicting transpiration (?ET) and evaporation (?EE) flux components of the terrestrial latent heat flux (?E), which has important implications for global climate change and water resource management. By physical integration of radiometric surface temperature (TR) into an integrated framework of the Penman?Monteith and Shuttleworth?Wallace models, we present a novel approach to directly quantify the canopy-scale biophysical controls on ?ET and ?EE over multiple plant functional types (PFTs) in the Amazon Basin. Combining data from six LBA (Large-scale Biosphere-Atmosphere Experiment in Amazonia) eddy covariance tower sites and a TR-driven physically based modeling approach, we identified the canopy-scale feedback-response mechanism between gC, ?ET, and atmospheric vapor pressure deficit (DA), without using any leaf-scale empirical parameterizations for the modeling. The TR-based model shows minor biophysical control on ?ET during the wet (rainy) seasons where ?ET becomes predominantly radiation driven and net radiation (RN) determines 75 to 80?% of the variances of ?ET. However, biophysical control on ?ET is dramatically increased during the dry seasons, and particularly the 2005 drought year, explaining 50 to 65?% of the variances of ?ET, and indicates ?ET to be substantially soil moisture driven during the rainfall deficit phase. Despite substantial differences in gA between forests and pastures, very similar canopy?atmosphere "coupling" was found in these two biomes due to soil moisture-induced decrease in gC in the pasture. This revealed the pragmatic aspect of the TR-driven model behavior that exhibits a high sensitivity of gC to per unit change in wetness as opposed to gA that is marginally sensitive to surface wetness variability. Our results reveal the occurrence of a significant hysteresis between ?ET and gC during the dry season for the pasture sites, which is attributed to relatively low soil water availability as compared to the rainforests, likely due to differences in rooting depth between the two systems. Evaporation was significantly influenced by gA for all the PFTs and across all wetness conditions. Our analytical framework logically captures the responses of gC and gA to changes in atmospheric radiation, DA, and surface radiometric temperature, and thus appears to be promising for the improvement of existing land?surface?atmosphere exchange parameterizations across a range of spatial scales.

Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (C) 2007 Elsevier B.V. All rights reserved.