Applications of finite element computer modelling and thermal infrared remote sensing to the study of geothermal anomalies

Buried heat sources can be investigated by examining thermal infrared images and comparing these with the results of theoretical models which predict the thermal anomaly a given heat source may generate. Key factors influencing surface temperature include the geometry and temperature of the heat source, the surface meteorological environment, and the thermal conductivity and anisotropy of the rock. In general, a geothermal heat flux of greater than 2% of solar insolation is required to produce a detectable thermal anomaly in a thermal infrared image. A heat source of, for example, 2-300K greater than the average surface temperature must be a t depth shallower than 50m for the detection of the anomaly in a thermal infrared image, for typical terrestrial conditions. Atmospheric factors are of critical importance. While the mean atmospheric temperature has little significance, the convection is a dominant factor, and can act to swamp the thermal signature entirely. Given a steady state heat source that produces a detectable thermal anomaly, it is possible to loosely constrain the physical properties of the heat source and surrounding rock, using the surface thermal anomaly as a basis. The success of this technique is highly dependent on the degree to which the physical properties of the host rock are known. Important parameters include the surface thermal properties and thermal conductivity of the rock. Modelling of transient thermal situations was carried out, to assess the effect of time dependant thermal fluxes. One-dimensional finite element models can be readily and accurately applied to the investigation of diurnal heat flow, as with thermal inertia models. Diurnal thermal models of environments on Earth, the Moon and Mars were carried out using finite elements and found to be consistent with published measurements. The heat flow from an injection of hot lava into a near surface lava tube was considered. While this approach was useful for study, and long term monitoring in inhospitable areas, it was found to have little hazard warning utility, as the time taken for the thermal energy to propagate to the surface in dry rock (several months) in very long. The resolution of the thermal infrared imaging system is an important factor. Presently available satellite based systems such as Landsat (resolution of 120m) are inadequate for detailed study of geothermal anomalies. Airborne systems, such as TIMS (variable resolution of 3-6m) are much more useful for discriminating small buried heat sources. Planned improvements in the resolution of satellite based systems will broaden the potential for application of the techniques developed in this thesis. It is important to note, however, that adequate spatial resolution is a necessary but not sufficient condition for successful application of these techniques.

Fabrication, characterisation and electroanalysis at 1-dimensional nanostructures

Integrated nanowire electrodes that permit direct, sensitive and rapid electrochemical based detection of chemical and biological species are a powerful emerging class of sensor devices. As critical dimensions of the electrodes enter the nanoscale, radial analyte diffusion profiles to the electrode dominate with a corresponding enhancement in mass transport, steady-state sigmoidal voltammograms, low depletion of target molecules and faster analysis. To optimise these sensors it is necessary to fully understand the factors that influence performance limits including: electrode geometry, electrode dimensions, electrode separation distances (within nanowire arrays) and diffusional mass transport. Therefore, in this thesis, theoretical simulations of analyte diffusion occurring at a variety of electrode designs were undertaken using Comsol Multiphysics®. Sensor devices were fabricated and corresponding experiments were performed to challenge simulation results. Two approaches for the fabrication and integration of metal nanowire electrodes are presented: Template Electrodeposition and Electron-Beam Lithography. These approaches allow for the fabrication of nanowires which may be subsequently integrated at silicon chip substrates to form fully functional electrochemical devices. Simulated and experimental results were found to be in excellent agreement validating the simulation model. The electrochemical characteristics exhibited by nanowire electrodes fabricated by electronbeam lithography were directly compared against electrochemical performance of a commercial ultra-microdisc electrode. Steady-state cyclic voltammograms in ferrocenemonocarboxylic acid at single ultra-microdisc electrodes were observed at low to medium scan rates (≤ 500 mV.s-1). At nanowires, steady-state responses were observed at ultra-high scan rates (up to 50,000 mV.s-1), thus allowing for much faster analysis (20 ms). Approaches for elucidating faradaic signal without the requirement for background subtraction were also developed. Furthermore, diffusional process occurring at arrays with increasing inter-electrode distance and increasing number of nanowires were explored. Diffusion profiles existing at nanowire arrays were simulated with Comsol Multiphysics®. A range of scan rates were modelled, and experiments were undertaken at 5,000 mV.s-1 since this allows rapid data capture required for, e.g., biomedical, environmental and pharmaceutical diagnostic applications.

