Spatio-temporal modelling of environmental data

Programa Doutoral em Matemática e Aplicações.

Spatio-temporal modelling of the duration of the cotton cycle in the State Of Goiás, Brazil

The objective of this work was to develop and validate a mathematical model to estimate the duration of cotton (Gossypium hirsutum L. r. latifolium hutch) cycle in the State of Goiás, Brazil, by applying the method of growing degree-days (GD), and considering, simultaneously, its time-space variation. The model was developed as a linear combination of elevation, latitude, longitude, and Fourier series of time variation. The model parameters were adjusted by using multiple-linear regression to the observed GD accumulated with air temperature in the range of 15°C to 40°C. The minimum and maximum temperature records used to calculate the GD were obtained from 21 meteorological stations, considering data varying from 8 to 20 years of observation. The coefficient of determination, resulting from the comparison between the estimated and calculated GD along the year was 0.84. Model validation was done by comparing estimated and measured crop cycle in the period from cotton germination to the stage when 90 percent of bolls were opened in commercial crop fields. Comparative results showed that the model performed very well, as indicated by the Pearson correlation coefficient of 0.90 and Willmott agreement index of 0.94, resulting in a performance index of 0.85.

Spatio-temporal analyses of the relationship between armed conflict and climate change in the eastern Africa

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

Spatio-temporal modeling of data imputation for daily rainfall series in homogeneous zones.

Spatio-temporal modelling is an area of increasing importance in which models and methods have often been developed to deal with specific applications. In this study, a spatio-temporal model was used to estimate daily rainfall data. Rainfall records from several weather stations, obtained from the Agritempo system for two climatic homogeneous zones, were used. Rainfall values obtained for two fixed dates (January 1 and May 1, 2012) using the spatio-temporal model were compared with the geostatisticals techniques of ordinary kriging and ordinary cokriging with altitude as auxiliary variable. The spatio-temporal model was more than 17% better at producing estimates of daily precipitation compared to kriging and cokriging in the first zone and more than 18% in the second zone. The spatio-temporal model proved to be a versatile technique, adapting to different seasons and dates.

Modelling the spatio-temporal pattern of primary dispersal in stone pine (Pinus pinea L.) stands in the Northern Plateau (Spain)

Natural regeneration in stone pine (Pinus pinea L.) managed forests in the Spanish Northern Plateau is not achieved successfully under current silviculture practices, constituting a main concern for forest managers. We modelled spatio-temporal features of primary dispersal to test whether (a) present low stand densities constrain natural regeneration success and (b) seed release is a climate-controlled process. The present study is based on data collected from a 6 years seed trap experiment considering different regeneration felling intensities. From a spatial perspective, we attempted alternate established kernels under different data distribution assumptions to fit a spatial model able to predict P. pinea seed rain. Due to P. pinea umbrella-like crown, models were adapted to account for crown effect through correction of distances between potential seed arrival locations and seed sources. In addition, individual tree fecundity was assessed independently from existing models, improving parameter estimation stability. Seed rain simulation enabled to calculate seed dispersal indexes for diverse silvicultural regeneration treatments. The selected spatial model of best fit (Weibull, Poisson assumption) predicted a highly clumped dispersal pattern that resulted in a proportion of gaps where no seed arrival is expected (dispersal limitation) between 0.25 and 0.30 for intermediate intensity regeneration fellings and over 0.50 for intense fellings. To describe the temporal pattern, the proportion of seeds released during monthly intervals was modelled as a function of climate variables – rainfall events – through a linear model that considered temporal autocorrelation, whereas cone opening took place over a temperature threshold. Our findings suggest the application of less intensive regeneration fellings, to be carried out after years of successful seedling establishment and, seasonally, subsequent to the main rainfall period (late fall). This schedule would avoid dispersal limitation and would allow for a complete seed release. These modifications in present silviculture practices would produce a more efficient seed shadow in managed stands.

SESAM - a new framework integrating macroecological and species distribution models for predicting spatio-temporal patterns of species assemblages

