Simultaneous EEG-fMRI at ultra-high field: Artifact prevention and safety assessment.

The simultaneous recording of scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can provide unique insights into the dynamics of human brain function, and the increased functional sensitivity offered by ultra-high field fMRI opens exciting perspectives for the future of this multimodal approach. However, simultaneous recordings are susceptible to various types of artifacts, many of which scale with magnetic field strength and can seriously compromise both EEG and fMRI data quality in recordings above 3T. The aim of the present study was to implement and characterize an optimized setup for simultaneous EEG-fMRI in humans at 7T. The effects of EEG cable length and geometry for signal transmission between the cap and amplifiers were assessed in a phantom model, with specific attention to noise contributions from the MR scanner coldheads. Cable shortening (down to 12cm from cap to amplifiers) and bundling effectively reduced environment noise by up to 84% in average power and 91% in inter-channel power variability. Subject safety was assessed and confirmed via numerical simulations of RF power distribution and temperature measurements on a phantom model, building on the limited existing literature at ultra-high field. MRI data degradation effects due to the EEG system were characterized via B0 and B1(+) field mapping on a human volunteer, demonstrating important, although not prohibitive, B1 disruption effects. With the optimized setup, simultaneous EEG-fMRI acquisitions were performed on 5 healthy volunteers undergoing two visual paradigms: an eyes-open/eyes-closed task, and a visual evoked potential (VEP) paradigm using reversing-checkerboard stimulation. EEG data exhibited clear occipital alpha modulation and average VEPs, respectively, with concomitant BOLD signal changes. On a single-trial level, alpha power variations could be observed with relative confidence on all trials; VEP detection was more limited, although statistically significant responses could be detected in more than 50% of trials for every subject. Overall, we conclude that the proposed setup is well suited for simultaneous EEG-fMRI at 7T.

From MuSIASEM theory to practice: report and reflections from field research in Kampot province, Cambodia

This document contains a report and summary of the field research activities in a rural community of rice farmers in Kampot province, Cambodia in 2011, which I conducted within the context of my PhD research at ICTA-UAB (Institute of Environmental Science and Technology, Autonomous University of Barcelona, Spain). The purpose of the field research was to gather data for a MuSIASEM analysis (Multi-Scale Integrated Analysis of Societal and Ecosystem Metabolism) at the village and household level, in order to analyze the multidimensional challenges that small farmers may face nowadays within the context of global rural change and declining access to land. While the literature on MuSIASEM offers a great variety of theoretical explanations and practical applications, there is little information available for students regarding the practical steps required for doing a MuSIASEM analysis at the local level. Within this context, this report offers not only a documentation of the field research design and data collection methods, but further provides a general overview on some organizational and preparative aspects, including some personal reflections, that one may face when preparing and conducting field research for MuSIASEM analysis. In summary, this document thus serves three objectives: (i) to assure methodological transparency for the future work, based on the collected data during field research, (ii) to share my personal experience on the preparative and practical steps required for field research and data collection for a MuSIASEM analysis at the local level, and (iii) to make available for the further interested reader some more detailed background information on the case study village.

Simple finite field method for calculation of static and dynamic vibrational hyperpolarizabilities: curvature contributions

In the static field limit, the vibrational hyperpolarizability consists of two contributions due to: (1) the shift in the equilibrium geometry (known as nuclear relaxation), and (2) the change in the shape of the potential energy surface (known as curvature). Simple finite field methods have previously been developed for evaluating these static field contributions and also for determining the effect of nuclear relaxation on dynamic vibrational hyperpolarizabilities in the infinite frequency approximation. In this paper the finite field approach is extended to include, within the infinite frequency approximation, the effect of curvature on the major dynamic nonlinear optical processes

Infestation and natural parasitism of aphids in single and mixed pastures of black oats and ryegrass

