Mycorrhizal preferences and fine spatial structure of the epiphytic orchid Epidendrum rhopalostele

• Premise of the study: The presence of compatible fungi is necessary for epiphytic orchid recruitment. Thus, identifying associated mycorrhizal fungi at the population level is essential for orchid conservation. Recruitment patterns may also be conditioned by factors such as seed dispersal range and specific environmental characteristics. • Methods: In a forest plot, all trees with a diameter at breast height >1 cm and all individuals of the epiphytic orchid Epidendrum rhopalostele were identified and mapped. Additionally, one flowering individual of E. rhopalostele per each host tree was randomly selected for root sampling and DNA extraction. • Key results: A total of 239 E. rhopalostele individuals were located in 25 of the 714 potential host trees. Light microscopy of sampled roots showed mycorrhizal fungi in 22 of the 25 sampled orchids. Phylogenetic analysis of ITS1-5.8S-ITS2 sequences yielded two Tulasnella clades. In four cases, plants were found to be associated with both clades. The difference between univariate and bivariate K functions was consistent with the random labeling null model at all spatial scales, indicating that trees hosting clades A and B of Tulasnella are not spatially segregated. The analysis of the inhomogenous K function showed that host trees are not clustered, suggesting no limitations to population-scale dispersal. χ2 analysis of contingency tables showed that E. rhopalostele is more frequent on dead trees than expected. • Conclusions: Epidendrum rhopalostele establishes mycorrhizal associations with at least two different Tulasnella species. The analysis of the distribution patterns of this orchid suggests a microsite preference for dead trees and no seed dispersal limitation.

Spatial distribution patterns of ascidians (Ascidiacea: Tunicata) in combination with information on bathymetry, oceanography, chlorophyll-a, and sea-ice on the continental shelves off the northern Antarctic Peninsula during POLARSTERN cruise ANT-XXIX/3

Ascidians (Ascidiacea: Tunicata) are sessile suspension feeders that represent dominant epifaunal components of the Southern Ocean shelf benthos and play a significant role in the pelagic-benthic coupling. Here, we report the results of a first study on the relationship between the distribution patterns of eight common and/or abundant (putative) ascidian species, and environmental drivers in the waters off the northern Antarctic Peninsula. During RV Polarstern cruise XXIX/3 (PS81) in January-March 2013, we used seabed imaging surveys along 28 photographic transects of 2 km length each at water depths from 70 to 770 m in three regions (northwestern Weddell Sea, southern Bransfield Strait and southern Drake Passage), differing in their general environmental setting, primarily oceanographic characteristics and sea-ice dynamics, to comparatively analyze the spatial patterns in the abundance of the selected ascidians, reliably to be identified in the photographs, at three nested spatial scales. At a regional (100-km) scale, the ascidian assemblages of the Weddell Sea differed significantly from those of the other two regions, whereas at an intermediate 10-km scale no such differences were detected among habitat types (bank, upper slope, slope, deep/canyon) on the shelf and at the shelf break within each region. These spatial patterns were superimposed by a marked small-scale (10-m) patchiness of ascidian distribution within the 2-km-long transects. Among the environmental variables considered in our study, a combination of water-mass characteristics, sea-ice dynamics (approximated by 5-year averages in sea-ice cover in the region of or surrounding the photographic stations), as well as the seabed ruggedness, was identified as explaining best the distribution patterns of the ascidians.

Linking riparian vegetation spatial structure in Australian tropical savannas to ecosystem health indicators: semi-variogram analysis of high spatial resolution satellite imagery

Government agencies responsible for riparian environments are assessing the combined utility of field survey and remote sensing for mapping and monitoring indicators of riparian zone health. The objective of this work was to determine if the structural attributes of savanna riparian zones in northern Australia can be detected from commercially available remotely sensed image data. Two QuickBird images and coincident field data covering sections of the Daly River and the South Alligator River - Barramundie Creek in the Northern Territory were used. Semi-variograms were calculated to determine the characteristic spatial scales of riparian zone features, both vegetative and landform. Interpretation of semi-variograms showed that structural dimensions of riparian environments could be detected and estimated from the QuickBird image data. The results also show that selecting the correct spatial resolution and spectral bands is essential to maximize the accuracy of mapping spatial characteristics of savanna riparian features. The distribution of foliage projective cover of riparian vegetation affected spectral reflectance variations in individual spectral bands differently. Pan-sharpened image data enabled small-scale information extraction (< 6 m) on riparian zone structural parameters. The semi-variogram analysis results provide the basis for an inversion approach using high spatial resolution satellite image data to map indicators of savanna riparian zone health.

