SPATIAL PATTERN OF FOUR FOREST COMMUNITIES BY RIPLEY'S K FUNCTION IN SAO PAULO STATE

Investigating tree's spatial patterns according to their size classes and according to their more abundant species can provide evidences about the structure of the vegetal community, since the spatial pattern is a key question for forestry ecology studies. The tree spatial organization patterns on the environment depend on several ecological processes and on the specific characteristics of each environment, so that the best understanding of this frame provides important elements for the knowledge on forestry formation. This paper aimed to study tree spatial patterns, according to the diameter classes and from four most abundant species in different forests, in order to provide evidences regarding to the ecology of each vegetal community. The spatial pattern description in each forestry formation was developed using Ripley's K function. The studied forestry formations were: Ombrophilous Forest, Cerradao, Seasonal Forest and Restinga Forest. In this work, a 10.24 ha plot was installed in each forestry formation, and every tree, with a circumference at breast height (CBH) larger than 15 cm were measured, georeferenced and identified. The obtained data highlights the aggregated character in tropical forests, as observed in every studied forest. The 'Cerraddo' and 'Restinga' forest trees showed close aggregate patterns. In the Ombrophilous forest, for all distance scales, the aggregate pattern was meaningful. In the seasonal forest, a random tendency was observed, although a meaningful aggregation was observed in all short distances. The spatial pattern by diameter classes was generally aggregated for trees smaller than 10 cm of diameter and between 10 and 20 cm and random for the others, proving the existence of a tendency which in young trees is more aggregated than in old ones. The spatial pattern of the dominant species is always strongly similar to the general pattern of each forestry formation. The differences between the spatial patterns of two or three coincident species, among the forestry formations, indicate that its pattern is influenced by each different environment. This stands out the importance of its self-ecology and of its ecological processes, intrinsic of each community that can explain the observed patterns.

Nuevas herramientas para el estudio de la interacción entre especies en el espacio y en el tiempo

Las funciones de segundo orden son cada vez más empleadas en el análisis de procesos ecológicos. En este trabajo presentamos dos funciones de 2º orden desarrolladas recientemente que permiten analizar la interacción espacio-temporal entre dos especies o tipos funcionales de individuos. Estas funciones han sido desarrolladas para el estudio de interacciones entre especies en masas forestales a partir de la actual distribución diamétrica de los árboles. La primera de ellas es la función bivariante para procesos de puntos con marca Krsmm, que permite analizar la correlación espacial de una variable entre los individuos pertenecientes a dos especies en función de la distancia. La segunda es la función de reemplazo , que permite analizar la asociación entre los individuos pertenecientes a dos especies en función de la diferencia entre sus diámetros u otra variable asociada a dichos individuos. Para mostrar el comportamiento de ambas funciones en el análisis de sistemas forestales en los que operan diferentes procesos ecológicos se presentan tres casos de estudio: una masa mixta de Pinus pinea L. y Pinus pinaster Ait. en la Meseta Norte, un bosque de niebla de la Región Tropical Andina y el ecotono entre las masas de Quercus pyrenaica Willd. y Pinus sylvestris L. en el Sistema Central, en los que tanto la función Krsmm como la función r se utilizan para analizar la dinámica forestal a partir de parcelas experimentales con todos los árboles localizados y de parcelas de inventario.

Spatial distribution of triatomines in domiciles of an urban area of the Brazilian Southeast Region