Interaction of additive and quantization noises in digital PLLs

Phase-locked loops (PLLs) are a crucial component in modern communications systems. Comprising of a phase-detector, linear filter, and controllable oscillator, they are widely used in radio receivers to retrieve the information content from remote signals. As such, they are capable of signal demodulation, phase and carrier recovery, frequency synthesis, and clock synchronization. Continuous-time PLLs are a mature area of study, and have been covered in the literature since the early classical work by Viterbi [1] in the 1950s. With the rise of computing in recent decades, discrete-time digital PLLs (DPLLs) are a more recent discipline; most of the literature published dates from the 1990s onwards. Gardner [2] is a pioneer in this area. It is our aim in this work to address the difficulties encountered by Gardner [3] in his investigation of the DPLL output phase-jitter where additive noise to the input signal is combined with frequency quantization in the local oscillator. The model we use in our novel analysis of the system is also applicable to another of the cases looked at by Gardner, that is the DPLL with a delay element integrated in the loop. This gives us the opportunity to look at this system in more detail, our analysis providing some unique insights into the variance `dip' seen by Gardner in [3]. We initially provide background on the probability theory and stochastic processes. These branches of mathematics are the basis for the study of noisy analogue and digital PLLs. We give an overview of the classical analogue PLL theory as well as the background on both the digital PLL and circle map, referencing the model proposed by Teplinsky et al. [4, 5]. For our novel work, the case of the combined frequency quantization and noisy input from [3] is investigated first numerically, and then analytically as a Markov chain via its Chapman-Kolmogorov equation. The resulting delay equation for the steady-state jitter distribution is treated using two separate asymptotic analyses to obtain approximate solutions. It is shown how the variance obtained in each case matches well to the numerical results. Other properties of the output jitter, such as the mean, are also investigated. In this way, we arrive at a more complete understanding of the interaction between quantization and input noise in the first order DPLL than is possible using simulation alone. We also do an asymptotic analysis of a particular case of the noisy first-order DPLL with delay, previously investigated by Gardner [3]. We show a unique feature of the simulation results, namely the variance `dip' seen for certain levels of input noise, is explained by this analysis. Finally, we look at the second-order DPLL with additive noise, using numerical simulations to see the effects of low levels of noise on the limit cycles. We show how these effects are similar to those seen in the noise-free loop with non-zero initial conditions.

Design of digital differentiators

A digital differentiator simply involves the derivation of an input signal. This work includes the presentation of first-degree and second-degree differentiators, which are designed as both infinite-impulse-response (IIR) filters and finite-impulse-response (FIR) filters. The proposed differentiators have low-pass magnitude response characteristics, thereby rejecting noise frequencies higher than the cut-off frequency. Both steady-state frequency-domain characteristics and Time-domain analyses are given for the proposed differentiators. It is shown that the proposed differentiators perform well when compared to previously proposed filters. When considering the time-domain characteristics of the differentiators, the processing of quantized signals proved especially enlightening, in terms of the filtering effects of the proposed differentiators. The coefficients of the proposed differentiators are obtained using an optimization algorithm, while the optimization objectives include magnitude and phase response. The low-pass characteristic of the proposed differentiators is achieved by minimizing the filter variance. The low-pass differentiators designed show the steep roll-off, as well as having highly accurate magnitude response in the pass-band. While having a history of over three hundred years, the design of fractional differentiator has become a ‘hot topic’ in recent decades. One challenging problem in this area is that there are many different definitions to describe the fractional model, such as the Riemann-Liouville and Caputo definitions. Through use of a feedback structure, based on the Riemann-Liouville definition. It is shown that the performance of the fractional differentiator can be improved in both the frequency-domain and time-domain. Two applications based on the proposed differentiators are described in the thesis. Specifically, the first of these involves the application of second degree differentiators in the estimation of the frequency components of a power system. The second example concerns for an image processing, edge detection application.

