Identifying appropriate sampling and modelling approaches for analysing distributional patterns of Antarctic terrestrial arthropods along the Victoria Land latitudinal gradient

Biotic communities in Antarctic terrestrial ecosystems are relatively simple and often lack higher trophic levels (e. g. predators); thus, it is often assumed that species' distributions are mainly affected by abiotic factors such as climatic conditions, which change with increasing latitude, altitude and/or distance from the coast. However, it is becoming increasingly apparent that factors other than geographical gradients affect the distribution of organisms with low dispersal capability such as the terrestrial arthropods. In Victoria Land (East Antarctica) the distribution of springtail (Collembola) and mite (Acari) species vary at scales that range from a few square centimetres to regional and continental. Different species show different scales of variation that relate to factors such as local geological and glaciological history, and biotic interactions, but only weakly with latitudinal/altitudinal gradients. Here, we review the relevant literature and outline more appropriate sampling designs as well as suitable modelling techniques (e. g. linear mixed models and eigenvector mapping), that will more adequately address and identify the range of factors responsible for the distribution of terrestrial arthropods in Antarctica.

Spatial patterns and autocorrelation in the response of microarthropods to soil pollutants: The example of oribatid mites in an abandoned mining and smelting area

Although exogenous factors such as pollutants can act on endogenous drivers (e.g. dispersion) of populations and create spatially autocorrelated distributions, most statistical techniques assume independence of error terms. As there are no studies on metal soil pollutants and microarthropods that explicitly analyse this key issue, we completed a field study of the correlation between Oribatida and metal concentrations in litter, organic matter and soil in an attempt to account for spatial patterns of both metals and mites. The 50-m wide study area had homogenous macroscopic features, steep Pb and Cu gradients and high levels of Zn and Cd. Spatial models failed to detect metal-oribatid relationships because the observed latitudinal and longitudinal gradients in oribatid assemblages were independent of the collinear gradients in the concentration of metals. It is therefore hypothesised that other spatially variable factors (e.g. fungi, reduced macrofauna) affect oribatid assemblages, which may be influenced by metals only indirectly. (C) 2009 Elsevier Ltd. All rights reserved.

Large-scale spatial patterns in the distribution of Collembola (Hexapoda) species in Antarctic terrestrial ecosystems

We tested whether the distribution of three common springtail species (Gressittacantha terranova, Gomphiocephalus hodgsoni and Friesea grisea) in Victoria Land (Antarctica) could be modelled as a function of latitude, longitude, altitude and distance from the sea.

Victoria Land, Ross Dependency, Antarctica.

Generalized linear models were constructed using species presence/absence data relative to geographical features (latitude, longitude, altitude, distance from sea) across the species' entire ranges. Model results were then integrated with the known phylogeography of each species and hypotheses were generated on the role of climate as a major driver of Antarctic springtail distribution.

Based on model selection using Akaike's information criterion, the species' distributions were: hump-shaped relative to longitude and monotonic with altitude for Gressittacantha terranova; hump-shaped relative to latitude and monotonic with altitude for Gomphiocephalus hodgsoni; and hump-shaped relative to longitude and monotonic with latitude, altitude and distance from the sea for Friesea grisea.

No single distributional pattern was shared by the three species. While distributions were partially a response to climatic spatial clines, the patterns observed strongly suggest that past geological events have influenced the observed distributions. Accordingly, present-day spatial patterns are likely to have arisen from the interaction of historical and environmental drivers. Future studies will need to integrate a range of spatial and temporal scales to further quantify their respective roles.

Modelling local-scale determinants and the probability of microarthropod species occurrence in Antarctic soils

Soil fauna in the extreme conditions of Antarctica consists of a few microinvertebrate species patchily distributed at different spatial scales. Populations of the prostigmatic mite Stereotydeus belli and the collembolan Gressittacantha terranova from northern Victoria Land (Antarctica) were used as models to study the effect of soil properties on microarthropod distributions. In agreement with the general assumption that the development and distribution of life in these ecosystems is mainly controlled by abiotic factors, we found that the probability of occurrence of S. belli depends on soil moisture and texture and on the sampling period (which affects the general availability of water); surprisingly, none of the analysed variables were significantly related to the G. terranova distribution. Based on our results and literature data, we propose a theoretical model that introduces biotic interactions among the major factors driving the local distribution of collembolans in Antarctic terrestrial ecosystems. (c) 2007 Elsevier Ltd. All rights reserved.

Understanding macroalgal dispersal in a complex hydrodynamic environment: a combined population genetic and physical modelling approach