Monitoring the condition of savanna riparian zones in Northern Australia using high spatial resolution imagery and a field based rapid assessment method

Government agencies responsible for riparian environments are assessing the combined utility of field survey and remote sensing for mapping and monitoring indicators of riparian zone condition. The objective of this work was to compare the Tropical Rapid Appraisal of Riparian Condition (TRARC) method to a satellite image based approach. TRARC was developed for rapid assessment of the environmental condition of savanna riparian zones. The comparison assessed mapping accuracy, representativeness of TRARC assessment, cost-effectiveness, and suitability for multi-temporal analysis. Two multi-spectral QuickBird images captured in 2004 and 2005 and coincident field data covering sections of the Daly River in the Northern Territory, Australia were used in this work. Both field and image data were processed to map riparian health indicators (RHIs) including percentage canopy cover, organic litter, canopy continuity, stream bank stability, and extent of tree clearing. Spectral vegetation indices, image segmentation and supervised classification were used to produce RHI maps. QuickBird image data were used to examine if the spatial distribution of TRARC transects provided a representative sample of ground based RHI measurements. Results showed that TRARC transects were required to cover at least 3% of the study area to obtain a representative sample. The mapping accuracy and costs of the image based approach were compared to those of the ground based TRARC approach. Results proved that TRARC was more cost-effective at smaller scales (1-100km), while image based assessment becomes more feasible at regional scales (100-1000km). Finally, the ability to use both the image and field based approaches for multi-temporal analysis of RHIs was assessed. Change detection analysis demonstrated that image data can provide detailed information on gradual change, while the TRARC method was only able to identify more gross scale changes. In conclusion, results from both methods were considered to complement each other if used at appropriate spatial scales.

Spatial and Temporal Variability of Dissolved Organic Matter Quantity and Composition in an Oilgotrophic Subtropical Coastal Wetland

Dissolved organic matter (DOM) is an essential component of the carbon cycle and a critical driver in controlling variety of biogeochemical and ecological processes in wetlands. The quality of this DOM as it relates to composition and reactivity is directly related to its sources and may vary on temporal and spatial scales. However, large scale, long-term studies of DOM dynamics in wetlands are still scarce in the literature. Here we present a multi-year DOM characterization study for monthly surface water samples collected at 14 sampling stations along two transects within the greater Everglades, a subtropical, oligotrophic, coastal freshwater wetland-mangrove-estuarine ecosystem. In an attempt to assess quantitative and qualitative variations of DOM on both spatial and temporal scales, we determined dissolved organic carbon (DOC) values and DOM optical properties, respectively. DOM quality was assessed using, excitation emission matrix (EEM) fluorescence coupled with parallel factor analysis (PARAFAC). Variations of the PARAFAC components abundance and composition were clearly observed on spatial and seasonal scales. Dry versus wet season DOC concentrations were affected by dry-down and re-wetting processes in the freshwater marshes, while DOM compositional features were controlled by soil and higher plant versus periphyton sources respectively. Peat-soil based freshwater marsh sites could be clearly differentiated from marl-soil based sites based on EEM–PARAFAC data. Freshwater marsh DOM was enriched in higher plant and soil-derived humic-like compounds, compared to estuarine sites which were more controlled by algae- and microbial-derived inputs. DOM from fringe mangrove sites could be differentiated between tidally influenced sites and sites exposed to long inundation periods. As such coastal estuarine sites were significantly controlled by hydrology, while DOM dynamics in Florida Bay were seasonally driven by both primary productivity and hydrology. This study exemplifies the application of long term optical properties monitoring as an effective technique to investigate DOM dynamics in aquatic ecosystems. The work presented here also serves as a pre-restoration condition dataset for DOM in the context of the Comprehensive Everglades Restoration Plan (CERP).