Quantitative analysis of the effect of derivatisation of [Ru(BPY) phen] with a quinoline moiety on the interaction with DNA

The bifunctional Ru(II) complex [Ru(BPY)2POQ-Nmet]2+ (1), in which the metallic unit is tethered by an aliphatic chain to an organic DNA binder, was designed in order to increase the affinity toward nucleic acids. The interaction of 1 with DNA was characterised from luminescence and absorption data and compared with the binding of its monofunctional metallic and organic analogues, [Ru(BPY)2(ac)phen]2+ (2) and Nmet-quinoline (3). The bifunctional complex has a binding affinity one order of magnitude higher than that of each of its separated moieties. Absorption changes induced upon addition of DNA at different pH indicate protonation of the organic sub-unit upon interaction with DNA under neutral conditions. The combination of the luminescence data under steady-state and time-resolved conditions shows that the attachment of the organic unit in 1 induces modifications of the association modes of the metallic unit, owing to the presence of the aliphatic chain which probably hinders the metallic moiety binding. The salt dependence of the binding constants was analysed in order to compare the thermodynamic parameters describing the association with DNA for each complex. This study demonstrates the interest of the derivatisation of a Ru(II) complex with an organic moiety (ia the bifunctional ligand POQ-Nmet) for the development of high affinity DNA probes or photoreactive agents.

Typical delay determines waiting time on periodic-food schedules: Static and dynamic tests.

Pigeons and other animals soon learn to wait (pause) after food delivery on periodic-food schedules before resuming the food-rewarded response. Under most conditions the steady-state duration of the average waiting time, t, is a linear function of the typical interfood interval. We describe three experiments designed to explore the limits of this process. In all experiments, t was associated with one key color and the subsequent food delay, T, with another. In the first experiment, we compared the relation between t (waiting time) and T (food delay) under two conditions: when T was held constant, and when T was an inverse function of t. The pigeons could maximize the rate of food delivery under the first condition by setting t to a consistently short value; optimal behavior under the second condition required a linear relation with unit slope between t and T. Despite this difference in optimal policy, the pigeons in both cases showed the same linear relation, with slope less than one, between t and T. This result was confirmed in a second parametric experiment that added a third condition, in which T + t was held constant. Linear waiting appears to be an obligatory rule for pigeons. In a third experiment we arranged for a multiplicative relation between t and T (positive feedback), and produced either very short or very long waiting times as predicted by a quasi-dynamic model in which waiting time is strongly determined by the just-preceding food delay.

Phase transfer catalysts drive diverse organic solvent solubility of single-walled carbon nanotubes helically wrapped by ionic, semiconducting polymers.

Use of phase transfer catalysts such as 18-crown-6 enables ionic, linear conjugated poly[2,6-{1,5-bis(3-propoxysulfonicacidsodiumsalt)}naphthylene]ethynylene (PNES) to efficiently disperse single-walled carbon nanotubes (SWNTs) in multiple organic solvents under standard ultrasonication methods. Steady-state electronic absorption spectroscopy, atomic force microscopy (AFM), and transmission electron microscopy (TEM) reveal that these SWNT suspensions are composed almost exclusively of individualized tubes. High-resolution TEM and AFM data show that the interaction of PNES with SWNTs in both protic and aprotic organic solvents provides a self-assembled superstructure in which a PNES monolayer helically wraps the nanotube surface with periodic and constant morphology (observed helical pitch length = 10 ± 2 nm); time-dependent examination of these suspensions indicates that these structures persist in solution over periods that span at least several months. Pump-probe transient absorption spectroscopy reveals that the excited state lifetimes and exciton binding energies of these well-defined nanotube-semiconducting polymer hybrid structures remain unchanged relative to analogous benchmark data acquired previously for standard sodium dodecylsulfate (SDS)-SWNT suspensions, regardless of solvent. These results demonstrate that the use of phase transfer catalysts with ionic semiconducting polymers that helically wrap SWNTs provide well-defined structures that solubulize SWNTs in a wide range of organic solvents while preserving critical nanotube semiconducting and conducting properties.

