Measuring the relative effect of factors affecting species distribution model predictions

1. Species distribution models are increasingly used to address conservation questions, so their predictive capacity requires careful evaluation. Previous studies have shown how individual factors used in model construction can affect prediction. Although some factors probably have negligible effects compared to others, their relative effects are largely unknown. 2. We introduce a general "virtual ecologist" framework to study the relative importance of factors involved in the construction of species distribution models. 3. We illustrate the framework by examining the relative importance of five key factors-a missing covariate, spatial autocorrelation due to a dispersal process in presences/absences, sample size, sampling design and modeling technique-in a real study framework based on plants in a mountain landscape at regional scale, and show that, for the parameter values considered here, most of the variation in prediction accuracy is due to sample size and modeling technique. Contrary to repeatedly reported concerns, spatial autocorrelation has only comparatively small effects. 4. This study shows the importance of using a nested statistical framework to evaluate the relative effects of factors that may affect species distribution models.

Cardiovascular magnetic resonance imaging techniques for motion-suppressed multi-dimensional coronary artery imaging with high spatial and temporal resolution

L'imagerie par résonance magnétique (IRM) peut fournir aux cardiologues des informations diagnostiques importantes sur l'état de la maladie de l'artère coronarienne dans les patients. Le défi majeur pour l'IRM cardiaque est de gérer toutes les sources de mouvement qui peuvent affecter la qualité des images en réduisant l'information diagnostique. Cette thèse a donc comme but de développer des nouvelles techniques d'acquisitions des images IRM, en changeant les techniques de compensation du mouvement, pour en augmenter l'efficacité, la flexibilité, la robustesse et pour obtenir plus d'information sur le tissu et plus d'information temporelle. Les techniques proposées favorisent donc l'avancement de l'imagerie des coronaires dans une direction plus maniable et multi-usage qui peut facilement être transférée dans l'environnement clinique. La première partie de la thèse s'est concentrée sur l'étude du mouvement des artères coronariennes sur des patients en utilisant la techniques d'imagerie standard (rayons x), pour mesurer la précision avec laquelle les artères coronariennes retournent dans la même position battement après battement (repositionnement des coronaires). Nous avons découvert qu'il y a des intervalles dans le cycle cardiaque, tôt dans la systole et à moitié de la diastole, où le repositionnement des coronaires est au minimum. En réponse nous avons développé une nouvelle séquence d'acquisition (T2-post) capable d'acquérir les données aussi tôt dans la systole. Cette séquence a été testée sur des volontaires sains et on a pu constater que la qualité de visualisation des artère coronariennes est égale à celle obtenue avec les techniques standard. De plus, le rapport signal sur bruit fourni par la séquence d'acquisition proposée est supérieur à celui obtenu avec les techniques d'imagerie standard. La deuxième partie de la thèse a exploré un paradigme d'acquisition des images cardiaques complètement nouveau pour l'imagerie du coeur entier. La technique proposée dans ce travail acquiert les données sans arrêt (free-running) au lieu d'être synchronisée avec le mouvement cardiaque. De cette façon, l'efficacité de la séquence d'acquisition est augmentée de manière significative et les images produites représentent le coeur entier dans toutes les phases cardiaques (quatre dimensions, 4D). Par ailleurs, l'auto-navigation de la respiration permet d'effectuer cette acquisition en respiration libre. Cette technologie rend possible de visualiser et évaluer l'anatomie du coeur et de ses vaisseaux ainsi que la fonction cardiaque en quatre dimensions et avec une très haute résolution spatiale et temporelle, sans la nécessité d'injecter un moyen de contraste. Le pas essentiel qui a permis le développement de cette technique est l'utilisation d'une trajectoire d'acquisition radiale 3D basée sur l'angle d'or. Avec cette trajectoire, il est possible d'acquérir continûment les données d'espace k, puis de réordonner les données et choisir les paramètres temporel des images 4D a posteriori. L'acquisition 4D a été aussi couplée avec un algorithme de reconstructions itératif (compressed sensing) qui permet d'augmenter la résolution temporelle tout en augmentant la qualité des images. Grâce aux images 4D, il est possible maintenant de visualiser les artères coronariennes entières dans chaque phase du cycle cardiaque et, avec les mêmes données, de visualiser et mesurer la fonction cardiaque. La qualité des artères coronariennes dans les images 4D est la même que dans les images obtenues avec une acquisition 3D standard, acquise en diastole Par ailleurs, les valeurs de fonction cardiaque mesurées au moyen des images 4D concorde avec les valeurs obtenues avec les images 2D standard. Finalement, dans la dernière partie de la thèse une technique d'acquisition a temps d'écho ultra-court (UTE) a été développée pour la visualisation in vivo des calcifications des artères coronariennes. Des études récentes ont démontré que les acquisitions UTE permettent de visualiser les calcifications dans des plaques athérosclérotiques ex vivo. Cepandent le mouvement du coeur a entravé jusqu'à maintenant l'utilisation des techniques UTE in vivo. Pour résoudre ce problème nous avons développé une séquence d'acquisition UTE avec trajectoire radiale 3D et l'avons testée sur des volontaires. La technique proposée utilise une auto-navigation 3D pour corriger le mouvement respiratoire et est synchronisée avec l'ECG. Trois échos sont acquis pour extraire le signal de la calcification avec des composants au T2 très court tout en permettant de séparer le signal de la graisse depuis le signal de l'eau. Les résultats sont encore préliminaires mais on peut affirmer que la technique développé peut potentiellement montrer les calcifications des artères coronariennes in vivo. En conclusion, ce travail de thèse présente trois nouvelles techniques pour l'IRM du coeur entier capables d'améliorer la visualisation et la caractérisation de la maladie athérosclérotique des coronaires. Ces techniques fournissent des informations anatomiques et fonctionnelles en quatre dimensions et des informations sur la composition du tissu auparavant indisponibles. CORONARY artery magnetic resonance imaging (MRI) has the potential to provide the cardiologist with relevant diagnostic information relative to coronary artery disease of patients. The major challenge of cardiac MRI, though, is dealing with all sources of motions that can corrupt the images affecting the diagnostic information provided. The current thesis, thus, focused on the development of new MRI techniques that change the standard approach to cardiac motion compensation in order to increase the efficiency of cardioavscular MRI, to provide more flexibility and robustness, new temporal information and new tissue information. The proposed approaches help in advancing coronary magnetic resonance angiography (MRA) in the direction of an easy-to-use and multipurpose tool that can be translated to the clinical environment. The first part of the thesis focused on the study of coronary artery motion through gold standard imaging techniques (x-ray angiography) in patients, in order to measure the precision with which the coronary arteries assume the same position beat after beat (coronary artery repositioning). We learned that intervals with minimal coronary artery repositioning occur in peak systole and in mid diastole and we responded with a new pulse sequence (T2~post) that is able to provide peak-systolic imaging. Such a sequence was tested in healthy volunteers and, from the image quality comparison, we learned that the proposed approach provides coronary artery visualization and contrast-to-noise ratio (CNR) comparable with the standard acquisition approach, but with increased signal-to-noise ratio (SNR). The second part of the thesis explored a completely new paradigm for whole- heart cardiovascular MRI. The proposed techniques acquires the data continuously (free-running), instead of being triggered, thus increasing the efficiency of the acquisition and providing four dimensional images of the whole heart, while respiratory self navigation allows for the scan to be performed in free breathing. This enabling technology allows for anatomical and functional evaluation in four dimensions, with high spatial and temporal resolution and without the need for contrast agent injection. The enabling step is the use of a golden-angle based 3D radial trajectory, which allows for a continuous sampling of the k-space and a retrospective selection of the timing parameters of the reconstructed dataset. The free-running 4D acquisition was then combined with a compressed sensing reconstruction algorithm that further increases the temporal resolution of the 4D dataset, while at the same time increasing the overall image quality by removing undersampling artifacts. The obtained 4D images provide visualization of the whole coronary artery tree in each phases of the cardiac cycle and, at the same time, allow for the assessment of the cardiac function with a single free- breathing scan. The quality of the coronary arteries provided by the frames of the free-running 4D acquisition is in line with the one obtained with the standard ECG-triggered one, and the cardiac function evaluation matched the one measured with gold-standard stack of 2D cine approaches. Finally, the last part of the thesis focused on the development of ultrashort echo time (UTE) acquisition scheme for in vivo detection of calcification in the coronary arteries. Recent studies showed that UTE imaging allows for the coronary artery plaque calcification ex vivo, since it is able to detect the short T2 components of the calcification. The heart motion, though, prevented this technique from being applied in vivo. An ECG-triggered self-navigated 3D radial triple- echo UTE acquisition has then been developed and tested in healthy volunteers. The proposed sequence combines a 3D self-navigation approach with a 3D radial UTE acquisition enabling data collection during free breathing. Three echoes are simultaneously acquired to extract the short T2 components of the calcification while a water and fat separation technique allows for proper visualization of the coronary arteries. Even though the results are still preliminary, the proposed sequence showed great potential for the in vivo visualization of coronary artery calcification. In conclusion, the thesis presents three novel MRI approaches aimed at improved characterization and assessment of atherosclerotic coronary artery disease. These approaches provide new anatomical and functional information in four dimensions, and support tissue characterization for coronary artery plaques.

