Water-soluble organic matter in incipient podzols: accumulation in B horizons or in fibres?

To elucidate the mechanism of podzolization in its first stages we studied the fate of the water-soluble organic matter (WSOM) in incipient podzols in sandy soils by comparing the composition of the WSOM from L, F and H horizons with that in the bulk of the Bh horizons and fibres of three profiles. The WSOM appeared to consist significantly of ligno-cellulose and proteins, but these biopolymers were hardly present in the Bh horizons. The material of the fibres, however, greatly resembled the WSOM composition, thereby suggesting that in these soils most of the WSOM is transported through the B horizon and accumulates hardly changed in thin bands where the water stops moving. This implies that in the early steps of podzolization, accumulation of organic matter in the B horizon is not likely to be caused by water-soluble material.     

Bioavailability of β-cryptoxanthin in the presence of phytosterols: in vitro and in vivo studies

Bioactive compounds are used in the design and development of new food products with potential health benefits, although little is known regarding their bioavailability and interactions. This study assessed the stability, in vitro bioaccessibility, and human bioavailability of β-cryptoxanthin from β-cryptoxanthin-rich drinks with and without added phytosterols developed for this purpose. The developed drinks showed no difference in the content of β-cryptoxanthin, and they were stable over 6 months. In vitro, hydrolysis of β-cryptoxanthin esters and the amount of free β-cryptoxanthin at duodenal and micellar phases were similar regardless of the presence of phytosterols. In the human study, the daily intake provoked significant increments of β-cryptoxanthin in serum regardless of the type of the drink. In conclusion, in vitro and in vivo human studies have shown that the bioavailability of β-cryptoxanthin is not significantly affected by the presence of phytosterols when they are simultaneously supplied in a drink.     

Difference in selenium accumulation in shoots of two rice cultivars

Two japonica rice （Oryza sativa L.） cultivars, Xiushui 48 and S. Andrea, differing in their ability to accumulate Se in the grain （as high as a three-fold difference）, were compared for selenium （Se） accumulation in their shoots when their growth media was supplied with different forms of Se. Results indicated that when treated with 0.25μmol L^-1 Na2SeO3, Xiushui 48 accumulation of Se in the shoots was significantly more rapid （P〈0.05） than S. Andrea, probably because of greater Se uptake and transport in Xiushui 48. Xiushui 48 rice seedlings had a higher shoot-Se accumulation rate and absorbed selenocysteine （Se-Cys） more rapidly than S. Andrea seedlings. However, when treated with Se as 0.25μmol L^-1 selenomethionine （Se-Met）, the S. Andrea seedlings＇ accumulation rate was significantly greater （P〈0.05） than that of Xiushui 48. Possibly, the high Se accumulation rate of Xiushui 48 seedling shoots compared to S. Andrea shoots was the result of a higher capacity of Xiushui 48 to transform selenite to organic Se compounds and a higher selenite uptake rate.     

Metal contamination in Nullah Dek water and accumulation in rice

Corn field experiments with two treatments, NP and NPK, where N in the form of urea, P in the form of calcium phosphate, and K in the form of KCl were applied at rates of 187.5, 33.3, and 125 kg ha^-1, respectively, on soils derived from Quaternary red clay were conducted in the hilly red soil region of Zhejiang Province, China. Plant grains and stalks were collected for determination of K content. Seven equations were used to describe the kinetics of K release from surface soil samples taken before the corn experiments under electric field strengths of 44.4 and 88.8 V cm^-1 by means of electro-ultrafiltration （EUF） and to determine if their parameters had a practical application. The second-order and Elovich equations excellently described K release; the first-order, power function, and parabolic diffusion equations also described K release well; but the zero-order and exponential equations were not so good at reflecting K release. Five reference standards from the field experiments, including relative grain yield （yield of the NP treatment/yield of the NPK treatment）, relative dry matter yield （dry matter of the NP treatment/dry matter of the NPK treatment）, quantity of K uptake in the NP treatment （no K application）, soil exchangeable K, and soil HNO3-soluble K, were used to test the effectiveness of equation parameters obtained from the slope or intercept of these equations. Correlations of the ymax （the maximum desorbable quantity of K） in the second-order equation and the constant b in the first-order and Elovich equations to all five reference standards were highly significant （P ≤ 0.01）. The constant a in the power function equation was highly significant （P 〈 0.01） for four of the five reference standards with the fifth being significant （P ≤ 0.05）. The constant b in the parabolic equation was also significantly correlated （P ≤ 0.05） to the relative grain yield and soil HNO3-soluble K. These suggested that all of these parameters could be used to estimate the soil K supplying capacity and the crop response to K fertilizer.     