The role of tectonic uplift, climate, and vegetation in the long-term terrestrial phosphorous cycle

Phosphorus (P) is a crucial element for life and therefore for maintaining ecosystem productivity. Its local availability to the terrestrial biosphere results from the interaction between climate, tectonic uplift, atmospheric transport, and biotic cycling. Here we present a mathematical model that describes the terrestrial P-cycle in a simple but comprehensive way. The resulting dynamical system can be solved analytically for steady-state conditions, allowing us to test the sensitivity of the P-availability to the key parameters and processes. Given constant inputs, we find that humid ecosystems exhibit lower P availability due to higher runoff and losses, and that tectonic uplift is a fundamental constraint. In particular, we find that in humid ecosystems the biotic cycling seem essential to maintain long-term P-availability. The time-dependent P dynamics for the Franz Josef and Hawaii chronosequences show how tectonic uplift is an important constraint on ecosystem productivity, while hydroclimatic conditions control the P-losses and speed towards steady-state. The model also helps describe how, with limited uplift and atmospheric input, as in the case of the Amazon Basin, ecosystems must rely on mechanisms that enhance P-availability and retention. Our novel model has a limited number of parameters and can be easily integrated into global climate models to provide a representation of the response of the terrestrial biosphere to global change. © 2010 Author(s).

Portable, Fiber-Based, Diffuse Reflection Spectroscopy (DRS) Systems for Estimating Tissue Optical Properties.

Steady-state diffuse reflection spectroscopy is a well-studied optical technique that can provide a noninvasive and quantitative method for characterizing the absorption and scattering properties of biological tissues. Here, we compare three fiber-based diffuse reflection spectroscopy systems that were assembled to create a light-weight, portable, and robust optical spectrometer that could be easily translated for repeated and reliable use in mobile settings. The three systems were built using a broadband light source and a compact, commercially available spectrograph. We tested two different light sources and two spectrographs (manufactured by two different vendors). The assembled systems were characterized by their signal-to-noise ratios, the source-intensity drifts, and detector linearity. We quantified the performance of these instruments in extracting optical properties from diffuse reflectance spectra in tissue-mimicking liquid phantoms with well-controlled optical absorption and scattering coefficients. We show that all assembled systems were able to extract the optical absorption and scattering properties with errors less than 10%, while providing greater than ten-fold decrease in footprint and cost (relative to a previously well-characterized and widely used commercial system). Finally, we demonstrate the use of these small systems to measure optical biomarkers in vivo in a small-animal model cancer therapy study. We show that optical measurements from the simple portable system provide estimates of tumor oxygen saturation similar to those detected using the commercial system in murine tumor models of head and neck cancer.

Tyrosine phosphorylation of G protein alpha subunits by pp60c-src.

A number of lines of evidence suggest that cross-talk exists between the cellular signal transduction pathways involving tyrosine phosphorylation catalyzed by members of the pp60c-src kinase family and those mediated by guanine nucleotide regulatory proteins (G proteins). In this study, we explore the possibility that direct interactions between pp60c-src and G proteins may occur with functional consequences. Preparations of pp60c-src isolated by immunoprecipitation phosphorylate on tyrosine residues the purified G-protein alpha subunits (G alpha) of several heterotrimeric G proteins. Phosphorylation is highly dependent on G-protein conformation, and G alpha(GDP) uncomplexed by beta gamma subunits appears to be the preferred substrate. In functional studies, phosphorylation of stimulatory G alpha (G alpha s) modestly increases the rate of binding of guanosine 5'-[gamma-[35S]thio]triphosphate to Gs as well as the receptor-stimulated steady-state rate of GTP hydrolysis by Gs. Heterotrimeric G proteins may represent a previously unappreciated class of potential substrates for pp60c-src.