Spatial radial maze procedures and setups to dissociate local and distal relational spatial frameworks in humans.

BACKGROUND: Radial maze tasks have been used to assess optimal foraging and spatial abilities in rodents. The spatial performance was based on a capacity to rely on a configuration of local and distant cues. We adapted maze procedures assessing the relative weight of local cues and distant landmarks for arm choice in humans. NEW METHOD: The procedure allowed testing memory of places in four experimental setups: a fingertip texture-groove maze, a tactile screen maze, a virtual radial maze and a walking size maze. During training, the four reinforced positions remained fixed relative to local and distal cues. During subsequent conflict trials, these frameworks were made conflictive in the prediction of reward locations. RESULTS: Three experiments showed that the relative weight of local and distal relational cues is affected by different factors such as cues' nature, visual access to the environment, real vs. virtual environment, and gender. A fourth experiment illustrated how a walking maze can be used with people suffering intellectual disability. COMPARISON WITH EXISTING METHODS: In our procedure, long-term (reference) and short-term (working) memory can be assessed. It is the first radial task adapted to human that enables dissociating local and distal cues, to provides an indication as to their relative salience. Our mazes are moveable and easily used in limited spaces. Tasks are performed with realistic and spontaneous though controlled exploratory movements. CONCLUSION: Our tasks enabled highlighting the use of different strategies. In a clinical perspective, considering the use of compensatory strategies should orient towards adapted behavioural rehabilitation.