Reports of triatomine infestation in urban areas have increased. We analysed the spatial distribution of infestation by triatomines in the urban area of Diamantina, in the state of Minas Gerais, Brazil. Triatomines were obtained by community-based entomological surveillance. Spatial patterns of infestation were analysed by Ripley’s K function and Kernel density estimator. Normalised difference vegetation index (NDVI) and land cover derived from satellite imagery were compared between infested and uninfested areas. A total of 140 adults of four species were captured (100 Triatoma vitticeps, 25 Panstrongylus geniculatus, 8 Panstrongylus megistus, and 7 Triatoma arthurneivai specimens). In total, 87.9% were captured within domiciles. Infection by trypanosomes was observed in 19.6% of 107 examined insects. The spatial distributions of T. vitticeps, P. geniculatus, T. arthurneivai, and trypanosome-positive triatomines were clustered, occurring mainly in peripheral areas. NDVI values were statistically higher in areas infested by T. vitticeps and P. geniculatus. Buildings infested by these species were located closer to open fields, whereas infestations of P. megistus and T. arthurneivai were closer to bare soil. Human occupation and modification of natural areas may be involved in triatomine invasion, exposing the population to these vectors.

Tree species community spatial structure in a terra firme Amazon forest, Brazil

All trees with diameter at breast height dbh >= 10.0 cm were stem-mapped in a "terra firme" tropical rainforest in the Brazilian Amazon, at the EMBRAPA Experimental Site, Manaus, Brazil. Specifically, the relationships of tree species with soil properties were determined by using canonical correspondence analyses based on nine soil variables and 68 tree species. From the canonical correspondence analyses, the species were grouped into two groups: one where species occur mainly in sandy sites, presenting low organic matter content; and another one where species occur mainly in dry and clayey sites. Hence, we used Ripley's K function to analyze the distribution of species in 32 plots ranging from 2,500 m(2) to 20,000 m(2) to determine whether each group presents some spatial aggregation as a soil variations result. Significant spatial aggregation for the two groups was found only at over 10,000 m(2) sampling units, particularly for those species found in clayey soils and drier environments, where the sampling units investigated seemed to meet the species requirements. Soil variables, mediated by topographic positions had influenced species spatial aggregation, mainly in an intermediate to large distances varied range (>= 20 m). Based on our findings, we conclude that environmental heterogeneity and 10,000 m(2) minimum sample unit sizes should be considered in forest dynamic studies in order to understand the spatial processes structuring the "terra firme" tropical rainforest in Brazilian Amazon.

Spatial clustering analysis of the foot-and-mouth disease outbreaks in Mato Grosso do Sul state, Brazil - 2005

In the southern region of Mato Grosso do Sul state, Brazil, a foot-and-mouth disease (FMD) epidemic started in September 2005. A total of 33 outbreaks were detected and 33,741 FMD-susceptible animals were slaughtered and destroyed. There were no reports of FMD cases in other species than bovines. Based on the data of this epidemic, it was carried out an analysis using the K-function and it was observed spatial clustering of outbreaks within a range of 25km. This observation may be related to the dynamics of foot-and-mouth disease spread and to the measures undertaken to control the disease dissemination. The control measures were effective once the disease did not spread to farms more than 47 km apart from the initial outbreaks.

Macropore sheath: quantification of plant root and soil macropore association

Plants require roots to supply water, nutrients and oxygen for growth. The spatial distribution of roots in relation to the macropore structure of the soil in which they are growing influences how effective they are at accessing these resources. A method for quantifying root-macropore associations from horizontal soil sections is illustrated using two black vertisols from the Darling Downs, Queensland, Australia. Two-dimensional digital images were obtained of the macropore structure and root distribution for an area 55 x 55 mm at a resolution of 64 mu m. The spatial distribution of roots was quantified over a range of distances using the K-function. In all specimens, roots were shown to be clustered at short distances (1-10 mm) becoming more random at longer distances. Root location in relation to macropores was estimated using the function describing the distance of each root to the nearest macropore. From this function, a summary variable, termed the macropore sheath, was defined. The macropore sheath is the distance from macropores within which 80% of roots are located. Measured root locations were compared to random simulations of root distribution to establish if there was a preferential association between roots and macropores. More roots were found in and around macropores than expected at random.

Spatial analysis and investigation of fire events occurrences in the Valencian Community, Spain

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

