Espacialização dos sedimentos produzidos na parte inicial da bacia hidrográfica do Rio Pardo, Botucatu/SP

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise multivariada na diferenciaçãoentre manejos de solos cultivados com café

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Classificação de imagens em áreas cultivadas com citrus por técnicas de sensoriamento remoto e geoestatística

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estimating the extent of spatial association of Mycobacterium bovis infection in badgers in Ireland

Mycobacterium bovis infects the wildlife species badgers Meles meles who are linked with the spread of the associated disease tuberculosis (TB) in cattle. Control of livestock infections depends in part on the spatial and social structure of the wildlife host. Here we describe spatial association of M. bovis infection in a badger population using data from the first year of the Four Area Project in Ireland. Using second-order intensity functions, we show there is strong evidence of clustering of TB cases in each the four areas, i.e. a global tendency for infected cases to occur near other infected cases. Using estimated intensity functions, we identify locations where particular strains of TB cluster. Generalized linear geostatistical models are used to assess the practical range at which spatial correlation occurs and is found to exceed 6 in all areas. The study is of relevance concerning the scale of localized badger culling in the control of the disease in cattle.

Spatial modelling and prediction on river networks: up model, down model or hybrid?

Preservation of rivers and water resources is crucial in most environmental policies and many efforts are made to assess water quality. Environmental monitoring of large river networks are based on measurement stations. Compared to the total length of river networks, their number is often limited and there is a need to extend environmental variables that are measured locally to the whole river network. The objective of this paper is to propose several relevant geostatistical models for river modeling. These models use river distance and are based on two contrasting assumptions about dependency along a river network. Inference using maximum likelihood, model selection criterion and prediction by kriging are then developed. We illustrate our approach on two variables that differ by their distributional and spatial characteristics: summer water temperature and nitrate concentration. The data come from 141 to 187 monitoring stations in a network on a large river located in the Northeast of France that is more than 5000 km long and includes Meuse and Moselle basins. We first evaluated different spatial models and then gave prediction maps and error variance maps for the whole stream network.

Geoelectrical mapping of contamination in the cemeteries: the case study in Piracicaba, Sao Paulo/Brazil

This paper presents the results of electrical resistivity methods in the area delineation that was potentially contaminated by liquefaction products, which are also called putrefactive liquids in Vila Rezende municipal cemetery, Piracicaba, So Paulo, Brazil. The results indicate a depth of water table between 3.1 and 5.1 m, with two groundwater direction flows, one to the SW and another to the SE. Due to the contamination plumes, which have the same groundwater direction flow, as well the conductive anomalies observed in the geoelectric sections, the contamination suspicions in the area were confirmed. The probable plume to the SE extends beyond the limits of the cemetery. The location of the conductive anomalies and the probable contamination plumes showed that the contamination is linked with the depth of the water table and the burial time. Mapping using the geostatistical method of ordinary kriging applied to the work drew structural characteristics of the regional phenomenon and spatial behavior of the electrical resistivity data, resulting in continued surfaces. Thus, this method has proved to be an important tool for mapping contamination plumes in cemeteries.

Integration of spatial and temporal models to predict extreme water table depths

The objective of this work was to evaluate extreme water table depths in a watershed, using methods for geographical spatial data analysis. Groundwater spatio-temporal dynamics was evaluated in an outcrop of the Guarani Aquifer System. Water table depths were estimated from monitoring of water levels in 23 piezometers and time series modeling available from April 2004 to April 2011. For generation of spatial scenarios, geostatistical techniques were used, which incorporated into the prediction ancillary information related to the geomorphological patterns of the watershed, using a digital elevation model. This procedure improved estimates, due to the high correlation between water levels and elevation, and aggregated physical sense to predictions. The scenarios showed differences regarding the extreme levels - too deep or too shallow ones - and can subsidize water planning, efficient water use, and sustainable water management in the watershed.

SPATIAL VARIABILITY OF THE RAINFALL EROSIVE POTENTIAL IN THE STATE OF MATO GROSSO DO SUL, BRAZIL