Potential impact of the global climate changes on the spatial distribution of areas of risk for the occurrence of eucalyptus rust in Brazil

Rust, caused by Puccinia psidii, is one of the most important diseases affecting eucalyptus in Brazil. This pathogen causes disease in mini-clonal garden and in young plants in the field, especially in leaves and juvenile shoots. Favorable climate conditions for infection by this pathogen in eucalyptus include temperature between 18 and 25 ºC, together with at least 6-hour leaf wetness periods, for 5 to 7 consecutive days. Considering the interaction between the environment and the pathogen, this study aimed to evaluate the potential impact of global climate changes on the spatial distribution of areas of risk for the occurrence of eucalyptus rust in Brazil. Thus, monthly maps of the areas of risk for the occurrence of this disease were elaborated, considering the current climate conditions, based on a historic series between 1961 and 1990, and the future scenarios A2 and B2, predicted by IPCC. The climate conditions were classified into three categories, according to the potential risk for the disease occurrence, considering temperature (T) and air relative humidity (RH): i) high risk (18 < T < 25 ºC and RH > 90%); ii) medium risk (18 < T < 25 ºC and RH < 90%; T< 18 or T > 25 ºC and RH > 90%); and iii) low risk (T < 18 or T > 25 ºC and RH < 90%). Data about the future climate scenarios were supplied by GCM Change Fields. In this study, the simulation model Hadley Centers for Climate Prediction and Research (HadCm3) was adopted, using the software Idrisi 32. The obtained results led to the conclusion that there will be a reduction in the area favorable to eucalyptus rust occurrence, and such a reduction will be gradual for the decades of 2020, 2050 and 2080 but more marked in scenario A2 than in B2. However, it is important to point out that extensive areas will still be favorable to the disease development, especially in the coldest months of the year, i.e., June and July. Therefore, the zoning of areas and periods of higher occurrence risk, considering the global climate changes, becomes important knowledge for the elaboration of predicting models and an alert for the integrated management of this disease.

GPS receivers for georeferencing of spatial variability of soil attributes