Metal Contamination in Nullah Dek Water and Accumulation in Rice

A research study was carried out to determine the electrical conductivity （EC）, residual sodium carbonate （RSC）, sodium adsorption ratio （SAR）, pH and metals in metal-polluted irrigation water from a nullah and those in soils over a period of time, and the effect of metals on rice yield and metal concentrations in rice grain and straw. Two sites （I and If） were selected on the bank of Nullah Dek at Thatta Wasiran in Sheikhupura District （Pakistan）, with two rice varieties, Super Basmati and Basmati 385, at both sites. Water samples were collected during rice crop growth at 15-day intervals from August 3 to November 1, 2002. The results showed that Nullah Dek water had an EC 〉1.0 dS m^-1 and RSC of 2.78-4.11 mmolc L^-1, which was hazardous for crops, but the SAR was within the safe limit. Cu, Mn Cd and Sr were also within safe limits. The soil analysis showed that Site Ⅱ was free from salinity/sodicity, whereas Site Ⅰ was saline sodic. Among metals, Zn was sometimes deficient, Cu, Mn and Fe were adequate, and St, Ni and Cd were within safe limits in the soil at both the sites. After the rice crop harvest, concentrations of all metals tested were usually slightly increased, being higher in the upper soil layer than the lower. In addition, Basmati 385 produced higher rice grain and straw yield than Super Basmati. Chemical analysis of rice grain indicated the presence of Zn, Cu, Fe, Mn, Pb and Sr, whereas rice straw contained Zn, Cu, Fe, Mn and Sr, with Cd and Ni both being found in minute quantities.     

Increases in phosphatase and β‐glucosidase activities in wheat seedlings in response to phosphorus‐deficient growth

The ability of plants to utilize P efficiently is important for crops growing in P‐deficient soils or on soils with a high P‐fixing capacity. The purpose of this work was to investigate early physiological changes which occur when wheat (Triticum aestivum L.) seedlings were grown under P‐deficient conditions. Wheat plants were grown in a greenhouse and watered with nutrient solution containing or lacking P. During the interval 12 to 18 days after planting, the dry weight of wheat seedlings was similar regardless of P treatment, although the P‐deficient plants had a greater proportion of the total plant weight in the roots. Sixteen days after planting, the roots and leaves of P‐deficient plants had only 20 to 30% the P content of P‐sufficient plants. After 16 days, plants grown under P stress had 41% more p‐nitrophenol phosphatase activity and 70% more β‐glucosidase activity in shoot homogenates than was found in P‐sufficient plants. Changes in both enzyme activities may be involved in the mobilization of plant resources during the early stages of P‐deficient growth.     

Manganese‐iron relationship in soybean grown in calcareous soils 1

Although a positive response to iron (Fe) is, usually, expected in calcareous soils; this has not been always the case; and in some instances a depressing effect has been observed. An induced micronutrient imbalance is suspected. This experiment was designed to study the effect of Fe fertilizer on the plant micronutrients. Twenty three highly calcareous soils (18–46% calcium carbonate equivalent; pH 7.7–8.4; and a wide range of extractable Fe) from southern Iran were used in an eight‐week greenhouse experiment to study the effect of Fe fertilizers on soybean [Glycine max (L.) Merr.] growth and chemical composition. The statistical design was a 23 × 3 factorial arranged in a completely randomized block with three replications. Treatments consisted of 23 soils and three levels of applied Fe (0, 10, and 20 mg Fe/kg as FeEDDHA). Uniform doses of nitrogen (N), phosphorus (P), copper (Cu), manganese (Mn), and zinc (Zn) were applied to all pots. Dry matter (DM) and micronutrients concentrations and uptakes of plant tops were determined and used as the plant responses. Application of Fe either had no significant effect on DM or even decreased it. The plant concentration and uptake of Fe increased significantly in all soils. The concentrations and uptakes of Cu and Zn did not change but those of Mn decreased significantly. The negative effect of Fe application was, therefore, attributed to the interference of Fe with Mn nutrition. The mechanism involved appears to be the restriction in Mn translocation from soil to root and/or from root to the plant tops.     

Diatoms in peat - Dominant producers in a changing environment?

