A study on the performance of flexible pavements on mature soil subgrades

The country has witnessed tremendous increase in the vehicle population and increased axle loading pattern during the last decade, leaving its road network overstressed and leading to premature failure. The type of deterioration present in the pavement should be considered for determining whether it has a functional or structural deficiency, so that appropriate overlay type and design can be developed. Structural failure arises from the conditions that adversely affect the load carrying capability of the pavement structure. Inadequate thickness, cracking, distortion and disintegration cause structural deficiency. Functional deficiency arises when the pavement does not provide a smooth riding surface and comfort to the user. This can be due to poor surface friction and texture, hydro planning and splash from wheel path, rutting and excess surface distortion such as potholes, corrugation, faulting, blow up, settlement, heaves etc. Functional condition determines the level of service provided by the facility to its users at a particular time and also the Vehicle Operating Costs (VOC), thus influencing the national economy. Prediction of the pavement deterioration is helpful to assess the remaining effective service life (RSL) of the pavement structure on the basis of reduction in performance levels, and apply various alternative designs and rehabilitation strategies with a long range funding requirement for pavement preservation. In addition, they can predict the impact of treatment on the condition of the sections. The infrastructure prediction models can thus be classified into four groups, namely primary response models, structural performance models, functional performance models and damage models. The factors affecting the deterioration of the roads are very complex in nature and vary from place to place. Hence there is need to have a thorough study of the deterioration mechanism under varied climatic zones and soil conditions before arriving at a definite strategy of road improvement. Realizing the need for a detailed study involving all types of roads in the state with varying traffic and soil conditions, the present study has been attempted. This study attempts to identify the parameters that affect the performance of roads and to develop performance models suitable to Kerala conditions. A critical review of the various factors that contribute to the pavement performance has been presented based on the data collected from selected road stretches and also from five corporations of Kerala. These roads represent the urban conditions as well as National Highways, State Highways and Major District Roads in the sub urban and rural conditions. This research work is a pursuit towards a study of the road condition of Kerala with respect to varying soil, traffic and climatic conditions, periodic performance evaluation of selected roads of representative types and development of distress prediction models for roads of Kerala. In order to achieve this aim, the study is focused into 2 parts. The first part deals with the study of the pavement condition and subgrade soil properties of urban roads distributed in 5 Corporations of Kerala; namely Thiruvananthapuram, Kollam, Kochi, Thrissur and Kozhikode. From selected 44 roads, 68 homogeneous sections were studied. The data collected on the functional and structural condition of the surface include pavement distress in terms of cracks, potholes, rutting, raveling and pothole patching. The structural strength of the pavement was measured as rebound deflection using Benkelman Beam deflection studies. In order to collect the details of the pavement layers and find out the subgrade soil properties, trial pits were dug and the in-situ field density was found using the Sand Replacement Method. Laboratory investigations were carried out to find out the subgrade soil properties, soil classification, Atterberg limits, Optimum Moisture Content, Field Moisture Content and 4 days soaked CBR. The relative compaction in the field was also determined. The traffic details were also collected by conducting traffic volume count survey and axle load survey. From the data thus collected, the strength of the pavement was calculated which is a function of the layer coefficient and thickness and is represented as Structural Number (SN). This was further related to the CBR value of the soil and the Modified Structural Number (MSN) was found out. The condition of the pavement was represented in terms of the Pavement Condition Index (PCI) which is a function of the distress of the surface at the time of the investigation and calculated in the present study using deduct value method developed by U S Army Corps of Engineers. The influence of subgrade soil type and pavement condition on the relationship between MSN and rebound deflection was studied using appropriate plots for predominant types of soil and for classified value of Pavement Condition Index. The relationship will be helpful for practicing engineers to design the overlay thickness required for the pavement, without conducting the BBD test. Regression analysis using SPSS was done with various trials to find out the best fit relationship between the rebound deflection and CBR, and other soil properties for Gravel, Sand, Silt & Clay fractions. The second part of the study deals with periodic performance evaluation of selected road stretches representing National Highway (NH), State Highway (SH) and Major District Road (MDR), located in different geographical conditions and with varying traffic. 8 road sections divided into 15 homogeneous sections were selected for the study and 6 sets of continuous periodic data were collected. The periodic data collected include the functional and structural condition in terms of distress (pothole, pothole patch, cracks, rutting and raveling), skid resistance using a portable skid resistance pendulum, surface unevenness using Bump Integrator, texture depth using sand patch method and rebound deflection using Benkelman Beam. Baseline data of the study stretches were collected as one time data. Pavement history was obtained as secondary data. Pavement drainage characteristics were collected in terms of camber or cross slope using camber board (slope meter) for the carriage way and shoulders, availability of longitudinal side drain, presence of valley, terrain condition, soil moisture content, water table data, High Flood Level, rainfall data, land use and cross slope of the adjoining land. These data were used for finding out the drainage condition of the study stretches. Traffic studies were conducted, including classified volume count and axle load studies. From the field data thus collected, the progression of each parameter was plotted for all the study roads; and validated for their accuracy. Structural Number (SN) and Modified Structural Number (MSN) were calculated for the study stretches. Progression of the deflection, distress, unevenness, skid resistance and macro texture of the study roads were evaluated. Since the deterioration of the pavement is a complex phenomena contributed by all the above factors, pavement deterioration models were developed as non linear regression models, using SPSS with the periodic data collected for all the above road stretches. General models were developed for cracking progression, raveling progression, pothole progression and roughness progression using SPSS. A model for construction quality was also developed. Calibration of HDM–4 pavement deterioration models for local conditions was done using the data for Cracking, Raveling, Pothole and Roughness. Validation was done using the data collected in 2013. The application of HDM-4 to compare different maintenance and rehabilitation options were studied considering the deterioration parameters like cracking, pothole and raveling. The alternatives considered for analysis were base alternative with crack sealing and patching, overlay with 40 mm BC using ordinary bitumen, overlay with 40 mm BC using Natural Rubber Modified Bitumen and an overlay of Ultra Thin White Topping. Economic analysis of these options was done considering the Life Cycle Cost (LCC). The average speed that can be obtained by applying these options were also compared. The results were in favour of Ultra Thin White Topping over flexible pavements. Hence, Design Charts were also plotted for estimation of maximum wheel load stresses for different slab thickness under different soil conditions. The design charts showed the maximum stress for a particular slab thickness and different soil conditions incorporating different k values. These charts can be handy for a design engineer. Fuzzy rule based models developed for site specific conditions were compared with regression models developed using SPSS. The Riding Comfort Index (RCI) was calculated and correlated with unevenness to develop a relationship. Relationships were developed between Skid Number and Macro Texture of the pavement. The effort made through this research work will be helpful to highway engineers in understanding the behaviour of flexible pavements in Kerala conditions and for arriving at suitable maintenance and rehabilitation strategies. Key Words: Flexible Pavements – Performance Evaluation – Urban Roads – NH – SH and other roads – Performance Models – Deflection – Riding Comfort Index – Skid Resistance – Texture Depth – Unevenness – Ultra Thin White Topping

