Semiparametric Latent Variable Regression Models for Spatio-temporal Modeling of Mobile Source Particles in the Greater Boston Area

Traffic particle concentrations show considerable spatial variability within a metropolitan area. We consider latent variable semiparametric regression models for modeling the spatial and temporal variability of black carbon and elemental carbon concentrations in the greater Boston area. Measurements of these pollutants, which are markers of traffic particles, were obtained from several individual exposure studies conducted at specific household locations as well as 15 ambient monitoring sites in the city. The models allow for both flexible, nonlinear effects of covariates and for unexplained spatial and temporal variability in exposure. In addition, the different individual exposure studies recorded different surrogates of traffic particles, with some recording only outdoor concentrations of black or elemental carbon, some recording indoor concentrations of black carbon, and others recording both indoor and outdoor concentrations of black carbon. A joint model for outdoor and indoor exposure that specifies a spatially varying latent variable provides greater spatial coverage in the area of interest. We propose a penalised spline formation of the model that relates to generalised kringing of the latent traffic pollution variable and leads to a natural Bayesian Markov Chain Monte Carlo algorithm for model fitting. We propose methods that allow us to control the degress of freedom of the smoother in a Bayesian framework. Finally, we present results from an analysis that applies the model to data from summer and winter separately

U/Th systematics of authigenic carbonates and d18O records from Hydrate Ridge, off the coast of Oregon, USA

Uranium (U) concentrations and activity ratios (d234U) of authigenic carbonates are sensitive recorders of different fluid compositions at submarine seeps of hydrocarbon-rich fluids ("cold seeps") at Hydrate Ridge, off the coast of Oregon, USA. The low U concentrations (mean: 1.3 ± 0.4 µg/g) and high 234U values (165-317 per mil) of gas hydrate carbonates reflect the influence of sedimentary pore water indicating that these carbonates were formed under reducing conditions below or at the seafloor. Their 230Th/234U ages span a time interval from 0.8 to 6.4 ka and cluster around 1.2 and 4.7 ka. In contrast, chemoherm carbonates precipitate from marine bottom water marked by relatively high U concentrations (mean: 5.2 ± 0.8 µg/g) and a mean d234U ratio of 166 ± 3 per mil. Their U isotopes reflect the d234U ratios of the bottom water being enriched in 234U relative to normal seawater. Simple mass balance calculations based on U concentrations and their corresponding d234U ratios reveal a contribution of about 11% of sedimentary pore water to the bottom water. From the U pore water flux and the reconstructed U pore water concentration a mean flow rate of about 147 ± 68 cm/a can be estimated. 230Th/234U ages of chemoherm carbonates range from 7.3 to 267.6 ka. 230Th/234U ages of two chemoherms (Alvin and SE-Knoll chemoherm) correspond to time intervals of low sealevel stands in marine isotope stages (MIS) 2, 4, 5, 6, 7 and 8. This observation indicates that fluid flow at cold seep sites sensitively reflects pressure changes of the hydraulic head in the sediments. The d18OPDB ratios of the chemoherm carbonates support the hypothesis of precipitation during glacial times. Deviations of the chemoherm d18O values from the marine d18O record can be interpreted as to reflect temporally and spatially varying bottom water and/or vent fluid temperatures during carbonate precipitation between 2.6 and 8.6°C.

Seasonally resolved chemical records of six ice cores from the Dronning Maud Land, Antarctica

The ammonium, calcium, and sodium concentrations from three intermediate depth ice cores drilled in the area of Dronning Maud Land, East Antarctica, have been investigated. Since all measurements were performed by a high-resolution Continuous Flow Analysis system, for the first time seasonal signals of chemical trace species could be obtained from the interior of central Antarctica over a period of approximately 2 millennia. Although the elevation as well as the accumulation rate differ between the drilling sites, similar values were obtained by comparing mean concentrations spanning the last 900 years. However, a distinct lack of intersite correlation was found on decadal timescales. Despite a noticeable accumulation change, apparent in one core, no significant concentration change of all three species has occurred. All the measured ions show clear seasonal signals over the whole records. While the sea-salt-related component sodium peaks simultaneously with calcium, the maximum ammonium concentration occurs in the snow with a time lag of 2 months after the sea-salt peak. More than 60% of the calcium concentration can be attributed to an ocean source. Elevated sodium concentrations were found within this millennium compared to mean values of the whole records, but the spatially varying shape of the increase suggests that a possible climatic signal is biased by local deposition effects.

