Convecção mista em cavidades com fontes de calor aletadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Critical points on growth curves in autoregressive and mixed models

Adjusting autoregressive and mixed models to growth data fits discontinuous functions, which makes it difficult to determine critical points. In this study we propose a new approach to determine the critical stability point of cattle growth using a first-order autoregressive model and a mixed model with random asymptote, using the deterministic portion of the models. Three functions were compared: logistic, Gompertz, and Richards. The Richards autoregressive model yielded the best fit, but the critical growth values were adjusted very early, and for this purpose the Gompertz model was more appropriate.

Nonlocal dissipative tunneling for high-fidelity quantum-state transfer between distant parties

In this work, we propose the nonlocal tunneling mechanism for high-fidelity state transfer between distant parties. The nonlocal tunneling follows from the overlap between the distant sending and receiving wave functions, which is indirectlymediated by the off-resonant normal modes of a quantum channel. This channel is made up of a network of dissipative quantum systems exhibiting the same bosonic or fermionic statistical nature as the sender and receiver. We demonstrate that the incoherence arising from quantum channel nonidealities is almost completely circumvented by the tunneling mechanism, which thus affords a high-fidelity transfer process.

Microelectrode array in mixed alkanethiol self-assembled monolayers: Electrochemical studies

We report an efficient alternative to obtain recessed microelectrodes device on gold electrode surface, in which mixed self-assembled monolayer of long and short carbon alkanethiol chains was used for this purpose. Development of the modified electrodes included the chemical adsorption of 11-mercaptoundecanoic acid and 2-mercaptoethanol solution, as well as their mixtures, on gold surface, resulting in the final mixed self-assembled monolayer configuration. For comparison, the electrochemical performance of self-assembled monolayer of 11-mercaptoundecanoic acid. 3-mercaptopropionic acid, 4-mercapto-1-butanol and 6-mercapto-1-hexanol modified electrodes was also investigated. It was verified that, in the mixed self-assembled monolayer, the 11-mercaptoundecanoic acid acts as a barrier for electron transfer while the short alkanethiol chair is deposited in an island-like shape through which electrons can be freely transferred to ions in solution, allowing electrochemical reactions to occur. The performance of the modified electrodes toward microelectrode behavior was investigated via cyclic voltammetry and electrochemical impedance spectroscopy measurements using [Fe(CN)(6)](3-/4-) redox couple as a probe. In this case, sigmoidal voltammetric responses were obtained, very similar to those observed for microelectrodes. Such behavior reinforces the proposition of electron transfer through the short alkanethiol chain layer and surface blockage by the long chain one. Electrochemical impedance results allowed calculated the mean radius value of each microelectrode disks of 3.8 mu m with about 22 mu m interval between them. The microelectrode environment provided by the mixed self-assembled monolayer can be conveniently used to provide an efficient catalytic conversion in biosensing applications. (C) 2012 Elsevier Ltd. All rights reserved.

Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous cofactor participating in numerous redox reactions. It is also a substrate for regulatory modifications of proteins and nucleic acids via the addition of ADP-ribose moieties or removal of acyl groups by transfer to ADP-ribose. In this study, we use in-depth sequence, structure and genomic context analysis to uncover new enzymes and substrate-binding proteins in NAD-utilizing metabolic and macromolecular modification systems. We predict that Escherichia coli YbiA and related families of domains from diverse bacteria, eukaryotes, large DNA viruses and single strand RNA viruses are previously unrecognized components of NAD-utilizing pathways that probably operate on ADP-ribose derivatives. Using contextual analysis we show that some of these proteins potentially act in RNA repair, where NAD is used to remove 2'-3' cyclic phosphodiester linkages. Likewise, we predict that another family of YbiA-related enzymes is likely to comprise a novel NAD-dependent ADP-ribosylation system for proteins, in conjunction with a previously unrecognized ADP-ribosyltransferase. A similar ADP-ribosyltransferase is also coupled with MACRO or ADP-ribosylglycohydrolase domain proteins in other related systems, suggesting that all these novel systems are likely to comprise pairs of ADP-ribosylation and ribosylglycohydrolase enzymes analogous to the DraG-DraT system, and a novel group of bacterial polymorphic toxins. We present evidence that some of these coupled ADP-ribosyltransferases/ribosylglycohydrolases are likely to regulate certain restriction modification enzymes in bacteria. The ADP-ribosyltransferases found in these, the bacterial polymorphic toxin and host-directed toxin systems of bacteria such as Waddlia also throw light on the evolution of this fold and the origin of eukaryotic polyADP-ribosyltransferases and NEURL4-like ARTs, which might be involved in centrosomal assembly. We also infer a novel biosynthetic pathway that might be involved in the synthesis of a nicotinate-derived compound in conjunction with an asparagine synthetase and AMPylating peptide ligase. We use the data derived from this analysis to understand the origin and early evolutionary trajectories of key NAD-utilizing enzymes and present targets for future biochemical investigations.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Retrieval of trace gases vertical profile in the lower atmosphere combining. Differential Optical Absorption Spectroscopy with radiative transfer models

