Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Advanced multivariate analysis to assess remediation of hydrocarbons in soils

Accurate monitoring of degradation levels in soils is essential in order to understand and achieve complete degradation of petroleum hydrocarbons in contaminated soils. We aimed to develop the use of multivariate methods for the monitoring of biodegradation of diesel in soils and to determine if diesel contaminated soils could be remediated to a chemical composition similar to that of an uncontaminated soil. An incubation experiment was set up with three contrasting soil types. Each soil was exposed to diesel at varying stages of degradation and then analysed for key hydrocarbons throughout 161 days of incubation. Hydrocarbon distributions were analysed by Principal Coordinate Analysis and similar samples grouped by cluster analysis. Variation and differences between samples were determined using permutational multivariate analysis of variance. It was found that all soils followed trajectories approaching the chemical composition of the unpolluted soil. Some contaminated soils were no longer significantly different to that of uncontaminated soil after 161 days of incubation. The use of cluster analysis allows the assignment of a percentage chemical similarity of a diesel contaminated soil to an uncontaminated soil sample. This will aid in the monitoring of hydrocarbon contaminated sites and the establishment of potential endpoints for successful remediation.

Multivariate calibrations in gas chromatographic profiles for prediction of several physicochemical parameters of Brazilian commercial gasoline

Several Brazilian commercial gasoline physicochemical parameters, such as relative density, distillation curve (temperatures related to 10%, 50% and 90% of distilled volume, final boiling point and residue), octane numbers (motor and research octane number and anti-knock index), hydrocarbon compositions (olefins, aromatics and saturates) and anhydrous ethanol and benzene content was predicted from chromatographic profiles obtained by flame ionization detection (GC-FID) and using partial least square regression (PLS). GC-FID is a technique intensively used for fuel quality control due to its convenience, speed, accuracy and simplicity and its profiles are much easier to interpret and understand than results produced by other techniques. Another advantage is that it permits association with multivariate methods of analysis, such as PLS. The chromatogram profiles were recorded and used to deploy PLS models for each property. The standard error of prediction (SEP) has been the main parameter considered to select the "best model". Most of GC-FID-PLS results, when compared to those obtained by the Brazilian Government Petroleum, Natural Gas and Biofuels Agency - ANP Regulation 309 specification methods, were very good. In general, all PLS models developed in these work provide unbiased predictions with lows standard error of prediction and percentage average relative error (below 11.5 and 5.0, respectively). (C) 2007 Elsevier B.V. All rights reserved.

Multivariate Techniques in Biomechanical Analysis

This work outlines the theoretical advantages of multivariate methods in biomechanical data, validates the proposed methods and outlines new clinical findings relating to knee osteoarthritis that were made possible by this approach. New techniques were based on existing multivariate approaches, Partial Least Squares (PLS) and Non-negative Matrix Factorization (NMF) and validated using existing data sets. The new techniques developed, PCA-PLS-LDA (Principal Component Analysis – Partial Least Squares – Linear Discriminant Analysis), PCA-PLS-MLR (Principal Component Analysis – Partial Least Squares –Multiple Linear Regression) and Waveform Similarity (based on NMF) were developed to address the challenging characteristics of biomechanical data, variability and correlation. As a result, these new structure-seeking technique revealed new clinical findings. The first new clinical finding relates to the relationship between pain, radiographic severity and mechanics. Simultaneous analysis of pain and radiographic severity outcomes, a first in biomechanics, revealed that the knee adduction moment’s relationship to radiographic features is mediated by pain in subjects with moderate osteoarthritis. The second clinical finding was quantifying the importance of neuromuscular patterns in brace effectiveness for patients with knee osteoarthritis. I found that brace effectiveness was more related to the patient’s unbraced neuromuscular patterns than it was to mechanics, and that these neuromuscular patterns were more complicated than simply increased overall muscle activity, as previously thought.

Purification of pharmaceuticals and nutraceutical compounds by sub- and supercritical extraction and chromatography