Solar radiation over and under sea ice at ROV station PS86/081-1 during POLARSTERN cruise PS86 (AURORA) in July 2014

The observed changes in physical properties of sea ice such as decreased thickness and increased melt pond cover severely impact the energy budget of Arctic sea ice. Increased light transmission leads to increased deposition of solar energy in the upper ocean and thus plays a crucial role for amount and timing of sea-ice-melt and under-ice primary production. Recent developments in underwater technology provide new opportunities to study light transmission below the largely inaccessible underside of sea ice. We measured spectral under-ice radiance and irradiance using the new Nereid Under-Ice (NUI) underwater robotic vehicle, during a cruise of the R/V Polarstern to 83°N 6°W in the Arctic Ocean in July 2014. NUI is a next generation hybrid remotely operated vehicle (H-ROV) designed for both remotely piloted and autonomous surveys underneath land-fast and moving sea ice. Here we present results from one of the first comprehensive scientific dives of NUI employing its interdisciplinary sensor suite. We combine under-ice optical measurements with three dimensional under-ice topography (multibeam sonar) and aerial images of the surface conditions. We investigate the influence of spatially varying ice-thickness and surface properties on the spatial variability of light transmittance during summer. Our results show that surface properties such as melt ponds dominate the spatial distribution of the under-ice light field on small scales (<1000 m**2), while sea ice-thickness is the most important predictor for light transmission on larger scales. In addition, we propose the use of an algorithm to obtain histograms of light transmission from distributions of sea ice thickness and surface albedo.

Recurrent excitatory postsynaptic potentials induced by synchronized fast cortical oscillations

Gamma frequency (about 20–70 Hz) oscillations occur during novel sensory stimulation, with tight synchrony over distances of at least 7 mm. Synchronization in the visual system has been proposed to reflect coactivation of different parts of the visual field by a single spatially extended object. We have shown that intracortical mechanisms, including spike doublet firing by interneurons, can account for tight long-range synchrony. Here we show that synchronous gamma oscillations in two sites also can cause long-lasting (>1 hr) potentiation of recurrent excitatory synapses. Synchronous oscillations lasting >400 ms in hippocampal area CA1 are associated with an increase in both excitatory postsynaptic potential (EPSP) amplitude and action potential afterhyperpolarization size. The resulting EPSPs stabilize and synchronize a prolonged beta frequency (about 10–25 Hz) oscillation. The changes in EPSP size are not expressed during non-oscillatory behavior but reappear during subsequent gamma-oscillatory events. We propose that oscillation-induced EPSPs serve as a substrate for memory, whose expression either enhances or blocks synchronization of spatially separated sites. This phenomenon thus provides a dynamical mechanism for storage and retrieval of stimulus-specific neuronal assemblies.

Recursive filtering of images with symmetric extension

Recursive filters are widely used in image analysis due to their efficiency and simple implementation. However these filters have an initialisation problem which either produces unusable results near the image boundaries or requires costly approximate solutions such as extending the boundary manually. In this paper, we describe a method for the recursive filtering of symmetrically extended images for filters with symmetric denominator. We begin with an analysis of symmetric extensions and their effect on non-recursive filtering operators. Based on the non-recursive case, we derive a formulation of recursive filtering on symmetric domains as a linear but spatially varying implicit operator. We then give an efficient method for decomposing and solving the linear implicit system, along with a proof that this decomposition always exists. This decomposition needs to be performed only once for each dimension of the image. This yields a filtering which is both stable and consistent with the ideal infinite extension. The filter is efficient, requiring less computation than the standard recursive filtering. We give experimental evidence to verify these claims. (c) 2005 Elsevier B.V. All rights reserved.

Correlations between the clustering patterns of the pathological changes in sporadic Creutzfeldt-Jakob disease

Correlations between the clustering patterns of the vacuolation ('spongiform change'), prion protein (PrP) deposits, and surviving neurons were studied in the cerebral cortex, hippocampus, and cerebellum in 11 cases of sporadic Creutzfeldt-Jakob disease (sCJD). Differences in the sizes of the clusters of vacuoles were observed between brain regions and in the cerebral cortex, between the upper and lower laminae. With the exception of the parietal cortex, mean cluster size of the vacuoles was similar to that of the PrP deposits in each brain area. Clusters of the vacuoles were spatially correlated with the density of surviving neurons and with the clusters of PrP deposits in 47% and 53% of cortical areas analysed respectively but there were few spatial correlation between the PrP deposits and the density of surviving neurons. The data suggest that the pathology of sCJD may spread through the brain via specific anatomical pathways. Development of the clusters of vacuoles is spatially related to surviving neurons while the appearance of clusters of PrP deposits is related to the development of the vacuolation.