Gene flow in macroalgal populations can be strongly influenced by spore or gamete dispersal. This, in turn, is influenced by a convolution of the effects of current flow and specific plant reproductive strategies. Although several studies have demonstrated genetic variability in macroalgal populations over a wide range of spatial scales, the associated current data have generally been poorly resolved spatially and temporally. In this study, we used a combination of population genetic analyses and high-resolution hydrodynamic modelling to investigate potential connectivity between populations of the kelp Laminaria digitata in the Strangford Narrows, a narrow channel characterized by strong currents linking the large semi-enclosed sea lough, Strangford Lough, to the Irish Sea. Levels of genetic structuring based on six microsatellite markers were very low, indicating high levels of gene flow and a pattern of isolation-by-distance, where populations are more likely to exchange migrants with geographically proximal populations, but with occasional long-distance dispersal. This was confirmed by the particle tracking model, which showed that, while the majority of spores settle near the release site, there is potential for dispersal over several kilometres. This combined population genetic and modelling approach suggests that the complex hydrodynamic environment at the entrance to Strangford Lough can facilitate dispersal on a scale exceeding that proposed for L. digitata in particular, and the majority of macroalgae in general. The study demonstrates the potential of integrated physical–biological approaches for the prediction of ecological changes resulting from factors such as anthropogenically induced coastal zone changes.

Is Mortality Spatial or Social?

Mortality modelling for the purposes of demographic forecasting and actuarial pricing is generally done at an aggregate level using national data. Modelling at this level fails to capture the variation in mortality within country and potentially leads to a mis-specification of mortality forecasts for a subset of the population. This can have detrimental effects for pricing and reserving in the actuarial context. In this paper we consider mortality rates at a regional level and analyse the variation in those rates. We consider whether variation in mortality rates within a country can be explained using local economic and social variables. Using Northern Ireland data on mortality and measures of deprivation we identify the variables explaining mortality variation. We create a population polarisation variable and find that this variable is significant in explaining some of the variation in mortality rates. Further, we consider whether spatial and non-spatial models have a part to play in explaining mortality differentials.

Numerical modelling of low temperature laser produced magnesium plasma

The two-dimensional laser-plasma-interaction hydrodynamic code POLLUX has been used to simulate the ablation of a magnesium target by a 30-ns, 248-nm KrF excimer laser at low laser fluences of ≤10 J cm2. This code, originally written for much higher laser intensities, has been recently extended to include a detailed description of the equation of state in order to treat changes of phase within the target material, and also includes a Thomas Fermi description of the electrons. The simulated temporal and spatial evolution of the plasma plume in the early phase of the expansion (≤100 ns) is compared with experimental interferometric measurements of electron density. The expansion dynamics are in good agreement, although the simulated electron number density is about 2.5 times higher than the experimental values.

Exploiting Human Resource Requirements to Infer Human Movement Patterns for Use in Modelling Disease Transmission Systems: An Example from Eastern Province, Zambia

In this research, an agent-based model (ABM) was developed to generate human movement routes between homes and water resources in a rural setting, given commonly available geospatial datasets on population distribution, land cover and landscape resources. ABMs are an object-oriented computational approach to modelling a system, focusing on the interactions of autonomous agents, and aiming to assess the impact of these agents and their interactions on the system as a whole. An A* pathfinding algorithm was implemented to produce walking routes, given data on the terrain in the area. A* is an extension of Dijkstra's algorithm with an enhanced time performance through the use of heuristics. In this example, it was possible to impute daily activity movement patterns to the water resource for all villages in a 75 km long study transect across the Luangwa Valley, Zambia, and the simulated human movements were statistically similar to empirical observations on travel times to the water resource (Chi-squared, 95% confidence interval). This indicates that it is possible to produce realistic data regarding human movements without costly measurement as is commonly achieved, for example, through GPS, or retrospective or real-time diaries. The approach is transferable between different geographical locations, and the product can be useful in providing an insight into human movement patterns, and therefore has use in many human exposure-related applications, specifically epidemiological research in rural areas, where spatial heterogeneity in the disease landscape, and space-time proximity of individuals, can play a crucial role in disease spread.

Numerical modelling of radionuclide migration in near-surface systems