Changes in hydrology and temperature can induce rapid changes in boreal wetland ecosystems. Factors such as hydrosere, permafrost, climate and human interference may disturb the prevailing mire vegetation, whereby a new dominant assemblage can develop. At the transition from one vegetation type to another, the old vegetation may be suppressed, die out or start to decay, and some time may pass until a new mire vegetation is fully established. Here, we demonstrate that diatoms may thrive during such transitions, creating isolated and shallow peat layers with significantly elevated biogenic silica content. Biogenic silica and other nutrients that would otherwise be lost during mineralization in runoff are in this way retained in the ecosystem. Our results imply that silica storage originating from diatoms can be expected to increase in today's rapidly changing boreal wetlands. The impacts on transport of Si through boreal watersheds are currently unknown.     

Some advances in important research on soil in a changing world

This special issue, for celebrating the 30th Anniversary of Pedosphere (founded in 1991), contains 15 articles reviewing recent advances and the state of the art of important topics, including some emerging topics, in soil science. The first article in this issue, by Du et al., is a meta-analysis systematically reviewing the effect of elevated carbon dioxide (CO₂) on nitrogen oxide (N₂O) emission from global agricultural soils, depending on other soil and experimental conditions.     

Genetic variability in apple fruit polyphenol composition in Malus × domestica and Malus sieversii germplasm grown in New Zealand

Variations in the concentrations of flavan-3-ol, oligomeric procyanidin, chlorogenic acid, dihydrochalcone, flavonol, and anthocyanin polyphenol groups and total polyphenols were examined in the fruit peel and cortical flesh of 93 (80 Malus × domestica and 13 Malus sieversii) apple genotypes in at least 1 year between 2003 and 2005 grown at one site in New Zealand (NZ). Differences among genotypes accounted for 46-97% of the total variation in the concentrations of total polyphenols and each of the individual phenol groups in the flesh and peel in both species, whereas effects of year and genotype × year were minimal, except for peel flavonols in M. × domestica and flesh flavonols in both species. In these cases, differences among genotypes accounted for less than 30% of the total variation, which was less than the variation found for the interaction between genotype and year. Total polyphenol concentrations among genotypes were spread over a 7- and 9-fold range in the flesh and a 4- and 3-fold range in the peel of M. sieversii and M. × domestica, respectively, with the spread in concentrations of individual polyphenol groups in each tissue and within each species varying from a 2-fold to over a 500-fold range. Higher concentrations were generally found in M. sieversii. In M. × domestica, cultivars and breeding selections originating in NZ had lower average flesh and peel total polyphenols and chlorogenic acid than older cultivars previously imported into NZ from overseas countries.     

Tree species effects on microbial respiration from decomposing leaf and fine root litter

Tree species have an impact on decomposition processes of woody litter, but the effects of different tree species on microbial heterotrophic respiration derived from decomposing litter are still unclear. Here we used leaf and fine root litter of six tree species differing in chemical and morphological traits in a temperate forest and elucidated the effects of tree species on the relationships between litter-derived microbial respiration rates and decomposition rates and morphological traits, including specific leaf area (cm² g⁻¹) and specific root length (m g⁻¹) of litter at the same site. Litterbags set in forest soil were sequentially collected five times over the course of 18 months. During litter decomposition, microbial respiration from leaf and fine root litter differed among the six tree species. Temporal changes in the remaining mass and morphology (specific leaf area and specific root length) were observed, and the magnitude of these changes differed among species. Positive correlations were observed between respiration and mass loss or morphology across species. These results revealed that litter mass loss and morphological dynamics during decomposition jointly enhanced microbial respiration, and these carbon-based litter traits explained species differences in decomposition of leaves and fine roots. In conclusion, tree species influenced the magnitude and direction of microbial respiration during leaf/fine root litter decomposition. Tree species also affected the relationship between microbial respiration and litter decomposition through direct effects of litter traits and indirect effects mediated by regulation of heterotroph requirements.     

Litter-type specific microbial responses to the transformation of leaf litter into millipede feces