Information about rainfall erosivity is important during soil and water conservation planning. Thus, the spatial variability of rainfall erosivity of the state Mato Grosso do Sul was analyzed using ordinary kriging interpolation. For this, three pluviograph stations were used to obtain the regression equations between the erosivity index and the rainfall coefficient EI30. The equations obtained were applied to 109 pluviometric stations, resulting in EI30 values. These values were analyzed from geostatistical technique, which can be divided into: descriptive statistics, adjust to semivariogram, cross-validation process and implementation of ordinary kriging to generate the erosivity map. Highest erosivity values were found in central and northeast regions of the State, while the lowest values were observed in the southern region. In addition, high annual precipitation values not necessarily produce higher erosivity values.

Influence diagnostics in Gaussian spatial linear models

Spatial linear models have been applied in numerous fields such as agriculture, geoscience and environmental sciences, among many others. Spatial dependence structure modelling, using a geostatistical approach, is an indispensable tool to estimate the parameters that define this structure. However, this estimation may be greatly affected by the presence of atypical observations in the sampled data. The purpose of this paper is to use diagnostic techniques to assess the sensitivity of the maximum-likelihood estimators, covariance functions and linear predictor to small perturbations in the data and/or the spatial linear model assumptions. The methodology is illustrated with two real data sets. The results allowed us to conclude that the presence of atypical values in the sample data have a strong influence on thematic maps, changing the spatial dependence structure.

Yield mapping, soil fertility and tree gaps in an orange orchard

The current high competition on Citrus industry demands from growers new management technologies for superior efficiency and sustainability. In this context, precision agriculture (PA) has developed techniques based on yield mapping and management systems that recognize field spatial variability, which contribute to increase profitability of commercial crops. Because spatial variability is often not perceived the orange orchards are still managed as uniform and adoption of PA technology on citrus farms is low. Thus, the objective of the present study was to characterize the spatial variability of three factors: fruit yield, soil fertility and occurrence of plant gaps caused by either citrus blight or huanglongbing (HLB) in a commercial Valencia orchard in Brotas, São Paulo State, Brazil. Data from volume, geographic coordinates and representative area of the bags used on harvest were recorded to generate yield points that were then interpolated to produce the yield map. Soil chemical characteristics were studied by analyzing samples collected along planting rows and inter-rows in 24 points distributed in the field. A map of density of tree gaps was produced by georeferencing individual gaps and later by counting the number of gaps within 500 m² cells. Data were submitted to statistical and geostatistical analyses. A t test was used to compare means of soil chemical characteristics between sampling regions. High variation on yield and density of tree gaps was observed from the maps. It was also demonstrated overlapping regions of high density of plant absence and low fruit yield. Soil fertility varied depending on the sampling region in the orchard. The spatial variability found on yield, soil fertility and on disease occurrence demonstrated the importance to adopt site specific nutrient management and disease control as tools to guarantee efficiency of fruit production.

Precision horticulture: application on apple orchards

Precision horticulture and spatial analysis applied to orchards are a growing and evolving part of precision agriculture technology. The aim of this discipline is to reduce production costs by monitoring and analysing orchard-derived information to improve crop performance in an environmentally sound manner. Georeferencing and geostatistical analysis coupled to point-specific data mining allow to devise and implement management decisions tailored within the single orchard. Potential applications range from the opportunity to verify in real time along the season the effectiveness of cultural practices to achieve the production targets in terms of fruit size, number, yield and, in a near future, fruit quality traits. These data will impact not only the pre-harvest but their effect will extend to the post-harvest sector of the fruit chain. Chapter 1 provides an updated overview on precision horticulture , while in Chapter 2 a preliminary spatial statistic analysis of the variability in apple orchards is provided before and after manual thinning; an interpretation of this variability and how it can be managed to maximize orchard performance is offered. Then in Chapter 3 a stratification of spatial data into management classes to interpret and manage spatial variation on the orchard is undertaken. An inverse model approach is also applied to verify whether the crop production explains environmental variation. In Chapter 4 an integration of the techniques adopted before is presented. A new key for reading the information gathered within the field is offered. The overall goal of this Dissertation was to probe into the feasibility, the desirability and the effectiveness of a precision approach to fruit growing, following the lines of other areas of agriculture that already adopt this management tool. As existing applications of precision horticulture already had shown, crop specificity is an important factor to be accounted for. This work focused on apple because of its importance in the area where the work was carried out, and worldwide.

Metodi avanzati di indagine di ammassi rocciosi a blocchi in matrice (bimrock).