Two different approaches currently prevail for predicting spatial patterns of species assemblages. The first approach (macroecological modelling, MEM) focuses directly on realised properties of species assemblages, whereas the second approach (stacked species distribution modelling, S-SDM) starts with constituent species to approximate assemblage properties. Here, we propose to unify the two approaches in a single 'spatially-explicit species assemblage modelling' (SESAM) framework. This framework uses relevant species source pool designations, macroecological factors, and ecological assembly rules to constrain predictions of the richness and composition of species assemblages obtained by stacking predictions of individual species distributions. We believe that such a framework could prove useful in many theoretical and applied disciplines of ecology and evolution, both for improving our basic understanding of species assembly across spatio-temporal scales and for anticipating expected consequences of local, regional or global environmental changes. In this paper, we propose such a framework and call for further developments and testing across a broad range of community types in a variety of environments.

The residual circulation of the Southern Ocean: which spatio-temporal scales are needed?

The Southern Ocean circulation consists of a complicated mixture of processes and phenomena that arise at different time and spatial scales which need to be parametrized in the state-of-the-art climate models. The temporal and spatial scales that give rise to the present-day residual mean circulation are here investigated by calculating the Meridional Overturning Circulation (MOC) in density coordinates from an eddy-permitting global model. The region sensitive to the temporal decomposition is located between 38°S and 63°S, associated with the eddy-induced transport. The ‘‘Bolus’’ component of the residual circulation corresponds to the eddy-induced transport. It is dominated by timescales between 1 month and 1 year. The temporal behavior of the transient eddies is examined in splitting the ‘‘Bolus’’ component into a ‘‘Seasonal’’, an ‘‘Eddy’’ and an ‘‘Inter-monthly’’ component, respectively representing the correlation between density and velocity fluctuations due to the average seasonal cycle, due to mesoscale eddies and due to large-scale motion on timescales longer than one month that is not due to the seasonal cycle. The ‘‘Seasonal’’ bolus cell is important at all latitudes near the surface. The ‘‘Eddy’’ bolus cell is dominant in the thermocline between 50°S and 35°S and over the whole ocean depth at the latitude of the Drake Passage. The ‘‘Inter-monthly’’ bolus cell is important in all density classes and is maximal in the Brazil–Malvinas Confluence and the Agulhas Return Current. The spatial decomposition indicates that a large part of the Eulerian mean circulation is recovered for spatial scales larger than 11.25°, implying that small-scale meanders in the Antarctic Circumpolar Current (ACC), near the Subantarctic and Polar Fronts, and near the Subtropical Front are important in the compensation of the Eulerian mean flow.

In-situ spatio-temporal measurements of the detailed azimuthal substructure of the substorm current wedge

The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100~km at altitudes of 4,000-7,000~km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240~s after Cluster 4 at 1,300-2,000~km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction, in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven ``wedgelets". Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

DYPTOP: a cost-efficient TOPMODEL implementation to simulate sub-grid spatio-temporal dynamics of global wetlands and peatlands

Simulating the spatio-temporal dynamics of inundation is key to understanding the role of wetlands under past and future climate change. Earlier modelling studies have mostly relied on fixed prescribed peatland maps and inundation time series of limited temporal coverage. Here, we describe and assess the the Dynamical Peatland Model Based on TOPMODEL (DYPTOP), which predicts the extent of inundation based on a computationally efficient TOPMODEL implementation. This approach rests on an empirical, grid-cell-specific relationship between the mean soil water balance and the flooded area. DYPTOP combines the simulated inundation extent and its temporal persistency with criteria for the ecosystem water balance and the modelled peatland-specific soil carbon balance to predict the global distribution of peatlands. We apply DYPTOP in combination with the LPX-Bern DGVM and benchmark the global-scale distribution, extent, and seasonality of inundation against satellite data. DYPTOP successfully predicts the spatial distribution and extent of wetlands and major boreal and tropical peatland complexes and reveals the governing limitations to peatland occurrence across the globe. Peatlands covering large boreal lowlands are reproduced only when accounting for a positive feedback induced by the enhanced mean soil water holding capacity in peatland-dominated regions. DYPTOP is designed to minimize input data requirements, optimizes computational efficiency and allows for a modular adoption in Earth system models.

