Laser beam control, combining, and propagation in atmospheric turbulence

This dissertation is concerned with the control, combining, and propagation of laser beams through a turbulent atmosphere. In the first part we consider adaptive optics: the process of controlling the beam based on information of the current state of the turbulence. If the target is cooperative and provides a coherent return beam, the phase measured near the beam transmitter and adaptive optics can, in principle, correct these fluctuations. However, for many applications, the target is uncooperative. In this case, we show that an incoherent return from the target can be used instead. Using the principle of reciprocity, we derive a novel relation between the field at the target and the scattered field at a detector. We then demonstrate through simulation that an adaptive optics system can utilize this relation to focus a beam through atmospheric turbulence onto a rough surface. In the second part we consider beam combining. To achieve the power levels needed for directed energy applications it is necessary to combine a large number of lasers into a single beam. The large linewidths inherent in high-power fiber and slab lasers cause random phase and intensity fluctuations occurring on sub-nanosecond time scales. We demonstrate that this presents a challenging problem when attempting to phase-lock high-power lasers. Furthermore, we show that even if instruments are developed that can precisely control the phase of high-power lasers; coherent combining is problematic for DE applications. The dephasing effects of atmospheric turbulence typically encountered in DE applications will degrade the coherent properties of the beam before it reaches the target. Finally, we investigate the propagation of Bessel and Airy beams through atmospheric turbulence. It has been proposed that these quasi-non-diffracting beams could be resistant to the effects of atmospheric turbulence. However, we find that atmospheric turbulence disrupts the quasi-non-diffracting nature of Bessel and Airy beams when the transverse coherence length nears the initial aperture diameter or diagonal respectively. The turbulence induced transverse phase distortion limits the effectiveness of Bessel and Airy beams for applications requiring propagation over long distances in the turbulent atmosphere.

THE PHYSICS OF IDEAS: INFERRING THE MECHANICS OF OPINION FORMATION FROM MACROSCOPIC STATISTICAL PATTERNS

In a microscopic setting, humans behave in rich and unexpected ways. In a macroscopic setting, however, distinctive patterns of group behavior emerge, leading statistical physicists to search for an underlying mechanism. The aim of this dissertation is to analyze the macroscopic patterns of competing ideas in order to discern the mechanics of how group opinions form at the microscopic level. First, we explore the competition of answers in online Q&A (question and answer) boards. We find that a simple individual-level model can capture important features of user behavior, especially as the number of answers to a question grows. Our model further suggests that the wisdom of crowds may be constrained by information overload, in which users are unable to thoroughly evaluate each answer and therefore tend to use heuristics to pick what they believe is the best answer. Next, we explore models of opinion spread among voters to explain observed universal statistical patterns such as rescaled vote distributions and logarithmic vote correlations. We introduce a simple model that can explain both properties, as well as why it takes so long for large groups to reach consensus. An important feature of the model that facilitates agreement with data is that individuals become more stubborn (unwilling to change their opinion) over time. Finally, we explore potential underlying mechanisms for opinion formation in juries, by comparing data to various types of models. We find that different null hypotheses in which jurors do not interact when reaching a decision are in strong disagreement with data compared to a simple interaction model. These findings provide conceptual and mechanistic support for previous work that has found mutual influence can play a large role in group decisions. In addition, by matching our models to data, we are able to infer the time scales over which individuals change their opinions for different jury contexts. We find that these values increase as a function of the trial time, suggesting that jurors and judicial panels exhibit a kind of stubbornness similar to what we include in our model of voting behavior.

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

Mixed Layer formation and restratification in presence of mesoscale and submesoscale turbulence

Recent realistic high resolution modeling studies show a net increase of submesoscale activity in fall and winter when the mixed layer depth is at its maximum. This submesoscale activity increase is associated with a reduced deepening of the mixed layer. Both phenomena can be related to the development of mixed layer instabilities, which convert available potential energy into submesoscale eddy kinetic energy and contribute to a fast restratification by slumping the horizontal density gradient in the mixed layer. In the present work, the mixed layer formation and restratification was studied by uniformly cooling a fully turbulent zonal jet in a periodic channel at different resolutions, from eddy resolving (10 km) to submesoscale permitting (2 km). The effect of the submesoscale activity, highlighted by these different horizontal resolutions, was quantified in terms of mixed layer depth, restratification rate and buoyancy fluxes. Contrary to many idealized studies focusing on the restratification phase only, this study addresses a continuous event of mixed layer formation followed by its complete restratification. The robustness of the present results was established by ensemble simulations. The results show that, at higher resolution, when submesoscale starts to be resolved, the mixed layer formed during the surface cooling is significantly shallower and the total restratification almost three times faster. Such differences between coarse and fine resolution models are consistent with the submesoscale upward buoyancy flux, which balances the convection during the formation phase and accelerates the restratification once the surface cooling is stopped. This submesoscale buoyancy flux is active even below the mixed layer. Our simulations show that mesoscale dynamics also cause restratification, but on longer time scales. Finally, the spatial distribution of the mixed layer depth is highly heterogeneous in the presence of submesoscale activity, prompting the question of whether it is possible to parameterize submesoscale effects and their effects on the marine biology as a function of a spatially-averaged mixed layer depth.

