A priori parameterisation of the CERES soil-crop models and tests against several European data sets

Mechanistic soil-crop models have become indispensable tools to investigate the effect of management practices on the productivity or environmental impacts of arable crops. Ideally these models may claim to be universally applicable because they simulate the major processes governing the fate of inputs such as fertiliser nitrogen or pesticides. However, because they deal with complex systems and uncertain phenomena, site-specific calibration is usually a prerequisite to ensure their predictions are realistic. This statement implies that some experimental knowledge on the system to be simulated should be available prior to any modelling attempt, and raises a tremendous limitation to practical applications of models. Because the demand for more general simulation results is high, modellers have nevertheless taken the bold step of extrapolating a model tested within a limited sample of real conditions to a much larger domain. While methodological questions are often disregarded in this extrapolation process, they are specifically addressed in this paper, and in particular the issue of models a priori parameterisation. We thus implemented and tested a standard procedure to parameterize the soil components of a modified version of the CERES models. The procedure converts routinely-available soil properties into functional characteristics by means of pedo-transfer functions. The resulting predictions of soil water and nitrogen dynamics, as well as crop biomass, nitrogen content and leaf area index were compared to observations from trials conducted in five locations across Europe (southern Italy, northern Spain, northern France and northern Germany). In three cases, the model’s performance was judged acceptable when compared to experimental errors on the measurements, based on a test of the model’s root mean squared error (RMSE). Significant deviations between observations and model outputs were however noted in all sites, and could be ascribed to various model routines. In decreasing importance, these were: water balance, the turnover of soil organic matter, and crop N uptake. A better match to field observations could therefore be achieved by visually adjusting related parameters, such as field-capacity water content or the size of soil microbial biomass. As a result, model predictions fell within the measurement errors in all sites for most variables, and the model’s RMSE was within the range of published values for similar tests. We conclude that the proposed a priori method yields acceptable simulations with only a 50% probability, a figure which may be greatly increased through a posteriori calibration. Modellers should thus exercise caution when extrapolating their models to a large sample of pedo-climatic conditions for which they have only limited information.

Soil quality assessment through a multi-approach analysis in soils of abandoned terraced land in NE Spain

The abandonment of agricultural land in mountainous areas has been an outstanding problem along the last century and has captured the attention of scientists, technicians and administrations, for the dramatic consequences sometimes occurred due to soil instability, steep slopes, rainfall regimes and wildfires. Hidromorfological and pedological alterations causing exceptional floods and accelerated erosion processes has therefore been studied, identifying the cause in the loss of landscape heterogeneity. Through the disappearance of agricultural works and drainage maintenance, slope stability has resulted severely affected. The mechanization of agriculture has caused the displacement of vines, olives and corks trees cultivation in terraced areas along the Mediterranean catchment towards more economically suitable areas. On the one hand, land use and management changes have implicated sociological changes as well, transforming areas inhabited by agricultural communities into deserted areas where the colonization of disorganized spontaneous vegetation has buried a valuable rural patrimony. On the other hand, lacking of planning and management of the abandoned areas has produced badlands and infertile soils due to wildfire and high erosion rates strongly degrading the whole ecosystems. In other cases, after land abandonment a process of soil regeneration has been recorded. Investigations have been conducted in a part of NE Spain where extended areas of terraced soils previously cultivated have been abandoned in the last century. The selected environments were semi-abandoned vineyards, semi-abandoned olive groves, abandoned stands of cork trees, abandoned stands of pine trees, scrubland of Cistaceaea, scrubland of Ericaceaea, and pasture. The research work was focused on the study of most relevant physical, chemical and biological soil properties, as well as runoff and erosion under soils with different plant cover to establish the abandonment effect on soil quality, due to the peculiarity and vulnerability of these soils with a much reduced depth. The period of observation was carried out from autumn 2009 to autumn 2010. The sediment concentration of soil erosion under vines was recorded as 34.52 g/l while under pasture it was 4.66 g/l. In addition, the soil under vines showed the least amount of organic matter, which was 12 times lower than all other soil environments. The carbon dioxide (CO2) and total glomalin (TG) ratio to soil organic carbon (SOC) in this soil was 0.11 and 0.31 respectively. However, the soil under pasture contained a higher amount of organic matter and showed that the CO2 and TG ratio to SOC was 0.02 and 0.11 respectively indicating that the soil under pasture better preserves the soil carbon pool. A similar trend was found in the intermediate soils in the sequence of land use change and abandonment. Soil structural stability increased in the two soil fractions investigated (0.25-2.00 mm, 2.0-5.6 mm) especially in those soils that did not undergo periodical perturbations like wildfires. Soil quality indexes were obtained by using relevant physical and chemical soil parameters. Factor analysis carried out to study the relationship between all soil parameters allowed to related variables and environments and identify those areas that better contribute to soil quality towards others that may need more attention to avoid further degradation processes