Em todo o mundo são usados, hoje em dia, modelos numéricos hidrogeoquímicos para simular fenómenos naturais e fenómenos decorrentes de actividades antrópicas. Estes modelos ajudam-nos a compreender o ambiente envolvente, a sua variabilidade espacial e evolução temporal. No presente trabalho apresenta-se o desenvolvimento de modelos numéricos hidrogeoquímicos aplicados no contexto do repositório geológico profundo para resíduos nucleares de elevada actividade. A avaliação da performance de um repositório geológico profundo inclui o estudo da evolução geoquímica do repositório, bem como a análise dos cenários de mau funcionamento do repositório, e respectivas consequências ambientais. Se se escaparem acidentalmente radionuclídeos de um repositório, estes poderão atravessar as barreiras de engenharia e barreiras naturais que constituem o repositório, atingindo eventualmente, os ecosistemas superficiais. Neste caso, os sedimentos subsuperficiais constituem a última barreira natural antes dos ecosistemas superficiais. No presente trabalho foram desenvolvidos modelos numéricos que integram processos biogeoquímicos, geoquímicos, hidrodinâmicos e de transporte de solutos, para entender e quantificar a influência destes processos na mobilidade de radionuclídeos em sistemas subsuperficiais. Os resultados alcançados reflectem a robustez dos instrumentos numéricos utilizados para desenvolver simulações descritivas e predictivas de processos hidrogeoquímicos que influenciam a mobilidade de radionuclídeos. A simulação (descritiva) de uma experiência laboratorial revela que a actividade microbiana induz a diminuição do potencial redox da água subterrânea que, por sua vez, favorece a retenção de radionuclídeos sensíveis ao potencial redox, como o urânio. As simulações predictivas indicam que processos de co-precipitação com minerais de elementos maioritários, precipitação de fases puras, intercâmbio catiónico e adsorção à superfície de minerais favorecem a retenção de U, Cs, Sr e Ra na fase sólida de uma argila glaciar e uma moreia rica em calcite. A etiquetagem dos radionuclídeos nas simulações numéricas permitiu concluir que a diluição isotópica joga um papel importante no potencial impacte dos radionuclídeos nos sistemas subsuperficiais. A partir dos resultados das simulações numéricas é possivel calcular coeficientes de distribuição efectivos. Esta metodologia proporciona a simulação de ensaios de traçadores de longa duração que não seriam exequíveis à escala da vida humana. A partir destas simulações podem ser obtidos coeficientes de retardamento que são úteis no contexto da avaliação da performance de repositórios geológicos profundos.

Exposure modelling to traffic-related air toxic pollutants

Atualmente, a poluição atmosférica constitui uma das principais causas ambientais de mortalidade. Cerca de 30% da população europeia residente em áreas urbanas encontra-se exposta a níveis de poluição atmosférica superiores aos valores- limite de qualidade do ar legislados para proteção da saúde humana, representando o tráfego rodoviário uma das principais fontes de poluição urbana. Além dos poluentes tradicionais avaliados em áreas urbanas, os poluentes classificados como perigosos para a saúde (Hazard Air Pollutants - HAPs) têm particular relevância devido aos seus conhecidos efeitos tóxicos e cancerígenos. Neste sentido, a avaliação da exposição tornase primordial para a determinação da relação entre a poluição atmosférica urbana e efeitos na saúde. O presente estudo tem como principal objetivo o desenvolvimento e implementação de uma metodologia para avaliação da exposição individual à poluição atmosférica urbana relacionada com o tráfego rodoviário. De modo a atingir este objetivo, foram identificados os parâmetros relevantes para a quantificação de exposição e analisados os atuais e futuros potenciais impactos na saúde associados com a exposição à poluição urbana. Neste âmbito, o modelo ExPOSITION (EXPOSure model to traffIc-relaTed aIr pOllutioN) foi desenvolvido baseado numa abordagem inovadora que envolve a análise da trajetória dos indivíduos recolhidas por telemóveis com tecnologia GPS e processadas através da abordagem de data mining e análise geoespacial. O modelo ExPOSITION considera também uma abordagem probabilística para caracterizar a variabilidade dos parâmetros microambientais e a sua contribuição para exposição individual. Adicionalmente, de forma a atingir os objetivos do estudo foi desenvolvido um novo módulo de cálculo de emissões de HAPs provenientes do transporte rodoviário. Neste estudo, um sistema de modelação, incluindo os modelos de transporteemissões- dispersão-exposição, foi aplicado na área urbana de Leiria para quantificação de exposição individual a PM2.5 e benzeno. Os resultados de modelação foram validados com base em medições obtidas por monitorização pessoal e monitorização biológica verificando-se uma boa concordância entre os resultados do modelo e dados de medições. A metodologia desenvolvida e implementada no âmbito deste trabalho permite analisar e estimar a magnitude, frequência e inter e intra-variabilidade dos níveis de exposição individual, bem como a contribuição de diferentes microambientes, considerando claramente a sequência de eventos de exposição e relação fonte-recetor, que é fundamental para avaliação dos efeitos na saúde e estudos epidemiológicos. O presente trabalho contribui para uma melhor compreensão da exposição individual em áreas urbanas, proporcionando novas perspetivas sobre a exposição individual, essenciais na seleção de estratégias de redução da exposição à poluição atmosférica urbana, e consequentes efeitos na saúde.