Spatial distribution of triatomines in domiciles of an urban area of the Brazilian Southeast Region

Reports of triatomine infestation in urban areas have increased. We analysed the spatial distribution of infestation by triatomines in the urban area of Diamantina, in the state of Minas Gerais, Brazil. Triatomines were obtained by community-based entomological surveillance. Spatial patterns of infestation were analysed by Ripley’s K function and Kernel density estimator. Normalised difference vegetation index (NDVI) and land cover derived from satellite imagery were compared between infested and uninfested areas. A total of 140 adults of four species were captured (100 Triatoma vitticeps, 25Panstrongylus geniculatus, 8 Panstrongylus megistus, and 7 Triatoma arthurneivai specimens). In total, 87.9% were captured within domiciles. Infection by trypanosomes was observed in 19.6% of 107 examined insects. The spatial distributions ofT. vitticeps, P. geniculatus, T. arthurneivai, and trypanosome-positive triatomines were clustered, occurring mainly in peripheral areas. NDVI values were statistically higher in areas infested by T. vitticeps and P. geniculatus. Buildings infested by these species were located closer to open fields, whereas infestations of P. megistus andT. arthurneivai were closer to bare soil. Human occupation and modification of natural areas may be involved in triatomine invasion, exposing the population to these vectors.

Pharmacokinetic and pharmacodynamic evaluation of valganciclovir in solid organ transplant patients

Goal: To validate oral vatgancictovir (VGC) in the prophylaxis of CMV infection in Lung (Lu) and Liver (L) recipients and in the treatment of CMV infection/disease in solid organ transplant recipients, using pharmacokinetic and pharmacodynamic studies in comparison with i/v gancicLovir (GCV). Methods: patients undergoing organ transpLantation donor or recipient CMV-seropositive receiving VGC prophylaxis for a period of 3 months (D+/R- Lung recipients, 6 months) were enroLLed. Heart (H), Lu, and L recipients received 900mg VGC q.d., adjusted to kidney (K) function. No K recipients received more than 450mg of VGC q.d. GCV trough (Ctrough) and peak (Cpeak = 3 hours after drug administration) LeveLs, and CMV DNA were measured at 7, 30, and 60 days post-transpLant (prophyLactic study). Patients who developed CMV infection/disease after stopping prophylaxis were treated with VGC (1800mg per day adjusted to K function and GCV blood LeveLs). GCV trough and peak LeveLs, and CMV DNA were measured weekly for the first 3 weeks and biweekly thereafter, until therapy cessation (therapeutic study). PLasma concentration of GCV is measured by HPLC. Results: In the first 8 prophyLaxed patients (6 K, and 1 L and 1 H transplant recipient) of 450mg VGC q.d., the average GCV concentration was 0.5±0.3 mg/t at trough, and 3.9±l.0mg/t 3 hours after administration. Inter-patient variability was substantiaL, especiaLLy for Ctrough (63% of total variance), which correlated with the patient's estimated gtomerutar filtration rate (r square = 42%). No CMV DNA was detected during VGC prophy- Laxis. Two patients (1 H and 1 L) were treated for Late CMV disease. Average GCV Cpeak were 8.9±2.3 mg/L and 4.6±0.5 rag/L, and GCV Ctrough were 2.0±0.9 mg/t and 1.6±0.2 mg/t respectively in each patient during induction phase. VGC treatment afforded a decrease in CMV DNA from 5.2 and 4.4 Log copies/10E6 cettutes at week 0, to 3.9 and 3.0 at week 1, and 3.3 and 2.1 at week 3, respectively.

Spatial Scale of Clustering of Motor Vehicle Crash Types and Appropriate Countermeasures, 2009

