Estimativa de densidade amostral para elaboração de mapas de produtividade

Yield mapping represents the spatial variability concerning the features of a productive area and allows intervening on the next year production, for example, on a site-specific input application. The trial aimed at verifying the influence of a sampling density and the type of interpolator on yield mapping precision to be produced by a manual sampling of grains. This solution is usually adopted when a combine with yield monitor can not be used. An yield map was developed using data obtained from a combine equipped with yield monitor during corn harvesting. From this map, 84 sample grids were established and through three interpolators: inverse of square distance, inverse of distance and ordinary kriging, 252 yield maps were created. Then they were compared with the original one using the coefficient of relative deviation (CRD) and the kappa index. The loss regarding yield mapping information increased as the sampling density decreased. Besides, it was also dependent on the interpolation method used. A multiple regression model was adjusted to the variable CRD, according to the following variables: spatial variability index and sampling density. This model aimed at aiding the farmer to define the sampling density, thus, allowing to obtain the manual yield mapping, during eventual problems in the yield monitor.

Disproportionate single-species contribution to canopy-soil nutrient flux in an Amazonian rainforest

Rainfall, throughfall and stemflow were monitored on an event basis in an undisturbed open tropical rainforest with a large number of palm trees located in the southwestern Amazon basin of Brazil. Stemflow samples were collected from 24 trees with a diameter at breast height (DBH) > 5 cm, as well as eight young and four full-grown babassu palms (Attalea speciosa Mart.) for 5 weeks during the peak of the wet season. We calculated rainfall, throughfall and stemflow concentrations and fluxes of Na+, K+, Ca2+, Mg2+,, Cl-, SO42-, NO3- and H+ and stemflow volume-weighted mean concentrations and fluxes for three size classes of broadleaf trees and three size classes of palms. The concentrations of most solutes were higher in stemflow than in rainfall and increased with increasing tree and palm size. Concentration enrichments from rainfall to stemflow and throughfall were particularly high (81-fold) for NO3-. Stemflow fluxes of NO3- and H+ exceeded throughfall fluxes but stemflow fluxes of other solutes were less than throughfall fluxes. Stemflow solute fluxes to the forest soil were dominated by fluxes on babassu palms, which represented only 4% of total stem number and 10% of total basal area. For NO3-, stemflow contributed 51% of the total mass of nitrogen delivered to the forest floor (stemflow + throughfall) and represented more than a 2000-fold increase in NO3- flux compared what would have been delivered by rainfall alone on the equivalent area. Because these highly localized fluxes of both water and NO3- persist in time and space, they have the potential to affect patterns of soil moisture, microbial populations and other features of soil biogeochemistry conducive to the creation of hotspots for nitrogen leaching and denitrification, which could amount to an important fraction of total ecosystem fluxes. Because these hotspots occur over very small areas, they have likely gone undetected in previous studies and need to be considered as an important feature of the biogeochemistry of palm-rich tropical forest. (C) 2011 Elsevier B.V. All rights reserved.

Estabilidade temporal e variabilidade espacial da distribuição da armazenagem de água no solo numa sucessão feijão/aveia-preta