Modelling the Minho river plume intrusion into the rias baixas

The Minho River, situated 30 km south of the Rias Baixas is the most important freshwater source flowing into the Western Galician Coast (NW of the Iberian Peninsula). This discharge is important to determine the hydrological patterns adjacent to its mouth, particularly close to the Galician coastal region. The buoyancy generated by the Minho plume can flood the Rias Baixas for long periods, reversing the normal estuarine density gradients. Thus, it becomes important to analyse its dynamics as well as the thermohaline patterns of the areas affected by the freshwater spreading. Thus, the main aim of this work was to study the propagation of the Minho estuarine plume to the Rias Baixas, establishing the conditions in which this plume affects the circulation and hydrographic features of these coastal systems, through the development and application of the numerical model MOHID. For this purpose, the hydrographic features of the Rias Baixas mouths were studied. It was observed that at the northern mouths, due to their shallowness, the heat fluxes between the atmosphere and ocean are the major forcing, influencing the water temperature, while at the southern mouths the influence of the upwelling events and the Minho River discharge were more frequent. The salinity increases from south to north, revealing that the observed low values may be caused by the Minho River freshwater discharge. An assessment of wind data along the Galician coast was carried out, in order to evaluate the applicability of the study to the dispersal of the Minho estuarine plume. Firstly, a comparative analysis between winds obtained from land meteorological stations and offshore QuikSCAT satellite were performed. This comparison revealed that satellite data constitute a good approach to study wind induced coastal phenomena. However, since the numerical model MOHID requires wind data with high spatial and temporal resolution close to the coast, results of the forecasted model WRF were added to the previous study. The analyses revealed that the WRF model data is a consistent tool to obtain representative wind data near the coast, showing good results when comparing with in situ wind observations from oceanographic buoys. To study the influence of the Minho buoyant discharge influence on the Rias Baixas, a set of three one-way nested models was developed and implemented, using the numerical model MOHID. The first model domain is a barotropic model and includes the whole Iberian Peninsula coast. The second and third domains are baroclinic models, where the second domain is a coarse representation of the Rias Baixas and adjacent coastal area, while the third includes the same area with a higher resolution. A bi-dimensional model was also implemented in the Minho estuary, in order to quantify the flow (and its properties) that the estuary injects into the ocean. The chosen period for the Minho estuarine plume propagation validation was the spring of 1998, since a high Minho River discharge was reported, as well as favourable wind patterns to advect the estuarine plume towards the Rias Baixas, and there was field data available to compare with the model predictions. The obtained results show that the adopted nesting methodology was successful implemented. Model predictions reproduce accurately the hydrodynamics and thermohaline patterns on the Minho estuary and Rias Baixas. The importance of the Minho river discharge and the wind forcing in the event of May 1998 was also studied. The model results showed that a continuous moderate Minho River discharge combined with southerly winds is enough to reverse the Rias Baixas circulation pattern, reducing the importance of the occurrence of specific events of high runoff values. The conditions in which the estuarine plume Minho affects circulation and hydrography of the Rias Baixas were evaluated. The numerical results revealed that the Minho estuarine plume responds rapidly to wind variations and is also influenced by the bathymetry and morphology of the coastline. Without wind forcing, the plume expands offshore, creating a bulge in front of the river mouth. When the wind blows southwards, the main feature is the offshore extension of the plume. Otherwise, northward wind spreads the river plume towards the Rias Baixas. The plume is confined close to the coast, reaching the Rias Baixas after 1.5 days. However, for Minho River discharges higher than 800 m3 s-1, the Minho estuarine plume reverses the circulation patterns in the Rias Baixas. It was also observed that the wind stress and Minho River discharge are the most important factors influencing the size and shape of the Minho estuarine plume. Under the same conditions, the water exchange between Rias Baixas was analysed following the trajectories particles released close to the Minho River mouth. Over 5 days, under Minho River discharges higher than 2100 m3 s-1 combined with southerly winds of 6 m s-1, an intense water exchange between Rias was observed. However, only 20% of the particles found in Ria de Pontevedra come directly from the Minho River. In summary, the model application developed in this study contributed to the characterization and understanding of the influence of the Minho River on the Rias Baixas circulation and hydrography, highlighting that this methodology can be replicated to other coastal systems.