Coral reef habitat maps derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Coral reef maps at various spatial scales and extents are needed for mapping, monitoring, modelling, and management of these environments. High spatial resolution satellite imagery, pixel <10 m, integrated with field survey data and processed with various mapping approaches, can provide these maps. These approaches have been accurately applied to single reefs (10-100 km**2), covering one high spatial resolution scene from which a single thematic layer (e.g. benthic community) is mapped. This article demonstrates how a hierarchical mapping approach can be applied to coral reefs from individual reef to reef-system scales (10-1000 km**2) using object-based image classification of high spatial resolution images guided by ecological and geomorphological principles. The approach is demonstrated for three individual reefs (10-35 km**2) in Australia, Fiji, and Palau; and for three complex reef systems (300-600 km**2) one in the Solomon Islands and two in Fiji. Archived high spatial resolution images were pre-processed and mosaics were created for the reef systems. Georeferenced benthic photo transect surveys were used to acquire cover information. Field and image data were integrated using an object-based image analysis approach that resulted in a hierarchically structured classification. Objects were assigned class labels based on the dominant benthic cover type, or location-relevant ecological and geomorphological principles, or a combination thereof. This generated a hierarchical sequence of reef maps with an increasing complexity in benthic thematic information that included: 'reef', 'reef type', 'geomorphic zone', and 'benthic community'. The overall accuracy of the 'geomorphic zone' classification for each of the six study sites was 76-82% using 6-10 mapping categories. For 'benthic community' classification, the overall accuracy was 52-75% with individual reefs having 14-17 categories and reef systems 20-30 categories. We show that an object-based classification of high spatial resolution imagery, guided by field data and ecological and geomorphological principles, can produce consistent, accurate benthic maps at four hierarchical spatial scales for coral reefs of various sizes and complexities.

Spatial Relationships among Hydroacoustic, Hydrographic and Top Predator Patterns: Cetacean Distributions in the Mid-Atlantic Bight

Effective conservation and management of top predators requires a comprehensive understanding of their distributions and of the underlying biological and physical processes that affect these distributions. The Mid-Atlantic Bight shelf break system is a dynamic and productive region where at least 32 species of cetaceans have been recorded through various systematic and opportunistic marine mammal surveys from the 1970s through 2012. My dissertation characterizes the spatial distribution and habitat of cetaceans in the Mid-Atlantic Bight shelf break system by utilizing marine mammal line-transect survey data, synoptic multi-frequency active acoustic data, and fine-scale hydrographic data collected during the 2011 summer Atlantic Marine Assessment Program for Protected Species (AMAPPS) survey. Although studies describing cetacean habitat and distributions have been previously conducted in the Mid-Atlantic Bight, my research specifically focuses on the shelf break region to elucidate both the physical and biological processes that influence cetacean distribution patterns within this cetacean hotspot.

In Chapter One I review biologically important areas for cetaceans in the Atlantic waters of the United States. I describe the study area, the shelf break region of the Mid-Atlantic Bight, in terms of the general oceanography, productivity and biodiversity. According to recent habitat-based cetacean density models, the shelf break region is an area of high cetacean abundance and density, yet little research is directed at understanding the mechanisms that establish this region as a cetacean hotspot.

In Chapter Two I present the basic physical principles of sound in water and describe the methodology used to categorize opportunistically collected multi-frequency active acoustic data using frequency responses techniques. Frequency response classification methods are usually employed in conjunction with net-tow data, but the logistics of the 2011 AMAPPS survey did not allow for appropriate net-tow data to be collected. Biologically meaningful information can be extracted from acoustic scattering regions by comparing the frequency response curves of acoustic regions to theoretical curves of known scattering models. Using the five frequencies on the EK60 system (18, 38, 70, 120, and 200 kHz), three categories of scatterers were defined: fish-like (with swim bladder), nekton-like (e.g., euphausiids), and plankton-like (e.g., copepods). I also employed a multi-frequency acoustic categorization method using three frequencies (18, 38, and 120 kHz) that has been used in the Gulf of Maine and Georges Bank which is based the presence or absence of volume backscatter above a threshold. This method is more objective than the comparison of frequency response curves because it uses an established backscatter value for the threshold. By removing all data below the threshold, only strong scattering information is retained.