Some species of aphids are major pests on cereal crops and grass pastures. Usually these pests are not adequately controlled in pasture lands that become sources of aphid infestations to cereal crops. The dynamics of aphids and the incidence of natural enemies are less known in pasture systems than in cereal fields. The objective of this work was to assess the aphid infestation and natural aphid parasitism in different pasture composition. Three hypotheses were tested: 1- the aphid species composition in pastures may vary according to the cereal species in the field; 2- the mixture of two plant species can modify the amount and diversity of aphids; 3- the plant species composition of pasture fields influences the parasitism of aphids. Empirical data were obtained from three Poaceae fields: black oats (Avena strigosa L.), ryegrass (Lolium multiflorum L.), and a mixed field of black oats and ryegrass. The most abundant aphid species was Rhopalosiphum padi followed by Sitobion avenae. Plant species composition increases the amount and the parasitism rates of aphids. The mixture of heavily infested black oats with a poorly infested ryegrass resulted in reduction of aphid infestation in comparison with heavily-infested single plant species field. This is possible because the conditions are favorable for the development of parasitoid populations. Aphidius colemani was the main parasitoid found in these areas.

Biological aspects of Leucothyreus alvarengai Frey and Leucothyreus aff. semipruinosus Ohaus (Coleoptera, Melolonthidae, Rutelinae) in crop succession at central Brazil

Biological aspects of Leucothyreus alvarengai Frey and Leucothyreus aff. semipruinosus Ohaus (Coleoptera, Melolonthidae, Rutelinae) in crop succession at central Brazil. Beetles of the family Melolonthidae make up a large group and some species are considered pests of planted crops. Little information is available on the basic biological aspects of the genus Leucothyreus, such as association with cultivated crops and their occurrence periods. Therefore studies were developed in soybean and corn crops in Tangará da Serra, Mato Grosso, Brazil, with the objective of studying the occurrence and biological aspects of Leucothyreus alvarengai Frey and Leucothyreus aff. semipruinosus Ohaus. For acquisition of immature specimens of both species, in April 2011 sampling was performed in corn fields, in July and October in the fallow area, and in soybeans fields planted in December; in 2012 sampling was performed in January and February in soybean fields and in March in corn fields. In 2011 the total number of larvae obtained in April, July, October and December were 100, 6, 30 and 27, and in January, February and March of 2012 these quantities were 32, 52 and 65 larvae, respectively. In all sampling events the larvae of L. alvarengai were collected in greater quantity. At the beginning of the reproductive period of L. alvarengai and L. aff. semipruinosus, it was observed that the adults began to fly and soon after started oviposition in the field in September. The appearance of larvae coincides with the time of soybean planting in the field, thus the larvae feed on roots of soybean plants at the beginning of their development and the cycle from egg to adult of the two species was completed in one year.

Surface topography and mixed-pixel effects on the simulated L-band brightness temperatures.

The impact of topography and mixed pixels on L-band radiometric observations over land needs to be quantified to improve the accuracy of soil moisture retrievals. For this purpose, a series of simulations has been performed with an improved version of the soil moisture and ocean salinity (SMOS) end-to-end performance simulator (SEPS). The brightness temperature generator of SEPS has been modified to include a 100-m-resolution land cover map and a 30-m-resolution digital elevation map of Catalonia (northeast of Spain). This high-resolution generator allows the assessment of the errors in soil moisture retrieval algorithms due to limited spatial resolution and provides a basis for the development of pixel disaggregation techniques. Variation of the local incidence angle, shadowing, and atmospheric effects (up- and downwelling radiation) due to surface topography has been analyzed. Results are compared to brightness temperatures that are computed under the assumption of an ellipsoidal Earth.

Surface and subsurface decomposition of a desiccated grass pasture biomass related to erosion and its prediction with RUSLE