The motivation for the work presented in this thesis is to retrieve profile information for the atmospheric trace constituents nitrogen dioxide (NO2) and ozone (O3) in the lower troposphere from remote sensing measurements. The remote sensing technique used, referred to as Multiple AXis Differential Optical Absorption Spectroscopy (MAX-DOAS), is a recent technique that represents a significant advance on the well-established DOAS, especially for what it concerns the study of tropospheric trace consituents. NO2 is an important trace gas in the lower troposphere due to the fact that it is involved in the production of tropospheric ozone; ozone and nitrogen dioxide are key factors in determining the quality of air with consequences, for example, on human health and the growth of vegetation. To understand the NO2 and ozone chemistry in more detail not only the concentrations at ground but also the acquisition of the vertical distribution is necessary. In fact, the budget of nitrogen oxides and ozone in the atmosphere is determined both by local emissions and non-local chemical and dynamical processes (i.e. diffusion and transport at various scales) that greatly impact on their vertical and temporal distribution: thus a tool to resolve the vertical profile information is really important. Useful measurement techniques for atmospheric trace species should fulfill at least two main requirements. First, they must be sufficiently sensitive to detect the species under consideration at their ambient concentration levels. Second, they must be specific, which means that the results of the measurement of a particular species must be neither positively nor negatively influenced by any other trace species simultaneously present in the probed volume of air. Air monitoring by spectroscopic techniques has proven to be a very useful tool to fulfill these desirable requirements as well as a number of other important properties. During the last decades, many such instruments have been developed which are based on the absorption properties of the constituents in various regions of the electromagnetic spectrum, ranging from the far infrared to the ultraviolet. Among them, Differential Optical Absorption Spectroscopy (DOAS) has played an important role. DOAS is an established remote sensing technique for atmospheric trace gases probing, which identifies and quantifies the trace gases in the atmosphere taking advantage of their molecular absorption structures in the near UV and visible wavelengths of the electromagnetic spectrum (from 0.25 μm to 0.75 μm). Passive DOAS, in particular, can detect the presence of a trace gas in terms of its integrated concentration over the atmospheric path from the sun to the receiver (the so called slant column density). The receiver can be located at ground, as well as on board an aircraft or a satellite platform. Passive DOAS has, therefore, a flexible measurement configuration that allows multiple applications. The ability to properly interpret passive DOAS measurements of atmospheric constituents depends crucially on how well the optical path of light collected by the system is understood. This is because the final product of DOAS is the concentration of a particular species integrated along the path that radiation covers in the atmosphere. This path is not known a priori and can only be evaluated by Radiative Transfer Models (RTMs). These models are used to calculate the so called vertical column density of a given trace gas, which is obtained by dividing the measured slant column density to the so called air mass factor, which is used to quantify the enhancement of the light path length within the absorber layers. In the case of the standard DOAS set-up, in which radiation is collected along the vertical direction (zenith-sky DOAS), calculations of the air mass factor have been made using “simple” single scattering radiative transfer models. This configuration has its highest sensitivity in the stratosphere, in particular during twilight. This is the result of the large enhancement in stratospheric light path at dawn and dusk combined with a relatively short tropospheric path. In order to increase the sensitivity of the instrument towards tropospheric signals, measurements with the telescope pointing the horizon (offaxis DOAS) have to be performed. In this circumstances, the light path in the lower layers can become very long and necessitate the use of radiative transfer models including multiple scattering, the full treatment of atmospheric sphericity and refraction. In this thesis, a recent development in the well-established DOAS technique is described, referred to as Multiple AXis Differential Optical Absorption Spectroscopy (MAX-DOAS). The MAX-DOAS consists in the simultaneous use of several off-axis directions near the horizon: using this configuration, not only the sensitivity to tropospheric trace gases is greatly improved, but vertical profile information can also be retrieved by combining the simultaneous off-axis measurements with sophisticated RTM calculations and inversion techniques. In particular there is a need for a RTM which is capable of dealing with all the processes intervening along the light path, supporting all DOAS geometries used, and treating multiple scattering events with varying phase functions involved. To achieve these multiple goals a statistical approach based on the Monte Carlo technique should be used. A Monte Carlo RTM generates an ensemble of random photon paths between the light source and the detector, and uses these paths to reconstruct a remote sensing measurement. Within the present study, the Monte Carlo radiative transfer model PROMSAR (PROcessing of Multi-Scattered Atmospheric Radiation) has been developed and used to correctly interpret the slant column densities obtained from MAX-DOAS measurements. In order to derive the vertical concentration profile of a trace gas from its slant column measurement, the AMF is only one part in the quantitative retrieval process. One indispensable requirement is a robust approach to invert the measurements and obtain the unknown concentrations, the air mass factors being known. For this purpose, in the present thesis, we have used the Chahine relaxation method. Ground-based Multiple AXis DOAS, combined with appropriate radiative transfer models and inversion techniques, is a promising tool for atmospheric studies in the lower troposphere and boundary layer, including the retrieval of profile information with a good degree of vertical resolution. This thesis has presented an application of this powerful comprehensive tool for the study of a preserved natural Mediterranean area (the Castel Porziano Estate, located 20 km South-West of Rome) where pollution is transported from remote sources. Application of this tool in densely populated or industrial areas is beginning to look particularly fruitful and represents an important subject for future studies.