cAMP stimulates transcription of the beta 2-adrenergic receptor gene in response to short-term agonist exposure.

In addition to conveying cellular responses to an effector molecule, receptors are often themselves regulated by their effectors. We have demonstrated that epinephrine modulates both the rate of transcription of the beta 2-adrenergic receptor (beta 2AR) gene and the steady-state level of beta 2AR mRNA in DDT1MF-2 cells. Short-term (30 min) exposure to epinephrine (100 nM) stimulates the rate of beta 2AR gene transcription, resulting in a 3- to 4-fold increase in steady-state beta 2AR mRNA levels. These effects are mimicked by 1 mM N6,O2'-dibutyryladenosine 3',5'-cyclic monophosphate (Bt2cAMP) or foskolin but not by phorbol esters. The half-life of the beta 2AR mRNA after addition of actinomycin D (46.7 +/- 10.2 min; mean +/- SEM; n = 5) remained unchanged after 30 min of epinephrine treatment (46.8 +/- 10.6 min; mean +/- SEM; n = 4), indicating that a change in transcription rate is the predominant factor responsible for the increase of beta 2AR mRNA. Whereas brief exposure to epinephrine or Bt2cAMP does not significantly affect the total number of cellular beta 2ARs (assessed by ligand binding), continued exposure results in a gradual decline in beta 2AR number to approximately 20% (epinephrine) or approximately 45% (Bt2cAMP) of the levels in control cells by 24 hr. Similar decreases in agonist-stimulated adenylyl cyclase activity are observed. This loss of receptors with prolonged agonist exposure is accompanied by a 50% reduction in beta 2AR mRNA. Transfection of the beta 2AR promoter region cloned onto a reporter gene (bacterial chloramphenicol acetyltransferase) allowed demonstration of a 2- to 4-fold induction of transcription by agents that elevate cAMP levels, such as forskolin or phosphodiesterase inhibitors. These results establish the presence of elements within the proximal promoter region of the beta 2AR gene responsible for the transcriptional enhancing activity of cAMP and demonstrate that beta 2AR gene expression is regulated by a type of feedback mechanism involving the second messenger cAMP.

Autophagy enhances NFκB activity in specific tissue macrophages by sequestering A20 to boost antifungal immunity.

Immune responses must be well restrained in a steady state to avoid excessive inflammation. However, such restraints are quickly removed to exert antimicrobial responses. Here we report a role of autophagy in an early host antifungal response by enhancing NFκB activity through A20 sequestration. Enhancement of NFκB activation is achieved by autophagic depletion of A20, an NFκB inhibitor, in F4/80(hi) macrophages in the spleen, peritoneum and kidney. We show that p62, an autophagic adaptor protein, captures A20 to sequester it in the autophagosome. This allows the macrophages to release chemokines to recruit neutrophils. Indeed, mice lacking autophagy in myeloid cells show higher susceptibility to Candida albicans infection due to impairment in neutrophil recruitment. Thus, at least in the specific aforementioned tissues, autophagy appears to break A20-dependent suppression in F4/80(hi) macrophages, which express abundant A20 and contribute to the initiation of efficient innate immune responses.

An evaluation of remifentanil-sevoflurane response surface models in patients emerging from anesthesia: model improvement using effect-site sevoflurane concentrations.