A silviculture-oriented spatio-temporal model for germination in Pinus pinea L. in the Spanish Northern Plateau based on a direct seeding experiment

Natural regeneration in Pinus pinea stands commonly fails throughout the Spanish Northern Plateau under current intensive regeneration treatments. As a result, extensive direct seeding is commonly conducted to guarantee regeneration occurrence. In a period of rationalization of the resources devoted to forest management, this kind of techniques may become unaffordable. Given that the climatic and stand factors driving germination remain unknown, tools are required to understand the process and temper the use of direct seeding. In this study, the spatio-temporal pattern of germination of P. pinea was modelled with those purposes. The resulting findings will allow us to (1) determine the main ecological variables involved in germination in the species and (2) infer adequate silvicultural alternatives. The modelling approach focuses on covariates which are readily available to forest managers. A two-step nonlinear mixed model was fitted to predict germination occurrence and abundance in P. pinea under varying climatic, environmental and stand conditions, based on a germination data set covering a 5-year period. The results obtained reveal that the process is primarily driven by climate variables. Favourable conditions for germination commonly occur in fall although the optimum window is often narrow and may not occur at all in some years. At spatial level, it would appear that germination is facilitated by high stand densities, suggesting that current felling intensity should be reduced. In accordance with other studies on P. pinea dispersal, it seems that denser stands during the regeneration period will reduce the present dependence on direct seeding.

Spatio-Temporal Tracking and Phylodynamics of an Urban Dengue 3 Outbreak in Sao Paulo, Brazil

The dengue virus has a single-stranded positive-sense RNA genome of similar to 10.700 nucleotides with a single open reading frame that encodes three structural (C, prM, and E) and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) proteins. It possesses four antigenically distinct serotypes (DENV 1-4). Many phylogenetic studies address particularities of the different serotypes using convenience samples that are not conducive to a spatio-temporal analysis in a single urban setting. We describe the pattern of spread of distinct lineages of DENV-3 circulating in Sao Jose do Rio Preto, Brazil, during 2006. Blood samples from patients presenting dengue-like symptoms were collected for DENV testing. We performed M-N-PCR using primers based on NS5 for virus detection and identification. The fragments were purified from PCR mixtures and sequenced. The positive dengue cases were geo-coded. To type the sequenced samples, 52 reference sequences were aligned. The dataset generated was used for iterative phylogenetic reconstruction with the maximum likelihood criterion. The best demographic model, the rate of growth, rate of evolutionary change, and Time to Most Recent Common Ancestor (TMRCA) were estimated. The basic reproductive rate during the epidemics was estimated. We obtained sequences from 82 patients among 174 blood samples. We were able to geo-code 46 sequences. The alignment generated a 399-nucleotide-long dataset with 134 taxa. The phylogenetic analysis indicated that all samples were of DENV-3 and related to strains circulating on the isle of Martinique in 2000-2001. Sixty DENV-3 from Sao Jose do Rio Preto formed a monophyletic group (lineage 1), closely related to the remaining 22 isolates (lineage 2). We assumed that these lineages appeared before 2006 in different occasions. By transforming the inferred exponential growth rates into the basic reproductive rate, we obtained values for lineage 1 of R(0) = 1.53 and values for lineage 2 of R(0) = 1.13. Under the exponential model, TMRCA of lineage 1 dated 1 year and lineage 2 dated 3.4 years before the last sampling. The possibility of inferring the spatio-temporal dynamics from genetic data has been generally little explored, and it may shed light on DENV circulation. The use of both geographic and temporally structured phylogenetic data provided a detailed view on the spread of at least two dengue viral strains in a populated urban area.

Spatio-temporal patterns of throughfall and solute deposition in an open tropical rain forest