Magnetohydrodynamic Effects On Mixed Convection Flows In Channels And Ducts

This work focuses on magnetohydrodynamic (MHD) mixed convection flow of electrically conducting fluids enclosed in simple 1D and 2D geometries in steady periodic regime. In particular, in Chapter one a short overview is given about the history of MHD, with reference to papers available in literature, and a listing of some of its most common technological applications, whereas Chapter two deals with the analytical formulation of the MHD problem, starting from the fluid dynamic and energy equations and adding the effects of an external imposed magnetic field using the Ohm's law and the definition of the Lorentz force. Moreover a description of the various kinds of boundary conditions is given, with particular emphasis given to their practical realization. Chapter three, four and five describe the solution procedure of mixed convective flows with MHD effects. In all cases a uniform parallel magnetic field is supposed to be present in the whole fluid domain transverse with respect to the velocity field. The steady-periodic regime will be analyzed, where the periodicity is induced by wall temperature boundary conditions, which vary in time with a sinusoidal law. Local balance equations of momentum, energy and charge will be solved analytically and numerically using as parameters either geometrical ratios or material properties. In particular, in Chapter three the solution method for the mixed convective flow in a 1D vertical parallel channel with MHD effects is illustrated. The influence of a transverse magnetic field will be studied in the steady periodic regime induced by an oscillating wall temperature. Analytical and numerical solutions will be provided in terms of velocity and temperature profiles, wall friction factors and average heat fluxes for several values of the governing parameters. In Chapter four the 2D problem of the mixed convective flow in a vertical round pipe with MHD effects is analyzed. Again, a transverse magnetic field influences the steady periodic regime induced by the oscillating wall temperature of the wall. A numerical solution is presented, obtained using a finite element approach, and as a result velocity and temperature profiles, wall friction factors and average heat fluxes are derived for several values of the Hartmann and Prandtl numbers. In Chapter five the 2D problem of the mixed convective flow in a vertical rectangular duct with MHD effects is discussed. As seen in the previous chapters, a transverse magnetic field influences the steady periodic regime induced by the oscillating wall temperature of the four walls. The numerical solution obtained using a finite element approach is presented, and a collection of results, including velocity and temperature profiles, wall friction factors and average heat fluxes, is provided for several values of, among other parameters, the duct aspect ratio. A comparison with analytical solutions is also provided, as a proof of the validity of the numerical method. Chapter six is the concluding chapter, where some reflections on the MHD effects on mixed convection flow will be made, in agreement with the experience and the results gathered in the analyses presented in the previous chapters. In the appendices special auxiliary functions and FORTRAN program listings are reported, to support the formulations used in the solution chapters.