The characterization of the spatial variability of soil attributes is essential to support agricultural practices in a sustainable manner. The use of geostatistics to characterize spatial variability of these attributes, such as soil resistance to penetration (RP) and gravimetric soil moisture (GM) is now usual practice in precision agriculture. The result of geostatistical analysis is dependent on the sample density and other factors according to the georeferencing methodology used. Thus, this study aimed to compare two methods of georeferencing to characterize the spatial variability of RP and GM as well as the spatial correlation of these variables. Sampling grid of 60 points spaced 20 m was used. For RP measurements, an electronic penetrometer was used and to determine the GM, a Dutch auger (0.0-0.1 m depth) was used. The samples were georeferenced using a GPS navigation receiver, Simple Point Positioning (SPP) with navigation GPS receiver, and Semi-Kinematic Relative Positioning (SKRP) with an L1 geodetic GPS receiver. The results indicated that the georeferencing conducted by PPS did not affect the characterization of spatial variability of RP or GM, neither the spatial structure relationship of these attributes.

Characteristics of urban heat island (UHI) in a high latitude coastal city - a case study of Turku, SW Finland

The term urban heat island (UHI) refers to the common situation in which the city is warmer than its rural surroundings. In this dissertation, the local climate, and especially the UHI, of the coastal city of Turku (182,000 inh.), SW Finland, was studied in different spatial and temporal scales. The crucial aim was to sort out the urban, topographical and water body impact on temperatures at different seasons and times of the day. In addition, the impact of weather on spatiotemporal temperature differences was studied. The relative importance of environmental factors was estimated with different modelling approaches and a large number of explanatory variables with various spatial scales. The city centre is the warmest place in the Turku area. Temperature excess relative to the coldest sites, i.e. rural areas about 10 kilometers to the NE from the centre, is on average 2 °C. Occasionally, the UHI intensity can be even 10 °C. The UHI does not prevail continuously in the Turku area, but occasionally the city centre can be colder than its surroundings. Then the term urban cool island or urban cold island (UCI) is used. The UCI is most common in daytime in spring and in summer, whereas during winter the UHI prevails throughout the day. On average, the spatial temperature differences are largest in summer, whereas the single extreme values are often observed in winter. The seasonally varying sea temperature causes the shift of relatively warm areas towards the coast in autumn and inland in spring. In the long term, urban land use was concluded to be the most important factor causing spatial temperature differences in the Turku area. The impact was mainly a warming one. The impact of water bodies was emphasised in spring and autumn, when the water temperature was relatively cold and warm, respectively. The impact of topography was on average the weakest, and was seen mainly in proneness of relatively low-lying places for cold air drainage during night-time. During inversions, however, the impact of topography was emphasised, occasionally outperforming those of urban land use and water bodies.

Spatial and temporal variation in fish populations and assemblages in coastal waters of the northern Baltic Proper

The distribution and traits of fish are of interest both ecologically and socio-economically. In this thesis, phenotypic and structural variation in fish populations and assemblages was studied on multiple spatial and temporal scales in shallow coastal areas in the archipelago of the northern Baltic Proper. In Lumparn basin in Åland Islands, the fish assemblage displayed significant seasonal variation in depth zone distribution. The results indicate that investigating both spatial and temporal variation in small scale is crucial for understanding patterns in fish distribution and community structure in large scale. The local population of Eurasian perch Perca fluviatilis L displayed habitat-specific morphological and dietary variation. Perch in the pelagic zone were on average deeper in their body shape than the littoral ones and fed on fish and benthic invertebrates. The results differ from previous studies conducted in freshwater habitats, where the pelagic perch typically are streamlined in body shape and zooplanktivorous. Stable isotopes of carbon and nitrogen differed between perch with different stomach contents, suggesting differentiation of individual diet preferences. In the study areas Lumparn and Ivarskärsfjärden in Åland Islands and Galtfjärden in Swedish east coast, the development in fish assemblages during the 2000’s indicated a general shift towards higher abundances of small-bodied lower-order consumers, especially cyprinids. For European pikeperch Sander lucioperca L., recent declines in adult fish abundances and high mortalities (Z = 1.06–1.16) were observed, which suggests unsustainably high fishing pressure on pikeperch. Based on the results it can be hypothesized that fishing has reduced the abundances of large predatory fish, which together with bottom-up forcing by eutrophication has allowed the lower-order consumer species to increase in abundances. This thesis contributes to the scientific understanding of aquatic ecosystems with new descriptions on morphological and dietary adaptations in perch in brackish water, and on the seasonal variation in small-scale spatial fish distribution. The results also demonstrate anthropogenic effects on coastal fish communities and underline the urgency of further reducing nutrient inputs and regulating fisheries in the Baltic Sea region.