Root effects on the turnover of grain legume residues in soil

Das Ziel dieser Arbeit war, die Einflüsse von Wurzeln und Rhizodeposition auf den Umsatz von Körnerleguminosenresiduen und damit verknüpfte mikrobielle Prozesse zu untersuchen. In einem integrierten Versuch wurden Ackerbohne (Vicia faba L.), Erbse (Pisum sativum L.) und Weiße Lupine (Lupinus albus L.) untersucht. Der Versuch bestand aus drei Teilen, zwei Gefäß-Experimenten und einem Inkubationsexperiment, in denen ausgehend von einem Gefäß-Experiment derselbe Boden und dasselbe Pflanzenmaterial verwendet wurden. In Experiment I wurde die Stickstoff-Rhizodeposition der Körnerleguminosenarten, definiert als wurzelbürtiger N nach dem Entfernen aller sichtbaren Wurzeln im Boden, gemessen und der Verbleib des Rhizodepositions-N in verschiednenen Bodenpools untersucht. Dazu wurden die Leguminosen in einem Gefäßversuch unter Verwendung einer in situ 15N-Docht-Methode mit einer 15N Harnstofflösung pulsmarkiert. In Experiment II wurde der Umsatz der N-Rhizodeposition der Körnerleguminosen und der Einfluss der Rhizodeposition auf den anschließenden C- und N-Umsatz der Körnerleguminosenresiduen in einem Inkubationsexperiment untersucht. In Experiment III wurde der N-Transfer aus den Körnerleguminosenresiduen einschließlich N-Rhizodeposition in die mikrobielle Biomasse und die Folgefrüchte Weizen (Triticum aestivum L.) und Raps (Brassica napus L.) in einem Gewächshaus-Gefäßversuch ermittelt. Die in situ 15N Docht-Markierungs-Methode wies hohe 15N Wiederfindungsraten von ungefähr 84 Prozent für alle drei Leguminosenarten auf und zeigte eine vergleichsweise homogene 15N Verteilung zwischen verschiedenen Pflanzenteilen zur Reife. Die Wurzeln zeigten deutliche Effekte auf die N-Dynamik nach dem Anbau von Körnerleguminosen. Die Effekte konnten auf die N-Rhizodeposition und deren anschließenden Umsatz, Einflüsse der Rhizodeposition von Körnerleguminosen auf den anschließenden Umsatz ihrer Residuen (Stängel, Blätter, erfassbare Wurzeln) und die Wirkungen nachfolgender Nichtleguminosen auf den Umsatzprozess der Residuen zurückgeführt werden: Die N-Rhizodeposition betrug zur Reife der Pflanzen bezogen auf die Gesamt-N- Aufnahme 13 Prozent bei Ackerbohne und Erbse und 16 Prozent bei Weißer Lupine. Bezogen auf den Residual N nach Ernte der Körner erhöhte sich der relative Anteil auf 35 - 44 Prozent. Die N-Rhizodeposition ist daher ein wesentlicher Pool für die N-Bilanz von Körnerleguminosen und trägt wesentlich zur Erklärung positiver Fruchtfolgeeffekte nach Körnerleguminosen bei. 7 - 21 Prozent des Rhizodepositions-N wurden als Feinwurzeln nach Nasssiebung (200 µm) wiedergefunden. Nur 14 - 18 Prozent des Rhizodepositions-N wurde in der mikrobiellen Biomasse und ein sehr kleiner Anteil von 3 - 7 Prozent in der mineralischen N Fraktion gefunden. 48 bis 72 Prozent der N-Rhizodeposition konnte in keinem der untersuchten Pools nachgewiesen werden. Dieser Teil dürfte als mikrobielle Residualmasse immobilisiert worden sein. Nach 168 Tagen Inkubation wurden 21 bis 27 Prozent des Rhizodepositions-N in den mineralisiert. Der mineralisierte N stammte im wesentlichen aus zwei Pools: Zwischen 30 Prozent und 55 Prozent wurde aus der mikrobiellen Residualmasse mineralisiert und eine kleinere Menge stammte aus der mikrobielle Biomasse. Der Einfluss der Rhizodeposition auf den Umsatz der Residuen war indifferent. Durch Rhizodeposition wurde die C Mineralisierung der Leguminosenresiduen nur in der Lupinenvariante erhöht, wobei der mikrobielle N und die Bildung von mikrobieller Residualmasse aus den Leguminosenresiduen in allen Varianten durch Rhizodepositionseinflüsse erhöht waren. Das Potential des residualen Körnerleguminosen-N für die N Ernährung von Folgefrüchten war gering. Nur 8 - 12 Prozent des residualen N wurden in den Folgenfrüchten Weizen und Raps wiedergefunden. Durch die Berücksichtigung des Rhizodepositions-N war der relative Anteil des Residual-N bezogen auf die Gesamt-N-Aufnahme der Folgefrucht hoch und betrug zwischen 18 und 46 Prozent. Dies lässt auf einen höheren N-Beitrag der Körnerleguminosen schließen als bisher angenommen wurde. Die residuale N-Aufnahme von Weizen von der Blüte bis zur Reife wurde durch den Residual-N gespeist, der zur Blüte in der mikrobiellen Biomasse immobilisiert worden war. Die gesamte Poolgröße, Residual-N in der mikrobiellen Biomasse und in Weizen, veränderte sich von der Blüte bis zur Reife nicht. Jedoch konnte ein Rest von 80 Prozent des Residual-N in keinem der untersuchten Pools nachgewiesen werden und dürfte als mikrobielle Residualmasse immobilisiert worden sein oder ist noch nicht abgebaut worden. Die zwei unterschiedlichen Folgefrüchte - Weizen und Raps - zeigten sehr ähnliche Muster bei der N-Aufnahme, der Residual-N Wiederfindung und bei mikrobiellen Parametern für die Residuen der drei Körnerleguminosenarten. Ein differenzierender Effekt auf den Umsatz der Residuen bzw. auf das Residual-N-Aneignungsvermögen der Folgefrüchte konnte nicht beobachtet werden.