This thesis discusses the use of sub- and supercritical fluids as the medium in extraction and chromatography. Super- and subcritical extraction was used to separate essential oils from herbal plant Angelica archangelica. The effect of extraction parameters was studied and sensory analyses of the extracts were done by an expert panel. The results of the sensory analyses were compared to the analytically determined contents of the extracts. Sub- and supercritical fluid chromatography (SFC) was used to separate and purify high-value pharmaceuticals. Chiral SFC was used to separate the enantiomers of racemic mixtures of pharmaceutical compounds. Very low (cryogenic) temperatures were applied to substantially enhance the separation efficiency of chiral SFC. The thermodynamic aspects affecting the resolving ability of chiral stationary phases are briefly reviewed. The process production rate which is a key factor in industrial chromatography was optimized by empirical multivariate methods. General linear model was used to optimize the separation of omega-3 fatty acid ethyl esters from esterized fish oil by using reversed-phase SFC. Chiral separation of racemic mixtures of guaifenesin and ferulic acid dimer ethyl ester was optimized by using response surface method with three variables per time. It was found that by optimizing four variables (temperature, load, flowate and modifier content) the production rate of the chiral resolution of racemic guaifenesin by cryogenic SFC could be increased severalfold compared to published results of similar application. A novel pressure-compensated design of industrial high pressure chromatographic column was introduced, using the technology developed in building the deep-sea submersibles (Mir 1 and 2). A demonstration SFC plant was built and the immunosuppressant drug cyclosporine A was purified to meet the requirements of US Pharmacopoeia. A smaller semi-pilot size column with similar design was used for cryogenic chiral separation of aromatase inhibitor Finrozole for use in its development phase 2.

Psychological resources, their social antecedents, and association with well-being and health behaviour in early adulthood

The success of entering work life, young people s psychological resources and self-reported well-being were studied in a longitudinal setting from a life-span developmental-contextual perspective in early adulthood. The aim was to analyse how psychosocial characteristics in early childhood and adolescence predict successful entrance into work life, how this is associated with well-being, and to assess the level of psychological resources such as dispositional optimism, personal meaning of work and coping in early adulthood. The role of these and social support, in the relationship between regional factors (such as place of residence and migration), self-reported health and life satisfaction was studied. The association between a specific coping strategy, i.e. eating and drinking in a stressful situation and eating habits, was studied to demonstrate how coping is associated with health behaviour. Multivariate methods, including binary logistic regression analyses and ANOVA, were used for statistical analyses. The subjects were members of the Northern Finland 1966 Birth Cohort, which consists of all women and men born in 1966 in the two northernmost provinces of Finland (n= 12,058). The most recent follow-up, at the age of 31 years when 11,637 subjects were alive, took place in 1997-1998. The results show, first, that social resources in the childhood family and adolescence school achievement predict entrance into the labour market. Secondly, psychosocial resources were found to mediate the relationship between migration from rural to urban areas, and subjective well-being. Thirdly, psychological resources at entrance into the labour market were found to develop from early infancy on. They are, however, influenced later by work history. Fourthly, stress-related eating and drinking, as a way of coping, was found to be directly associated with unhealthy eating habits and alcohol use. Gender differences were found in psychosocial resources predicting, and being associated with success in entering the labour market. For men, the role of attitudinal and psychological factors seems to be especially important in entrance into work life and in the development of psychological resources. For women, academic attainment was more important for successfully entering work life, and lack of emotional social support was a risk factor for stress-related eating only among women. Stress-related eating and drinking habits were predicted by a long history of unemployment as well as a low level of education among both genders, but not excluding an academic degree among men. The results emphasize the role of childhood psychosocial factors in preventing long-term unemployment and in enhancing psychological well-being in early adulthood. Success in entering work life, in terms of continuous work history, plays a crucial role for well-being and the amount of psychological resources in early adulthood. The results emphasize the crucial role of enhancing psychological resources for promoting positive health behaviour and diminishing regional differences in subjective well-being.

Estimating crop area using seasonal time series of enhanced vegetation index from MODIS satellite imagery