This project analyzes the characteristics and spatial distributions of motor vehicle crash types in order to evaluate the degree and scale of their spatial clustering. Crashes occur as the result of a variety of vehicle, roadway, and human factors and thus vary in their clustering behavior. Clustering can occur at a variety of scales, from the intersection level, to the corridor level, to the area level. Conversely, other crash types are less linked to geographic factors and are more spatially “random.” The degree and scale of clustering have implications for the use of strategies to promote transportation safety. In this project, Iowa's crash database, geographic information systems, and recent advances in spatial statistics methodologies and software tools were used to analyze the degree and spatial scale of clustering for several crash types within the counties of the Iowa Northland Regional Council of Governments. A statistical measure called the K function was used to analyze the clustering behavior of crashes. Several methodological issues, related to the application of this spatial statistical technique in the context of motor vehicle crashes on a road network, were identified and addressed. These methods facilitated the identification of crash clusters at appropriate scales of analysis for each crash type. This clustering information is useful for improving transportation safety through focused countermeasures directly linked to crash causes and the spatial extent of identified problem locations, as well as through the identification of less location-based crash types better suited to non-spatial countermeasures. The results of the K function analysis point to the usefulness of the procedure in identifying the degree and scale at which crashes cluster, or do not cluster, relative to each other. Moreover, for many individual crash types, different patterns and processes and potentially different countermeasures appeared at different scales of analysis. This finding highlights the importance of scale considerations in problem identification and countermeasure formulation.

Rockfall detection from LIDAR point clouds: A clustering approach using R

The analysis of rockfall characteristics and spatial distribution is fundamental to understand and model the main factors that predispose to failure. In our study we analysed LiDAR point clouds aiming to: (1) detect and characterise single rockfalls; (2) investigate their spatial distribution. To this end, different cluster algorithms were applied: 1a) Nearest Neighbour Clutter Removal (NNCR) in combination with the Expectation?Maximization (EM) in order to separate feature points from clutter; 1b) a density based algorithm (DBSCAN) was applied to isolate the single clusters (i.e. the rockfall events); 2) finally we computed the Ripley's K-function to investigate the global spatial pattern of the extracted rockfalls. The method allowed proper identification and characterization of more than 600 rockfalls occurred on a cliff located in Puigcercos (Catalonia, Spain) during a time span of six months. The spatial distribution of these events proved that rockfall were clustered distributed at a welldefined distance-range. Computations were carried out using R free software for statistical computing and graphics. The understanding of the spatial distribution of precursory rockfalls may shed light on the forecasting of future failures.

Spatial pattern of landslides in Swiss Rhone valley

The present study analyses the spatial pattern of quaternary gravitational slope deformations (GSD) and historical/present-day instabilities (HPI) inventoried in the Swiss Rhone Valley. The main objective is to test if these events are clustered (spatial attraction) or randomly distributed (spatial independency). Moreover, analogies with the cluster behaviour of earthquakes inventoried in the same area were examined. The Ripley's K-function was applied to measure and test for randomness. This indicator allows describing the spatial pattern of a point process at increasing distance values. To account for the non-constant intensity of the geological phenomena, a modification of the K-function for inhomogeneous point processes was adopted. The specific goal is to explore the spatial attraction (i.e. cluster behaviour) among landslide events and between gravitational slope deformations and earthquakes. To discover if the two classes of instabilities (GSD and HPI) are spatially independently distributed, the cross K-function was computed. The results show that all the geological events under study are spatially clustered at a well-defined distance range. GSD and HPI show a similar pattern distribution with clusters in the range 0.75?9 km. The cross K-function reveals an attraction between the two classes of instabilities in the range 0?4 km confirming that HPI are more prone to occur within large-scale slope deformations. The K-function computed for GSD and earthquakes indicates that both present a cluster tendency in the range 0?10 km, suggesting that earthquakes could represent a potential predisposing factor which could influence the GSD distribution.

Analysis of space-time environmental data in complex regions