In Chapter Three I analyze the distribution of the categorized acoustic regions of interest during the daytime cross shelf transects. Over all transects, plankton-like acoustic regions of interest were detected most frequently, followed by fish-like acoustic regions and then nekton-like acoustic regions. Plankton-like detections were the only significantly different acoustic detections per kilometer, although nekton-like detections were only slightly not significant. Using the threshold categorization method by Jech and Michaels (2006) provides a more conservative and discrete detection of acoustic scatterers and allows me to retrieve backscatter values along transects in areas that have been categorized. This provides continuous data values that can be integrated at discrete spatial increments for wavelet analysis. Wavelet analysis indicates significant spatial scales of interest for fish-like and nekton-like acoustic backscatter range from one to four kilometers and vary among transects.

In Chapter Four I analyze the fine scale distribution of cetaceans in the shelf break system of the Mid-Atlantic Bight using corrected sightings per trackline region, classification trees, multidimensional scaling, and random forest analysis. I describe habitat for common dolphins, Risso’s dolphins and sperm whales. From the distribution of cetacean sightings, patterns of habitat start to emerge: within the shelf break region of the Mid-Atlantic Bight, common dolphins were sighted more prevalently over the shelf while sperm whales were more frequently found in the deep waters offshore and Risso’s dolphins were most prevalent at the shelf break. Multidimensional scaling presents clear environmental separation among common dolphins and Risso’s dolphins and sperm whales. The sperm whale random forest habitat model had the lowest misclassification error (0.30) and the Risso’s dolphin random forest habitat model had the greatest misclassification error (0.37). Shallow water depth (less than 148 meters) was the primary variable selected in the classification model for common dolphin habitat. Distance to surface density fronts and surface temperature fronts were the primary variables selected in the classification models to describe Risso’s dolphin habitat and sperm whale habitat respectively. When mapped back into geographic space, these three cetacean species occupy different fine-scale habitats within the dynamic Mid-Atlantic Bight shelf break system.

In Chapter Five I present a summary of the previous chapters and present potential analytical steps to address ecological questions pertaining the dynamic shelf break region. Taken together, the results of my dissertation demonstrate the use of opportunistically collected data in ecosystem studies; emphasize the need to incorporate middle trophic level data and oceanographic features into cetacean habitat models; and emphasize the importance of developing more mechanistic understanding of dynamic ecosystems.

Integrated Optimal Power Flow for Distribution Networks in Local and Urban Scales

This paper develops an integrated optimal power flow (OPF) tool for distribution networks in two spatial scales. In the local scale, the distribution network, the natural gas network, and the heat system are coordinated as a microgrid. In the urban scale, the impact of natural gas network is considered as constraints for the distribution network operation. The proposed approach incorporates unbalance three-phase electrical systems, natural gas systems, and combined cooling, heating, and power systems. The interactions among the above three energy systems are described by energy hub model combined with components capacity constraints. In order to efficiently accommodate the nonlinear constraint optimization problem, particle swarm optimization algorithm is employed to set the control variables in the OPF problem. Numerical studies indicate that by using the OPF method, the distribution network can be economically operated. Also, the tie-line power can be effectively managed.

Seep communities from two mud volcanoes in the deep eastern Mediterranean Sea: faunal composition, spatial patterns and environmental control