Cereal grain is one of the main export commodities of Australian agriculture. Over the past decade, crop yield forecasts for wheat and sorghum have shown appreciable utility for industry planning at shire, state, and national scales. There is now an increasing drive from industry for more accurate and cost-effective crop production forecasts. In order to generate production estimates, accurate crop area estimates are needed by the end of the cropping season. Multivariate methods for analysing remotely sensed Enhanced Vegetation Index (EVI) from 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery within the cropping period (i.e. April-November) were investigated to estimate crop area for wheat, barley, chickpea, and total winter cropped area for a case study region in NE Australia. Each pixel classification method was trained on ground truth data collected from the study region. Three approaches to pixel classification were examined: (i) cluster analysis of trajectories of EVI values from consecutive multi-date imagery during the crop growth period; (ii) harmonic analysis of the time series (HANTS) of the EVI values; and (iii) principal component analysis (PCA) of the time series of EVI values. Images classified using these three approaches were compared with each other, and with a classification based on the single MODIS image taken at peak EVI. Imagery for the 2003 and 2004 seasons was used to assess the ability of the methods to determine wheat, barley, chickpea, and total cropped area estimates. The accuracy at pixel scale was determined by the percent correct classification metric by contrasting all pixel scale samples with independent pixel observations. At a shire level, aggregated total crop area estimates were compared with surveyed estimates. All multi-temporal methods showed significant overall capability to estimate total winter crop area. There was high accuracy at pixel scale (>98% correct classification) for identifying overall winter cropping. However, discrimination among crops was less accurate. Although the use of single-date EVI data produced high accuracy for estimates of wheat area at shire scale, the result contradicted the poor pixel-scale accuracy associated with this approach, due to fortuitous compensating errors. Further studies are needed to extrapolate the multi-temporal approaches to other geographical areas and to improve the lead time for deriving cropped-area estimates before harvest.

The role of fruit traits of bird-dispersed plants in invasiveness and weed risk assessment

Aim: Birds play a major role in the dispersal of seeds of many fleshy-fruited invasive plants. The fruits that birds choose to consume are influenced by fruit traits. However, little is known of how the traits of invasive plant fruits contribute to invasiveness or to their use by frugivores. We aim to gain a greater understanding of these relationships to improve invasive plant management. Location: South-east Queensland, Australia. Methods: We measure a variety of fruit morphology, pulp nutrient and phenology traits of a suite of bird-dispersed alien plants. Frugivore richness of these aliens was derived from the literature. Using regressions and multivariate methods, we investigate relationships between fruit traits, frugivore richness and invasiveness. Results: Plant invasiveness was negatively correlated to fruit size, and all highly invasive species had quite similar fruit morphology [smaller fruits, seeds of intermediate size and few (<10) seeds per fruit]. Lower pulp water was the only pulp nutrient trait associated with invasiveness. There were strong positive relationships between the diversity of bird frugivores and plant invasiveness, and in the diversity of bird frugivores in the study region and another part of the plants' alien range. Main conclusions: Our results suggest that weed risk assessments (WRA) and predictions of invasive success for bird-dispersed plants can be improved. Scoring criteria for WRA regarding fruit size would need to be system-specific, depending on the fruit-processing capabilities of local frugivores. Frugivore richness could be quantified in the plant's natural range, its invasive range elsewhere, or predictions made based on functionally similar fruits.

Distribution of non-biting midges (Diptera, Chironomidae) in subarctic lakes of Finnish Lapland applications in lake classification and palaeolimnology

Climate change contributes directly or indirectly to changes in species distributions, and there is very high confidence that recent climate warming is already affecting ecosystems. The Arctic has already experienced the greatest regional warming in recent decades, and the trend is continuing. However, studies on the northern ecosystems are scarce compared to more southerly regions. Better understanding of the past and present environmental change is needed to be able to forecast the future. Multivariate methods were used to explore the distributional patterns of chironomids in 50 shallow (≤ 10m) lakes in relation to 24 variables determined in northern Fennoscandia at the ecotonal area from the boreal forest in the south to the orohemiarctic zone in the north. Highest taxon richness was noted at middle elevations around 400 m a.s.l. Significantly lower values were observed from cold lakes situated in the tundra zone. Lake water alkalinity had the strongest positive correlation with the taxon richness. Many taxa had preference for lakes either on tundra area or forested area. The variation in the chironomid abundance data was best correlated with sediment organic content (LOI), lake water total organic carbon content, pH and air temperature, with LOI being the strongest variable. Three major lake groups were separated on the basis of their chironomid assemblages: (i) small and shallow organic-rich lakes, (ii) large and base-rich lakes, and (iii) cold and clear oligotrophic tundra lakes. Environmental variables best discriminating the lake groups were LOI, taxon richness, and Mg. When repeated, this kind of an approach could be useful and efficient in monitoring the effects of global change on species ranges. Many species of fast spreading insects, including chironomids, show a remarkable ability to track environmental changes. Based on this ability, past environmental conditions have been reconstructed using their chitinous remains in the lake sediment profiles. In order to study the Holocene environmental history of subarctic aquatic systems, and quantitatively reconstruct the past temperatures at or near the treeline, long sediment cores covering the last 10000 years (the Holocene) were collected from three lakes. Lower temperature values than expected based on the presence of pine in the catchment during the mid-Holocene were reconstructed from a lake with great water volume and depth. The lake provided thermal refuge for profundal, cold adapted taxa during the warm period. In a shallow lake, the decrease in the reconstructed temperatures during the late Holocene may reflect the indirect response of the midges to climate change through, e.g., pH change. The results from three lakes indicated that the response of chironomids to climate have been more or less indirect. However, concurrent shifts in assemblages of chironomids and vegetation in two lakes during the Holocene time period indicated that the midges together with the terrestrial vegetation had responded to the same ultimate cause, which most likely was the Holocene climate change. This was also supported by the similarity in the long-term trends in faunal succession for the chironomid assemblages in several lakes in the area. In northern Finnish Lapland the distribution of chironomids were significantly correlated with physical and limnological factors that are most likely to change as a result of future climate change. The indirect and individualistic response of aquatic systems, as reconstructed using the chironomid assemblages, to the climate change in the past suggests that in the future, the lake ecosystems in the north do not respond in one predictable way to the global climate change. Lakes in the north may respond to global climate change in various ways that are dependent on the initial characters of the catchment area and the lake.