This thesis develops a comprehensive and a flexible statistical framework for the analysis and detection of space, time and space-time clusters of environmental point data. The developed clustering methods were applied in both simulated datasets and real-world environmental phenomena; however, only the cases of forest fires in Canton of Ticino (Switzerland) and in Portugal are expounded in this document. Normally, environmental phenomena can be modelled as stochastic point processes where each event, e.g. the forest fire ignition point, is characterised by its spatial location and occurrence in time. Additionally, information such as burned area, ignition causes, landuse, topographic, climatic and meteorological features, etc., can also be used to characterise the studied phenomenon. Thereby, the space-time pattern characterisa- tion represents a powerful tool to understand the distribution and behaviour of the events and their correlation with underlying processes, for instance, socio-economic, environmental and meteorological factors. Consequently, we propose a methodology based on the adaptation and application of statistical and fractal point process measures for both global (e.g. the Morisita Index, the Box-counting fractal method, the multifractal formalism and the Ripley's K-function) and local (e.g. Scan Statistics) analysis. Many measures describing the space-time distribution of environmental phenomena have been proposed in a wide variety of disciplines; nevertheless, most of these measures are of global character and do not consider complex spatial constraints, high variability and multivariate nature of the events. Therefore, we proposed an statistical framework that takes into account the complexities of the geographical space, where phenomena take place, by introducing the Validity Domain concept and carrying out clustering analyses in data with different constrained geographical spaces, hence, assessing the relative degree of clustering of the real distribution. Moreover, exclusively to the forest fire case, this research proposes two new methodologies to defining and mapping both the Wildland-Urban Interface (WUI) described as the interaction zone between burnable vegetation and anthropogenic infrastructures, and the prediction of fire ignition susceptibility. In this regard, the main objective of this Thesis was to carry out a basic statistical/- geospatial research with a strong application part to analyse and to describe complex phenomena as well as to overcome unsolved methodological problems in the characterisation of space-time patterns, in particular, the forest fire occurrences. Thus, this Thesis provides a response to the increasing demand for both environmental monitoring and management tools for the assessment of natural and anthropogenic hazards and risks, sustainable development, retrospective success analysis, etc. The major contributions of this work were presented at national and international conferences and published in 5 scientific journals. National and international collaborations were also established and successfully accomplished. -- Cette thèse développe une méthodologie statistique complète et flexible pour l'analyse et la détection des structures spatiales, temporelles et spatio-temporelles de données environnementales représentées comme de semis de points. Les méthodes ici développées ont été appliquées aux jeux de données simulées autant qu'A des phénomènes environnementaux réels; nonobstant, seulement le cas des feux forestiers dans le Canton du Tessin (la Suisse) et celui de Portugal sont expliqués dans ce document. Normalement, les phénomènes environnementaux peuvent être modélisés comme des processus ponctuels stochastiques ou chaque événement, par ex. les point d'ignition des feux forestiers, est déterminé par son emplacement spatial et son occurrence dans le temps. De plus, des informations tels que la surface bru^lée, les causes d'ignition, l'utilisation du sol, les caractéristiques topographiques, climatiques et météorologiques, etc., peuvent aussi être utilisées pour caractériser le phénomène étudié. Par conséquent, la définition de la structure spatio-temporelle représente un outil puissant pour compren- dre la distribution du phénomène et sa corrélation avec des processus sous-jacents tels que les facteurs socio-économiques, environnementaux et météorologiques. De ce fait, nous proposons une méthodologie basée sur l'adaptation et l'application de mesures statistiques et fractales des processus ponctuels d'analyse global (par ex. l'indice de Morisita, la dimension fractale par comptage de boîtes, le formalisme multifractal et la fonction K de Ripley) et local (par ex. la statistique de scan). Des nombreuses mesures décrivant les structures spatio-temporelles de phénomènes environnementaux peuvent être trouvées dans la littérature. Néanmoins, la plupart de ces mesures sont de caractère global et ne considèrent pas de contraintes spatiales com- plexes, ainsi que la haute variabilité et la nature multivariée des événements. A cet effet, la méthodologie ici proposée prend en compte les complexités de l'espace géographique ou le phénomène a lieu, à travers de l'introduction du concept de Domaine de Validité et l'application des mesures d'analyse spatiale dans des données en présentant différentes contraintes géographiques. Cela permet l'évaluation du degré relatif d'agrégation spatiale/temporelle des structures du phénomène observé. En plus, exclusif au cas de feux forestiers, cette recherche propose aussi deux nouvelles méthodologies pour la définition et la cartographie des zones périurbaines, décrites comme des espaces anthropogéniques à proximité de la végétation sauvage ou de la forêt, et de la prédiction de la susceptibilité à l'ignition de feu. A cet égard, l'objectif principal de cette Thèse a été d'effectuer une recherche statistique/géospatiale avec une forte application dans des cas réels, pour analyser et décrire des phénomènes environnementaux complexes aussi bien que surmonter des problèmes méthodologiques non résolus relatifs à la caractérisation des structures spatio-temporelles, particulièrement, celles des occurrences de feux forestières. Ainsi, cette Thèse fournit une réponse à la demande croissante de la gestion et du monitoring environnemental pour le déploiement d'outils d'évaluation des risques et des dangers naturels et anthro- pogéniques. Les majeures contributions de ce travail ont été présentées aux conférences nationales et internationales, et ont été aussi publiées dans 5 revues internationales avec comité de lecture. Des collaborations nationales et internationales ont été aussi établies et accomplies avec succès.