A study of electrochemical transport and diffuse charge dynamics using a Langevin equation

This thesis aims to develop new numerical and computational tools to study electrochemical transport and diffuse charge dynamics at small scales. Previous efforts at modeling electrokinetic phenomena at scales where the noncontinuum effects become significant have included continuum models based on the Poisson-Nernst-Planck equations and atomic simulations using molecular dynamics algorithms. Neither of them is easy to use or conducive to electrokinetic transport modeling in strong confinement or over long time scales. This work introduces a new approach based on a Langevin equation for diffuse charge dynamics in nanofluidic devices, which incorporates features from both continuum and atomistic methods. The model is then extended to include steric effects resulting from finite ion size, and applied to the phenomenon of double layer charging in a symmetric binary electrolyte between parallel-plate blocking electrodes, between which a voltage is applied. Finally, the results of this approach are compared to those of the continuum model based on the Poisson-Nernst-Planck equations.

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

Context dependency of trait repeatability and its relevance for management and conservation of fish populations

Repeatability of behavioural and physiological traits is increasingly a focus for animal researchers, for which fish have become important models. Almost all of this work has been done in the context of evolutionary ecology, with few explicit attempts to apply repeatability and context dependency of trait variation toward understanding conservation-related issues. Here, we review work examining the degree to which repeatability of traits (such as boldness, swimming performance, metabolic rate and stress responsiveness) is context dependent. We review methods for quantifying repeatability (distinguishing between within-context and across-context repeatability) and confounding factors that may be especially problematic when attempting to measure repeatability in wild fish. Environmental factors such temperature, food availability, oxygen availability, hypercapnia, flow regime and pollutants all appear to alter trait repeatability in fishes. This suggests that anthropogenic environmental change could alter evolutionary trajectories by changing which individuals achieve the greatest fitness in a given set of conditions. Gaining a greater understanding of these effects will be crucial for our ability to forecast the effects of gradual environmental change, such as climate change and ocean acidification, the study of which is currently limited by our ability to examine trait changes over relatively short time scales. Also discussed are situations in which recent advances in technologies associated with electronic tags (biotelemetry and biologging) and respirometry will help to facilitate increased quantification of repeatability for physiological and integrative traits, which so far lag behind measures of repeatability of behavioural traits.

First metabolomic approach of the epiphytic bacteria-marine diatom Haslea ostrearia relationships

The marine diatom Haslea ostrearia [1] produces a water-soluble blue-pigment named marennine [2] of economic interest. But the lack of knowledge of the ecological conditions, under which this microalga develops in its natural ecosystem, more especially bacteria H. ostrearia interactions, prevents any optimization of its culture in well-controlled conditions. The structure of the bacterial community was analyzed by PCR-TTGE before and after the isolation of H. ostrearia cells recovered from 4 localities, to distinguish the relative part of the biotope and the biocenose and eventually to describe the temporal dynamic of the structure of the bacterial community at two time-scales. The differences in genetic fingerprints, more especially high between two H. ostrearia isolates (HO-R and HO-BM) showed also the highest differences in the bacterial structure [3] as the result of specific metabolomics profiles. The non-targeted metabolomic investigation showed that these profiles were more distinct in case of bacteria-alga associations than for the H. ostrearia monoculture Here we present a Q-TOF LC/MS metabolomic fingerprinting approach [3]: - to investigate differential metabolites of axenic versus non axenic H. ostrearia cultures. - to focus on the specific metabolites of a bacterial surrounding associated with the activation or inhibition of the microalga growing. The Agilent suite of data processing software makes feature finding, statistical analysis, and identification easier. This enables rapid transformation of complex raw data into biologically relevant metabolite information.