Flächenhafte und funktionale Analyse kleinräumiger standortbedingter Variationen der Grünlandvegetation und -produktivität

Zur Abbildung heterogener Standorteigenschaften und Ertragspotenziale werden zunehmend flächenhafte Daten nachgefragt. Insbesondere für Grünland, das häufig durch ausgeprägte Standortheterogenität gekennzeichnet ist, ergeben sich hohe Anforderungen an die Wiedergabequalität, denn die realen Verhältnisse sollen in praktikabler Weise möglichst exakt abgebildet werden. Außerdem können flächenhafte Daten genutzt werden, um Zusammenhänge zwischen teilflächenspezifischen Standorteigenschaften und Grünlandaspekten detaillierter zu analysieren und bisher nicht erkannte Wechselbeziehungen nachzuweisen. Für mitteleuropäisches Grünland lagen zu Beginn dieser Arbeit derartige räumliche Untersuchungen nicht oder nur in Teilaspekten vor. Diese Arbeit befasste sich mit der Analyse von Wirkungsbeziehungen zwischen Standort- und Grünlandmerkmalen auf einer im Nordhessischen Hügelland (Deutschland) weitgehend praxisüblicher bewirtschafteten 20 ha großen Weidefläche. Erhoben wurden als Standortfaktoren die Geländemorphologie, die Bodentextur, die Grundnährstoffgehalten sowie als Parameter des Grünlandbestandes die botanische Zusammensetzung, der Ertrag und die Qualitätsparameter. Sie wurden sowohl in einem 50 m-Raster ganzflächig, als auch auf drei 50x50 m großen Teilflächen in erhöhter Beprobungsdichte (6,25 m-Rasterweite) aufgenommen. Die relevanten Fragestellungen zielen auf die räumliche und zeitliche Variabilität von Grünlandbestandesparametern innerhalb von Grünlandflächen sowie deren Abhängigkeit von den Standortfaktoren. Ein weiterer Schwerpunkt war die Überprüfung der Frage, ob die reale Variabilität der Zielvariablen durch die Interpolierung der punktuell erfassten Daten wiedergegeben werden kann. Die Beziehungen zwischen Standort- und Grünlandmerkmalen wurden mit monokausalen und multivariaten Ansätzen untersucht. Die Ergebnisse ließen, unabhängig vom Jahreseinfluss, bereits bestimmte Zusammenhänge zwischen botanischer Zusammensetzung und Standort, auch auf dem untersuchten kleinen Maßstab innerhalb der Grünlandfläche, finden. Demzufolge können unterschiedliche Areale abgegrenzt und charakterisiert werden, die als Grundlage für Empfehlungen zur Ausweisung von Arealen zur teilspezifischen Bewirtschaftung erarbeitet wurden. Die Validierung der interpolierten Daten zeigte, dass die 50-m Rasterbeprobung nur eine begrenzte Wiedergabe der räumlichen Variabilität ermöglicht. Inwieweit derartige Beziehungen quantitativ genauer beschreibbar sind, bleibt auf Grund der verbliebenen unerklärten Varianz im Datensatz dieser Studie offen.