Reinforcing and expanding the predictions of the disturbance vicariance hypothesis in Amazonian harlequin frogs: A molecular phylogenetic and climate envelope modelling approach

The disturbance vicariance hypothesis (DV) has been proposed to explain speciation in Amazonia, especially its edge regions, e. g. in eastern Guiana Shield harlequin frogs (Atelopus) which are suggested to have derived from a cool-adapted Andean ancestor. In concordance with DV predictions we studied that (i) these amphibians display a natural distribution gap in central Amazonia; (ii) east of this gap they constitute a monophyletic lineage which is nested within a pre-Andean/western clade; (iii) climate envelopes of Atelopus west and east of the distribution gap show some macroclimatic divergence due to a regional climate envelope shift; (iv) geographic distributions of climate envelopes of western and eastern Atelopus range into central Amazonia but with limited spatial overlap. We tested if presence and apparent absence data points of Atelopus were homogenously distributed with Ripley's K function. A molecular phylogeny (mitochondrial 16S rRNA gene) was reconstructed using Maximum Likelihood and Bayesian Inference to study if Guianan Atelopus constitute a clade nested within a larger genus phylogeny. We focused on climate envelope divergence and geographic distribution by computing climatic envelope models with MaxEnt based on macroscale bioclimatic parameters and testing them by using Schoener's index and modified Hellinger distance. We corroborated existing DV predictions and, for the first time, formulated new DV predictions aiming on species' climate envelope change. Our results suggest that cool-adapted Andean Atelopus ancestors had dispersed into the Amazon basin and further onto the eastern Guiana Shield where, under warm conditions, they were forced to change climate envelopes. © 2010 The Author(s).

The spatial distribution of Copernicia alba (Morong) in the district of Bahia Negra, Paraguay

The Humid Chaco of Northeast Paraguay harbors monoculture palm savannas in which Copernicia alba is the only dominant overstory species. The study’s objective was to provide the complete spatial distribution of a simple ecosystem lacking confounding factors of overstory competition and changes in slope. Palms within six, 50 x 50m plots were marked by their GPS location and measured for dbh and total stem height. The spatial distribution was individually analyzed for each plot at the local scale up to 12 m using Ripley’s K test. For the total population including juvenile and adult plants, the sample plots contained both random and clustered distribution patterns. In each of the six plots, the juvenile populations exhibited more clustered patterns than the adult population of each plot.