The Mediterranean Sea constitutes a unique environment to study cold-seep ecosystems due to the presence of different geodynamic settings, from an active margin along the Mediterranean Ridge (MR) to a passive margin in the Nile Deep-Sea Fan (NDSF). We attempted to identify the structure of benthic communities associated with the Napoli and Amsterdam mud volcanoes (MVs) located on the MR and to establish the links between faunal distribution and environmental conditions at different spatial scales. Comparison between the 2 MVs revealed that the faunal distribution seemed to be mainly controlled by the characteristics of the microhabitats. On both geological structures, the variability between the different microhabitats was higher than the variability observed between replicates of the same microhabitat, and the distribution of macro-fauna was apparently linked to gradients in physico-chemical conditions. The peripheral sites from Napoli were generally more oxygenated and harboured lower species richness than the active sites. The reduced sediment microhabitat from Amsterdam presented the highest methane concentrations and was mainly colonised by symbiont-bearing vesicomyid bivalves and heterotrophic dorvilleid polychaetes. Overall, a higher taxonomic diversity was observed on Napoli. Sub-stratum type was hypothesised to be the second factor influencing faunal distribution. The results of this study highlight the high heterogeneity of faunal communities associated with seep ecosystems within this region and the need to pursue investigations at various spatial and temporal scales.

Environmental and grazing influence on spatial variability of intertidal biofilm on subtropical rocky shores

Epilithic biofilm on rocky shores is regulated by physico-chemical and biological factors and is important as a source of food for benthic organisms. The influences of environmental and grazing pressure on spatial variability of biomass of biofilm were evaluated on shores on the north coast of São Paulo State (SE Brazil). A general trend of greater abundance of microalgae was observed lower on the shore, but neither of the environmental factors evaluated (wave exposure and shore level) showed consistent effects, and differences were found among specific shores or times (September 2007 and March 2008). The abundance of slow-moving grazers (limpets and littorinids) showed a negative correlation with chlorophyll a concentration on shores. However, experimental exclusion of these grazers failed to show consistent results at small spatial scales. Observations of divergent abundances of the isopod Ligia exotica and biomass of biofilm on isolated boulders on shores led to a short exclusion experiment, where the grazing pressure by L. exotica significantly decreased microalgal biomass. The result suggests that grazing activities of this fast-moving consumer probably mask the influence of slow-moving grazers at small spatial scales, while both have an additive effect at larger scales that masks environmental influences. This is the first evaluation of the impact of the fast-moving herbivore L. exotica on microalgal biomass on rocky shores and opens an interesting discussion about the role of these organisms in subtropical coastal environments.

Spatial patterns of genetic diversity in Mediterranean eagle owl Bubo bubo populations

Little information is available on the patterns of genetic connectivity in owls. We studied the genetic structure of the eagle owl Bubo bubo (Linnaeus, 1758) in southeastern Spain at two different spatial scales. Seven microsatellites previously described for this species were used, although only six loci amplified correctly. The observed low genetic variation could be explained by the short dispersal distance, high mortality rate and high degree of monogamy shown by this large nocturnal predator. As expected, the highest genetic isolation was detected in the geographically most isolated population. Significant genetic differentiation was found among study units separated by less than 50 km. The territorial analysis showed interesting connectivity patterns related with the gene flow and turnover rate of the breeding individuals. The lowest genetic diversity was found in the region with the largest population, which could imply incipient inbreeding.

Genetic diversity across geographical scales in marine coastal ecosystems: Holothuria arguinensis a model species

Coastal lagoons are considered one of the most productive areas of our planet harboring a large variety of habitats. Their transitional character, between terrestrial and marine environments, creates a very particular ecosystem with important variations of its environmental conditions. The organisms that are able to survive on these ecosystems frequently experience strong selective pressures and constrictions to gene flowwith marine populations, which could contribute to genetic divergence among populations inhabiting coastal lagoon and marine environments. Therefore, the main aims of this study are to asses the genetic diversity and population structure of Holothuria arguinensis across geographical ranges, to test the hypothesis of coastal lagoons as hotspots of genetic diversity in the Ria Formosa lagoon, and to determine the role of exporting standing genetic variation from the lagoon to open sea and their implications to recent geographical expansion events. To reach these objectives, we investigate the genetic structure of H. arguinensis using two mitochondrial DNA markers (COI and 16S) at different spatial scales: i) small, inside Ria Formosa coastal lagoon, South Portugal; 2) large, including most of the geographical distribution of this species (South and Western Portuguese coast and Canary islands); these results will allow us to compare the genetic diversity of lagoonal and marine populations of H. arguinensis. On this framework, its recent geographical expansion events, recorded by Rodrigues (2012) and González-Wangüemert and Borrero-Pérez (2012), will be analyzed considering the potential contribution from lagoonal genetic pool. Non-significant genetic structure and high haplotypic diversity were found inside the Ria Formosa coastal lagoon. Both genes were unable to detect significant genetic differentiation among lagoonal and marine localities, suggesting a high rate of gene flow. The results supported our hypotheses that coastal lagoons are not only acting as hotspots of genetic diversity, but also contributing for the genetic variability of the species, working as a source of new haplotypes and enhancing adaptation to the high variable conditions. Different genetic patterns of colonization were found on H. arguinensis, but they must be studied more deeply.