Effects of fertilizer type and rate on partitioning of soil organic matter into pools of different stability

Type and rate of fertilizers influence the level of soil organic carbon (Corg) and total nitrogen (Nt) markedly, but the effect on C and N partitioning into different pools is open to question. The objectives of the present work were to: (i) quantify the impact of fertilizer type and rate on labile, intermediate and passive C and N pools by using a combination of biological, chemical and mathematical methods; (ii) explain previously reported differences in the soil organic matter (SOM) levels between soils receiving farmyard manure with or without biodynamic preparations by using Corg time series and information on SOM partitioning; and (iii) quantify the long-term and short-term dynamics of SOM in density fractions and microbial biomass as affected by fertilizer type and rate and determine the incorporation of crop residues into labile SOM fractions. Samples were taken from a sandy Cambisol from the long-term fertilization trial in Darmstadt, Germany, founded in 1980. The nine treatments (four field replicates) were: straw incorporation plus application of mineral fertilizer (MSI) and application of rotted farmyard manure with (DYN) or without (FYM) addition of biodynamic preparations, each at high (140 – 150 kg N ha-1 year-1; MSIH, DYNH, FYMH), medium (100 kg N ha-1 year-1; MSIM, DYNM, FYMM) and low (50 – 60 kg N ha-1 year-1; MSIL, DYNL, FYML) rates. The main findings were: (i) The stocks of Corg (t ha-1) were affected by fertilizer type and rate and increased in the order MSIL (23.6), MSIM (23.7), MSIH (24.2) < FYML (25.3) < FYMM (28.1), FYMH (28.1). Stocks of Nt were affected in the same way (C/N ratio: 11). Storage of C and N in the modelled labile pools (turnover times: 462 and 153 days for C and N, respectively) were not influenced by the type of fertilizer (FYM and MSI) but depended significantly (p ≤ 0.05) on the application rate and ranged from 1.8 to 3.2 t C ha 1 (7 – 13% of Corg) and from 90 to 140 kg N ha-1 (4-5% of Nt). In the calculated intermediate pool (C/N ratio 7), stocks of C were markedly higher in FYM treatments (15-18 t ha-1) compared to MSI treatments (12-14 t ha-1). This showed that differences in SOM stocks in the sandy Cambisol induced by fertilizer rate may be short-lived in case of changing management, but differences induced by fertilizer type may persist for decades. (ii) Crop yields, estimated C inputs (1.5 t ha-1 year-1) with crop residue, microbial bio¬mass C (Cmic, 118 – 150 mg kg-1), microbial biomass N (17 – 20 mg kg-1) and labile C and N pools did not differ significantly between FYM and DYN treatments. However, labile C increased linearly with application rate (R2 = 0.53) from 7 to 11% of Corg. This also applied for labile N (3.5 to 4.9% of Nt). The higher contents of Corg in DYN treatments existed since 1982, when the first sampling was conducted for all individual treatments. Contents of Corg between DYN and FYM treatments con-verged slightly since then. Furthermore, at least 30% of the difference in Corg was located in the passive pool where a treatment effect could be excluded. Therefore, the reported differences in Corg contents existed most likely since the beginning of the experiment and, as a single factor of biodynamic agriculture, application of bio-dynamic preparations had no effect on SOM stocks. (iii) Stocks of SOM, light fraction organic C (LFOC, ρ ≤ 2.0 g cm-3), light fraction organic N and Cmic decreased in the order FYMH > FYML > MSIH, MSIL for all sampling dates in 2008 (March, May, September, December). However, statistical significance of treatment effects differed between the dates, probably due to dif-ferences in the spatial variation throughout the year. The high proportion of LFOC on total Corg stocks (45 – 55%) highlighted the importance of selective preservation of OM as a stabilization mechanism in this sandy Cambisol. The apparent turnover time of LFOC was between 21 and 32 years, which agreed very well with studies with substantially longer vegetation change compared to our study. Overall, both approaches; (I) the combination of incubation, chemical fractionation and simple modelling and (II) the density fractionation; provided complementary information on the partitioning of SOM into pools of different stability. The density fractionation showed that differences in Corg stocks between FYM and MSI treatments were mainly located in the light fraction, i.e. induced by higher recalcitrance of the organic input in the FYM treatments. Moreover, the use of the combination of biological, chemical and mathematical methods indicated that effects of fertilizer rate on total Corg and Nt stocks may be short-lived, but that the effect of fertilizer type may persist for longer time spans in the sandy Cambisol.