Atmospheric radiative transfer in multiple scattering conditions. Application to NWP models.

High spectral resolution radiative transfer (RT) codes are essential tools in the study of the radiative energy transfer in the Earth atmosphere and a support for the development of parameterizations for fast RT codes used in climate and weather prediction models. Cirrus clouds cover permanently 30% of the Earth's surface, representing an important contribution to the Earth-atmosphere radiation balance. The work has been focussed on the development of the RT model LBLMS. The model, widely tested in the infra-red spectral range, has been extended to the short wave spectrum and it has been used in comparison with airborne and satellite measurements to study the optical properties of cirrus clouds. A new database of single scattering properties has been developed for mid latitude cirrus clouds. Ice clouds are treated as a mixture of ice crystals with various habits. The optical properties of the mixture are tested in comparison to radiometric measurements in selected case studies. Finally, a parameterization of the mixture for application to weather prediction and global circulation models has been developed. The bulk optical properties of ice crystals are parameterized as functions of the effective dimension of measured particle size distributions that are representative of mid latitude cirrus clouds. Tests with the Limited Area Weather Prediction model COSMO have shown the impact of the new parameterization with respect to cirrus cloud optical properties based on ice spheres.

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Factorization method and inclusions of mixed type in an inverse elliptic boundary value problem

In various imaging problems the task is to use the Cauchy data of the solutions to an elliptic boundary value problem to reconstruct the coefficients of the corresponding partial differential equation. Often the examined object has known background properties but is contaminated by inhomogeneities that cause perturbations of the coefficient functions. The factorization method of Kirsch provides a tool for locating such inclusions. In this paper, the factorization technique is studied in the framework of coercive elliptic partial differential equations of the divergence type: Earlier it has been demonstrated that the factorization algorithm can reconstruct the support of a strictly positive (or negative) definite perturbation of the leading order coefficient, or if that remains unperturbed, the support of a strictly positive (or negative) perturbation of the zeroth order coefficient. In this work we show that these two types of inhomogeneities can, in fact, be located simultaneously. Unlike in the earlier articles on the factorization method, our inclusions may have disconnected complements and we also weaken some other a priori assumptions of the method. Our theoretical findings are complemented by two-dimensional numerical experiments that are presented in the framework of the diffusion approximation of optical tomography.

Study of the mechanisms of shear transfer in fiber-reinforced concrete elements

matlab functions for the validation of push-off tests results

The impact of ice crystals on radiative forcing and remote sensing of arctic boundary-layer mixed-phase clouds