Avaluació de la capacitat de segrest de carboni d'un sòl restaurat amb fangs de depuradora: paper de les propietats físiques del sòl

S’ha estudiat el grau de protecció física del carboni dins d’agregats de diferents mides en un sòl d’una pedrera restaurada fa 17 anys amb terres adobades amb fangs de depuradora, i s’ha interpretat el seu paper en el context del segrest de carboni en el sòl. La metodologia aplicada es basa en la d’humitejament ràpid dels agregats del sòl per immersió en aigua (Le Bissonnais, 1996), que simula l’estabilitat d'un sòl sec que es veu sotmès a processos tals com la inundació local o saturació ràpida. També s’ha determinat la quantitat de carboni oxidable present en els agregats del sòl amb el mètode de Nelson i Sommers (1982). Els resultats han mostrat que l’adobat amb fangs de depuradora contribueix a augmentar el contingut de carboni orgànic en els agregats del sòl i n’estimula el segrest a mitjà termini (unes 10 tones ha-1 en 17 anys), aportant estabilitat al sòl i protegint físicament el carboni orgànic dins dels agregats de mida major (5-2 mm). A més s’ha constatat que per determinar l’estabilitat del carboni segrestat en el sòl cal conèixer com es distribueix entre les diferents mides d’agregats. Finalment, l’augment del segrest de carboni en el sòl propiciat per l’aplicació dels fangs de depuradora li dóna més capacitat per fixar CO2 atmosfèric.

Environmental Kuznets Curves for Carbon Emissions: A Critical Survey

The empirical finding of an inverse U-shaped relationship between per capita income and pollution, the so-called Environmental Kuznets Curve (EKC), suggests that as countries experience economic growth, environmental deterioration decelerates and thus becomes less of an issue. Focusing on the prime example of carbon emissions, the present article provides a critical review of the new econometric techniques that have questioned the baseline polynomial specification in the EKC literature. We discuss issues related to the functional form, heterogeneity, “spurious” regressions and spatial dependence to address whether and to what extent the EKC can be observed. Despite these new approaches, there is still no clear-cut evidence supporting the existence of the EKC for carbon emissions. JEL classifications: C20; Q32; Q50; O13 Keywords: Environmental Kuznets Curve; Carbon emissions; Functional form; Heterogeneity; “Spurious” regressions; Spatial dependence.Residential satisfaction is often used as a barometer to assess the performance of public policy and programmes designed to raise individuals' well-being. However, the fact that responses elicited from residents might be biased by subjective, non-observable factors casts doubt on whether these responses can be taken as trustable indicators of the individuals' housing situation. Emotional factors such as aspirations or expectations might affect individuals' cognitions of their true residential situation. To disentangle this puzzle, we investigated whether identical residential attributes can be perceived differently depending on tenure status. Our results indicate that tenure status is crucial not only in determining the level of housing satisfaction, but also regarding how dwellers perceive their housing characteristics. Keywords: Housing satisfaction, subjective well-being, homeownership. JEL classification: D1, R2.

Why are slip lengths so large in carbon nanotubes?