GIS linked environmental process models

In recent years there has been a great effort to combine the technologies and techniques of GIS and process models. This project examines the issues of linking a standard current generation 2½d GIS with several existing model codes. The focus for the project has been the Shropshire Groundwater Scheme, which is being developed to augment flow in the River Severn during drought periods by pumping water from the Shropshire Aquifer. Previous authors have demonstrated that under certain circumstances pumping could reduce the soil moisture available for crops. This project follows earlier work at Aston in which the effects of drawdown were delineated and quantified through the development of a software package that implemented a technique which brought together the significant spatially varying parameters. This technique is repeated here, but using a standard GIS called GRASS. The GIS proved adequate for the task and the added functionality provided by the general purpose GIS - the data capture, manipulation and visualisation facilities - were of great benefit. The bulk of the project is concerned with examining the issues of the linkage of GIS and environmental process models. To this end a groundwater model (Modflow) and a soil moisture model (SWMS2D) were linked to the GIS and a crop model was implemented within the GIS. A loose-linked approach was adopted and secondary and surrogate data were used wherever possible. The implications of which relate to; justification of a loose-linked versus a closely integrated approach; how, technically, to achieve the linkage; how to reconcile the different data models used by the GIS and the process models; control of the movement of data between models of environmental subsystems, to model the total system; the advantages and disadvantages of using a current generation GIS as a medium for linking environmental process models; generation of input data, including the use of geostatistic, stochastic simulation, remote sensing, regression equations and mapped data; issues of accuracy, uncertainty and simply providing adequate data for the complex models; how such a modelling system fits into an organisational framework.

A generative model for online depth fusion

We present a probabilistic, online, depth map fusion framework, whose generative model for the sensor measurement process accurately incorporates both long-range visibility constraints and a spatially varying, probabilistic outlier model. In addition, we propose an inference algorithm that updates the state variables of this model in linear time each frame. Our detailed evaluation compares our approach against several others, demonstrating and explaining the improvements that this model offers, as well as highlighting a problem with all current methods: systemic bias. © 2012 Springer-Verlag.

Laser cooling of a metastable argon atomic beam

A production of low velocity and monoenergetic atomic beams would increase the resolution in spectroscopic studies and many other experiments in atomic physics. Laser Cooling uses the radiation pressure to decelerate and cool atoms. The effusing from a glow discharge metastable argon atomic beam is affected by a counterpropagating laser light tuned to the cycling transition in argon. The Zeeman shift caused by a spatially varying magnetic field compensates for the changing Doppler shift that takes the atoms out of resonance as they decelerated. Deceleration and velocity bunching of atoms to a final velocity that depends on the detuning of the laser relative to a frequency of the transition have been observed. Time-of-Flight (TOF) spectroscopy is used to examine the velocity distribution of the cooled atomic beam. These TOF studies of the laser cooled atomic beam demonstrate the utility of laser deceleration for atomic-beam "velocity selection".

Standing stocks and carbon isotopes of live benthic foraminifera from the South Atlantic