Channel modelling for MIMO systems

Systems equipped with multiple antennas at the transmitter and at the receiver, known as MIMO (Multiple Input Multiple Output) systems, offer higher capacities, allowing an efficient exploitation of the available spectrum and/or the employment of more demanding applications. It is well known that the radio channel is characterized by multipath propagation, a phenomenon deemed problematic and whose mitigation has been achieved through techniques such as diversity, beamforming or adaptive antennas. By exploring conveniently the spatial domain MIMO systems turn the characteristics of the multipath channel into an advantage and allow creating multiple parallel and independent virtual channels. However, the achievable benefits are constrained by the propagation channel’s characteristics, which may not always be ideal. This work focuses on the characterization of the MIMO radio channel. It begins with the presentation of the fundamental results from information theory that triggered the interest on these systems, including the discussion of some of their potential benefits and a review of the existing channel models for MIMO systems. The characterization of the MIMO channel developed in this work is based on experimental measurements of the double-directional channel. The measurement system is based on a vector network analyzer and a two-dimensional positioning platform, both controlled by a computer, allowing the measurement of the channel’s frequency response at the locations of a synthetic array. Data is then processed using the SAGE (Space-Alternating Expectation-Maximization) algorithm to obtain the parameters (delay, direction of arrival and complex amplitude) of the channel’s most relevant multipath components. Afterwards, using a clustering algorithm these data are grouped into clusters. Finally, statistical information is extracted allowing the characterization of the channel’s multipath components. The information about the multipath characteristics of the channel, induced by existing scatterers in the propagation scenario, enables the characterization of MIMO channel and thus to evaluate its performance. The method was finally validated using MIMO measurements.

Combining sampling, profiling, remote sensing and modelling to efficiently monitoring estuarine and coastal water quality

Monitoring of coastal and estuarine water quality has been traditionally performed by sampling with subsequent laboratory analysis. This has the disadvantages of low spatial and temporal resolution and high cost. In the last decades two alternative techniques have emerged to overcome this drawback: profiling and remote sensing. Profiling using multi-parameter sensors is now in a commercial stage. It can be used, tied to a boat, to obtain a quick “picture” of the system. The spatial resolution thus increases from single points to a line coincident with the boat track. The temporal resolution however remains unchanged since campaigns and resources involved are basically the same. The need for laboratory analysis was reduced but not eliminated because parameters like nutrients, microbiology or metals are still difficult to obtain with sensors and validation measurements are still needed. In the last years the improvement in satellite resolution has enabled its use for coastal and estuarine water monitoring. Although spatial coverage and resolution of satellite images in the present is already suitable to coastal and estuarine monitoring, temporal resolution is naturally limited to satellite passages and cloud cover. With this panorama the best approach to water monitoring is to integrate and combine data from all these sources. The natural tools to perform this integration are numerical models. Models benefit from the different sources of data to obtain a better calibration. After calibration they can be used to extend spatially and temporally the methods resolution. In Algarve (South of Portugal) a monitoring effort using this approach is being undertaken. The monitoring effort comprises five different locations including coastal waters, estuaries and coastal lagoons. The objective is to establish the base line situation to evaluate the impact of Waste Water Treatment Plants design and retrofitting. The field campaigns include monthly synoptic profiling, using an YSI 6600 multi-parameter system, laboratory analysis and fixed stations. The remote sensing uses ENVISAT\MERIS Level 2 Full Resolution data. This data is combined and used with the MOHID modelling system to obtain an integrate description of the systems. The results show the limitations of each method and the ability of the modelling system to integrate the results and to produce a comprehensive picture of the system.

Modelling underwater acoustic noise as a tool for coastal management

Dissertação de Mestrado, Gestão da Água e da Costa, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2009

Defining Spatial and Temporal Hydromorphological Sampling Strategies for the Leigh Brook River Site

Detailed surveys of depth and velocity are undertaken to describe hydro-ecological status of rivers. Fieldwork for these surveys is time consuming and expensive. This paper aims to describe the methodology applied in order to determine the most suitable depth sampling strategy for effective field data collection and river representation in time and space at the Leigh Brook river site, Worcester, UK. The accuracy of three different sampling strategies for predicting depth at non-measured points has been compared and the mesohabitats that better characterise depth changes due to variations in discharge have been identified. The results show that depth changes due to discharge change are mainly located at shallow and deep glide mesohabitat types. The analysis for the comparison of sampling strategies indicates that grid sampling strategies give better results than regular transects. Since the results also show that higher errors in predictions are obtained in the deepest areas, higher sampling densities should be applied in these locations.