Genetic variation of and environmental effects on inducibility of resistance in tomatoes (Solanum lycopersicum L.) to Phytophthora infestans (Mont) de Bary

Many plant strengtheners are promoted for their supposed effects on nutrient uptake and/or resistance induction (IR). In addition, many organic fertilizers are supposed to enhance plant health and several studies have shown that tomatoes grown organically are more resistant to late blight, caused by Phytophthora infestans to tomatoes grown conventionally. Much is known about the mechanisms underlying IR. In contrast, there is no systematic knowledge about genetic variation for IR. Therefore, the following questions were addressed in the presented dissertation: (i) Is there genetic variation among tomato genotypes for inducibility of resistance to P. infestans? (ii) How do different PS compare with the chemical inducer BABA in their ability to IR? (iii) Does IR interact with the inducer used and different organic fertilizers? A varietal screening showed that contrary to the commonly held belief IR in tomatoes is genotype and isolate specific. These results indicate that it should be possible to select for inducibility of resistance in tomato breeding. However, isolate specificity also suggests that there could be pathogen adaptation. The three tested PS as well as two of the three tested organic fertilisers all induced resistance in the tomatoes. Depending on PS or BABA variety and isolate effects varied. In contrast, there were no variety and isolate specific effects of the fertilisers and no interactions with the PS and fertilisers. This suggests that the different PS should work independent of the soil substrate used. In contrast the results were markedly different when isolate mixtures were used for challenge inoculations. Plants were generally less susceptible to isolate mixtures than to single isolates. In addition, the effectiveness of the PS was greater and more similar to BABA when isolate mixtures were used. The fact that the different PS and BABA differed in their ability to induce resistance in different host genotype -pathogen isolate combinations puts the usefulness of IR as a breeding goal in question. This would result in varieties depending on specific inducers. The results with the isolate mixtures are highly relevant. On the one hand they increase the effectiveness of the resistance inducers. On the other hand, measures that increase the pathogen diversity such as the use of diversified host populations will also increase the overall resistance of the hosts. For organic tomato production the results indicate that it is possible to enhance the tomato growing system with respect to plant health management by using optimal fertilisers, plant strengtheners and any measures that increase system diversity.

Soil fertility and nutrient status of traditional Gayo coffee agroforestry systems in the Takengon region, Aceh Province, Indonesia