Bird depredation of grapes in Niagara vineyards : a novel approach to identifying spatial and temporal trends /

Although local grape growers view bird depredation as a significant economic issue, the most recent research on the problem in the Niagara Peninsula is three decades old. Peer-reviewed publications on the subject are rare, and researchers have struggled to develop bird-damage assessment techniques useful for facilitating management programmes. I used a variation of Stevenson and Virgo's (1971) visual estimation procedure to quantify spatial and temporal trends in bird damage to grapes within single vineyard plots at two locations near St. Catharines, Ontario. I present a novel approach to managing the rank-data from visual estimates, which is unprecedented in its sensitivity to spatial trends in bird damage. I also review its valid use in comparative statistical analysis. Spatial trends in 3 out of 4 study plots confirmed a priori predictions about localisation in bird damage based on optimal foraging from a central location (staging area). Damage to grape clusters was: (1) greater near the edges of vineyard plots and decreased with distance towards the center, (2) greater in areas adjacent to staging areas for birds, and (3) vertically stratified, with upper-tier clusters sustaining more damage than lower-tier clusters. From a management perspective, this predictive approach provides vineyard owners with the ability to identify the portions of plots likely to be most susceptible to bird damage, and thus the opportunity to focus deterrent measures in these areas. Other management considerations at Henry of Pelham were: (1) wind damage to ice-wine Riesling and Vidal was much higher than bird damage, (2) plastic netting with narrow mesh provided more effective protection agsiinst birds than nylon netting with wider mesh, and (3) no trends in relative susceptibility of varietals by colour (red vs green) were evident.

Aspects of spatial and habitat ecology of multiple Anopheles species (Diptera: Culicidae): malaria vectors in the highlands and foothills of Ecuador

The resurgence of malaria in highland regions of Africa, Oceania and recently in South America underlines the importance of the study of the ecology of highland mosquito vectors of malaria. Since the incidence of malaria is limited by the distribution of its vectors, the purpose of this PhD thesis was to examine aspects of the ecology of Anopheles mosquitoes in the Andes of Ecuador, South America. A historical literature and archival data review (Chapter 2) indicated that Anopheles pseudopunctipennis transmitted malaria in highland valleys of Ecuador prior to 1950, although it was eliminated through habitat removal and the use of chemical insecticides. Other anopheline species were previously limited to low-altitude regions, except in a few unconfirmed cases. A thorough larval collection effort (n=438 attempted collection sites) in all road-accessible parts of Ecuador except for the lowland Amazon basin was undertaken between 2008 - 2010 (Chapter 3). Larvae were identified morphologically and using molecular techniques (mitochondrial COl gene), and distribution maps indicated that all five species collected (Anopheles albimanus, An. pseudopunctipennis, Anopheles punctimacula, Anopheles oswaldoi s.l. and Anopheles eiseni) were more widespread throughout highland regions than previously recorded during the 1940s, with higher maximum altitudes for all except An. pseudopunctipennis (1541 m, 1930 m, 1906 m, 1233 m and 1873 m, respectively). During larval collections, to characterize species-specific larval habitat, a variety of abiotic and biotic habitat parameters were measured and compared between species-present and species-absent sites using chi-square tests and stepwise binary logistic regression analyses (Chapter 4). An. albimanus was significantly associated with permanent pools with sand substrates and An. pseudopunctipennis with gravel and boulder substrates. Both species were significantly associated with floating cyanobacterial mats and warmer temperatures, which may limit their presence in cooler highland regions. Anopheles punctimacula was collected more often than expected from algae-free, shaded pools with higher-than-average calculated dissolved oxygen. Anopheles oswaldoi s.l., the species occurring on the Amazonian side of the Andes, was associated with permanent, anthropogenic habitats such as roadside ditches and ponds. To address the hypothesis that human land use change is responsible for the emergence of multiple highland Anopheles species by creating larval habitat, common land uses in the western Andes were surveyed for standing water and potential larval habitat suitability (Chapter 5). Rivers and road edges provided large amounts of potentially suitable anopheline habitat in the western Andes, while cattle pasture also created potentially suitable habitat in irrigation canals and watering ponds. Other common land uses surveyed (banana farms, sugarcane plantations, mixed tree plantations, and empty lots) were usually established on steep slopes and had very little standing water present. Using distribution and larval habitat data, a GIS-based larval habitat distribution model for the common western species was constructed in ArcGIS v.l 0 (ESRI 2010) using derived data layers from field measurements and other sources (Chapter 6). The additive model predicted 76.4 - 97.9% of the field-observed collection localities of An. albimanus, An. pseudopunctipennis and An. punctimacula, although it could not accurately distinguish between species-absent and speciespresent sites due to its coarse scale. The model predicted distributional expansion and/or shift of one or more anopheline species into the following highland valleys with climate warming: Mira/Chota, Imbabura province, Tumbaco, Pichincha province, Pallatanga and Sibambe, Chimborazo province, and Yungilla, Azuay province. These valleys may serve as targeted sites of future monitoring to prevent highland epidemics of malaria. The human perceptions of malaria and mosquitoes in relation to land management practices were assessed through an interview-based survey (n=262) in both highlands and lowlands, of male and female land owners and managers of five property types (Chapter 7). Although respondents had a strong understanding of where the disease occurs in their own country and of the basic relationship among standing water, mosquitoes and malaria, about half of respondents in potential risk areas denied the current possibility of malaria infection on their own property. As well, about half of respondents with potential anopheline larval habitat did not report its presence, likely due to a highly specific definition of suitable mosquito habitat. Most respondents who are considered at risk of malaria currently use at least one type of mosquito bite prevention, most commonly bed nets. In conclusion, this interdisciplinary thesis examines the occurrence of Anopheles species in the lowland transition area and highlands in Ecuador, from a historic, geographic, ecological and sociological perspective.