Preliminary metabolomic approach on the bacterial structure community of Haslea ostrearia

The marine diatom Haslea ostrearia [1] produces a water-soluble blue-pigment named marennine [2] of economic interest. But the lack of knowledge of the ecological conditions, under which this microalga develops in its natural ecosystem, more especially bacteria H. ostrearia interactions, prevents any optimization of its culture in well-controlled conditions. The structure of the bacterial community was analyzed by PCR-TTGE before and after the isolation of H. ostrearia cells recovered from 4 localities, to distinguish the relative part of the biotope and the biocenose and eventually to describe the temporal dynamic of the structure of the bacterial community at two time-scales. The differences in genetic fingerprints, more especially high between two H. ostrearia isolates (HO-R and HO-BM) showed also the highest differences in the bacterial structure [3] as the result of specific metabolomics profiles. The non-targeted metabolomic investigation showed that these profiles were more distinct in case of bacteria-alga associations than for the H. ostrearia monoculture Here we present a Q-TOF LC/MS metabolomic fingerprinting approach [3]: - to investigate differential metabolites of axenic versus non axenic H. ostrearia cultures. - to focus on the specific metabolites of a bacterial surrounding associated with the activation or inhibition of the microalga growing. The Agilent suite of data processing software makes feature finding, statistical analysis, and identification easier. This enables rapid transformation of complex raw data into biologically relevant metabolite information.

Engineering a Quantum Many-body Hamiltonian with Trapped Ions

While fault-tolerant quantum computation might still be years away, analog quantum simulators offer a way to leverage current quantum technologies to study classically intractable quantum systems. Cutting edge quantum simulators such as those utilizing ultracold atoms are beginning to study physics which surpass what is classically tractable. As the system sizes of these quantum simulators increase, there are also concurrent gains in the complexity and types of Hamiltonians which can be simulated. In this work, I describe advances toward the realization of an adaptable, tunable quantum simulator capable of surpassing classical computation. We simulate long-ranged Ising and XY spin models which can have global arbitrary transverse and longitudinal fields in addition to individual transverse fields using a linear chain of up to 24 Yb+ 171 ions confined in a linear rf Paul trap. Each qubit is encoded in the ground state hyperfine levels of an ion. Spin-spin interactions are engineered by the application of spin-dependent forces from laser fields, coupling spin to motion. Each spin can be read independently using state-dependent fluorescence. The results here add yet more tools to an ever growing quantum simulation toolbox. One of many challenges has been the coherent manipulation of individual qubits. By using a surprisingly large fourth-order Stark shifts in a clock-state qubit, we demonstrate an ability to individually manipulate spins and apply independent Hamiltonian terms, greatly increasing the range of quantum simulations which can be implemented. As quantum systems grow beyond the capability of classical numerics, a constant question is how to verify a quantum simulation. Here, I present measurements which may provide useful metrics for large system sizes and demonstrate them in a system of up to 24 ions during a classically intractable simulation. The observed values are consistent with extremely large entangled states, as much as ~95% of the system entangled. Finally, we use many of these techniques in order to generate a spin Hamiltonian which fails to thermalize during experimental time scales due to a meta-stable state which is often called prethermal. The observed prethermal state is a new form of prethermalization which arises due to long-range interactions and open boundary conditions, even in the thermodynamic limit. This prethermalization is observed in a system of up to 22 spins. We expect that system sizes can be extended up to 30 spins with only minor upgrades to the current apparatus. These results emphasize that as the technology improves, the techniques and tools developed here can potentially be used to perform simulations which will surpass the capability of even the most sophisticated classical techniques, enabling the study of a whole new regime of quantum many-body physics.

Guided Migration and Collective Behavior: Cell Dynamics during Cancer Progression