Live (Rose Bengal stained) and dead benthic foraminifera of surface and subsurface sediments from 25 stations in the eastern South Atlantic Ocean and the Atlantic sector of the Southern Ocean were analyzed to decipher a potential influence of seasonally and spatially varying high primary productivity on the stable carbon isotopic composition of foraminiferal tests. Therefore, stations were chosen so that productivity strongly varied, whereas conservative water mass properties changed only little. To define the stable carbon isotopic composition of dissolved inorganic carbon (d13CDIC) in ambient water masses, we compiled new and previously published d13CDIC data in a section running from Antarctica through Agulhas, Cape and Angola Basins, via the Guinea Abyssal Plain to the Equator. We found that intraspecific d13C variability of all species at a single site is constantly low throughout their distribution within the sediments, i.e. species specific and site dependent mean values calculated from all subbottom depths on average only varied by +/-0.09 per mil. This is important because it makes the stable carbon isotopic signal of species independent of the particular microhabitat of each single specimen measured and thus more constant and reliable than has been previously assumed. So-called vital and/or microhabitat effects were further quantified: (1) d13C values of endobenthic Globobulimina affinis, Fursenkoina mexicana, and Bulimina mexicana consistently are by between -1.5 and -1.0 per mil VPDB more depleted than d13C values of preferentially epibenthic Fontbotia wuellerstorfi, Cibicidoides pachyderma, and Lobatula lobatula. (2) In contrast to the Antarctic Polar Front region, at all stations except one on the African continental slope Fontbotia wuellerstorfi records bottom water d13CDIC values without significant offset, whereas L. lobatula and C. pachyderma values deviate from bottom water values by about -0.4 per mil and -0.6 per mil, respectively. This adds to the growing amount of data on contrasting cibicid d13C values which on the one hand support the original 1:1-calibration of F. wuellerstorfi and bottom water d13CDIC, and on the other hand document severe depletions of taxonomically close relatives such as L. lobatula and C. pachyderma. At one station close to Bouvet Island at the western rim of Agulhas Basin, we interpret the offset of -1.5 per mil between bottom water d13CDIC and d13C values of infaunal living Bulimina aculeata in contrast to about -0.6 +/- 0.1 per mil measured at eight stations close-by, as a direct reflection of locally increased organic matter fluxes and sedimentation rates. Alternatively, we speculate that methane locally released from gas vents and related to hydrothermal venting at the mid-ocean ridge might have caused this strong depletion of 13C in the benthic foraminiferal carbon isotopic composition. Along the African continental margin, offsets between deep infaunal Globobulimina affinis and epibenthic Fontbotia wuellerstorfi as well as between shallow infaunal Uvigerina peregrina and F. wuellerstorfi, d13C values tend to increase with generally increasing organic matter decomposition rates. Although clearly more data are needed, these offsets between species might be used for quantification of biogeochemical paleogradients within the sediment and thus paleocarbon flux estimates. Furthermore, our data suggest that in high-productivity areas where sedimentary carbonate contents are lower than 15 weight %, epibenthic and endobenthic foraminiferal d13C values are strongly influenced by 13C enrichment probably due to carbonate-ion undersaturation, whereas above this sedimentary carbonate threshold endobenthic d13C values reflect depleted pore water d13CDIC values.

Surface wave and storm surge prediction during Hurricane Sandy

Hurricane Sandy was the largest storm on historical record in the Atlantic Ocean basin with extensive coastal damage caused by large waves and high storm surge. The primary objectives of this thesis are to compare and evaluate three different spatially-varying surface wind fields of Hurricane Sandy to investigate the impact of the differences between the complex wind fields on predictions of the sea surface evolution, and to evaluate the impact of the storm on the hydrodynamics in Great South Bay (GSB) and the discharge of ocean water into the back-barrier bay from overwash over Fire Island. Three different spatially-varying surface wind fields were evaluated and compared to wind observations, including the parametric Holland (1980) model (H80), the parametric Generalized Asymmetric Holland Model (GAHM), and results from the WeatherFlow Regional Atmospheric Modelling System (WRAMS). The winds were used to drive the coupled Delft3D-SWAN hydrodynamic and ocean wave models on a regional grid. The results indicate that the WRAMS wind field produces wave model predictions in the best agreement with significant wave height observations, followed by the GAHM and H80 wind fields and that a regional atmospheric wind model is best for hindcasting hurricane waves and water levels when detailed observations are available, while a parametric vortex model is best for forecasting hurricane sea surface conditions. Using a series of four connected Delft3D-SWAN grids to achieve finer resolution over Fire Island and GSB, a higher resolution WRAMS was used to predict waves and storm surge. The results indicate that strong local winds have the largest influence on water level fluctuations in GSB. Three numerical solutions were conducted with varying extents of barrier island overwash. The simulations allowing for minor and major overwash indicated good agreement with observations in the east end of GSB and suggest that island overwash provided a significant contribution of ocean water to GSB during the storm. Limiting the overwash in the numerical model directly impacts the total discharge into GSB from the ocean through existing inlets. The results of this study indicate that barrier island overwash had a significant impact on the water levels in eastern GSB.

Lesion detection performance: comparative analysis of low-dose CT data of the chest on two hybrid imaging systems