Little is known about the traditional coffee cultivation systems in Central Aceh, Indonesia, where coffee production is a major source of income for local Gayo people. Based on field observations and farmer interviews, 14 representative agroforestry coffee plantations of different age classes (60-70 years, 30-40 years, and 20 years) as well as seven adjacent grassland and native forest sites were selected for this study, and soil and coffee leaf samples collected for nutrient analysis. Significant differences in soil and coffee leaf parameters were found between former native forest and Sumatran pine (Pinus merkusii) forest as previous land cover indicating the importance of the land use history for today’s coffee cultivation. Soil pH as well as exchangeable Na and Ca concentrations were significantly lower on coffee plantations compared to grassland and forest sites. Soil C, N, plant available P, exchangeable K, and Mg concentrations showed no consistent differences between land use groups. Nitrogen (N), phosphorus (P), and potassium (K) concentrations of coffee leaves were in the sufficiency range, whereas zinc (Zn) contents were found to be consistently below the sufficiency threshold and significantly lower in coffee plantations of previous pine forest cover compared to those of previous native forest cover. While the results of this study provided insights into the nutrient status of coffee plantations in Central Aceh, the heterogeneity of site conditions, limited sampling size, and scarcity of reliable data about the land use history and initial soil conditions of sampled sites preclude more definitive conclusions about the sustainability of the studied systems.

Yield gap analysis and assessment of climate-induced yield trends of irrigated rice in selected provinces of the Philippines

This study describes a combined empirical/modeling approach to assess the possible impact of climate variability on rice production in the Philippines. We collated climate data of the last two decades (1985-2002) as well as yield statistics of six provinces of the Philippines, selected along a North-South gradient. Data from the climate information system of NASA were used as input parameters of the model ORYZA2000 to determine potential yields and, in the next steps, the yield gaps defined as the difference between potential and actual yields. Both simulated and actual yields of irrigated rice varied strongly between years. However, no climate-driven trends were apparent and the variability in actual yields showed no correlation with climatic parameters. The observed variation in simulated yields was attributable to seasonal variations in climate (dry/wet season) and to climatic differences between provinces and agro-ecological zones. The actual yield variation between provinces was not related to differences in the climatic yield potential but rather to soil and management factors. The resulting yield gap was largest in remote and infrastructurally disfavored provinces (low external input use) with a high production potential (high solar radiation and day-night temperature differences). In turn, the yield gap was lowest in central provinces with good market access but with a relatively low climatic yield potential. We conclude that neither long-term trends nor the variability of the climate can explain current rice yield trends and that agroecological, seasonal, and management effects are over-riding any possible climatic variations. On the other hand the lack of a climate-driven trend in the present situation may be superseded by ongoing climate change in the future.

Interrelatedness of grazing livestock with vegetation parameters and farmers’ livelihoods in the Mahafaly region, southwestern Madagascar

The study aims to get deeper insight into the highly extensive system of animal husbandry in the Mahafaly region of southwestern Madagascar. It tries to understand the major drivers for pastoral dynamics, land and resource use along a gradient in altitude and vegetation to consider the area’s high spatial and temporal heterogeneity. The study also analyzes the reproductive performance of local livestock as well as the owners’ culling strategies to determine herd dynamics, opportunities for economic growth, and future potential for rural development. Across seasons, plateau herds from both livestock species covered longer distances (cattle 13.6±3.02 km, goats 12.3±3.48 km) and were found further away from the settlements (cattle 3.1±0.96 km, goats 2.8±0.98 km) than those from the coastal plain (walking_dist: cattle 9.5±3.25 km, goats 9.2±2.57 km; max_dist: cattle 2.6±1.28 km, goats 1.8±0.61 km). Transhumant cattle were detected more vulnerable through limited access to pasture land and water resources compared to local herds. Seasonal water shortage has been confirmed as a key constraint on the plateau while livestock keeping along the coast is more limited by dry season forage availability. However, recent security issues and land use conflicts with local crop farmers are gaining importance and force livestock owners to adapt their traditional grazing management, resulting in spatio-temporal variation of livestock numbers and in the impending risk of local overgrazing and degradation of rangelands. Among the 133 plant species consumed by livestock, 13 were determined of major importance for the animals’ nutrition. The nutritive value and digestibility of the natural forage, as well as its abundance in the coastal zone, substantially decreased over the course of the dry season and emphasized the importance of supplementary forage plants, in particular Euphorbia stenoclada. At the same time, an unsustainable utilization and overexploitation of its wild stocks may raise the pressure on the vegetation and pasture resources within the nearby Tsimanampetsotsa National Park. Age at first parturition was 40.5±0.59 months for cattle and 21.3±0.63 months for goats. Both species showed long parturition intervals (cattle 24.2±0.48 months, goats 12.4±0.30 months), mostly due to the maintenance of poorly performing breeding females within the herds. Reported offspring mortality, however, was low with 2.5% of cattle and 18.8% of goats dying before reaching maturity. The analysis of economic information revealed higher than expected market dynamics, especially for zebus, resulting in annual contribution margins of 33 € per cattle unit and 11 € per goat unit. The application of the PRY Herd Life model to simulate herd development for present management and two alternate scenarios confirmed the economic profitability of the current livestock system and showed potential for further productive and economic development. However, this might be clearly limited by the region’s restricted carrying capacity. Summarizing, this study illustrates the highly extensive and resources-driven character of the livestock system in the Mahafaly region, with herd mobility being a central element to cope with seasonal shortages in forage and water. But additional key drivers and external factors are gaining importance and increasingly affect migration decisions and grazing management. This leads to an increased risk of local overgrazing and overexploitation of natural pasture resources and intensifies the tension between pastoral and conservation interests. At the same time, it hampers the region’s agronomic development, which has not yet been fully exploited. The situation therefore demonstrates the need for practical improvement suggestions and implication measures, such as the systematic forestation of supplemental forage plant species in the coastal zone or a stronger integration of animal husbandry and crop production, to sustain the traditional livestock system without compromising peoples’ livelihoods while at the same time minimizing the pastoral impact on the area’s unique nature and environment.