Suokasvillisuuden vaste muuttuvaan ilmastoon: kokeellinen tutkimus kuivumisen vaikutuksesta

Soilla on merkittävä rooli ilmastonmuutoksen hillitsemisessä suuren hiilivarastonsa sekä ekosysteemin ja ilmakehän välisen kaasunvaihdon ansiosta. Ilmastonmuutoksen ennustetaan vaikuttavan suokasvillisuuteen ja suon toimintaan epäsuorasti. Vedenpinnan ennustetaan laskevan 14–21 cm johtuen kasveista ja avoimilta pinnoilta tapahtuvan haihdunnan lisääntymisestä lämpötilan noustessa, mikäli sadanta ei lisäänny. Aiemmat vedenpinnan laskun jälkeistä kasvillisuutta seuranneet tutkimukset ovat osoittaneet, että putkilokasvit hyötyvät alhaisemmasta vedenpinnan tasosta ja että kuljuun sopeutuneet rahkasammalet kärsivät kuivuneista oloista. Kasvillisuuden runsaussuhteiden muuttumisen lisäksi kasviyhteisöjen monimuotoisuus vähenee. Erityisen herkkiä vedenpinnan laskulle ovat olleet välipinta- ja kuljurahkasammalet ja sarat. Funktionaalisten kasviryhmien vasteiden selvittämiksesi käytettiin BACI (before-after-control-impact) –tutkimusotetta. Tutkimuksessa oli kolme verrokkialaa ja kolme käsittelyalaa, joissa vedenpintaa oli laskettu 14–21 senttimetriin. Lisäksi vertailukohdaksi tutkimuksessa oli mukana kolme alaa, joissa oli tehty metsäojitus n. 50 vuotta sitten. Nämä toistot sijaitsivat meso-, oligo ombrotrofisilla suotyypeillä Oriveden Lakkasuolla. Kasvillisuus kartoitettiin ja vedenpinnat mitattiin aloilta ennen käsittelyä vuonna 2000 sekä vuosina 2001–2003 ja 2009. Aineisto analysoitiin TWINSPAN- (PC-Ord), PRC ja DCA (CANOCO)-monimuuttujamenetelmillä. Tulokset osoittivat, että verrokki- ja käsittelyalat olivat samanlaisia lähtökohdiltaan, joten niitä voitiin käsittelyn jälkeen verrata toisiinsa. Kasvillisuuden rakenne vaihteli vuosien välillä myös verrokkialoilla, mikä osoittaa kasvien sopeutumiskyvyn muuttuviin sääoloihin (lämpötila, sademäärä). Vuosi 2003 erottui tutkimuksessa alhaisella vedenpinnantasolla, mutta toisaalta myös ainavihantien varpujen suuren peittävyyden osalta. Vuoteen 2009 mennessä kasvillisuuden erityisesti sarojen peittävyys väheni. Ravinteikkaimmilla toistoilla kasvillisuuden vasteet vaikuttivat olevan vahvemmat kuin vähäravinteisilla toistoilla. Kasviryhmistä kulju- ja välipintasammalilla oli vahvimmat vasteetvedenpinnan laskuun ja mätäslajeilla heikoimmat. Tulosten mukaan kasviryhmien vasteet vaihtelevat riippuen tarkasteltavasta aikajaksosta: ensimmäiset kolme vuotta käsittelyn jälkeen suo oli häiriötilassa ja vasta sen jälkeen kasvillisuus sopeutui muuttuneisiin oloihin.