Analysis of spatial variation of tree heights in an urban area using LiDAR data.

Estimating with greater precision and accuracy the height of plants has been a challenge for the scientific community. The objective this study is to evaluate the spatial variation of tree heights at different spatial scales in areas of the city of Recife, Brazil, using LiDAR remote sensing data. The LiDAR data were processed in the QT Modeler (Quick Terrain Modeler v. 8.0.2) software from Applied Imagery. The TreeVaW software was utilized to estimate the heights and crown diameters of trees. The results obtained for tree height were consistent with field measurements.

Sexual selection in a tropical fruit fly : role of a plant derived chemical in mate choice

The specific mechanisms by which selective pressures affect individuals are often difficult to resolve. In tephritid fruit flies, males respond strongly and positively to certain plant derived chemicals. Sexual selection by female choice has been hypothesized as the mechanism driving this behaviour in certain species, as females preferentially mate with males that have fed on these chemicals. This hypothesis is, to date, based on studies of only very few species and its generality is largely untested. We tested the hypothesis on different spatial scales (small cage and seminatural field-cage) using the monophagous fruit fly, Bactrocera cacuminata. This species is known to respond to methyl eugenol (ME), a chemical found in many plant species and one upon which previous studies have focused. Contrary to expectation, no obvious female choice was apparent in selecting ME-fed males over unfed males as measured by the number of matings achieved over time, copulation duration, or time of copulation initiation. However, the number of matings achieved by ME-fed males was significantly greater than unfed males 16 and 32 days after exposure to ME in small cages (but not in a field-cage). This delayed advantage suggests that ME may not influence the pheromone system of B. cacuminata but may have other consequences, acting on some other fitness consequence (e.g., enhancement of physiology or survival) of male exposure to these chemicals. We discuss the ecological and evolutionary implications of our findings to explore alternate hypotheses to explain the patterns of response of dacine fruit flies to specific plant-derived chemicals.

Socio-environmental factors and suicide in Queensland, Australia