The transformation of leaf litter into fecal pellets by saprophagous macroarthropods has long been suggested to play an important role in litter decomposition by altering microbial processes. However, conflicting results are reported in the literature, and it is currently not clear to what extent varying initial litter quality contributes to distinct microbial responses to the transformation of litter into feces. Here we performed a screening test using a wide range of distinct leaf litter from 26 tree species. We fed these litters to the macroarthropod species Glomeris marginata during one week under controlled conditions, and compared microbial responses in uningested leaf litter with that of feces produced from the 26 different leaf litter types. We assessed substrate induced respiration (SIR) as an integrative measure of microbial responses. We found that litter SIR was highly variable across species and well related to initial litter quality. However, variability in feces SIR was strongly reduced and only weakly related to initial litter quality. Moreover, the difference between feces and litter SIR decreased with increasing litter SIR as a result of higher microbial stimulation in litter with low associated litter SIR. Our data clearly showed that the direction and magnitude of microbial stimulation in feces depend strongly on the litter type. Therefore, the consequence of litter transformation into macroarthropod fecal pellets for microbial decomposers and possibly for subsequent decomposition of feces is specific to litter species.     

Spatial variability of the functional stability of microbial respiration process: a microcosm study using tropical forest soil

Purpose

Understanding the ability of ecosystem processes to resist to and to recover from disturbances is critical to sustainable land use. However, the spatial variability of the stability has rarely been addressed. Here, we investigated the functional stability of a soil microbial process for 24 soils collected from adjacent locations from a 0.3?ha tropical rainforest plot in Paracou, French Guiana.

Materials and methods

The 24 locations were characterized regarding soil chemical and biological (microbial diversity) parameters and forest structure. The corresponding soils were submitted to an experimental transient heat disturbance during a microcosm experiment. The response of the respiration process was followed using substrate-induced respiration (SIR).

Results and discussion

The response of soil SIR to heat disturbance varied widely between samples. The variability of the SIR response increased just after the disturbance, and a global rather homogeneous decrease in SIR rates was observed 15 and 30?days after. The stability of SIR in response to heat disturbance could not be related to either the genetic or the metabolic diversity of the microbial community. The initial level of SIR before the disturbance was the soil variable that best correlated with the impact of the disturbance: the soil locations with the highest initial SIR rates were the most affected 15 and 30?days after the heat disturbance.

Conclusions

Such a heterogeneous response suggests that the response of soil processes to a disturbance will be difficult to assess from only local-scale analyses and highlights the need for spatial explicitness in understanding biogeochemical processes.     

Land use change and soil nutrient transformations in the Los Haitises region of the Dominican Republic

We characterized soil cation, carbon (C) and nitrogen (N) transformations within a variety of land use types in the karst region of the northeastern Dominican Republic. We examined a range of soil pools and fluxes during the wet and dry seasons in undisturbed forest, regenerating forest and active agricultural sites within and directly adjacent to Los Haitises National Park. Soil moisture, soil organic matter (SOM), soil cations, leaf litter C and pH were significantly greater in regenerating forest sites than agricultural sites, while bulk density was greater in active agricultural sites. Potential denitrification, microbial biomass C and N, and microbial respiration g⁻¹ dry soil were significantly greater in the regenerating forest sites than in the active agricultural sites. However, net mineralization, net nitrification, microbial biomass C, and microbial respiration were all significantly greater in the agricultural sites on g⁻¹ SOM basis. These results suggest that land use is indirectly affecting microbial activity and C storage through its effect on SOM quality and quantity. While agriculture can significantly decrease soil fertility, it appears that the trend can begin to rapidly reverse with the abandonment of agriculture and the subsequent regeneration of forest. The regenerating forest soils were taken out of agricultural use only 5-7 years before our study and already have soil properties and processes similar to an undisturbed old forest site. Compared to undisturbed mogote forest sites, regenerating sites had smaller amounts of SOM and microbial biomass N, as well as lower rates of microbial respiration, mineralization and nitrification g⁻¹ SOM. Initial recovery of soil pools and processes appeared to be rapid, but additional research must be done to address the long-term rate of recovery in these forest stands.     

Inter- and intra-specific interactions of Lumbricus rubellus (Hoffmeister, 1843) and Octolasion lacteum (Örley, 1881) (Lumbricidae) and the implication for C cycling