Erosion is deleterious because it reduces the soil's productivity capacity for growing crops and causes sedimentation and water pollution problems. Surface and buried crop residue, as well as live and dead plant roots, play an important role in erosion control. An efficient way to assess the effectiveness of such materials in erosion reduction is by means of decomposition constants as used within the Revised Universal Soil Loss Equation - RUSLE's prior-land-use subfactor - PLU. This was investigated using simulated rainfall on a 0.12 m m-1 slope, sandy loam Paleudult soil, at the Agriculture Experimental Station of the Federal University of Rio Grande do Sul, in Eldorado do Sul, State of Rio Grande do Sul, Brazil. The study area had been covered by native grass pasture for about fifteen years. By the middle of March 1996, the sod was mechanically mowed and the crop residue removed from the field. Late in April 1996, the sod was chemically desiccated with herbicide and, about one month later, the following treatments were established and evaluated for sod biomass decomposition and soil erosion, from June 1996 to May 1998, on duplicated 3.5 x 11.0 m erosion plots: (a) and (b) soil without tillage, with surface residue and dead roots; (c) soil without tillage, with dead roots only; (d) soil tilled conventionally every two-and-half months, with dead roots plus incorporated residue; and (e) soil tilled conventionally every six months, with dead roots plus incorporated residue. Simulated rainfall was applied with a rotating-boom rainfall simulator, at an intensity of 63.5 mm h-1 for 90 min, eight to nine times during the experimental period (about every two-and-half months). Surface and subsurface sod biomass amounts were measured before each rainfall test along with the erosion measurements of runoff rate, sediment concentration in runoff, soil loss rate, and total soil loss. Non-linear regression analysis was performed using an exponential and a power model. Surface sod biomass decomposition was better depicted by the exponential model, while subsurface sod biomass was by the power model. Subsurface sod biomass decomposed faster and more than surface sod biomass, with dead roots in untilled soil without residue on the surface decomposing more than dead roots in untilled soil with surface residue. Tillage type and frequency did not appreciably influence subsurface sod biomass decomposition. Soil loss rates increased greatly with both surface sod biomass decomposition and decomposition of subsurface sod biomass in the conventionally tilled soil, but they were minimally affected by subsurface sod biomass decomposition in the untilled soil. Runoff rates were little affected by the studied treatments. Dead roots plus incorporated residues were effective in reducing erosion in the conventionally tilled soil, while consolidation of the soil surface was important in no-till. The residual effect of the turned soil on erosion diminished gradually with time and ceased after two years.

Efecto de los usos del suelo en la emisión de dióxido de carbono del suelo a la atmósfera en un agroecosistema semiárido del Valle del Ebro

Los cambios en los usos del suelo han contribuido de manera importante al incremento de gases de efecto invernadero en la atmósfera, especialmente de dióxido de carbono, aumentando sus emisiones desde 1970 en un 80%. Estos cambios causan la alteración de los suelos provocando un impacto sobre el ciclo del carbono, aumentando las tasas de descomposición de la fracción orgánica creando así un flujo de CO2 a la atmosfera. Entre las recomendaciones del Panel Intergubernamental de expertos sobre el Cambio Climático (IPCC, en inglés), y contemplado en el Protocolo de Kyoto, se encuentra el proceso de secuestro de carbono en suelos, que implica la eliminación del CO2 atmosférico por parte de las plantas y su almacenamiento como materia orgánica del suelo. Para poder favorecer dicho proceso, en un determinado tipo de ecosistema, es fundamental conocer cuáles son los factores que gobiernan la respiración del suelo y el impacto que tienen los diferentes usos en la emisión de CO2. En el presente trabajo se han estudiado 4 usos del suelo representativos del secano aragonés: un monocultivo de cebada en siembra directa (NT), un suelo abandonado labrado (AC), un suelo abandonado no alterado (AU) y un suelo forestal (FR) con el objetivo de conocer sus tasas de respiración, la influencia de diferentes parámetros edáficos en ellas, y proponer cambios en el uso del suelo que ayuden a mitigar estas emisiones. Además, se ha dedicado un apartado para conocer cómo influyen diferentes técnicas de fertilización nitrogenada (mineral y orgánica) en la respiración de un monocultivo de cebada en siembra directa. En cuanto a los usos, los resultados obtenidos tanto in situ como en laboratorio muestran una mayor respiración en AC, siendo los valores más bajos los de NT y FR. Una de las principales conclusiones es que la supresión del laboreo y del periodo de barbecho largo, así como la conversión de tierras abandonadas y marginales a cultivos y zonas forestales se presentan en este tipo de ecosistemas como prácticas de secuestro de carbono. En el estudio de aplicación de fertilizantes, no se observó ningún cambio en la respiración del suelo después de la aplicación de nitrógeno mineral. En cambio, el suelo fertilizado con purín sí que mostró picos de emisión durante las siguientes horas a la incorporación de éste, debido fundamentalmente a su alto contenido de carbono lábil.

The combined use of reflectance, emissivity and elevation Aster/Terra data for tropical soil studies

