Soil water content (theta) and bulk density (rhos) greatly influence important soil and plant processes, such as water movement, soil compaction, soil aeration, and plant root system development. Spatial and temporal variability of theta and rhos during different periods of the year and different phases of crops are of fundamental interest. This work involves the characterization of spatial and temporal patterns of theta and rhos during different climatic periods of year, aiming to verify whether there are significant temporal changes in rhos at the soil surface layer when submitted to wetting and drying cycles. The field experiment was carried out in a coffee plantation, Rhodic Kandiudalf soil, clayey texture. Using a neutron/gamma surface probe, theta and rhos were measured meter by meter along a 200 m spatial transect, along an interrow contour line. During the wet period there was no difference of spatial patterns of theta while during the dry period differences were observed, and can be associated to precipitation events. It was also observed that there are rhos temporal changes at the soil surface along the studied period as a consequence of the in situ wetting and drying cycles.
neutron/gamma surface probe; wetting/drying cycles; spatial and temporal variability
Umidade (teta) e densidade do solo (ros) influenciam importantes processos no solo e planta tais como: movimento de água, compactação do solo, aeração do solo e desenvolvimento radicular. Baseado neste fato, questões referentes à variabilidade espacial e temporal de teta e ros para diferentes períodos do ano e diferentes fases de desenvolvimento da cultura do café tornam-se de extremo interesse. Este trabalho teve como objetivo caracterizar padrões espaciais e temporais de teta e ros durante diferentes períodos do ano e verificar se existem mudanças temporais significativas de ros na superfície do solo quando submetida a ciclos de umedecimento/secagem. O experimento foi conduzido em campo cultivado com café em um solo classificado como Nitossolo Vermelho Eutrófico. O experimento iniciou-se em maio de 2001 com espaçamento de 1,75 m entre linhas e 0,75 m entre plantas. Medidas de teta e ros foram feitas metro a metro ao longo de 200 m usando uma sonda de superfície nêutron/gama. Durante o período úmido não houve padrões de diferença espacial para teta, enquanto que para o período seco o contrário foi observado. Estes padrões podem ser associados a eventos de precipitação alterando a estrutura de correlação espacial para observações adjacentes de teta. Para ros existem mudanças temporais na superfície do solo ao longo do período estudado como uma conseqüência dos ciclos de umedecimento/secagem no campo.
sonda de superfície nêutron/gama; ciclos de umedecimento/secagem; variabilidade espacial e temporal
SOILS AND PLANT NUTRITION
Field spatial and temporal patterns of soil water content and bulk density changes
Padrões espacial e temporal de mudanças da umidade e densidade do solo no campo
Luís Carlos TimmI; Luiz Fernando PiresII; Renato RoverattiII; Robson Clayton Jacques ArthurII; Klaus ReichardtII, * * Corresponding author < email@example.com> ; Júlio Cesar Martins de OliveiraIII; Osny Oliveira Santos BacchiII
IUFPel/FAEM - Depto de Engenharia Rural, C.P. 354 - 96001-970 - Pelotas, RS - Brasil
IIUSP/CENA - Lab. de Física de Solo, C.P. 96 - 13416-000 - Piracicaba, SP - Brasil
IIIEEP - Depto. de Física, Av. Monsenhor Martinho Salgot, 560 - 13414-040 - Piracicaba, SP - Brasil
Soil water content (q) and bulk density (rs) greatly influence important soil and plant processes, such as water movement, soil compaction, soil aeration, and plant root system development. Spatial and temporal variability of q and rs during different periods of the year and different phases of crops are of fundamental interest. This work involves the characterization of spatial and temporal patterns of q and rs during different climatic periods of year, aiming to verify whether there are significant temporal changes in rs at the soil surface layer when submitted to wetting and drying cycles. The field experiment was carried out in a coffee plantation, Rhodic Kandiudalf soil, clayey texture. Using a neutron/gamma surface probe, q and rs were measured meter by meter along a 200 m spatial transect, along an interrow contour line. During the wet period there was no difference of spatial patterns of q while during the dry period differences were observed, and can be associated to precipitation events. It was also observed that there are rs temporal changes at the soil surface along the studied period as a consequence of the in situ wetting and drying cycles.
Key words: neutron/gamma surface probe, wetting/drying cycles, spatial and temporal variability
Umidade (q) e densidade do solo (rs) influenciam importantes processos no solo e planta tais como: movimento de água, compactação do solo, aeração do solo e desenvolvimento radicular. Baseado neste fato, questões referentes à variabilidade espacial e temporal de q e rs para diferentes períodos do ano e diferentes fases de desenvolvimento da cultura do café tornam-se de extremo interesse. Este trabalho teve como objetivo caracterizar padrões espaciais e temporais de q e rs durante diferentes períodos do ano e verificar se existem mudanças temporais significativas de rs na superfície do solo quando submetida a ciclos de umedecimento/secagem. O experimento foi conduzido em campo cultivado com café em um solo classificado como Nitossolo Vermelho Eutrófico. O experimento iniciou-se em maio de 2001 com espaçamento de 1,75 m entre linhas e 0,75 m entre plantas. Medidas de q e rs foram feitas metro a metro ao longo de 200 m usando uma sonda de superfície nêutron/gama. Durante o período úmido não houve padrões de diferença espacial para q, enquanto que para o período seco o contrário foi observado. Estes padrões podem ser associados a eventos de precipitação alterando a estrutura de correlação espacial para observações adjacentes de q. Para rs existem mudanças temporais na superfície do solo ao longo do período estudado como uma conseqüência dos ciclos de umedecimento/secagem no campo.
Palavras-chave: sonda de superfície nêutron/gama, ciclos de umedecimento/secagem, variabilidade espacial e temporal
Several factors affect crop productivity, especially the edaphic, climatic and management factors. Among edaphic factors, physico-chemical attributes of farmed soils stand out because they affect directly plants dry matter accumulation in different periods of the year and in all development phases of the crop. Tominaga et al. (2002) showed that on tropical soils, physical attributes such as soil water content and bulk density, greatly influence important soil and plant processes like water movement (Reichardt & Timm, 2004), soil compaction (Logsdon & Karlen, 2004), soil aeration (Stepniewski et al., 1994), and plant root system development (Boone & Veen, 1994).
An important question is whether soil bulk density can be considered a static soil property or if small temporal changes are important (Logsdon & Karlen, 2004). In this sense, Logsdon & Cambardella (2000) detected significant temporal changes in near-surface soil bulk densities for tillage systems in a sub-humid climate.
At soil surface, soil water content influences the partitioning of precipitation into infiltration and runoff, and controls evapotranspiration by controlling evaporation and water availability to plants, therefore affecting also the partitioning of latent and sensible heat (Grayson et al., 1997). An improved understanding of the spatial distribution of soil water content is thus useful for a range of hydrologic factors affecting crop production. However, an important question for hydrology is "how does the temporal variation in soil water content affect its spatial patterns evaluated through geostatistical tools". Ahuja et al. (1993) reported that temporal changes of soil surface moisture two days after wetting revealed a lot about the spatial distribution of drainage and hydraulic properties within the soil profile.
Wendroth et al. (1999) monitored the soil water pressure heads in surface horizons at two field sites, a sandy loam and a heavy clay soil, between April and November 1995 in north-east Germany, and concluded that with decreasing soil water pressure head, the variance of log10(-h) decreased to a critical value, for which a spatial correlation structure disappeared. On the other hand, with further drying, its variance increased again, and a spatial range of correlation appeared.
Soil spatial variability occurs at different scales and is related to variations of the parent material, climate, relief, organisms and time, i.e., related to the processes of soil formation and/or effects of adopted management practices for each agricultural use (McGraw, 1994). Large areas have been managed as homogeneous, although presenting a considerable spatial variability inherent to the soil, causing the appearance of zones of distinct soil physical properties and distinct soil fertility (Wendroth et al., 2003).
Several statistical tools have been used to evaluate the variabilities of soil attributes, which can contribute to improve management and the understanding of atmosphere-plant-soil interactions (Vieira et al., 1981; Cassel, 1983; Morkoc et al., 1985; Wendroth et al., 2001; 2003; Dourado-Neto et al., 1999; Timm et al., 2003; 2004; Tominaga et al., 2002; Martínez & Zinck, 2004). This work studies the importance of soil water content and bulk density changes, using some of these tools to characterize spatial and temporal patterns of changes during different climatic periods of the year, with the objective of examining possible differences that might occur when a soil is submitted to wetting or drying processes.
MATERIAL AND METHODS
The field experiment was carried out along interrows of a coffee plantation in Piracicaba, SP, Brazil (22º42' S; 47º38' W; altitude 580 m), on a Rhodic Kandiudalf, clayey texture. According to Köppen's classification, the region's climate is the Cwa type, tropical highland, mesothermal with dry winter. Average values for air temperature, rainfall, and relative humidity are 21.2ºC; 1,253 mm/year; and 74%, respectively. The dry season covers June-August, July being the driest month. During spring-summer, October to March, very high intensity rainfall events are common, several of them reaching 50 mm h-1 or more. Soil water content q (m3 m-3) and soil bulk density rs (Mg m-3) were evalueted in a coffee crop field (variety Catuai Red IAC44) planted in May 2001, along rows 1.75 m apart, with 0.75 m between plants. This narrow spacing between coffee plants (0.75 m) resulted in high density plant ranks. During the development of the crop, according to Rena & Maestri (1987), the ranked plants undergo less water stress due to their adaptation to a higher intra-specific competition, i.e., plant roots grow deeper into the soil seeking water and nutrients. The ranked plant system also reduces leaf and soil temperatures. Mendes (2001) noted that although the progenies of the Catuai variety are vigorous, they should not be submitted to stress conditions such as high periods of hydric deficiency, among others.
Monitoring of q and rs at the soil surface layer (0 0.15 m) was made along the bare interrow, when the coffee crop was two years old (April to August 2003), during the end of its vegetative cycle (Marin, 2003). At the first q and rs evaluation (April 2003) coffee plants were all green along the spatial transect, presenting large quantity of leaves, being about 1m tall, shading a 1 m strip of interrow, leaving a bare central strip of about 0.75 m on which measurements were taken, in a period without hydric stress problems. During the July to August 2003 period, hydric stress problems could visually be observed on the plants, causing leaf loss and a consequent higher exposition of the soil surface to solar radiation. This period corresponds to the winter season in the Southern Hemisphere, in which events of precipitation are rare and of low intensity in Piracicaba.
This study focuses q and rs measurements carried out every meter along a 200 meter spatial transect. Measurements were made using a neutron-gamma surface gauge (SG), which yields average values of the 0-0.15 m soil layer, calibrated for this particular soil according to Cássaro et al. (2000). The advantage of this methodology is that the "sampled" layer remains undisturbed for several measurements made at the very same point. Two distinct climatic periods were chosen: i. April to May 2003, a wet period - 1st measurement MW1: 04/24/2003; 2nd measurement MW2: 05/07/2003; 3rd measurement MW3: 05/23/2003; and ii. July to August 2003, a dry period - 1st measurement MD1: 07/04/2003; 2nd measurement MD2: 07/18/2003; 3rd measurement MD3: 08/02/2003. These 6 dates (measurements) were chosen according to the possibility of obtaining a good contrast between wetting and drying processes. In any other season of the year this would not be possible.
Daily variation of precipitation and reference evapotranspiration ETo of the experimental field are presented for the period under study in Figure 1. MW1 was measured five days after a precipitation event of 12.2 mm (04/19/2003); MW2 after two precipitation events of 26.1 mm and 9.4 mm on days 05/02 and 05/05/2003, shown on Figure 1. The third spatial series MW3 was measured after a period of sixteen days without precipitation events, the crop being submitted to a long period of evapotranspiration only. Changes in q between consecutive measurements (Dq) were analysed to examine wetting and drying cycles during both seasons.
Soil texture (clay, sand and silt contents) was evaluated using disturbed soil samples, collected at the end of the experiment, also every meter along the 200 meter spatial transect, in the 0-0.15 m soil layer at the same points where q and rs were taken. Details about the soil mechanical analysis can be found in EMBRAPA (1997). Relations of soil texture variability to water content and bulk density were also investigated.
Statistical analysis was performed using the classical concepts of mean, standard deviation, coefficients of variation and the t test (a = 0.05) through the GLM procedure of SAS® (SAS Institute, 1996). The significance level used for the variance analysis was of 5% through the Duncan test, using the LSMEANS procedure. To analyse and characterize spatial data distributions, geostatistical tools such as autho- and cross-correlations (Nielsen & Wendroth, 2003) were used:
where n is the number of observations; j is the number of lags h between observations used for correlation; xi the variable at position i and xi+j at position i+j; r(j) are correlation coefficients for j lags; s2 the variance; and the cross-correlations rc (j): xi and yi+j are the variables under study, with variances and .
To determine whether precipitation events change the spatial correlation structure among the adjacent q observations, autocorrelograms of each spatial series were calculated, since according to Nielsen & Wendroth (2003) the autocorrelation function ACF is a primary diagnostic tool which indicates if there is a spatial or temporal pattern of on-site sampled data. The number of points considered as replicates for calculating the rs mean values for each data set was also calculated. After this, the Duncan statistical test (a = 0.05) was used to compare the mean values of each data set using SAS software (SAS Institute, 1996).
RESULTS AND DISCUSSION
Wet period MW1, MW2, MW3 (April 24 to May 23, 2003): precipitation 56.4 mm
To visualize the spatial distribution of the soil water content data sets, Figure 2 is presented as a sample, corresponding to: MW1 04/24/2003. The statistical parameters for all six dates are shown in Table 1.
For the MW series, qave is inversely related to the CV showing a clear tendency of increased variability with soil drying, i.e., the dispersion of soil water contents around their mean is lower when the soil is wet. The MW2 measurement had the highest qave because of recent precipitation, whereas the MW3 had the lowest qave because of the previous dry period. The wetting spatial patterns of Dq presented higher dispersion of data in relation to its mean (CV = 32.3%) as compared to the drying spatial pattern (CV = 22.1%). This fact may be explained by the hysteresis effect, which can strongly affect the water redistribution in the soil profile (Kutilek & Nielsen, 1994). Figure 3 presents the spatial dependence of q data sets in relation to their neighborhood for the three different dates.
Adjacent q observations for MW1 are spatially dependent up to 2 m (2 lags), i.e., based on the t test (a = 0.05), observations separated by a distance of less than 2 m would not have different soil water content. It is important to note that 1 lag = 1 m so that the number of lags becomes equal to distances. For MW2 in which the highest q values were measured at any location is representative for a nearly 10 m domain, i.e., adjacent q observations are spatially dependent up to 5 m using the t test (a = 0.05). Adjacent q observations on MW3 are spatially dependent up to 4 m using the same statistical test. It was not possible to identify for the wet period, a relationship between the precipitation event and spatial pattern of q along the transect, since the spatial correlation structure presented in Figures 3B (MW2) and 3C (MW3) were similar. Grayson et al. (1997) proposed two hypotheses in respect to surface q spatial patterns in temperate regions of Australia (and similar climates elsewhere): i. in periods when precipitation continually exceeds evapotranspiration, it would be expected that the soil water content data set would present a spatial correlation structure; and ii. in periods when evapotranspiration continually exceeds precipitation, it would be expected that the soil water content data set would present a more random appearance, i.e., a data set without a spatial correlation structure among observations. These hypotheses were formulated for temperate regions, different of the tropical highland climate of the Piracicaba region, which has higher rainfall and a clayey loam soil. From all soil bulk density (rs) data sets (MW1, MW2 and MW3) a sample (MW1) is presented in Figure 4. Table 2 shows all classical statistics parameters for these data sets.
The CV ranged from 3.8 (MW2) to 4.5% (MW3), lower than in the case of q, an indication that observations of rs present a certain spatial homogeneity along the spatial transect for the three different dates. For the date at which the highest q values were observed (MW2) it was registered that rs presented the lowest dispersion in relation to their mean (CV = 3.8%), i.e., CVs for q and rs data sets are lower when the soil is wet, indicating higher spatial homogeneity of data in this condition. The autocorrelation function ACF for the rs data sets are shown in Figure 5.
Adjacent rs observations are spatially correlated only up to 1 m for MW1 (Figure 5A) and MW2 (Figure 5B), with no spatial correlation for MW3 (Figure 5C), presenting a random behavior at this date. Crosscorrelation functions CCF for the differences Dq and Drs for both wet periods, are shown in Figure 6.
For the wet soil (MW1 to MW2, Figure 6A) as well as for the dry soil (MW2 to MW3, Figure 6B) there is only a significant spatial correlation between q and rs variations at the same spatial point of transect (lag = 0), showing that they are not spatially correlated, following a random distribution. The spatial distribution of clay, sand and silt contents of the 0-0.15 m layer along the transect can be observed in Figure 7. The soil, as mentioned before, is clayey with predominance of 1:1 clay minerals (mainly kaolinite) with visible swelling and shrinking movements causing cracks down to the depth of 0.20 m during dry periods, as observed in situ.
The spatial degree of the linkage between Dq and soil texture for both wet periods was also made through crosscorrelograms (Figure 8). A significant correlation was found at the same spatial point between Dq in the MW1 to MW2 period and silt content (Figure 8A), and no significant correlation with clay or sand. From Figures 8 B and C , a significant correlation was found at the same spatial point between Dq and clay content as well as silt content in the next MW2 to MW3 period, indicating that during the drying of the soil, Dq at the same point is correlated with the Drs (soil structure) and with clay and silt contents which are the soil particles of lowest diameter and more important for soil water retention when the soil is drying. Nielsen & Wendroth (2003) comment that near soil water saturation, retention of water is not dominated by the total soil particle surface area, but is influenced more by soil aggregation and the internal architecture of soil pores (soil structure).
Dry period MD1, MD2, MD3 (July 04 to August 02, 2003): 18.6 mm
The classical statistics parameters for these data sets are also presented in Table 1. MD3 presented mean of 0.147 m3 m-3 which corresponds to the lowest qave measured along the spatial transect in all periods. MD2 period was taken six days after two precipitation events and this is the reason why we observed the highest qave of the dry period at this date. After 16 days without precipitation, when the highest evapotranspiration was observed, measurements of MD3 were taken with the lowest q observations along the transect.
The wetting spatial pattern of the differences Dq presented a higher dispersion in relation to its mean (CV of 91.4%) as compared to the drying spatial pattern (CV of 60.3%), the same behavior as in the wet period. Figure 9 presents the autocorrelograms of these data sets.
The spatial dependency of q up to one lag for MD1 (Figure 9A) means that q measured at any location represents the soil water content of a 2 m domain, observations being spatially independent among them beyond this. For MD2 (Figure 9B), which includes the highest q values in this period, adjacent q observations are spatially dependent up to eight lags (9 m), and for MD3 (Figure 9C), up to 2 m. A joint analysis indicates that the precipitation event (07/10/2003) modified the spatial correlation structure of soil water content data sets for the dry period, i.e., the spatial correlation range for the soil water content data set measured at date of MD2 (Figure 9B) was higher in comparison to the other two dates. This fact can be associated to the presence of cracks along the spatial transect due to wetting and drying cycles, causing a soil mass contraction and increasing rs of the soil matrix between cracks. This increase in rs reduces vertical water infiltration and facilitates lateral water movement (cracks form preferential channels), increasing the similarity between adjacent q observations.
Simultaneously measured rs data sets (Table 2), show again that the observations present a certain spatial homogeneity along the transect for the three different dates. According to the CV classification proposed by Wilding & Drees (1983), these rs measurements can be classified as low variability (CV < 15%). As an example, Grossman & Reinsch (2002) reported rs spatial variability of 10% and Warrick & Nielsen (1980) found a rs of 7% for rs in three different experimental areas. Autocorrelograms for our data sets are shown in Figure 10.
Only for MD3 (Figure 10C) there was a spatial correlation for 2 m distance. Figures 4 and 10 reveal that adjacent rs observations present a higher spatial dependence range in the dry period as compared to the wet period. This fact was also observed for the q measurements. Crosscorrelation functions CCF between the differences Dq and Drs and soil texture, are shown in Figure 11.
For the MD1 to MD2 period there is no significant spatial correlation between Dq and: (i) Drs (Figure 11A); (ii) clay content (Figure 11B); (iii) sand content (Figure 11C); and (iv) silt content (Figure 11D); i.e., in this period, for lag = 0, there is no correlation. Same results (not presented) were found for the MD2 to MD3 period.
From the soil water content autocorrelograms shown for wet (Figures 3A - C) and dry (Figures 9A - C) periods it is verified that the spatial correlation structure of q measurements is not unique along the transect for the different dates, i.e., it is time-dependent, not being unique within the field. The same behavior was observed for the crosscorrelograms shown in Figures 6A-B, 8A-C and 11A-D with a different spatial correlation structure among q and rs variations, and q variations and soil texture. This can be attributed to the fact that q is rarely in hydraulic equilibrium in the soil profile, constantly being under redistribution within the profile as a result of plant root extraction, evaporation and infiltration at the soil surface. Therefore, autocorrelograms and crosscorrelograms of soil water properties and attributes are not expected to be unique within a field (Nielsen & Wendroth, 2003).
Comparisons made among the rs mean values for each date, to verify if there were significant temporal changes of the rs mean value along the studied period (April to August 2003) due to the natural mechanisms of consolidation of soil parent material in the 0-0.15 m soil layer of the experimental field, presented interesting results (Table 3). As mentioned before, the soil presents swelling and shrinking movements during the wetting and drying cycles along the different periods of the year. Based on the rs autocorrelograms in the wet (Figures 5A-C) and dry (Figures 10A-C) periods, and on the low CV values found for each soil bulk density data set (CVs ranged from 3.8% to 4.9%), mean values of the rs data sets measured at MW1 and MW2 dates can be obtained considering one hundred points as replicates; 200 points for MW3, MD1 and MD2; and sixty six points for MD3.
Soil bulk density temporal changes here presented (Table 3) indicate the rs cannot be considered as a static soil property. It can also be observed that there is an increase of its mean value during the dry period as a consequence of shrinking soil movements. There are no significant differences among mean values of rs for MW3, MD1 and MD2 and between rs mean values for MW1 and MW2. Although MW3 was considered to belong to the wet period, its mean value along the spatial transect was not statistically different from the mean values for MD1 and MD2 which belong to the dry period.
To FAPESP, Capes and CNPq for financial support.
Received February 03, 2005
Accepted December 19, 2005
- AHUJA, L.R.; WENDROTH, O.; NIELSEN, D.R. Relationship between initial drainage of surface soil and average profile saturated conductivity. Soil Science Society of America Journal, v.57, p.19-25, 1993.
- BOONE, F.R.; VEEN, D.E. Mechanisms of crop responses to soil compaction. In: SOANE, B.D.; VAN OUWERKERK, C. (Ed.) Soil compaction in crop production New York: Elsevier, 1994. p.237-264.
- CÁSSARO, F.A.M.; TOMINAGA, T.T.; BACCHI, O.O.S.; REICHARDT, K.; OLIVEIRA, J.C.M.; TIMM, L.C. Improved laboratory calibration of a single-probe surface gamma-neutron gauge. Australian Journal Soil Research, v.38, p.937-946, 2000.
- CASSEL, D.K. Spatial and temporal variability of soil physical properties following tillage of Norfolk loamy sand. Soil Science Society of America Journal, v.47, p.196-201, 1983.
- DOURADO-NETO, D.; TIMM, L.C.; OLIVEIRA, J.C.M.; REICHARDT, K.; BACCHI, O.O.S.; TOMINAGA, T.T.; CASSARO, F.A.M. State-space approach for the analysis of soil water content and temperature in a sugarcane crop. Scientia Agricola, v.56, p.1215-1221, 1999.
- EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA-EMBRAPA. Centro Nacional de pesquisa de Solo. Manual de métodos de análise de solos 2.ed. Rio de Janeiro: EMBRAPA, CNPS, 1997. 212p.
- GRAYSON, R.B.; WESTERN, A.W.; CHIEW, F.H.S. Preferred states in spatial soil moisture patterns: local and nonlocal controls. Water Resources Research, v.33, p.2897-2908, 1997.
- GROSSMAN, R.B.; REINSCH, T.G. The solid phase: 2.1. Bulk density and linear extensibility. In: DANE, J.H.; TOPP, G.C. (Ed.) Methods of soil analysis Part 4. Physical methods. Madison: SSSA, 2002. ch.2, p.201-228. (Book Series, 5).
- KUTILEK, M.; NIELSEN, D.R. Soil hydrology Cremlingen-Destedt: Catena Verlag, 1994. 370p.
- LOGSDON, S.D.; CAMBARDELLA, C.A. Temporal changes in small depth incremental soil bulk density. Soil Science Society of America Journal, v.64, p.710-714, 2000.
- LOGSDON, S.D.; KARLEN, D.L. Bulk density as a soil quality indicator during conversion to no-tillage. Soil & Tillage Research, v.78, p.143-149, 2004.
- MARIN, F.R. Evapotranspiration and maximum transpiration in a high density coffee plantation. Piracicaba: USP/ESALQ, 2003. 118p. (Dissertação - Mestrado).
- MARTÍNEZ, L.J.; ZINCK, J.A. Temporal variation of soil compaction and deterioration of soil quality in pasture areas of Colombian Amazonia. Soil & Tillage Research, v.75, p.3-17, 2004.
- MCGRAW, T. Soil test level variability in Southern Minnesota. Better Crops with Plant Foods, v.78, p.24-25, 1994.
- MENDES, A.N.G. Cultivares com potencialidade para lavouras irrigadas. In: SANTOS, C.M.; MENDONÇA, F.C.; MELO, B.; TEODORO, R.E.F.; SANTOS, V.L.M. (Ed.) Irrigação da cafeicultura no cerrado Uberlândia: Universidade Federal de Uberlândia, 2001. p.125-135.
- MORKOC, F.; BIGGAR, J.W.; NIELSEN, D.R.; ROLSTON, D.E. Analysis of soil water content and temperature using State-space approach. Soil Science Society of America Journal, v.49, p.798-803, 1985.
- NIELSEN, D.R.; WENDROTH, O. Spatial and temporal statistics: sampling field soils and their vegetation Reiskirchen: Catena Verlag, 2003. 398p.
- REICHARDT, K.; TIMM, L.C. Solo, planta e atmosfera: conceitos, processos e aplicações. Barueri: Manole, 2004. 478p.
- RENA, A.B.; MAESTRI, M. Ecofisiologia do cafeeiro. In: CASTRO, P.R.C.; FERREIRA, S.O.; YAMADA, T. (Ed.) Ecofisiologia da produção agrícola. Piracicaba: Potafos, 1987. p.119-145.
- SAS INSTITUTE. System for information Version 6.11. Cary: SAS Institute Inc., 1996.
- STEPNIEWSKI, W.; GLINSKI, J.; BALL, B.C. Effects of compaction on soil aeration properties. In: SOANE, B.D.; VAN OUWERKERK, C. (Ed.) Soil compaction in crop production New York: Elsevier, 1994. p.167-169.
- TIMM, L.C.; REICHARDT, K.; OLIVEIRA, J.C.M.; CASSARO, F.A.M.; TOMINAGA, T.T.; BACCHI, O.O.S.; DOURADO-NETO, D. Sugarcane production evaluated by the state-space approach. Journal of Hydrology, v.272, p.226-237, 2003.
- TIMM, L.C.; REICHARDT, K.; OLIVEIRA, J.C.M.; CASSARO, F.A.M.; TOMINAGA, T.T.; BACCHI, O.O.S.; DOURADO-NETO, D.; NIELSEN, D.R. State-space approach to evaluate the relation between soil physical and chemical properties. Brazilian Journal of Soil Science, v.28, p.49-58, 2004.
- TOMINAGA, T.T.; CÁSSARO, F.A.M.; BACCHI, O.O.S.; REICHARDT, K.; OLIVEIRA, J.C.M.; Timm, L.C. Variability of soil water content and bulk density in a sugarcane field. Australian Journal of Soil Research, v.40, p.605-614, 2002.
- VIEIRA, S.R.; NIELSEN, D.R.; BIGGAR, J.W. Spatial variability of field-measured infiltration rate. Soil Science Society of America Journal, v.45, p.1040-1048, 1981.
- WARRICK, A.W.; NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D. (Ed.) Applications of soil physics New York: Academic Press, 1980. p.319-344.
- WENDROTH, O.; REUTER, H.I.; KERSEBAUM, K.C. Predicting yield of barley across a landscape: a state-space modeling approach. Journal of Hydrology, v.272, p.250-263, 2003.
- WENDROTH, O.; JÜRSCHIK, P.; KERSEBAUM, K.C.; REUTER, H.; VAN KESSEL, C.; Nielsen, D.R. Identifying, understanding, and describing spatial processes in agricultural landscapes four case studies. Soil & Tillage Research, v.58, p.113-127, 2001.
- WENDROTH, O.; POHL, W.; KOSZINSKI, S.; ROGASIK, H.; RITSEMA, C.J.; NIELSEN, D.R. Spatio-temporal patterns and covariance structures of soil water status in two Northeast-German field sites. Journal of Hydrology, v.215, p.38-58, 1999.
- WILDING, L.P.; DREES, L.R. Spatial variability and pedology. In: WILDING, L.P.; SMECK, N.E.; HALL, G.F. (Ed.) Pedogenesis and soil taxonomy: concepts and interactions. New York: Elsevier, 1983. p.83-116.
Publication in this collection
22 Feb 2006
Date of issue
19 Dec 2005
03 Feb 2005