Modeling, Estimation, and Control of Robot-Soil Interactions

This thesis presents the development of hardware, theory, and experimental methods to enable a robotic manipulator arm to interact with soils and estimate soil properties from interaction forces. Unlike the majority of robotic systems interacting with soil, our objective is parameter estimation, not excavation. To this end, we design our manipulator with a flat plate for easy modeling of interactions. By using a flat plate, we take advantage of the wealth of research on the similar problem of earth pressure on retaining walls. There are a number of existing earth pressure models. These models typically provide estimates of force which are in uncertain relation to the true force. A recent technique, known as numerical limit analysis, provides upper and lower bounds on the true force. Predictions from the numerical limit analysis technique are shown to be in good agreement with other accepted models. Experimental methods for plate insertion, soil-tool interface friction estimation, and control of applied forces on the soil are presented. In addition, a novel graphical technique for inverting the soil models is developed, which is an improvement over standard nonlinear optimization. This graphical technique utilizes the uncertainties associated with each set of force measurements to obtain all possible parameters which could have produced the measured forces. The system is tested on three cohesionless soils, two in a loose state and one in a loose and dense state. The results are compared with friction angles obtained from direct shear tests. The results highlight a number of key points. Common assumptions are made in soil modeling. Most notably, the Mohr-Coulomb failure law and perfectly plastic behavior. In the direct shear tests, a marked dependence of friction angle on the normal stress at low stresses is found. This has ramifications for any study of friction done at low stresses. In addition, gradual failures are often observed for vertical tools and tools inclined away from the direction of motion. After accounting for the change in friction angle at low stresses, the results show good agreement with the direct shear values.

Model Order Reduction for Determining Bubble Parameters to Attain a Desired Fluid Surface Shape

In this paper, a new methodology for predicting fluid free surface shape using Model Order Reduction (MOR) is presented. Proper Orthogonal Decomposition combined with a linear interpolation procedure for its coefficient is applied to a problem involving bubble dynamics near to a free surface. A model is developed to accurately and efficiently capture the variation of the free surface shape with different bubble parameters. In addition, a systematic approach is developed within the MOR framework to find the best initial locations and pressures for a set of bubbles beneath the quiescent free surface such that the resultant free surface attained is close to a desired shape. Predictions of the free surface in two-dimensions and three-dimensions are presented.

The relationship between diffuse spectral reflectance of the soil and its cation exchange capacity is scale-dependent

Diffuse reflectance spectroscopy (DRS) is increasingly being used to predict numerous soil physical, chemical and biochemical properties. However, soil properties and processes vary at different scales and, as a result, relationships between soil properties often depend on scale. In this paper we report on how the relationship between one such property, cation exchange capacity (CEC), and the DRS of the soil depends on spatial scale. We show this by means of a nested analysis of covariance of soils sampled on a balanced nested design in a 16 km × 16 km area in eastern England. We used principal components analysis on the DRS to obtain a reduced number of variables while retaining key variation. The first principal component accounted for 99.8% of the total variance, the second for 0.14%. Nested analysis of the variation in the CEC and the two principal components showed that the substantial variance components are at the > 2000-m scale. This is probably the result of differences in soil composition due to parent material. We then developed a model to predict CEC from the DRS and used partial least squares (PLS) regression do to so. Leave-one-out cross-validation results suggested a reasonable predictive capability (R2 = 0.71 and RMSE = 0.048 molc kg− 1). However, the results from the independent validation were not as good, with R2 = 0.27, RMSE = 0.056 molc kg− 1 and an overall correlation of 0.52. This would indicate that DRS may not be useful for predictions of CEC. When we applied the analysis of covariance between predicted and observed we found significant scale-dependent correlations at scales of 50 and 500 m (0.82 and 0.73 respectively). DRS measurements can therefore be useful to predict CEC if predictions are required, for example, at the field scale (50 m). This study illustrates that the relationship between DRS and soil properties is scale-dependent and that this scale dependency has important consequences for prediction of soil properties from DRS data