Reflectance, emissivity and elevation data of the sensor ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer)/Terra were used to characterize soil composition variations according to the toposequence position. Normalized data of SWIR (shortwave infrared) reflectance and TIR (thermal infrared) emissivity, coupled to a soil-fraction image from a spectral mixture model, were evaluated to separate bare soils from nonphotosynthetic vegetation. Regression relationships of some soil properties with reflectance and emissivity data were then applied on the exposed soil pixels. The resulting estimated values were plotted on the ASTER-derived digital elevation model. Results showed that the SWIR bands 5 and 6 and the TIR bands 10 and 14 measured the clay mineral absorption band and the quartz emissivity feature, respectively. These bands improved also the discrimination between nonphotosynthetic vegetation and soils. Despite the differences in pixel size and field sampling size, some soil properties were correlated with reflectance (R² of 0.65 for Al2O3 in band 6; 0.61 for Fe2O3 in band 3) and emissivity (R² of 0.65 for total sand fraction in the 10/14 band ratio). The combined use of reflectance, emissivity and elevation data revealed variations in soil composition with topography in specific parts of the landscape. From higher to lower slope positions, a general decrease in Al2O3 and increase in total sand fraction was observed, due to the prevalence of Rhodic Acrustox at the top and its gradual transition to Typic Acrustox at the bottom.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Improvement of a testing apparatus for dynamometry: procedures for penetrometry and influence of strain rate to quantify the tensile strength of soil aggregates

Soil penetration resistance (PR) and the tensile strength of aggregates (TS) are commonly used to characterize the physical and structural conditions of agricultural soils. This study aimed to assess the functionality of a dynamometry apparatus by linear speed and position control automation of its mobile base to measure PR and TS. The proposed equipment was used for PR measurement in undisturbed samples of a clayey "Nitossolo Vermelho eutroférrico" (Kandiudalfic Eutrudox) under rubber trees sampled in two positions (within and between rows). These samples were also used to measure the volumetric soil water content and bulk density, and determine the soil resistance to penetration curve (SRPC). The TS was measured in a sandy loam "Latossolo Vermelho distrófico" (LVd) - Typic Haplustox - and in a very clayey "Nitossolo Vermelho distroférrico" (NVdf) - Typic Paleudalf - under different uses: LVd under "annual crops" and "native forest", NVdf under "annual crops" and "eucalyptus plantation" (> 30 years old). To measure TS, different strain rates were applied using two dynamometry testing devices: a reference machine (0.03 mm s-1), which has been widely used in other studies, and the proposed equipment (1.55 mm s-1). The determination coefficient values of the SRPC were high (R² > 0.9), regardless of the sampling position. Mean TS values in LVd and NVdf obtained with the proposed equipment did not differ (p > 0.05) from those of the reference testing apparatus, regardless of land use and soil type. Results indicate that PR and TS can be measured faster and accurately by the proposed procedure.

Effect of compaction on microbial activity and carbon and nitrogen transformations in two oxisols with different mineralogy

The use of machinery in agricultural and forest management activities frequently increases soil compaction, resulting in greater soil density and microporosity, which in turn reduces hydraulic conductivity and O2 and CO2 diffusion rates, among other negative effects. Thus, soil compaction has the potential to affect soil microbial activity and the processes involved in organic matter decomposition and nutrient cycling. This study was carried out under controlled conditions to evaluate the effect of soil compaction on microbial activity and carbon (C) and nitrogen (N) mineralization. Two Oxisols with different mineralogy were utilized: a clayey oxidic-gibbsitic Typic Acrustox and a clayey kaolinitic Xantic Haplustox (Latossolo Vermelho-Amarelo ácrico - LVA, and Latossolo Amarelo distrófico - LA, respectively, in the Brazil Soil Classification System). Eight treatments (compaction levels) were assessed for each soil type in a complete block design, with six repetitions. The experimental unit consisted of PVC rings (height 6 cm, internal diameter 4.55 cm, volume 97.6 cm³). The PVC rings were filled with enough soil mass to reach a final density of 1.05 and 1.10 kg dm-3, respectively, in the LVA and LA. Then the soil samples were wetted (0.20 kg kg-1 = 80 % of field capacity) and compacted by a hydraulic press at pressures of 0, 60, 120, 240, 360, 540, 720 and 900 kPa. After soil compression the new bulk density was calculated according to the new volume occupied by the soil. Subsequently each PVC ring was placed within a 1 L plastic pot which was then tightly closed. The soils were incubated under aerobic conditions for 35 days and the basal respiration rate (CO2-C production) was estimated in the last two weeks. After the incubation period, the following soil chemical and microbiological properties were detremined: soil microbial biomass C (C MIC), total soil organic C (TOC), total N, and mineral N (NH4+-N and NO3--N). After that, mineral N, organic N and the rate of net N mineralization was calculated. Soil compaction increased NH4+-N and net N mineralization in both, LVA and LA, and NO3--N in the LVA; diminished the rate of TOC loss in both soils and the concentration of NO3--N in the LA and CO2-C in the LVA. It also decreased the C MIC at higher compaction levels in the LA. Thus, soil compaction decreases the TOC turnover probably due to increased physical protection of soil organic matter and lower aerobic microbial activity. Therefore, it is possible to conclude that under controlled conditions, the oxidic-gibbsitic Oxisol (LVA) was more susceptible to the effects of high compaction than the kaolinitic (LA) as far as organic matter cycling is concerned; and compaction pressures above 540 kPa reduced the total and organic nitrogen in the kaolinitic soil (LA), which was attributed to gaseous N losses.