Suicide has drawn much attention from both the scientific community and the public. Examining the impact of socio-environmental factors on suicide is essential in developing suicide prevention strategies and interventions, because it will provide health authorities with important information for their decision-making. However, previous studies did not examine the impact of socio-environmental factors on suicide using a spatial analysis approach. The purpose of this study was to identify the patterns of suicide and to examine how socio-environmental factors impact on suicide over time and space at the Local Governmental Area (LGA) level in Queensland. The suicide data between 1999 and 2003 were collected from the Australian Bureau of Statistics (ABS). Socio-environmental variables at the LGA level included climate (rainfall, maximum and minimum temperature), Socioeconomic Indexes for Areas (SEIFA) and demographic variables (proportion of Indigenous population, unemployment rate, proportion of population with low income and low education level). Climate data were obtained from Australian Bureau of Meteorology. SEIFA and demographic variables were acquired from ABS. A series of statistical and geographical information system (GIS) approaches were applied in the analysis. This study included two stages. The first stage used average annual data to view the spatial pattern of suicide and to examine the association between socio-environmental factors and suicide over space. The second stage examined the spatiotemporal pattern of suicide and assessed the socio-environmental determinants of suicide, using more detailed seasonal data. In this research, 2,445 suicide cases were included, with 1,957 males (80.0%) and 488 females (20.0%). In the first stage, we examined the spatial pattern and the determinants of suicide using 5-year aggregated data. Spearman correlations were used to assess associations between variables. Then a Poisson regression model was applied in the multivariable analysis, as the occurrence of suicide is a small probability event and this model fitted the data quite well. Suicide mortality varied across LGAs and was associated with a range of socio-environmental factors. The multivariable analysis showed that maximum temperature was significantly and positively associated with male suicide (relative risk [RR] = 1.03, 95% CI: 1.00 to 1.07). Higher proportion of Indigenous population was accompanied with more suicide in male population (male: RR = 1.02, 95% CI: 1.01 to 1.03). There was a positive association between unemployment rate and suicide in both genders (male: RR = 1.04, 95% CI: 1.02 to 1.06; female: RR = 1.07, 95% CI: 1.00 to 1.16). No significant association was observed for rainfall, minimum temperature, SEIFA, proportion of population with low individual income and low educational attainment. In the second stage of this study, we undertook a preliminary spatiotemporal analysis of suicide using seasonal data. Firstly, we assessed the interrelations between variables. Secondly, a generalised estimating equations (GEE) model was used to examine the socio-environmental impact on suicide over time and space, as this model is well suited to analyze repeated longitudinal data (e.g., seasonal suicide mortality in a certain LGA) and it fitted the data better than other models (e.g., Poisson model). The suicide pattern varied with season and LGA. The north of Queensland had the highest suicide mortality rate in all the seasons, while there was no suicide case occurred in the southwest. Northwest had consistently higher suicide mortality in spring, autumn and winter. In other areas, suicide mortality varied between seasons. This analysis showed that maximum temperature was positively associated with suicide among male population (RR = 1.24, 95% CI: 1.04 to 1.47) and total population (RR = 1.15, 95% CI: 1.00 to 1.32). Higher proportion of Indigenous population was accompanied with more suicide among total population (RR = 1.16, 95% CI: 1.13 to 1.19) and by gender (male: RR = 1.07, 95% CI: 1.01 to 1.13; female: RR = 1.23, 95% CI: 1.03 to 1.48). Unemployment rate was positively associated with total (RR = 1.40, 95% CI: 1.24 to 1.59) and female (RR=1.09, 95% CI: 1.01 to 1.18) suicide. There was also a positive association between proportion of population with low individual income and suicide in total (RR = 1.28, 95% CI: 1.10 to 1.48) and male (RR = 1.45, 95% CI: 1.23 to 1.72) population. Rainfall was only positively associated with suicide in total population (RR = 1.11, 95% CI: 1.04 to 1.19). There was no significant association for rainfall, minimum temperature, SEIFA, proportion of population with low educational attainment. The second stage is the extension of the first stage. Different spatial scales of dataset were used between the two stages (i.e., mean yearly data in the first stage, and seasonal data in the second stage), but the results are generally consistent with each other. Compared with other studies, this research explored the variety of the impact of a wide range of socio-environmental factors on suicide in different geographical units. Maximum temperature, proportion of Indigenous population, unemployment rate and proportion of population with low individual income were among the major determinants of suicide in Queensland. However, the influence from other factors (e.g. socio-culture background, alcohol and drug use) influencing suicide cannot be ignored. An in-depth understanding of these factors is vital in planning and implementing suicide prevention strategies. Five recommendations for future research are derived from this study: (1) It is vital to acquire detailed personal information on each suicide case and relevant information among the population in assessing the key socio-environmental determinants of suicide; (2) Bayesian model could be applied to compare mortality rates and their socio-environmental determinants across LGAs in future research; (3) In the LGAs with warm weather, high proportion of Indigenous population and/or unemployment rate, concerted efforts need to be made to control and prevent suicide and other mental health problems; (4) The current surveillance, forecasting and early warning system needs to be strengthened, to trace the climate and socioeconomic change over time and space and its impact on population health; (5) It is necessary to evaluate and improve the facilities of mental health care, psychological consultation, suicide prevention and control programs; especially in the areas with low socio-economic status, high unemployment rate, extreme weather events and natural disasters.