This dissertation focuses on gaining understanding of cell migration and collective behavior through a combination of experiment, analysis, and modeling techniques. Cell migration is a ubiquitous process that plays an important role during embryonic development and wound healing as well as in diseases like cancer, which is a particular focus of this work. As cancer cells become increasingly malignant, they acquire the ability to migrate away from the primary tumor and spread throughout the body to form metastatic tumors. During this process, changes in gene expression and the surrounding tumor environment can lead to changes in cell migration characteristics. In this thesis, I analyze how cells are guided by the texture of their environment and how cells cooperate with their neighbors to move collectively. The emergent properties of collectively moving groups are a particular focus of this work as collective cell dynamics are known to change in diseases such as cancer. The internal machinery for cell migration involves polymerization of the actin cytoskeleton to create protrusions that---in coordination with retraction of the rear of the cell---lead to cell motion. This actin machinery has been previously shown to respond to the topography of the surrounding surface, leading to guided migration of amoeboid cells. Here we show that epithelial cells on nanoscale ridge structures also show changes in the morphology of their cytoskeletons; actin is found to align with the ridge structures. The migration of the cells is also guided preferentially along the ridge length. These ridge structures are on length scales similar to those found in tumor microenvironments and as such provide a system for studying the response of the cells' internal migration machinery to physiologically relevant topographical cues. In addition to sensing surface topography, individual cells can also be influenced by the pushing and pulling of neighboring cells. The emergent properties of collectively migrating cells show interesting dynamics and are relevant for cancer progression, but have been less studied than the motion of individual cells. We use Particle Image Velocimetry (PIV) to extract the motion of a collectively migrating cell sheet from time lapse images. The resulting flow fields allow us to analyze collective behavior over multiple length and time scales. To analyze the connection between individual cell properties and collective migration behavior, we compare experimental flow fields with the migration of simulated cell groups. Our collective migration metrics allow for a quantitative comparison between experimental and simulated results. This comparison shows that tissue-scale decreases in collective behavior can result from changes in individual cell activity without the need to postulate the existence of subpopulations of leader cells or global gradients. In addition to tissue-scale trends in collective behavior, the migration of cell groups includes localized dynamic features such as cell rearrangements. An individual cell may smoothly follow the motion of its neighbors (affine motion) or move in a more individualistic manner (non-affine motion). By decomposing individual motion into both affine and non-affine components, we measure cell rearrangements within a collective sheet. Finally, finite-time Lyapunov exponent (FTLE) values capture the stretching of the flow field and reflect its chaotic character. Applying collective migration analysis techniques to experimental data on both malignant and non-malignant human breast epithelial cells reveals differences in collective behavior that are not found from analyzing migration speeds alone. Non-malignant cells show increased cooperative motion on long time scales whereas malignant cells remain uncooperative as time progresses. Combining multiple analysis techniques also shows that these two cell types differ in their response to a perturbation of cell-cell adhesion through the molecule E-cadherin. Non-malignant MCF10A cells use E-cadherin for short time coordination of collective motion, yet even with decreased E-cadherin expression, the cells remain coordinated over long time scales. In contrast, the migration behavior of malignant and invasive MCF10CA1a cells, which already shows decreased collective dynamics on both time scales, is insensitive to the change in E-cadherin expression.

The Westerly Index as complementary indicator of the North Atlantic oscillation in explaining drought variability across Europe

This paper analyses the influence of different atmospheric circulation indices on the multi-scalar drought variability across Europe by using the Standardized Precipitation Evapotranspiration Index (SPEI). The monthly circulation indices used in this study include the North Atlantic oscillation (NAO), the East Atlantic (EA), the Scandinavian (SCAN) and the East Atlantic-Western Russia (EA-WR) patterns, as well as the recently published Westerly Index (WI), defined as the persistence of westerly winds over the eastern north Atlantic region. The results indicate that European drought variability is better explained by the station-based NAO index and the WI than by any other combination of circulation indices. In northern and central Europe the variability of drought severity for different seasons and time-scales is strongly associated with the WI. On the contrary, the influence of the NAO on southern Europe droughts is stronger than that exerted by the WI. The correlation patterns of the NAO and WI with the SPEI show a spatial complementarity in shaping drought variability across Europe. Lagged correlations of the NAO and WI with the SPEI also indicate enough skill of both indices to anticipate drought severity several months in advance. As long as instrumental series of the NAO and WI are available, their combined use would allow inferring European drought variability for the last two centuries and improve the calibration and interpretation of paleoclimatic proxies associated with drought.

Investigation of capacitive discharge heating of metallic glasses

In recent years, the discovery of bulk metallic glasses with exceptional properties has generated much interest. One of their most intriguing features is their capacity for viscous flow above the glass transition temperature. This characteristic allows metallic glasses to be formed like plastics at modest temperatures. However, crystallization of supercooled metallic liquids in the best bulk metallic glass-formers is much more rapid than in most polymers and silicate glass-forming liquids. The short times to crystallization impairs experimentation on and processing of supercooled glass-forming metallic liquids. A technique to rapidly and uniformly heat metallic glasses at rates of 10⁵ to 10⁶ kelvin per second is presented. A capacitive discharge is used to ohmically heat metallic glasses to temperatures in the super cooled liquid region in millisecond time-scales. By heating samples rapidly, the most time-consuming step in experiments on supercooled metallic liquids is reduced orders of magnitude in length. This allows for experimentation on and processing of metallic liquids in temperature ranges that were previously inaccessible because of crystallization.