This PhD thesis is embedded into the Arctic Study of Tropospheric Aerosol, Clouds and Radiation (ASTAR) and investigates the radiative transfer through Arctic boundary-layer mixed-phase (ABM) clouds. For this purpose airborne spectral solar radiation measurements and simulations of the solar and thermal infrared radiative transfer have been performed. This work reports on measurements with the Spectral Modular Airborne Radiation measurement sysTem (SMART-Albedometer) conducted in the framework of ASTAR in April 2007 close to Svalbard. For ASTAR the SMART-Albedometer was extended to measure spectral radiance. The development and calibration of the radiance measurements are described in this work. In combination with in situ measurements of cloud particle properties provided by the Laboratoire de M¶et¶eorologie Physique (LaMP) and simultaneous airborne lidar measurements by the Alfred Wegener Institute for Polar and Marine Research (AWI) ABM clouds were sampled. The SMART-Albedometer measurements were used to retrieve the cloud thermodynamic phase by three different approaches. A comparison of these results with the in situ and lidar measurements is presented in two case studies. Beside the dominating mixed-phase clouds pure ice clouds were found in cloud gaps and at the edge of a large cloud field. Furthermore the vertical distribution of ice crystals within ABM clouds was investigated. It was found that ice crystals at cloud top are necessary to describe the observed SMART-Albedometer measurements. The impact of ice crystals on the radiative forcing of ABM clouds is in vestigated by extensive radiative transfer simulations. The solar and net radiative forcing was found to depend on the ice crystal size, shape and the mixing ratio of ice crystals and liquid water droplets.

3D-ultrastructure, functions and stress responses of rhogocytes from the freshwater gastropods Biomphalaria glabrata and Lymnaea stagnalis

Rhogocytes, also termed ‘pore cells’, exist free in the hemolymph or embedded in the connective tissue of different body parts of molluscs, notably gastropods. These unique cells can be round, elongated or irregularly shaped, and up to 30 μm in diameter. Their hallmark is the so-called slit apparatus: i.e. pocket-like invaginations of the plasma membrane creating extracellular lacunae, bridged by cytoplasmic bars. These bars form distinctive slits of ca. 20 nm width. A slit diaphragm composed of proteins establishes a molecular sieve with holes of 20 x 20 nm. Different functions have been assigned to this special molluscan cell type, notably biosynthesis of the hemolymph respiratory protein hemocyanin. It has further been proposed, but not proven, that in the case of red-blooded snail species rhogocytes might synthesize the hemoglobin. However, the secretion pathway of these hemolymph proteins, and the functional role of the enigmatic slit apparatus remained unclear. Additionally proposed functions of rhogocytes, such as heavy metal detoxification or hemolymph protein degradation, are also not well studied. This work provides more detailed electron microscopical, histological and immunobiochemical information on the structure and function of rhogocytes of the freshwater snails Biomphalaria glabrata and Lymnaea stagnalis. By in situ hybridization on mantle tissues, it proves that B. glabrata rhogocytes synthesize hemoglobin and L. stagnalis rhogocytes synthesize hemocyanin. Hemocyanin is present, in endoplasmic reticulum lacunae and in vesicles, as individual molecules or pseudo-crystalline arrays. The first 3D reconstructions of rhogocytes are provided by means of electron tomography and show unprecedented details of the slit apparatus. A highly dense material in the cytoplasmic bars close to the diaphragmatic slits was shown, by immunogold labeling, to contain actin. By immunofluorescence microscopy, the protein nephrin was localized at the periphery of rhogocytes. The presence of both proteins in the slit apparatus supports the previous hypothesis, hitherto solely based on similarities of the ultrastructure, that the molluscan rhogocytes are phylogenetically related to mammalian podocytes and insect nephrocytes. A possible secretion pathway of respiratory proteins that includes a transfer mechanism of vesicles through the diaphragmatic slits is proposed and discussed. We also studied, by electron microscopy, the reaction of rhogocytes in situ to two forms of animal stress: deprivation of food and cadmium contamination of the tank water. Significant cellular reactions to both stressors were observed and documented. Notably, the slit apparatus surface and the number of electron-dense cytoplasmic vesicles increased in response to cadmium stress. Food deprivation led to an increase in hemocyanin production. These observations are also discussed in the framework of using such animals as potential environmental biomarkers.