Incidental findings on low-dose CT images obtained during hybrid imaging are an increasing phenomenon as CT technology advances. Understanding the diagnostic value of incidental findings along with the technical limitations is important when reporting image results and recommending follow-up, which may result in an additional radiation dose from further diagnostic imaging and an increase in patient anxiety. This study assessed lesions incidentally detected on CT images acquired for attenuation correction on two SPECT/CT systems. Methods: An anthropomorphic chest phantom containing simulated lesions of varying size and density was imaged on an Infinia Hawkeye 4 and a Symbia T6 using the low-dose CT settings applied for attenuation correction acquisitions in myocardial perfusion imaging. Twenty-two interpreters assessed 46 images from each SPECT/CT system (15 normal images and 31 abnormal images; 41 lesions). Data were evaluated using a jackknife alternative free-response receiver-operating-characteristic analysis (JAFROC). Results: JAFROC analysis showed a significant difference (P < 0.0001) in lesion detection, with the figures of merit being 0.599 (95% confidence interval, 0.568, 0.631) and 0.810 (95% confidence interval, 0.781, 0.839) for the Infinia Hawkeye 4 and Symbia T6, respectively. Lesion detection on the Infinia Hawkeye 4 was generally limited to larger, higher-density lesions. The Symbia T6 allowed improved detection rates for midsized lesions and some lower-density lesions. However, interpreters struggled to detect small (5 mm) lesions on both image sets, irrespective of density. Conclusion: Lesion detection is more reliable on low-dose CT images from the Symbia T6 than from the Infinia Hawkeye 4. This phantom-based study gives an indication of potential lesion detection in the clinical context as shown by two commonly used SPECT/CT systems, which may assist the clinician in determining whether further diagnostic imaging is justified.

Superconducting Transition-Edge Sensor Physics

Despite record-setting performance demonstrated by superconducting Transition Edge Sensors (TESs) and growing utilization of the technology, a theoretical model of the physics governing TES devices superconducting phase transition has proven elusive. Earlier attempts to describe TESs assumed them to be uniform superconductors. Sadleir et al. 2010 shows that TESs are weak links and that the superconducting order parameter strength has significant spatial variation. Measurements are presented of the temperature T and magnetic field B dependence of the critical current Ic measured over 7 orders of magnitude on square Mo/Au bilayers ranging in length from 8 to 290 microns. We find our measurements have a natural explanation in terms of a spatially varying order parameter that is enhanced in proximity to the higher transition temperature superconducting leads (the longitudinal proximity effect) and suppressed in proximity to the added normal metal structures (the lateral inverse proximity effect). These in-plane proximity effects and scaling relations are observed over unprecedentedly long lengths (in excess of 1000 times the mean free path) and explained in terms of a Ginzburg-Landau model. Our low temperature Ic(B) measurements are found to agree with a general derivation of a superconducting strip with an edge or geometric barrier to vortex entry and we also derive two conditions that lead to Ic rectification. At high temperatures the Ic(B) exhibits distinct Josephson effect behavior over long length scales and following functional dependences not previously reported. We also investigate how film stress changes the transition, explain some transition features in terms of a nonequilibrium superconductivity effect, and show that our measurements of the resistive transition are not consistent with a percolating resistor network model.

Laser cooling and trapping of argon metastable atomic beam

The high velocity of free atoms associated with the thermal motion, together with the velocity distribution of atoms has imposed the ultimate limitation on the precision of ultrahigh resolution spectroscopy. A sample consisting of low velocity atoms would provide a substantial improvement in spectroscopy resolution. To overcome the problem of thermal motion, atomic physicists have pursued two goals; first, the reduction of the thermal motion (cooling); and second, the confinement of the atoms by means of electromagnetic fields (trapping). Cooling carried sufficiently far, eliminates the motional problems, whereas trapping allows for long observation times. In this work the laser cooling and trapping of an argon atomic beam will be discussed. The experiments involve a time-of-flight spectroscopy on metastable argon atoms. Laser deceleration or cooling of atoms is achieved by counter propagating a photon against an atomic beam of metastable atoms. The solution to the Doppler shift problem is achieved using spatially varying magnetic field along the beam path to Zeeman shift the atomic resonance frequency so as to keep the atoms in resonance with a fixed frequency cooling laser. For trapping experiments a Magnetooptical trap (MOT) will be used. The MOT is formed by three pairs of counter-propagating laser beams with mutual opposite circular polarization and a frequency tuned slightly below the center of the atomic resonance and superimposed on a magnetic quadrupole field.