Analyse des réponses neuronales du cortex visuel primaire du chat à la fréquence spatiale suite à des adaptations répétées

Les neurones du cortex visuel primaire (aire 17) du chat adulte répondent de manière sélective à différentes propriétés d’une image comme l’orientation, le contraste ou la fréquence spatiale. Cette sélectivité se manifeste par une réponse sous forme de potentiels d’action dans les neurones visuels lors de la présentation d’une barre lumineuse de forme allongée dans les champs récepteurs de ces neurones. La fréquence spatiale (FS) se mesure en cycles par degré (cyc./deg.) et se définit par la quantité de barres lumineuses claires et sombres présentées à une distance précise des yeux. Par ailleurs, jusqu’à récemment, l’organisation corticale chez l’adulte était considérée immuable suite à la période critique post-natale. Or, lors de l'imposition d'un stimulus non préféré, nous avons observé un phénomène d'entrainement sous forme d'un déplacement de la courbe de sélectivité à la suite de l'imposition d'une FS non-préférée différente de la fréquence spatiale optimale du neurone. Une deuxième adaptation à la même FS non-préférée induit une réponse neuronale différente par rapport à la première imposition. Ce phénomène de "gain cortical" avait déjà été observé dans le cortex visuel primaire pour ce qui est de la sélectivité à l'orientation des barres lumineuses, mais non pour la fréquence spatiale. Une telle plasticité à court terme pourrait être le corrélat neuronal d'une modulation de la pondération relative du poids des afférences synaptiques.

Spatial and temporal characteristics of rain intensity in the peninsular Malaysia using TRMM rain rate

The present study is focused on the intensity distribution of rainfall in different classes and their contribution to the total seasonal rainfall. In addition, we studied the spatial and diurnal variation of the rainfall in the study areas. For the present study, we retrieved data from TRMM (Tropical Rain Measuring Mission) rain rate available in every 3 h temporal and 25 km spatial resolutions. Moreover, station rainfall data is used to validate the TRMM rain rate and found significant correlation between them (linear correlation coefficients are 0.96, 0.85, 0.75 and 0.63 for the stations Kota Bharu, Senai, Cameron highlands and KLIA, respectively). We selected four areas in the Peninsular Malaysia and they are south coastal, east coastal, west coastal and highland regions. Diurnal variation of frequency of rain occurrence is different for different locations. We noticed bimodal variation in the coastal areas in most of the seasons and unimodal variation in the highland/inland area. During the southwest monsoon period in the west coastal stations, there is no distinct diurnal variation. The distribution of different intensity classes during different seasons are explained in detail in the results