A variety of forming techniques, including injection molding and forging, were coupled with capacitive discharge heating to produce near net-shaped metallic glass parts. In addition, a new forming technique, which combines a magnetic field with the heating current to produce a forming force, was developed. Viscosities were measured in previously inaccessible temperature ranges using parallel plate rheometry combined with capacitive discharge heating. Lastly, a rapid pulse calorimeter was developed with this technique to investigate the thermophysical behavior of metallic glasses at these rapid heating rates.

Guided Migration and Collective Behavior: Cell Dynamics during Cancer Progression

This dissertation focuses on gaining understanding of cell migration and collective behavior through a combination of experiment, analysis, and modeling techniques. Cell migration is a ubiquitous process that plays an important role during embryonic development and wound healing as well as in diseases like cancer, which is a particular focus of this work. As cancer cells become increasingly malignant, they acquire the ability to migrate away from the primary tumor and spread throughout the body to form metastatic tumors. During this process, changes in gene expression and the surrounding tumor environment can lead to changes in cell migration characteristics. In this thesis, I analyze how cells are guided by the texture of their environment and how cells cooperate with their neighbors to move collectively. The emergent properties of collectively moving groups are a particular focus of this work as collective cell dynamics are known to change in diseases such as cancer. The internal machinery for cell migration involves polymerization of the actin cytoskeleton to create protrusions that---in coordination with retraction of the rear of the cell---lead to cell motion. This actin machinery has been previously shown to respond to the topography of the surrounding surface, leading to guided migration of amoeboid cells. Here we show that epithelial cells on nanoscale ridge structures also show changes in the morphology of their cytoskeletons; actin is found to align with the ridge structures. The migration of the cells is also guided preferentially along the ridge length. These ridge structures are on length scales similar to those found in tumor microenvironments and as such provide a system for studying the response of the cells' internal migration machinery to physiologically relevant topographical cues. In addition to sensing surface topography, individual cells can also be influenced by the pushing and pulling of neighboring cells. The emergent properties of collectively migrating cells show interesting dynamics and are relevant for cancer progression, but have been less studied than the motion of individual cells. We use Particle Image Velocimetry (PIV) to extract the motion of a collectively migrating cell sheet from time lapse images. The resulting flow fields allow us to analyze collective behavior over multiple length and time scales. To analyze the connection between individual cell properties and collective migration behavior, we compare experimental flow fields with the migration of simulated cell groups. Our collective migration metrics allow for a quantitative comparison between experimental and simulated results. This comparison shows that tissue-scale decreases in collective behavior can result from changes in individual cell activity without the need to postulate the existence of subpopulations of leader cells or global gradients. In addition to tissue-scale trends in collective behavior, the migration of cell groups includes localized dynamic features such as cell rearrangements. An individual cell may smoothly follow the motion of its neighbors (affine motion) or move in a more individualistic manner (non-affine motion). By decomposing individual motion into both affine and non-affine components, we measure cell rearrangements within a collective sheet. Finally, finite-time Lyapunov exponent (FTLE) values capture the stretching of the flow field and reflect its chaotic character. Applying collective migration analysis techniques to experimental data on both malignant and non-malignant human breast epithelial cells reveals differences in collective behavior that are not found from analyzing migration speeds alone. Non-malignant cells show increased cooperative motion on long time scales whereas malignant cells remain uncooperative as time progresses. Combining multiple analysis techniques also shows that these two cell types differ in their response to a perturbation of cell-cell adhesion through the molecule E-cadherin. Non-malignant MCF10A cells use E-cadherin for short time coordination of collective motion, yet even with decreased E-cadherin expression, the cells remain coordinated over long time scales. In contrast, the migration behavior of malignant and invasive MCF10CA1a cells, which already shows decreased collective dynamics on both time scales, is insensitive to the change in E-cadherin expression.

Dinâmica dos nutrientes na Ria Formosa: efeitos da interacçäo da laguna com as suas interfaces na reciclagem do azoto, fósforo e sílica

Tese de dout. em Ciências do Mar, Instituto de Investigação das Pescas e do Mar, Univ. do Algarve, 1996