INTRODUCTION: We previously reported models that characterized the synergistic interaction between remifentanil and sevoflurane in blunting responses to verbal and painful stimuli. This preliminary study evaluated the ability of these models to predict a return of responsiveness during emergence from anesthesia and a response to tibial pressure when patients required analgesics in the recovery room. We hypothesized that model predictions would be consistent with observed responses. We also hypothesized that under non-steady-state conditions, accounting for the lag time between sevoflurane effect-site concentration (Ce) and end-tidal (ET) concentration would improve predictions. METHODS: Twenty patients received a sevoflurane, remifentanil, and fentanyl anesthetic. Two model predictions of responsiveness were recorded at emergence: an ET-based and a Ce-based prediction. Similarly, 2 predictions of a response to noxious stimuli were recorded when patients first required analgesics in the recovery room. Model predictions were compared with observations with graphical and temporal analyses. RESULTS: While patients were anesthetized, model predictions indicated a high likelihood that patients would be unresponsive (> or = 99%). However, after termination of the anesthetic, models exhibited a wide range of predictions at emergence (1%-97%). Although wide, the Ce-based predictions of responsiveness were better distributed over a percentage ranking of observations than the ET-based predictions. For the ET-based model, 45% of the patients awoke within 2 min of the 50% model predicted probability of unresponsiveness and 65% awoke within 4 min. For the Ce-based model, 45% of the patients awoke within 1 min of the 50% model predicted probability of unresponsiveness and 85% awoke within 3.2 min. Predictions of a response to a painful stimulus in the recovery room were similar for the Ce- and ET-based models. DISCUSSION: Results confirmed, in part, our study hypothesis; accounting for the lag time between Ce and ET sevoflurane concentrations improved model predictions of responsiveness but had no effect on predicting a response to a noxious stimulus in the recovery room. These models may be useful in predicting events of clinical interest but large-scale evaluations with numerous patients are needed to better characterize model performance.

Local kinetic interpretation of entropy production through reversed diffusion.

The time reversal of stochastic diffusion processes is revisited with emphasis on the physical meaning of the time-reversed drift and the noise prescription in the case of multiplicative noise. The local kinematics and mechanics of free diffusion are linked to the hydrodynamic description. These properties also provide an interpretation of the Pope-Ching formula for the steady-state probability density function along with a geometric interpretation of the fluctuation-dissipation relation. Finally, the statistics of the local entropy production rate of diffusion are discussed in the light of local diffusion properties, and a stochastic differential equation for entropy production is obtained using the Girsanov theorem for reversed diffusion. The results are illustrated for the Ornstein-Uhlenbeck process.

In vitro fluid dynamics of the Ahmed glaucoma valve modified with expanded polytetrafluoroethylene.

PURPOSE: Long-term intraocular pressure reduction by glaucoma drainage devices (GDDs) is often limited by the fibrotic capsule that forms around them. Prior work demonstrates that modifying a GDD with a porous membrane promotes a vascularized and more permeable capsule. This work examines the in vitro fluid dynamics of the Ahmed valve after enclosing the outflow tract with a porous membrane of expanded polytetrafluoroethylene (ePTFE). MATERIALS AND METHODS: The control and modified Ahmed implants (termed porous retrofitted implant with modified enclosure or PRIME-Ahmed) were submerged in saline and gelatin and perfused in a system that monitored flow (Q) and pressure (P). Flow rates of 1-50 μl/min were applied and steady state pressure recorded. Resistance was calculated by dividing pressure by flow. RESULTS: Modifying the Ahmed valve implant outflow with expanded ePTFE increased pressure and resistance. Pressure at a flow of 2 μl/min was increased in the PRIME-Ahmed (11.6 ± 1.5 mm Hg) relative to the control implant (6.5 ± 1.2 mm Hg). Resistance at a flow of 2 μl/min was increased in the PRIME-Ahmed (5.8 ± 0.8 mm Hg/μl/min) when compared to the control implant (3.2 ± 0.6 mm Hg/μl/min). CONCLUSIONS: Modifying the outflow tract of the Ahmed valve with a porous membrane adds resistance that decreases with increasing flow. The Ahmed valve implant behaves as a variable resistor. It is partially open at low pressures and provides reduced resistance at physiologic flow rates.