Quantifying rock fabrics: a test of autocorrelation of the spatial distribution of cristals

A novel test of spatial independence of the distribution of crystals or phases in rocks based on compositional statistics is introduced. It improves and generalizes the common joins-count statistics known from map analysis in geographic information systems. Assigning phases independently to objects in RD is modelled by a single-trial multinomial random function Z(x), where the probabilities of phases add to one and are explicitly modelled as compositions in the K-part simplex SK. Thus, apparent inconsistencies of the tests based on the conventional joins{count statistics and their possibly contradictory interpretations are avoided. In practical applications we assume that the probabilities of phases do not depend on the location but are identical everywhere in the domain of de nition. Thus, the model involves the sum of r independent identical multinomial distributed 1-trial random variables which is an r-trial multinomial distributed random variable. The probabilities of the distribution of the r counts can be considered as a composition in the Q-part simplex SQ. They span the so called Hardy-Weinberg manifold H that is proved to be a K-1-affine subspace of SQ. This is a generalisation of the well-known Hardy-Weinberg law of genetics. If the assignment of phases accounts for some kind of spatial dependence, then the r-trial probabilities do not remain on H. This suggests the use of the Aitchison distance between observed probabilities to H to test dependence. Moreover, when there is a spatial uctuation of the multinomial probabilities, the observed r-trial probabilities move on H. This shift can be used as to check for these uctuations. A practical procedure and an algorithm to perform the test have been developed. Some cases applied to simulated and real data are presented. Key words: Spatial distribution of crystals in rocks, spatial distribution of phases, joins-count statistics, multinomial distribution, Hardy-Weinberg law, Hardy-Weinberg manifold, Aitchison geometry

The spatial distribution and evolution characteristics of North Atlantic cyclones

A climatology of extratropical cyclones is produced using an objective method of identifying cyclones based on gradients of 1-km height wet-bulb potential temperature. Cyclone track and genesis density statistics are analyzed and this method is found to compare well with other cyclone identification methods. The North Atlantic storm track is reproduced along with the major regions of genesis. Cyclones are grouped according to their genesis location and the corresponding lysis regions are identified. Most of the cyclones that cross western Europe originate in the east Atlantic where the baroclinicity and the sea surface temperature gradients are weak compared to the west Atlantic. East Atlantic cyclones also have higher 1-km height relative vorticity and lower mean sea level pressure at their genesis point than west Atlantic cyclones. This is consistent with the hypothesis that they are secondary cyclones developing on the trailing fronts of preexisting “parent” cyclones. The evolution characteristics of composite west and east Atlantic cyclones have been compared. The ratio of their upper- to lower-level forcing indicates that type B cyclones are predominant in both the west and east Atlantic, with strong upper- and lower-level features. Among the remaining cyclones, there is a higher proportion of type C cyclones in the east Atlantic, whereas types A and C are equally frequent in the west Atlantic.

Tropospheric planetary wave dynamics and mixture modeling: Two preferred regimes and a regime shift

Investigation of preferred structures of planetary wave dynamics is addressed using multivariate Gaussian mixture models. The number of components in the mixture is obtained using order statistics of the mixing proportions, hence avoiding previous difficulties related to sample sizes and independence issues. The method is first applied to a few low-order stochastic dynamical systems and data from a general circulation model. The method is next applied to winter daily 500-hPa heights from 1949 to 2003 over the Northern Hemisphere. A spatial clustering algorithm is first applied to the leading two principal components (PCs) and shows significant clustering. The clustering is particularly robust for the first half of the record and less for the second half. The mixture model is then used to identify the clusters. Two highly significant extratropical planetary-scale preferred structures are obtained within the first two to four EOF state space. The first pattern shows a Pacific-North American (PNA) pattern and a negative North Atlantic Oscillation (NAO), and the second pattern is nearly opposite to the first one. It is also observed that some subspaces show multivariate Gaussianity, compatible with linearity, whereas others show multivariate non-Gaussianity. The same analysis is also applied to two subperiods, before and after 1978, and shows a similar regime behavior, with a slight stronger support for the first subperiod. In addition a significant regime shift is also observed between the two periods as well as a change in the shape of the distribution. The patterns associated with the regime shifts reflect essentially a PNA pattern and an NAO pattern consistent with the observed global warming effect on climate and the observed shift in sea surface temperature around the mid-1970s.