Earthworms are important engineering species of many terrestrial ecosystems as they play a significant role in regulating C turnover. The effects of earthworms on moderating C decomposition processes differ across species and with interactions between species, which is not fully understood. We carried out an experiment to study the interactions of Lumbricus rubellus and Octolasion lacteum, and their effects on soil respiration. Laboratory mesocosms were set up using tulip poplar (Liriodendron tulipifera) leaf litter and varying densities of earthworms in single and combined species treatments. CO₂ efflux rate was used as an indicator of C decomposition rates, and measured with CO₂ sensors every five days over one month. L. rubellus induced higher leaf consumption rate and higher CO₂ efflux than O. lacteum; meanwhile O. lacteum grew more than L. rubellus. Both litter consumption rate and growth rate of earthworms decreased with increasing earthworm density. Soil CO₂ efflux increased with increasing earthworm density (from ∼1-2 μg CO₂ g⁻¹ hr⁻¹ with no earthworms to ∼ 4 μg CO₂ g⁻¹ hr⁻¹ with 8 earthworms). Combining the two species had a synergistic effect on leaf litter consumption, and neutralizing effects on soil respiration. The data suggest that the strength of intra- and inter-specific interactions among earthworm ecological groups varies at different absolute and relative densities, leading to altered leaf litter decomposition and C cycling.     

Variation of soil respiration at three spatial scales: Components within measurements, intra-site variation and patterns on the landscape

Soil respiration is an important component of terrestrial carbon cycling and can be influenced by many factors that vary spatially. This research aims to determine the extent and causes of spatial variation of soil respiration, and to quantify the importance of scale on measuring and modeling soil respiration within and among common forests of Northern Wisconsin. The potential sources of variation were examined at three scales: [1] variation among the litter, root, and bulk soil respiration components within individual 0.1 m measurement collars, [2] variation between individual soil respiration measurements within a site (<1 m to 10 m), and [3] variation on the landscape caused by topographic influence (100 m to 1000 m). Soil respiration was measured over a two-year period at 12 plots that included four forest types. Root exclusion collars were installed at a subset of the sites, and periodic removal of the litter layer allowed litter and bulk soil contributions to be estimated by subtraction. Soil respiration was also measured at fixed locations in six northern hardwood sites and two aspen sites to examine the stability of variation between individual measurements. These study sites were added to an existing data set where soil respiration was measured in a random, rotating, systematic clustering which allowed the examination of spatial variability from scales of <1 m to 100+ m. The combined data set for this area was also used to examine the influence of topography on soil respiration at scales of over 1000 m by using a temperature and moisture driven soil respiration model and a 4 km² digital elevation model (DEM) to model soil moisture. Results indicate that, although variation of soil respiration and soil moisture is greatest at scales of 100 m or more, variation from locations 1 m or less can be large (standard deviation during summer period of 1.58 and 1.28 μmol CO₂ m⁻² s⁻¹, respectively). At the smallest of scales, the individual contributions of the bulk soil, the roots, and the litter mat changed greatly throughout the season and between forest types, although the data were highly variable within any given site. For scales of 1-10 m, variation between individual measurements could be explained by positive relationships between forest floor mass, root mass, carbon and nitrogen pools, or root nitrogen concentration. Lastly, topography strongly influenced soil moisture and soil properties, and created spatial patterns of soil respiration which changed greatly during a drought event. Integrating soil fluxes over a 4 km² region using an elevation dependent soil respiration model resulted in a drought induced reduction of peak summer flux rates by 37.5%, versus a 31.3% when only plot level data was used. The trends at these important scales may help explain some inter-annual and spatial variability of the net ecosystem exchange of carbon.     

Impact of land degradation on soil respiration in a steppe (Stipa tenacissima L.) semi-arid ecosystem in the SE of Spain

Climate change scenarios predict increases in temperature, changes in precipitation patterns, and longer drought periods in most semi-arid regions of the world. Ecosystems in these regions are prone to land degradation, which may be aggravated by climate change. Soil respiration is one of the main processes responsible for organic carbon losses from arid and semi-arid ecosystems. We measured soil respiration over one year in two steppe ecosystems having different degrees of land degradation under three ground-covers: with vegetation, bare soil, and an intermediate situation between plants and bare soil.The largest differences in soil respiration rates between the sites were observed in spring, coinciding with the highest level of plant activity. The degraded site had drier and hotter soils with less soil water availability and a longer drought period. As a result, vegetation on the degraded site did not respond to spring rainfall events. Soil respiration showed a strong seasonal variability, with average annual rates of 1.1 and 0.8 μmol CO₂ m⁻² s⁻¹ in the natural and degraded sites, respectively. We did not observe significant differences in soil respiration rates associated with ground-cover i.e., the temporal variation was much larger than the spatial variation. At both sites, soil moisture was the controlling driver of soil respiration for most of the year, when temperatures were above 20 °C and constrained the response to temperature for the few months when the temperature was below 20 °C. An empirical model based on soil temperature and soil moisture explained 90% and 72% of the seasonal variability of soil respiration on the natural and degraded sites, respectively. For the first time, this study suggests that land degradation may alter the carbon balance of these ecosystems through changes in the temporal dynamics of soil respiration and plant productivity, which have important negative consequences for ecosystem functioning and sustainability.     