Hydraulic conductivity and water retention in leptosols-regosols and saprolite derived from sandstone, Brazil

Leptosols and Regosols are soils with a series of restrictions for use, mainly related to the effective depth, which have been poorly studied in Brazil. These soils, when derived from sedimentary rocks should be treated with particular care to avoid environmental damage such as aquifer contamination. The purpose of this study was to verify the behavior of hydraulic conductivity and water retention capacity in profiles of Leptosols and Regosols derived from sandstone of the Caturrita formation in Rio Grande do Sul state. The morphology, particle size distribution, porosity, soil density (Ds), saturated hydraulic conductivity (Ks), basic water infiltration in the field (BI) and water retention were determined in soil and saprolite samples of six soil profiles. High Ds, low macroporosity and high microporosity were observed in the profiles, resulting in a low Ks and BI, even under conditions of sandy texture and a highly fractured saprolite layer. The variation coefficients of data of Ks and BI were high among the studied profiles and between replications of a same profile. Water retention of the studied soils was higher in Cr layers than in the A horizons and the volume of plant-available water greater and variable among A horizons and Cr layers.

N2O emissions from a cultivated mollisol: optimal time of day for sampling and the role of soil temperature

The correct use of closed field chambers to determine N2O emissions requires defining the time of day that best represents the daily mean N2O flux. A short-term field experiment was carried out on a Mollisol soil, on which annual crops were grown under no-till management in the Pampa Ondulada of Argentina. The N2O emission rates were measured every 3 h for three consecutive days. Fluxes ranged from 62.58 to 145.99 ∝g N-N2O m-2 h-1 (average of five field chambers) and were negatively related (R² = 0.34, p < 0.01) to topsoil temperature (14 - 20 ºC). N2O emission rates measured between 9:00 and 12:00 am presented a high relationship to daily mean N2O flux (R² = 0.87, p < 0.01), showing that, in the study region, sampling in the mornings is preferable for GHG.

Comparison between detailed digital and conventional soil maps of an area with complex geology

Since different pedologists will draw different soil maps of a same area, it is important to compare the differences between mapping by specialists and mapping techniques, as for example currently intensively discussed Digital Soil Mapping. Four detailed soil maps (scale 1:10.000) of a 182-ha sugarcane farm in the county of Rafard, São Paulo State, Brazil, were compared. The area has a large variation of soil formation factors. The maps were drawn independently by four soil scientists and compared with a fifth map obtained by a digital soil mapping technique. All pedologists were given the same set of information. As many field expeditions and soil pits as required by each surveyor were provided to define the mapping units (MUs). For the Digital Soil Map (DSM), spectral data were extracted from Landsat 5 Thematic Mapper (TM) imagery as well as six terrain attributes from the topographic map of the area. These data were summarized by principal component analysis to generate the map designs of groups through Fuzzy K-means clustering. Field observations were made to identify the soils in the MUs and classify them according to the Brazilian Soil Classification System (BSCS). To compare the conventional and digital (DSM) soil maps, they were crossed pairwise to generate confusion matrices that were mapped. The categorical analysis at each classification level of the BSCS showed that the agreement between the maps decreased towards the lower levels of classification and the great influence of the surveyor on both the mapping and definition of MUs in the soil map. The average correspondence between the conventional and DSM maps was similar. Therefore, the method used to obtain the DSM yielded similar results to those obtained by the conventional technique, while providing additional information about the landscape of each soil, useful for applications in future surveys of similar areas.