Assessing the sustainability of brackish water aquaculture systems in the Phillipines. Working Paper 1: A typology of brackish-water pond aquaculture systems in the Philippines

Aquaculture in the Philippines is a long-standing activity but has witnessed relatively recent, rapid, technical change with the introduction of hatchery technology and commercial feed-mills changing the production possibilities for a fishpond operator. We are confronted with a diversity of aquaculture practices in the coastal areas of the Philippines, with new technologies being incorporated into more traditional systems. As a first step to understanding the sector, we therefore present a typology of farming systems with the motivation of generating domains (farm “types”) over which we can compare performance on a number of indicators. Our typology, restricted to brackish-water pond systems, is constructed using multivariate methods (principal components analysis, cluster analysis). Eight variables are used relating to the management of the farm across all the major factors of production. A stratified net sample of 136 observations provides the data for the analysis, from a farm-level survey carried out between January and June 2003 in the two main brackish-water production regions in the Philippines. We define five farm types from this analysis. In later work we will show how the use of this typology can be used for comparative study of economic, social and ecological performance at the farm-level. [PDF contains 42 pages]

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

Deep-sea benthic footprint of the Deepwater Horizon blowout

The Deepwater Horizon (DWH) accident in the northern Gulf of Mexico occurred on April 20, 2010 at a water depth of 1525 meters, and a deep-sea plume was detected within one month. Oil contacted and persisted in parts of the bottom of the deep-sea in the Gulf of Mexico. As part of the response to the accident, monitoring cruises were deployed in fall 2010 to measure potential impacts on the two main soft-bottom benthic invertebrate groups: macrofauna and meiofauna. Sediment was collected using a multicorer so that samples for chemical, physical and biological analyses could be taken simultaneously and analyzed using multivariate methods. The footprint of the oil spill was identified by creating a new variable with principal components analysis where the first factor was indicative of the oil spill impacts and this new variable mapped in a geographic information system to identify the area of the oil spill footprint. The most severe relative reduction of faunal abundance and diversity extended to 3 km from the wellhead in all directions covering an area about 24 km2. Moderate impacts were observed up to 17 km towards the southwest and 8.5 km towards the northeast of the wellhead, covering an area 148 km2. Benthic effects were correlated to total petroleum hydrocarbon, polycyclic aromatic hydrocarbons and barium concentrations, and distance to the wellhead; but not distance to hydrocarbon seeps. Thus, benthic effects are more likely due to the oil spill, and not natural hydrocarbon seepage. Recovery rates in the deep sea are likely to be slow, on the order of decades or longer.

Compositional Data Analysis in Forensic Geoscience

Statistics are regularly used to make some form of comparison between trace evidence or deploy the exclusionary principle (Morgan and Bull, 2007) in forensic investigations. Trace evidence are routinely the results of particle size, chemical or modal analyses and as such constitute compositional data. The issue is that compositional data including percentages, parts per million etc. only carry relative information. This may be problematic where a comparison of percentages and other constraint/closed data is deemed a statistically valid and appropriate way to present trace evidence in a court of law. Notwithstanding an awareness of the existence of the constant sum problem since the seminal works of Pearson (1896) and Chayes (1960) and the introduction of the application of log-ratio techniques (Aitchison, 1986; Pawlowsky-Glahn and Egozcue, 2001; Pawlowsky-Glahn and Buccianti, 2011; Tolosana-Delgado and van den Boogaart, 2013) the problem that a constant sum destroys the potential independence of variances and covariances required for correlation regression analysis and empirical multivariate methods (principal component analysis, cluster analysis, discriminant analysis, canonical correlation) is all too often not acknowledged in the statistical treatment of trace evidence. Yet the need for a robust treatment of forensic trace evidence analyses is obvious. This research examines the issues and potential pitfalls for forensic investigators if the constant sum constraint is ignored in the analysis and presentation of forensic trace evidence. Forensic case studies involving particle size and mineral analyses as trace evidence are used to demonstrate the use of a compositional data approach using a centred log-ratio (clr) transformation and multivariate statistical analyses.

Assessment of reserve effect in a marine protected area: the case study of the Professor Luiz Saldanha Marine Park (Portugal)

Dissertação de mest., Biologia Marinha (Ecologia e Conservação Marinhas), Faculdade de Ciências e Tecnologia, Univ. do Algarve, 2011