Changes in soil chemical and microbial properties after a wildfire in a tropical rainforest in Sabah, Malaysia

Changes in soil caused by drought and wildfire in a Dipterocarp rainforest in Sabah, Malaysia were assessed by phosphorus fractionation, extractable nitrogen and nutrient limited respiration kinetics (after addition of glucose+N or P). Fire increased the concentration of total phosphorus (P) in the litter layer (per ha and per dry soil) by raising the 0.2 M NaOH extractable-P. In the soil organic layer, membrane exchangeable P was reduced by fire while 1.0 M HCl extractable-P, and 0.5 M NaHCO₃ extractable-P increased. Microbially available P increased after the fire and was most closely related to NaOH extractable-P that has been considered available to plants only over long time-scales. Total nitrogen (N) increased in the litter layer (per ha and per dry soil) due to post-fire litter fall, while the NO₃⁻ increased up to 10-fold down to the 10 cm mineral soil. In contrast, the microbially available N decreased by 50%. Basal respiration and substrate-induced respiration increased in the litter layer and decreased in the organic horizon (per dry soil and per organic matter). P limited microbial growth resulted in a slow and non-exponential increase in respiration, presumably reflecting the P-fixing nature of the soils, while N limitation resulted in a fast exponential increase. However, higher respiration rates were eventually achieved under P limitation than under N limitation.     

Evaluation of soil quality parameters in a tropical paddy soil amended with rice residues and tree litters

Laboratory and greenhouse experiments were conducted to study the effects of applications of rice residue and Pongamia pinnata and Azadirachta indica leaf litters on biochemical properties (extraction yield of humus, composition of humus, microbial biomass carbon, activities of urease and acid phosphatase) of a lowland rice soil under flooded conditions. Bulk soil sample collected from the Mandya paddy fields was used for the green house trials and the laboratory incubation studies. The organic materials were added at three rates – zero, 25.0 g carbon kg⁻¹ (2.5% C) and 50.0 g carbon kg⁻¹ dry soil (5.0% C). Results showed that tree leaf litter and rice residue at 5.0% C rate decreased instantaneous decay constant (k), there by retarded the rate of C mineralization. Carbon contents of HA increased with the rate of C added. Study of delta–log K values and C contents of humic acids revealed that greatest molecular weight of HA was in the pongamia litter treatment, followed by neem litter and rice residue. Grain and straw yields of rice crop in the pot culture study were statistically correlated to the soil quality parameters. Neem and pongamia tree litter incorporation at 2.5% C could be considered for improving soil health and crop yields of rice under flooded conditions; however, application at higher rates significantly (P ≤ 0.05) lowered total dry matter production in rice, despite favorable soil health parameters such as humic yields, microbial biomass – C content and acid phosphatase and urease activity. Among different soil health parameters, microbial quotient was found to be more sensitive indicator of decline in soil quality.     

Temporal change in spatial variability of soil respiration on a slope of Japanese cedar (Cryptomeria japonica D. Don) forest

Although information regarding the spatial variability of soil respiration is important for understanding carbon cycling and developing a suitable sampling design for estimating average soil respiration, it remains relatively understudied compared to temporal changes. In this study, soil respiration was measured at 35 locations by season on a slope of Japanese cedar forest in order to examine temporal changes in the spatial distribution of soil respiration. Spatial variability of soil respiration varied between seasons, with the highest coefficient variation in winter (42%) and lowest in summer (26%). Semivariogram analysis and kriged maps revealed different patterns of spatial distribution in each season. Factors affecting the spatial variability were relief index (autumn), soil hardness of the A layer (winter), soil hardness at 50 cm depth (spring) and the altitude and relief index (summer). Annual soil respiration (average: 39 mol m⁻² y⁻¹) varied from 26 mol m⁻² y⁻¹ to 55 mol m⁻² y⁻¹ between the 35 locations and was higher in the upper part of the slope and lower in the lower part. The average Q₁₀ value was 2.3, varying from 1.3 to 3.0 among the locations. These findings suggest that insufficient information on the spatial variability of soil respiration and imbalanced sampling could bias estimates of current and future carbon budgets.