The research focuses on the implementation and validation of advanced direct and indirect methods of investigation for the structural and mechanical characterisation of bimrocks. In particular, a non conventional in situ shear test has been develop in order to evaluate the strength parameters of bimrocks by properly taking into account the influence of blocks. Also, a geostatistical approach has been introduced for the investigation of block morphological and spatial properties from digital images, by means of a variographic analysis of the block Indicator Variable.

Distribution and state of fishery resources in the Northern and Central Adriatic Sea using trawl surveys data

The Ecosystem Approach to Fisheries represents the most recent research line in the international context, showing interest both towards the whole community and toward the identification and protection of all the “critical habitats” in which marine resources complete their life cycles. Using data coming from trawl surveys performed in the Northern and Central Adriatic from 1996 to 2010, this study provides the first attempt to appraise the status of the whole demersal community. It took into account not only fishery target species but also by-catch and discharge species by the use of a suite of biological indicators both at population and multi-specific level, allowing to have a global picture of the status of the demersal system. This study underlined the decline of extremely important species for the Adriatic fishery in recent years; adverse impact on catches is expected for these species in the coming years, since also minimum values of recruits recently were recorded. Both the excessive exploitation and environmental factors affected availability of resources. Moreover both distribution and nursery areas of the most important resources were pinpointed by means of geostatistical methods. The geospatial analysis also confirmed the presence of relevant recruitment areas in the North and Central Adriatic for several commercial species, as reported in the literature. The morphological and oceanographic features, the relevant rivers inflow together with the mosaic pattern of biocenoses with different food availability affected the location of the observed relevant nursery areas.

Footprint of recycled water subsidies downwind of Lake Michigan

Continental evaporation is a significant and dynamic flux within the atmospheric water budget, but few methods provide robust observational constraints on the large-scale hydroclimatological and hydroecological impacts of this ‘recycled-water' flux. We demonstrate a geospatial analysis that provides such information, using stable isotope data to map the distribution of recycled water in shallow aquifers downwind from Lake Michigan. The δ2H and δ18O values of groundwater in the study region decrease from south to north, as expected based on meridional gradients in climate and precipitation isotope ratios. In contrast, deuterium excess (d = δ2H − 8 × δ18O) values exhibit a significant zonal gradient and finer-scale spatially patterned variation. Local d maxima occur in the northwest and southwest corners of the Lower Peninsula of Michigan, where ‘lake-effect' precipitation events are abundant. We apply a published model that describes the effect of recycling from lakes on atmospheric vapor d values to estimate that up to 32% of recharge into individual aquifers may be derived from recycled Lake Michigan water. Applying the model to geostatistical surfaces representing mean d values, we estimate that between 10% and 18% of the vapor evaporated from Lake Michigan is re-precipitated within downwind areas of the Lake Michigan drainage basin. Our approach provides previously unavailable observational constraints on regional land-atmosphere water fluxes in the Great Lakes Basin and elucidates patterns in recycled-water fluxes that may influence the biogeography of the region. As new instruments and networks facilitate enhanced spatial monitoring of environmental water isotopes, similar analyses can be widely applied to calibrate and validate water cycle models and improve projections of regional hydroecological change involving the coupled lake-atmosphere-land system. Read More: http://www.esajournals.org/doi/abs/10.1890/ES12-00062.1

Functional error modeling for uncertainty quantification in hydrogeology

Approximate models (proxies) can be employed to reduce the computational costs of estimating uncertainty. The price to pay is that the approximations introduced by the proxy model can lead to a biased estimation. To avoid this problem and ensure a reliable uncertainty quantification, we propose to combine functional data analysis and machine learning to build error models that allow us to obtain an accurate prediction of the exact response without solving the exact model for all realizations. We build the relationship between proxy and exact model on a learning set of geostatistical realizations for which both exact and approximate solvers are run. Functional principal components analysis (FPCA) is used to investigate the variability in the two sets of curves and reduce the dimensionality of the problem while maximizing the retained information. Once obtained, the error model can be used to predict the exact response of any realization on the basis of the sole proxy response. This methodology is purpose-oriented as the error model is constructed directly for the quantity of interest, rather than for the state of the system. Also, the dimensionality reduction performed by FPCA allows a diagnostic of the quality of the error model to assess the informativeness of the learning set and the fidelity of the proxy to the exact model. The possibility of obtaining a prediction of the exact response for any newly generated realization suggests that the methodology can be effectively used beyond the context of uncertainty quantification, in particular for Bayesian inference and optimization.