The brief interaction of precipitation with a forest canopy can create a high spatial variability of both throughfall and solute deposition. We hypothesized that (i) the variability in natural forest systems is high but depends on system-inherent stability, (ii) the spatial variability of solute deposition shows seasonal dynamics depending on the increase in rainfall frequency, and (iii) spatial patterns persist only in the short-term. The study area in the north-western Brazilian state of Rondonia is subject to a climate with a distinct wet and dry season. We collected rain and throughfall on an event basis during the early wet season (n = 14) and peak of the wet season (n = 14) and analyzed the samples for pH and concentrations of NH4+, Na+, K+, Ca2+ Mg2+,, Cl-, NO3-, SO42- and DOC. The coefficient 3 4 cient of variation for throughfall based on both sampling intervals was 29%, which is at the lower end of values reported from other tropical forest sites, but which is higher than in most temperate forests. Coefficients of variation of solute deposition ranged from 29% to 52%. This heterogeneity of solute deposition is neither particularly high nor particularly tow compared with a range of tropical and temperate forest ecosystems. We observed an increase in solute deposition variability with the progressing wet season, which was explained by a negative correlation between heterogeneity of solute deposition and antecedent dry period. The temporal stability of throughfall. patterns was Low during the early wet season, but gained in stability as the wet season progressed. We suggest that rapid plant growth at the beginning of the rainy season is responsible for the lower stability, whereas less vegetative activity during the later rainy season might favor the higher persistence of ""hot"" and ""cold"" spots of throughfall. quantities. The relatively high stability of throughfall patterns during later stages of the wet season may influence processes at the forest floor and in the soil. Solute deposition patterns showed less clear trends but all patterns displayed a short-term stability only. The weak stability of those patterns is apt to impede the formation of solute deposition -induced biochemical microhabitats in the soil. (C) 2008 Elsevier B.V. All rights reserved.

Spatio-temporal variability behavior of land surface soil water content in shrub- and grass-land

The knowledge of the relationship between spatial variability of the surface soil water content (theta) and its mean across a spatial domain (theta(m)) is crucial for hydrological modeling and understanding soil water dynamics at different scales. With the aim to compare the soil moisture dynamics and variability between the two land uses and to explore the relationship between the spatial variability of theta and theta(m), this study analyzed sets of surface theta measurements performed with an impedance soil moisture probe, collected 136 times during a period of one year in two transects covering different land uses, i.e., korshinsk peashrub transect (KPT) and bunge needlegrass transect (BNT), in a watershed of the Loess Plateau, China. Results showed that the temporal pattern of theta behaved similarly for the two land uses, with both relative wetter soils during wet period and relative drier soils during dry period recognized in BNT. Soil moisture tended to be temporally stable among different dates, and more stable patterns could be observed for dates with more similar soil water conditions. The magnitude of the spatial variation of theta in KPT was greater than that in ENT. For both land uses, the standard deviation (SD) of theta in general increased as theta(m) increased, a behavior that could be well described with a natural logarithmic function. Convex relationship of CV and theta(m) and the maximum CV for both land uses (43.5% in KPT and 41.0% in BNT) can, therefore, be ascertained. Geostatistical analysis showed that the range in KPT (9.1 m) was shorter than that in BNT (15.1 m). The nugget effects, the structured variability, hence the total variability increased as theta(m) increased. For both land uses, the spatial dependency in general increased with increasing theta(m). 2011 Elsevier B.V. All rights reserved.

Spatio-temporal modeling of agricultural yield data with an application to pricing crop insurance contracts

This article presents a statistical model of agricultural yield data based on a set of hierarchical Bayesian models that allows joint modeling of temporal and spatial autocorrelation. This method captures a comprehensive range of the various uncertainties involved in predicting crop insurance premium rates as opposed to the more traditional ad hoc, two-stage methods that are typically based on independent estimation and prediction. A panel data set of county-average yield data was analyzed for 290 counties in the State of Parana (Brazil) for the period of 1990 through 2002. Posterior predictive criteria are used to evaluate different model specifications. This article provides substantial improvements in the statistical and actuarial methods often applied to the calculation of insurance premium rates. These improvements are especially relevant to situations where data are limited.