A armazenagem de água no solo é muito variável no tempo e no espaço devido à influência de vários fatores ambientais e topográficos. Objetivou-se neste estudo: quantificar a armazenagem de água no solo; avaliar sua estabilidade temporal e sua variabilidade espacial em um local próximo e outro distante do sistema radicular numa sucessão feijão/aveia-preta; e constatar sua variabilidade espacial em função do relevo. Em área experimental de 1.500 m², situada em Piracicaba-SP (latitude de 22° 42' 30" S, longitude de 47° 38' 00" W e 546 m de altitude), estabeleceram-se 60 pontos de amostragem, distanciados entre si de 5 m, numa grade de 10 por 6 pontos (50 x 30 m). Os valores da armazenagem de água no solo apresentaram comportamento-padrão para o solo estudado, com valores maiores em profundidade do que na camada mais superficial. Houve maior estabilidade temporal da armazenagem de água no solo durante o período chuvoso, porém na fase de secagem do solo a estabilidade temporal também foi constatada, mas com valores de coeficiente de correlação mais elevados na camada de 0,0-0,80 m, mostrando que essa estabilidade foi claramente devida à posição topográfica desses pontos, os quais estão localizados na porção mais baixa do relevo. Assim, o ponto 52 foi escolhido como representativo da média na fase de recarga de água no solo em ambas as camadas estudadas, e o ponto 46, na camada superior durante a fase de secagem do solo. A dependência espacial da variabilidade da armazenagem de água no solo foi detectada em ambas as camadas de solo, porém com maior alcance na camada de 0,0-0,40 m, quando houve ocorrência de precipitação. Quando as chuvas cessaram, o alcance foi maior na camada de 0,0-0,80 m. Durante o período chuvoso, o padrão de variabilidade espacial foi muito semelhante em ambas as camadas de solo, com armazenagens maiores na faixa de maior inclinação e depressões do terreno.

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.

Long-term effects of oil pollution in mangrove forests (Baixada Santista, Southeast Brazil) detected using a GIS-based multitemporal analysis of aerial photographs

Oil spills are potential threats to the integrity of highly productive coastal wetlands, such as mangrove forests. In October 1983, a mangrove area of nearly 300 ha located on the southeastern coast of Brazil was impacted by a 3.5 million liter crude oil spill released by a broken pipeline. In order to assess the long-term effects of oil pollution on mangrove vegetation, we carried out a GIS-based multitemporal analysis of aerial photographs of the years 1962, 1994, 2000 and 2003. Photointerpretation, visual classification, class quantification, ground-truth and vegetation structure data were combined to evaluate the oil impact. Before the spill, the mangroves exhibited a homogeneous canopy and well-developed stands. More than ten years after the spill, the mangrove vegetation exhibited three distinct zones reflecting the long-term effects of the oil pollution. The most impacted zone (10.5 ha) presented dead trees, exposed substrate and recovering stands with reduced structural development. We suggest that the distinct impact and recovery zones reflect the spatial variability of oil removal rates in the mangrove forest. This study identifies the multitemporal analysis of aerial photographs as a useful tool for assessing a system's capacity for recovery and monitoring the long-term residual effects of pollutants on vegetation dynamics, thus giving support to mangrove forest management and conservation.

Mapeamento do índice de vegetação da diferença normalizada em lavoura de algodão

Os dados de sensoriamento remoto em campo podem fornecer informações detalhadas sobre a variabilidade de parâmetros biofísicos ligados à produtividade em grandes áreas e apresentam potencial para o monitoramento destes parâmetros, ao longo de todo o ciclo de desenvolvimento da cultura. Este trabalho objetivou mapear a variabilidade espacial do índice de vegetação da diferença normalizada (NDVI) e seus componentes, em duas lavouras comerciais de algodão (Gossipium hirsutum L.), utilizando sensor óptico ativo, em nível terrestre. Os dados foram coletados utilizando-se sensor instalado em um pulverizador autopropelido agrícola. Um receptor GPS foi acoplado ao sensor, para a obtenção das coordenadas dos pontos de amostragem. As leituras foram realizadas em faixas espaçadas em 21,0 m, aproveitando-se as passadas do veículo no momento da pulverização de agroquímicos, e os dados submetidos à análise estatística clássica e geoestatística. Mapas de distribuição espacial das variáveis foram elaborados pela interpolação por krigagem. Observou-se maior variabilidade espacial do NDVI e da reflectância espectral da vegetação na região do infravermelho próximo (IVP) (880 nm) e do visível (590 nm) na lavoura com maior estresse fisiológico, devido ao ataque do percevejo castanho [Scaptocoris castanea (Hem.: Cydnidae)], em relação à lavoura sadia.

The aplication of land use and nighttime lights on the development of high-resolution emission inventories for South America

In this paper, nighttime light data are suggested as a proxy for spatial distribution of vehicles running in urban and nearby areas. Nighttime lights focus on human activities, in contrast to traditional Earth observing systems that focus on natural systems. It is the human activity being visible in the form of brightness of nocturnal lights. Two available nighttime lights dataset were used in this work. The first one was provided by the U.S. Air Force Defense Meteorological Satellite Program (DMSP) Operational Linescan System (OLS), henceforth, DMSO-OLS. The second one is the NASA-NOAA Suomi National Polar-orbiting Polar-orbiting Partnership (NPP) satellite, henceforth, Suomi-NPP. To validate the new proposed methodology, hundreds of urban areas of South America were analyzed in a high degree of resolution. The results of this study showed that night-time lights are very well correlated with vehicle fleet, population, and impervious surfaces but with strong spatial variability. The results of this study suggest a better understanding of the human activities in the context of a vehicular-based city conception.

Dynamics of the cold surface layer of polythermal Storglaciären, Sweden

Polythermal glaciers, i.e. glaciers with a combination of ice at and below the freezing point, are widespread in arctic and subarctic environments. The polythermal structure has major implications for glacier hydrology, ice flow and glacial erosion. However, the interplay of factors governing its spatial and temporal variations such as net mass balance, ice advection and water content in the ice is poorly investigated and as yet not fully understood. This study deals with a thorough investigation of the polythermal regime on Storglaciären, northern Sweden, a small valley glacier with a cold surface layer in the ablation area. Extensive field work was performed including mapping of the cold surface layer using ground-penetrating radar, ice temperature measurements, mass balance and ice velocity measurements. Analyses of these data combined with numerical modelling were used specifically to investigate the spatial and temporal variability of the cold surface layer, the spatial distribution of the water content just below the cold surface layer transition, the effect of radar frequency on the detection of the surface layer, and the sensitivity of the cold surface layer to changes in forcing. A comparison between direct temperature measurements in boreholes and ground-penetrating surveys shows that the radar-inferred cold-temperate transition depth is within ±1 m from the melting point of ice at frequencies above ~300 MHz. At frequencies below ~155 MHz, the accuracy degrades because of reduced scattering efficiency that occurs when the scatterers become much smaller compared to the wavelength. The mapped spatial pattern of the englacial cold-temperate transition boundary is complex. This pattern reflects the observed spatial variation in net loss of ice at the surface by ablation and vertical advection of ice, which is suggested to provide the predominant forcing of the cold surface layer thickness pattern. This is further supported by thermomechanical modeling of the cold surface layer, which indicates high sensitivity of the cold surface layer thickness to changes in vertical advection rates. The water content is the least investigated quantity that is relevant for the thermal regime of glaciers, but also the most difficult to assess. Spatial variability of absolute water content in the temperate ice immediately below the cold surface layer on Storglaciären was determined by combining relative estimates of water content from ground-penetrating radar data with absolute determination from temperature measurements and the thermal boundary condition at the freezing front. These measurements indicate large-scale spatial variability in the water content, which seems to arise from variations in entrapment of water at the firn-ice transition. However, this variability cannot alone explain the spatial pattern in the thermal regime on Storglaciären. Repeated surveys of the cold surface layer show a 22% average thinning of the cold surface layer on Storglaciären between 1989 and 2001. Transient thermomechanical modeling results suggest that the cold surface layer adapts to new equilibrium conditions in only a few decades after a perturbation in the forcing is introduced. An increased winter air temperature since mid-1980s seems to be the cause of the observed thinning of the cold surface layer. Over the last decades, mass balance measurements indicate that the glacier has been close to a steady state. The quasi-steady state situation is also reflected in the vertical advection, which shows no significant changes during the last decades. Increased winter temperatures at the ice surface would result in a slow-down of the formation of cold ice at the base of the cold surface layer and lead to a larger imbalance between net loss of ice at the surface and freezing of temperate ice at the cold-temperate transition.

Variabilidad del nivel del mar a partir de datos del altímetro del satélite ERS-2 al sur del Archipiélago Canario en 1998

En este trabajo se presentan algunos resultados obtenidos del análisis de la variabilidad de la altura de la superficie del mar a partir de las anomalías del nivel del mar proporcionadas por los datos del altímetro a bordo del satélite ERS-2. La finalidad del estudio has sido la determinación de la variación estacional que las estructuras oceanográficas mesoescalares presentan en las proximidades del archipiélago canario durante el año 1998. En esta zona, caracterizada por la generación de remolinos ciclónicos y anticiclónicos al sur de las islas debida a la perturbación que experimenta la corriente de Canarias a su paso por los canales entre las islas, y por los filamentos de agua fría procedente del afloramiento, el altímetro se muestra como una herramienta importante en la detección y posterior análisis de estas estructuras oceanográficas. Los resultados muestran que la variabilidad espacial y temporal del nivel del mar es máxima en el segundo semestre del año, y ésta se centra, fundamentalmente, en una estrecha banda situada al sudoeste del archipiélago. ABSTRACT: Some results obtained from the analysis of the sea surface height variability using sea level anomalies given by ERS-2 altimeter data are shown in this work. The aim of the study is to work out the seasonal variations of the mesoscale oceanographic features that appear in the vicinity of the Canary Archipelago during 1998 year. This area is characterized by cyclonic and anticyclonic eddies southward of the islands, which are generated by the interference suffered by the Canary Current through the canals between the islands, and also owing to cold water filaments coming from the Upwelling. The altimeter demonstrates to be an important tool in the detection and posterior analysis of these features. The results show that the temporal and spatial variability of the sea level is associated, fundamentally, to a narrow band located to the southwest of the archipelago, and which has been clearly seen with greater intensity during the periods of summer and autumn of 1998.

Scaling properties of rainfall and desert dust in the Canary Islands

[EN] Precipitation and desert dust event occurrence time series measured in the Canary Islands region are examined with the primary intention of exploring their scaling characteristics as well as their spatial variability in terms of the islands topography and geographical orientation. In particular, the desert dust intrusion regime in the islands is studied in terms of its relationship with visibility. Analysis of dust and rainfall events over the archipelago exhibits distributions in time that obey power laws. Results show that the rain process presents a high clustering and irregular pattern on short timescales and a more scattered structure for long ones. In contrast, dustiness presents a more uniform and dense structure and, consequently, a more persistent behaviour on short timescales. It was observed that the fractal dimension of rainfall events shows an important spatial variability, which increases with altitude, as well as towards northern latitudes and western longitudes.

Applicabilità degli idrogrammi sintetici di progetto nello studio del comportamento idraulico dei corsi d'acqua naturali

Nella protezione idraulica del territorio la previsione e il controllo delle piene sono di fondamentale importanza. I territori sono sempre più antropizzati, pertanto la riduzione dei rischi connessi a eventi idrometeorologici estremi è di notevole interesse. La previsione delle piene è resa difficile dall’innumerevole quantità di variabili che intervengono nel processo della loro formazione. Nelle attività di progettazione e nella verifica di opere idrauliche la identificazione dell’idrogramma di progetto spesso riveste un’importanza fondamentale. Un idrogramma di progetto è definito come un’onda di piena, realmente osservata o sintetica, associata ad un determinato livello di rischio, quantificato usualmente in termini di tempo di ritorno. Con il presente lavoro si cerca di verificare la possibilità di applicazione una metodologia per la stima degli idrogrammi di progetto associati ad un determinato tempo di ritorno, recentemente proposta dalla letteratura scientifica (Maione et al., 2001, Una metodologia per la stima indiretta degli idrogrammi sintetici per il progetto di opere di difesa idraulica del territorio). Il lavoro è riferito al Fiume Secchia, un affluente importante del Po che scorre tra le provincie di Modena e Reggio Emilia.

Effects of Paramuricea clavata (Risso, 1826) (Anthozoa: Plexauridae) forests on settlement of benthic species: an experimental approach

Rationale: Coralligenous habitat is considered the second most important subtidal “hot spot” of species diversity in the Mediterranean Sea after the Posidonia oceanica meadows. It can be defined as a typical Mediterranean biogenic hard bottom, mainly produced by the accumulation of calcareous encrusting algae that, together with other builder organisms, form a multidimensional framework with a high micro-spatial variability. The development of this habitat depends on physical factors (i.e. light, hydrodynamism, nutrients, etc.), but also biologic interactions can play a relevant role in structuring the benthic assemblages. This great environmental heterogeneity allows several different assemblages to coexist in a reduced space. One of the most beautiful is that characterised by the Mediterranean gorgonian Paramuricea clavata (Risso, 1826) that can contribute to above 40% of total biomass of the community and brings significant structural complexity into the coralligenous habitat. In sites moderately exposed to waves and currents, P. clavata can form high-density populations (up to 60 colonies m-2) between 20 – 70 m in depth. Being a suspension feeder, where it forms dense populations, P. clavata plays a significant role in transferring energy from planktonic to benthic system. The effects of the branched colonies of P. clavata could be comparable to those of the forests on land. They can affect the micro scale hydrodynamism and light, promoting or inhibiting the growth of other species. Unfortunately, gorgonians are threatened by several anthropogenic disturbance factors (i.e. fishing, pollution, tourism) and by climatic anomalies, linked to the global changes, that are responsible of thermal stress, development of mucilage and enhanced pathogens activity, leading to mass mortality events in last decades. Till now, the possible effects of gorgonian forest loss are largely unknown. Our goal was to analyse the ecological role of these sea fan forests on the coralligenous benthic assemblages. Experimental setup and main results: The influence of P. clavata in the settlement and recruitment of epibenthic organisms was analysed by a field experiment carried out in two randomly selected places: Tavolara island and Portofino promontory. The experiment consisted in recreate the presence and absence of the gorgonian forest on recruitment panels, arranged in four plots per type (forested and non-forested), interspersed each other, and deployed at the same depth. On every forested panel 3 gorgonian colonies about 20 cm height were grafted with the use of Eppendorf tubes and epoxy resin bicomponent simulating a density of 190 sea fans per m-2. This density corresponds to a mean biomass of 825 g DW m-2,3 which is of the same order of magnitude of the natural high-density populations. After about 4 months, the panels were collected and analysed in laboratory in order to estimate the percent cover of all the species that have colonized the substrata. The gorgonian forest effects were tested by multivariate and univariate permutational analyses of the variance (PERMANOVA). Recruited assemblages largely differed between the two study sites, probably due to different environmental conditions including water quality and turbidity. On overall, the presence of P. clavata reduced the settlement and recruitment of several algae: the shadow caused by the gorgonian might reduce light availability and therefore their growth. This effect might be greater in places where the waters are on average more clear, since at Portofino it is less visible and could be masked by the high turbidity of the water. The same pattern was registered for forams, more abundant outside gorgonian forest, probably linked with algal distribution, shadowing effect or alimentary competition. The last one hypothesis could be valid also for serpulids polychaetes that growth mainly on non-forested panels. An opposite trend, was showed by a species of bryozoan and by an hydroid that is facilitated by the presence of P. clavata, probably because it attenuates irradiance level and hydrodynamism. Species diversity was significantly reduced by the presence of P. clavata forests at both sites. This seems in contrast with what we expected, but the result may be influenced by the large algal component on non-forested panels. The analysis confirmed the presence of differences in the species diversity among plots and between sites respectively due to natural high variability of the coralligenous system and to different local environment conditions. The reduction of species diversity due to the presence of gorgonians appeared related to a worst evenness rather than to less species richness. With our experiment it is demonstrated that the presence of P. clavata forests can significantly alter local coralligenous assemblages patterns, promoting or inhibiting the recruitment of some species, modifying trophic relationships and adding heterogeneity and complexity to the habitat. Moreover, P. clavata could have a stabilising effect on the coralligenous assemblages.

Natural patterns and anthropogenic disturbance in north Adriatic marine benthic assemblages: descriptive and methodological studies

Marine soft bottom systems show a high variability across multiple spatial and temporal scales. Both natural and anthropogenic sources of disturbance act together in affecting benthic sedimentary characteristics and species distribution. The description of such spatial variability is required to understand the ecological processes behind them. However, in order to have a better estimate of spatial patterns, methods that take into account the complexity of the sedimentary system are required. This PhD thesis aims to give a significant contribution both in improving the methodological approaches to the study of biological variability in soft bottom habitats and in increasing the knowledge of the effect that different process (both natural and anthropogenic) could have on the benthic communities of a large area in the North Adriatic Sea. Beta diversity is a measure of the variability in species composition, and Whittaker’s index has become the most widely used measure of beta-diversity. However, application of the Whittaker index to soft bottom assemblages of the Adriatic Sea highlighted its sensitivity to rare species (species recorded in a single sample). This over-weighting of rare species induces biased estimates of the heterogeneity, thus it becomes difficult to compare assemblages containing a high proportion of rare species. In benthic communities, the unusual large number of rare species is frequently attributed to a combination of sampling errors and insufficient sampling effort. In order to reduce the influence of rare species on the measure of beta diversity, I have developed an alternative index based on simple probabilistic considerations. It turns out that this probability index is an ordinary Michaelis-Menten transformation of Whittaker's index but behaves more favourably when species heterogeneity increases. The suggested index therefore seems appropriate when comparing patterns of complexity in marine benthic assemblages. Although the new index makes an important contribution to the study of biodiversity in sedimentary environment, it remains to be seen which processes, and at what scales, influence benthic patterns. The ability to predict the effects of ecological phenomena on benthic fauna highly depends on both spatial and temporal scales of variation. Once defined, implicitly or explicitly, these scales influence the questions asked, the methodological approaches and the interpretation of results. Problem often arise when representative samples are not taken and results are over-generalized, as can happen when results from small-scale experiments are used for resource planning and management. Such issues, although globally recognized, are far from been resolved in the North Adriatic Sea. This area is potentially affected by both natural (e.g. river inflow, eutrophication) and anthropogenic (e.g. gas extraction, fish-trawling) sources of disturbance. Although few studies in this area aimed at understanding which of these processes mainly affect macrobenthos, these have been conducted at a small spatial scale, as they were designated to examine local changes in benthic communities or particular species. However, in order to better describe all the putative processes occurring in the entire area, a high sampling effort performed at a large spatial scale is required. The sedimentary environment of the western part of the Adriatic Sea was extensively studied in this thesis. I have described, in detail, spatial patterns both in terms of sedimentary characteristics and macrobenthic organisms and have suggested putative processes (natural or of human origin) that might affect the benthic environment of the entire area. In particular I have examined the effect of off shore gas platforms on benthic diversity and tested their effect over a background of natural spatial variability. The results obtained suggest that natural processes in the North Adriatic such as river outflow and euthrophication show an inter-annual variability that might have important consequences on benthic assemblages, affecting for example their spatial pattern moving away from the coast and along a North to South gradient. Depth-related factors, such as food supply, light, temperature and salinity play an important role in explaining large scale benthic spatial variability (i.e., affecting both the abundance patterns and beta diversity). Nonetheless, more locally, effects probably related to an organic enrichment or pollution from Po river input has been observed. All these processes, together with few human-induced sources of variability (e.g. fishing disturbance), have a higher effect on macrofauna distribution than any effect related to the presence of gas platforms. The main effect of gas platforms is restricted mainly to small spatial scales and related to a change in habitat complexity due to a natural dislodgement or structure cleaning of mussels that colonize their legs. The accumulation of mussels on the sediment reasonably affects benthic infauna composition. All the components of the study presented in this thesis highlight the need to carefully consider methodological aspects related to the study of sedimentary habitats. With particular regards to the North Adriatic Sea, a multi-scale analysis along natural and anthopogenic gradients was useful for detecting the influence of all the processes affecting the sedimentary environment. In the future, applying a similar approach may lead to an unambiguous assessment of the state of the benthic community in the North Adriatic Sea. Such assessment may be useful in understanding if any anthropogenic source of disturbance has a negative effect on the marine environment, and if so, planning sustainable strategies for a proper management of the affected area.

Modelling coupled physical-biogeochemical processes in ice-covered oceans

The last decades have seen a large effort of the scientific community to study and understand the physics of sea ice. We currently have a wide - even though still not exhaustive - knowledge of the sea ice dynamics and thermodynamics and of their temporal and spatial variability. Sea ice biogeochemistry is instead largely unknown. Sea ice algae production may account for up to 25% of overall primary production in ice-covered waters of the Southern Ocean. However, the influence of physical factors, such as the location of ice formation, the role of snow cover and light availability on sea ice primary production is poorly understood. There are only sparse localized observations and little knowledge of the functioning of sea ice biogeochemistry at larger scales. Modelling becomes then an auxiliary tool to help qualifying and quantifying the role of sea ice biogeochemistry in the ocean dynamics. In this thesis, a novel approach is used for the modelling and coupling of sea ice biogeochemistry - and in particular its primary production - to sea ice physics. Previous attempts were based on the coupling of rather complex sea ice physical models to empirical or relatively simple biological or biogeochemical models. The focus is moved here to a more biologically-oriented point of view. A simple, however comprehensive, physical model of the sea ice thermodynamics (ESIM) was developed and coupled to a novel sea ice implementation (BFM-SI) of the Biogeochemical Flux Model (BFM). The BFM is a comprehensive model, largely used and validated in the open ocean environment and in regional seas. The physical model has been developed having in mind the biogeochemical properties of sea ice and the physical inputs required to model sea ice biogeochemistry. The central concept of the coupling is the modelling of the Biologically-Active-Layer (BAL), which is the time-varying fraction of sea ice that is continuously connected to the ocean via brines pockets and channels and it acts as rich habitat for many microorganisms. The physical model provides the key physical properties of the BAL (e.g., brines volume, temperature and salinity), and the BFM-SI simulates the physiological and ecological response of the biological community to the physical enviroment. The new biogeochemical model is also coupled to the pelagic BFM through the exchange of organic and inorganic matter at the boundaries between the two systems . This is done by computing the entrapment of matter and gases when sea ice grows and release to the ocean when sea ice melts to ensure mass conservation. The model was tested in different ice-covered regions of the world ocean to test the generality of the parameterizations. The focus was particularly on the regions of landfast ice, where primary production is generally large. The implementation of the BFM in sea ice and the coupling structure in General Circulation Models will add a new component to the latters (and in general to Earth System Models), which will be able to provide adequate estimate of the role and importance of sea ice biogeochemistry in the global carbon cycle.

Characterization of geothermal reservoirs’ parameters by inverse problem resolution and geostatistical simulations

BTES (borehole thermal energy storage)systems exchange thermal energy by conduction with the surrounding ground through borehole materials. The spatial variability of the geological properties and the space-time variability of hydrogeological conditions affect the real power rate of heat exchangers and, consequently, the amount of energy extracted from / injected into the ground. For this reason, it is not an easy task to identify the underground thermal properties to use when designing. At the current state of technology, Thermal Response Test (TRT) is the in situ test for the characterization of ground thermal properties with the higher degree of accuracy, but it doesn’t fully solve the problem of characterizing the thermal properties of a shallow geothermal reservoir, simply because it characterizes only the neighborhood of the heat exchanger at hand and only for the test duration. Different analytical and numerical models exist for the characterization of shallow geothermal reservoir, but they are still inadequate and not exhaustive: more sophisticated models must be taken into account and a geostatistical approach is needed to tackle natural variability and estimates uncertainty. The approach adopted for reservoir characterization is the “inverse problem”, typical of oil&gas field analysis. Similarly, we create different realizations of thermal properties by direct sequential simulation and we find the best one fitting real production data (fluid temperature along time). The software used to develop heat production simulation is FEFLOW 5.4 (Finite Element subsurface FLOW system). A geostatistical reservoir model has been set up based on literature thermal properties data and spatial variability hypotheses, and a real TRT has been tested. Then we analyzed and used as well two other codes (SA-Geotherm and FV-Geotherm) which are two implementation of the same numerical model of FEFLOW (Al-Khoury model).