The enhanced flow in carbon nanotubes is explained using a mathematical model that includes a depletion layer with reduced viscosity near the wall. In the limit of large tubes the model predicts no noticeable enhancement. For smaller tubes the model predicts enhancement that increases as the radius decreases. An analogy between the reduced viscosity and slip-length models shows that the term slip-length is misleading and that on surfaces which are smooth at the nanoscale it may be thought of as a length-scale associated with the size of the depletion region and viscosity ratio. The model therefore provides a physical interpretation of the classical Navier slip condition and explains why `slip-lengths' may be greater than the tube radius.

Inclusion of soil erosion impacts in life cycle assessment on a global scale: application to energy crops in Spain

Purpose: Despite the fundamental role of ecosystem goods and services in sustaining human activities, there is no harmonized and internationally agreed method for including them in life cycle assessment (LCA). The main goal of this study was to develop a globally applicable and spatially resolved method for assessing land-use impacts on the erosion regulation ecosystem service.Methods: Soil erosion depends much on location. Thus, unlike conventional LCA, the endpoint method was regionalized at the grid-cell level (5 arc-minutes, approximately 10×10 km2) to reflect the spatial conditions of the site. Spatially explicit characterization factors were not further aggregated at broader spatial scales. Results and discussion: Life cycle inventory data of topsoil and topsoil organic carbon (SOC) losses were interpreted at the endpoint level in terms of the ultimate damage to soil resources and ecosystem quality. Human health damages were excluded from the assessment. The method was tested on a case study of five three-year agricultural rotations, two of them with energy crops, grown in several locations in Spain. A large variation in soil and SOC losses was recorded in the inventory step, depending on climatic and edaphic conditions. The importance of using a spatially explicit model and characterization factors is shown in the case study.Conclusions and outlook: The regionalized assessment takes into account the differences in soil erosion-related environmental impacts caused by the great variability of soils. Taking this regionalized framework as the starting point, further research should focus on testing the applicability of the method trough the complete life cycle of a product and on determining an appropriate spatial scale at which to aggregate characterization factors, in order to deal with data gaps on location of processes, especially in the background system. Additional research should also focus on improving reliability of the method by quantifying and, insofar as it is possible, reducing uncertainty.

Analysis of sensitivity with respect to a compositional parameter : A comment on 'Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling' by L. Loosvelt and co-authors

A comment about the article “Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling” writen by L. Loosvelt and co-authors. The present comment is centered in three specific points. The first one is related to the fact that the authors avoid the use of ilr-coordinates. The second one refers to some generalization of sensitivity analysis when input parameters are compositional. The third tries to show that the role of the Dirichlet distribution in the sensitivity analysis is irrelevant

Does soil moisture overrule temperature dependence of soil respiration in Mediterranean riparian forests?

Soil respiration (SR) is a major component of ecosystems' carbon cycles and represents the second largest CO2 flux in the terrestrial biosphere. Soil temperature is considered to be the primary abiotic control on SR, whereas soil moisture is the secondary control factor. However, soil moisture can become the dominant control on SR in very wet or dry conditions. Determining the trigger that makes soil moisture as the primary control factor of SR will provide a deeper understanding on how SR changes under the projected future increase in droughts. Specific objectives of this study were (1) to investigate the seasonal variations and the relationship between SR and both soil temperature and moisture in a Mediterranean riparian forest along a groundwater level gradient; (2) to determine soil moisture thresholds at which SR is controlled by soil moisture rather than by temperature; (3) to compare SR responses under different tree species present in a Mediterranean riparian forest (Alnus glutinosa, Populus nigra and Fraxinus excelsior). Results showed that the heterotrophic soil respiration rate, groundwater level and 30 cm integral soil moisture (SM30) decreased significantly from the riverside moving uphill and showed a pronounced seasonality. SR rates showed significant differences between tree species, with higher SR for P. nigra and lower SR for A. glutinosa. The lower threshold of soil moisture was 20 and 17% for heterotrophic and total SR, respectively. Daily mean SR rate was positively correlated with soil temperature when soil moisture exceeded the threshold, with Q10 values ranging from 1.19 to 2.14; nevertheless, SR became decoupled from soil temperature when soil moisture dropped below these thresholds.

Improved understanding of drought controls on seasonal variation in Mediterranean forest canopy CO2 and water fluxes through combined in situ measurements and ecosystem modelling

Water stress is a defining characteristic of Mediterranean ecosystems, and is likely to become more severe in the coming decades. Simulation models are key tools for making predictions, but our current understanding of how soil moisture controls ecosystem functioning is not sufficient to adequately constrain parameterisations. Canopy-scale flux data from four forest ecosystems with Mediterranean-type climates were used in order to analyse the physiological controls on carbon and water flues through the year. Significant non-stomatal limitations on photosynthesis were detected, along with lesser changes in the conductance-assimilation relationship. New model parameterisations were derived and implemented in two contrasting modelling approaches. The effectiveness of two models, one a dynamic global vegetation model ('ORCHIDEE'), and the other a forest growth model particularly developed for Mediterranean simulations ('GOTILWA+'), was assessed and modelled canopy responses to seasonal changes in soil moisture were analysed in comparison with in situ flux measurements. In contrast to commonly held assumptions, we find that changing the ratio of conductance to assimilation under natural, seasonally-developing, soil moisture stress is not sufficient to reproduce forest canopy CO2 and water fluxes. However, accurate predictions of both CO2 and water fluxes under all soil moisture levels encountered in the field are obtained if photosynthetic capacity is assumed to vary with soil moisture. This new parameterisation has important consequences for simulated responses of carbon and water fluxes to seasonal soil moisture stress, and should greatly improve our ability to anticipate future impacts of climate changes on the functioning of ecosystems in Mediterranean-type climates.

Rhizodeposition of organic carbon by plants with contrasting traits for resource acquisition: responses to different fertility regimes

Background and aims Rhizodeposition plays an important role in mediating soil nutrient availability in ecosystems. However, owing to methodological difficulties (i.e., narrow zone of soil around roots, rapid assimilation by soil microbes) fertility-induced changes in rhizodeposition remain mostly unknown. Methods We developed a novel long-term continuous 13C labelling method to address the effects of two levels of nitrogen (N) fertilization on rhizodeposited carbon (C) by species with different nutrient acquisition strategies. Results Fertility-induced changes in rhizodeposition were modulated by root responses to N availability rather than by changes in soil microbial biomass. Differences among species were mostly related to plant biomass: species with higher total leaf and root biomass also had higher total rhizodeposited C, whereas species with lower root biomass had higher specific rhizodeposited C (per gram root mass). Experimental controls demonstrated that most of the biases commonly associated with this type of experiment (i.e., long-term steady-state labelling) were avoided using our methodological approach. Conclusions These results suggest that the amount of rhizodeposited C from plants grown under different levels of N were driven mainly by plant biomass and root morphology rather than microbial biomass. They also underline the importance of plant characteristics (i.e., biomass allocation) as opposed to traits associated with plant resource acquisition strategies in predicting total C rhizodeposition.

Comparative performance of the stable isotope signatures of carbon, nitrogen and oxygen in assessing early vigour and grain yield in durum wheat

The present paper studied the performance of the stable isotope signatures of carbon (δ13C), nitrogen (δ15N) and oxygen (δ18O) in plants when used to assess early vigour and grain yield (GY) in durum wheat growing under mild and moderate Mediterranean stress conditions. A collection of 114 recombinant inbred lines was grown under rainfed (RF) and supplementary irrigation (IR) conditions. Broad sense heritabilities (H2) for GY and harvest index (HI) were higher under RF conditions than under IR. Broad sense heritabilities for δ13C were always above 0·60, regardless of the plant part studied, with similar values for IR and RF trials. Some of the largest genetic correlations with GY were those shown by the δ13C content of the flag leaf blade and mature grains. Under both water treatments, mature grains showed the highest negative correlations between δ13C and GY across genotypes. Flag leaf δ13C was negatively correlated with GY only under RF conditions. The δ13C in seedlings was negatively correlated, under IR conditions only, with GY but also with early vigour. The sources of variation in early vigour were studied by stepwise analysis using the stable isotope signatures measured in seedlings. The δ13C was able to explain almost 0·20 of this variation under RF, but up to 0·30 under IR. In addition, nitrogen concentration in seedlings accounted for another 0·05 of variation, increasing the amount explained to 0·35. The sources of variation in GY were also studied through stable isotope signatures and biomass of different plant parts: δ13C was always the first parameter to appear in the models for both water conditions, explaining c. 0·20 of the variation. The second parameter (δ15N or N concentration of grain, or biomass at maturity) depended on the water conditions and the plant tissue being analysed. Oxygen isotope composition (δ18O) was only able to explain a small amount of the variation in GY. In this regard, despite the known and previously described value of δ13C as a tool in breeding, δ15N is confirmed as an additional tool in the present study. Oxygen isotope composition does not seem to offer any potential, at least under the conditions of the present study.