Spatial covariation of Azotobacter abundance and soil properties: A case study using the wavelet transform

The soil microflora is very heterogeneous in its spatial distribution. The origins of this heterogeneity and its significance for soil function are not well understood. A problem for understanding spatial variation better is the assumption of statistical stationarity that is made in most of the statistical methods used to assess it. These assumptions are made explicit in geostatistical methods that have been increasingly used by soil biologists in recent years. Geostatistical methods are powerful, particularly for local prediction, but they require the assumption that the variability of a property of interest is spatially uniform, which is not always plausible given what is known about the complexity of the soil microflora and the soil environment. We have used the wavelet transform, a relatively new innovation in mathematical analysis, to investigate the spatial variation of abundance of Azotobacter in the soil of a typical agricultural landscape. The wavelet transform entails no assumptions of stationarity and is well suited to the analysis of variables that show intermittent or transient features at different spatial scales. In this study, we computed cross-variograms of Azotobacter abundance with the pH, water content and loss on ignition of the soil. These revealed scale-dependent covariation in all cases. The wavelet transform also showed that the correlation of Azotobacter abundance with all three soil properties depended on spatial scale, the correlation generally increased with spatial scale and was only significantly different from zero at some scales. However, the wavelet analysis also allowed us to show how the correlation changed across the landscape. For example, at one scale Azotobacter abundance was strongly correlated with pH in part of the transect, and not with soil water content, but this was reversed elsewhere on the transect. The results show how scale-dependent variation of potentially limiting environmental factors can induce a complex spatial pattern of abundance in a soil organism. The geostatistical methods that we used here make assumptions that are not consistent with the spatial changes in the covariation of these properties that our wavelet analysis has shown. This suggests that the wavelet transform is a powerful tool for future investigation of the spatial structure and function of soil biota. (c) 2006 Elsevier Ltd. All rights reserved.

Incorporating Legacy Soil pH Databases into Digital Soil Maps

Soil data and reliable soil maps are imperative for environmental management. conservation and policy. Data from historical point surveys, e.g. experiment site data and farmers fields can serve this purpose. However, legacy soil information is not necessarily collected for spatial analysis and mapping such that the data may not have immediately useful geo-references. Methods are required to utilise these historical soil databases so that we can produce quantitative maps of soil propel-ties to assess spatial and temporal trends but also to assess where future sampling is required. This paper discusses two such databases: the Representative Soil Sampling Scheme which has monitored the agricultural soil in England and Wales from 1969 to 2003 (between 400 and 900 bulked soil samples were taken annually from different agricultural fields); and the former State Chemistry Laboratory, Victoria, Australia where between 1973 and 1994 approximately 80,000 soil samples were submitted for analysis by farmers. Previous statistical analyses have been performed using administrative regions (with sharp boundaries) for both databases, which are largely unrelated to natural features. For a more detailed spatial analysis that call be linked to climate and terrain attributes, gradual variation of these soil properties should be described. Geostatistical techniques such as ordinary kriging are suited to this. This paper describes the format of the databases and initial approaches as to how they can be used for digital soil mapping. For this paper we have selected soil pH to illustrate the analyses for both databases.

The spatial prediction of soil mineral N and potentially available N using elevation

The precision farmer wants to manage the variation in soil nutrient status continuously, which requires reliable predictions at places between sampling sites. Ordinary kriging can be used for prediction if the data are spatially dependent and there is a suitable variogram model. However, even if data are spatially correlated, there are often few soil sampling sites in relation to the area to be managed. If intensive ancillary data are available and these are coregionalized with the sparse soil data, they could be used to increase the accuracy of predictions of the soil properties by methods such as cokriging, kriging with external drift and regression kriging. This paper compares the accuracy of predictions of the plant available N properties (mineral N and potentially available N) for two arable fields in Bedfordshire, United Kingdom, from ordinary kriging, cokriging, kriging with external drift and regression kriging. For the last three, intensive elevation data were used with the soil data. The mean squared errors of prediction from these methods of kriging were determined at validation sites where the values were known. Kriging with external drift resulted in the smallest mean squared error for two of the three properties examined, and cokriging for the other. The results suggest that the use of intensive ancillary data can increase the accuracy of predictions of soil properties in arable fields provided that the variables are related spatially. (c) 2005 Elsevier B.V. All rights reserved.