Measurement of the vascular input function in mice for DCE-MRI

DCE-MRI is an important technique in the study of small animal cancer models because its sensitivity to vascular changes opens the possibility of quantitative assessment of early therapeutic response. However, extraction of physiologically descriptive parameters from DCE-MRI data relies upon measurement of the vascular input function (VIF), which represents the contrast agent concentration time course in the blood plasma. This is difficult in small animal models due to artifacts associated with partial volume, inflow enhancement, and the limited temporal resolution achievable with MR imaging. In this work, the development of a suite of techniques for high temporal resolution, artifact resistant measurement of the VIF in mice is described. One obstacle in VIF measurement is inflow enhancement, which decreases the sensitivity of the MR signal to the presence of contrast agent. Because the traditional techniques used to suppress inflow enhancement degrade the achievable spatiotemporal resolution of the pulse sequence, improvements can be achieved by reducing the time required for the suppression. Thus, a novel RF pulse which provides spatial presaturation contemporaneously with the RF excitation was implemented and evaluated. This maximizes the achievable temporal resolution by removing the additional RF and gradient pulses typically required for suppression of inflow enhancement. A second challenge is achieving the temporal resolution required for accurate characterization of the VIF, which exceeds what can be achieved with conventional imaging techniques while maintaining adequate spatial resolution and tumor coverage. Thus, an anatomically constrained reconstruction strategy was developed that allows for sampling of the VIF at extremely high acceleration factors, permitting capture of the initial pass of the contrast agent in mice. Simulation, phantom, and in vivo validation of all components were performed. Finally, the two components were used to perform VIF measurement in the murine heart. An in vivo study of the VIF reproducibility was performed, and an improvement in the measured injection-to-injection variation was observed. This will lead to improvements in the reliability of quantitative DCE-MRI measurements and increase their sensitivity.

Sparse regularization for fiber ODF reconstruction: from the suboptimality of ℓ2 and ℓ1 priors to ℓ0.

Diffusion MRI is a well established imaging modality providing a powerful way to probe the structure of the white matter non-invasively. Despite its potential, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a large variety of methods have been recently proposed to shorten the acquisition times. Among them, spherical deconvolution approaches have gained a lot of interest for their ability to reliably recover the intra-voxel fiber configuration with a relatively small number of data samples. To overcome the intrinsic instabilities of deconvolution, these methods use regularization schemes generally based on the assumption that the fiber orientation distribution (FOD) to be recovered in each voxel is sparse. The well known Constrained Spherical Deconvolution (CSD) approach resorts to Tikhonov regularization, based on an ℓ(2)-norm prior, which promotes a weak version of sparsity. Also, in the last few years compressed sensing has been advocated to further accelerate the acquisitions and ℓ(1)-norm minimization is generally employed as a means to promote sparsity in the recovered FODs. In this paper, we provide evidence that the use of an ℓ(1)-norm prior to regularize this class of problems is somewhat inconsistent with the fact that the fiber compartments all sum up to unity. To overcome this ℓ(1) inconsistency while simultaneously exploiting sparsity more optimally than through an ℓ(2) prior, we reformulate the reconstruction problem as a constrained formulation between a data term and a sparsity prior consisting in an explicit bound on the ℓ(0)norm of the FOD, i.e. on the number of fibers. The method has been tested both on synthetic and real data. Experimental results show that the proposed ℓ(0) formulation significantly reduces modeling errors compared to the state-of-the-art ℓ(2) and ℓ(1) regularization approaches.

Mid- to Late Holocene shoreline reconstruction and human occupation in Ancient Eretria (South Central Euboea, Greece)

Few studies have aimed to reconstruct landscape change in the area of Eretria (South Central Euboea, Greece) during the last 6000 years. The aim of this paper is to partially fill in this gap by examining the interaction be- tween Mid- to Late Holocene shoreline evolution and human occupation, which is documented in the area from the Late Neolithic to the Late Roman period (with discontinuities). Evidence of shoreline displacements is derived from the study of five boreholes (maximum depth of 5.25 m below the surface) drilled in the lowlands of Eretria. Based on sedimentological analyses and micro/macrofaunal identifications, different facies have been identified in the cores and which reveal typical features of deltaic progradation with marine, lagoonal, fluvio- deltaic and fluvial environments. In addition, a chronostratigraphy has been obtained based on 20 AMS 14C radio- carbon dates performed on samples of plant remains and marine/lagoonal shells found in situ. The main sequences of landscape reconstruction in the plain of Eretria can be summarized as follows: a marine environ- ment predominated from ca. 4000 to 3200 cal. BC and a gradual transition to shallow marine conditions is ob- served ca. 3200-3000 cal. BC due to the general context of deltaic progradation west of the ancient city. Subsequently, from ca. 3000 to 2000 cal. BC, a lagoon occupied the area in the vicinity of the Temple of Apollo and the settlement's development was restricted to several fluvio-deltaic levees, thus severely limiting human activities in the plain. From ca. 2000 to 800 cal. BC, a phase of shallow marine presence prevailed and constrained settlement on higher ground, forcing abandonment of the major part of the plain. Finally, since the eighth century BC, the sea has regressed southward and created the modern landscape.

Structural analysis, clay mineralogy and K-Ar dating of fault gouges from Centovalli Line (Central Alps) for reconstruction of their recent activity

Between the cities of Domodossola and Locarno, the complex ``Centovalli Line'' tectonic zone of the Central Alps outlines deformation phases over a long period of time (probably starting similar to 30 Ma ago) and under variable P-T conditions. The last deformation phases developed gouge-bearing faults with a general E-W trend that crosscuts the roots of the Alpine Canavese zone and the Finero ultramafic body. Kinematic indicators show that the general motion was mainly dextral associated with back thrusting towards the S. The <2 mu m clay fractions of fault gouges from Centovalli Line consist mainly of illite, smectite and chlorite with varied illite-smectite, chlorite-smectite and chlorite-serpentine mixed-layers. Constrained with the illite crystallinity index, the thermal conditions induced by the tectonic activity show a gradual trend from anchizonal to diagenetic conditions. The <2 and <0.2 mu M clay fractions, and hydrothermal K-feldspar separates all provide K-Ar ages between 14.2 +/- 2.9 Ma and roughly 0 Ma, with major episodes at about 12,8, 6 and close to 0 Ma These ages set the recurrent tectonic activity and the associated fluid circulations between Upper Miocene and Recent. On the basis of the K-Ar ages and with a thermal gradient of 25-30 degrees C/km, the studied fault zones were located at a depth of 4-7 km. If they were active until now as observed in field, the exhumation was approximately 2.5-3.0 km for the last 12 Ma with a mean velocity of 0.4 mm/y. Comparison with available models on the recent Alpine evolution shows that the tectonic activity in the area relates to a continuum of the back-thrusting movements of the Canavese Line, and/or to several late-extensional phases of the Rhone-Simplon line. The Centovalli-Val Vigezzo zone therefore represents a major tectonic zone of the Central-Western Alps resulting from different interacting tectonic events. (C) 2011 Elsevier B.V. All rights reserved.

Transfer entropy reconstruction and labeling of neuronal connections from simulated calcium imaging

Neuronal dynamics are fundamentally constrained by the underlying structural network architecture, yet much of the details of this synaptic connectivity are still unknown even in neuronal cultures in vitro. Here we extend a previous approach based on information theory, the Generalized Transfer Entropy, to the reconstruction of connectivity of simulated neuronal networks of both excitatory and inhibitory neurons. We show that, due to the model-free nature of the developed measure, both kinds of connections can be reliably inferred if the average firing rate between synchronous burst events exceeds a small minimum frequency. Furthermore, we suggest, based on systematic simulations, that even lower spontaneous inter-burst rates could be raised to meet the requirements of our reconstruction algorithm by applying a weak spatially homogeneous stimulation to the entire network. By combining multiple recordings of the same in silico network before and after pharmacologically blocking inhibitory synaptic transmission, we show then how it becomes possible to infer with high confidence the excitatory or inhibitory nature of each individual neuron.

Reconstruction of gusty wind speed time series from autonomous data logger records

The collection of wind speed time series by means of digital data loggers occurs in many domains, including civil engineering, environmental sciences and wind turbine technology. Since averaging intervals are often significantly larger than typical system time scales, the information lost has to be recovered in order to reconstruct the true dynamics of the system. In the present work we present a simple algorithm capable of generating a real-time wind speed time series from data logger records containing the average, maximum, and minimum values of the wind speed in a fixed interval, as well as the standard deviation. The signal is generated from a generalized random Fourier series. The spectrum can be matched to any desired theoretical or measured frequency distribution. Extreme values are specified through a postprocessing step based on the concept of constrained simulation. Applications of the algorithm to 10-min wind speed records logged at a test site at 60 m height above the ground show that the recorded 10-min values can be reproduced by the simulated time series to a high degree of accuracy.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

A tephra lattice for Greenland and a reconstruction of volcanic events spanning 25–45 ka b2k

Tephra layers preserved within the Greenland ice-cores are crucial for the independent synchronisation of these high-resolution records to other palaeoclimatic archives. Here we present a new and detailed tephrochronological framework for the time period 25,000 e 45,000 a b2k that brings together results from 4 deep Greenland ice-cores. In total, 99 tephra deposits, the majority of which are preserved as cryptotephra, are described from the NGRIP, NEEM, GRIP and DYE-3 records. The major element signatures of single glass shards within these deposits indicate that 93 are basaltic in composition all originating from Iceland. Specifically, 43 originate from Grimsv € otn, 20 are thought to be sourced from the Katla volcanic system and 17 show affinity to the Kverkfj € oll system. Robust geochemical characterisations, independent ages derived from the GICC05 ice-core chronology, and the stratigraphic positions of these deposits relative to the Dansgaard-Oeschger climate events represent a key framework that provides new information on the frequency and nature of volcanic events in the North Atlantic region between GS-3 and GI-12. Of particular importance are 19 tephra deposits that lie on the rapid climatic transitions that punctuate the last glacial period. This framework of well-constrained, time-synchronous tie-lines represents an important step towards the independent synchronisation of marine, terrestrial and ice-core records from the North Atlantic region, in order to assess the phasing of rapid climatic changes during the last glacial period.

(Figure 2) Organic carbon content, stable isotope ratios, TEX86 index and sea surface temperature reconstruction for the early Aptian of ODP Hole 198-1207B

The Cretaceous has long been recognized as a time when greenhouse conditions were fueled by elevated atmospheric CO2 and accompanied by perturbations of the global carbon cycle described as oceanic anoxic events (OAEs). Yet, the magnitude and frequency of temperature change during this interval of warm and equable climate are poorly constrained. Here we present a high-resolution record of sea-surface temperatures (SSTs) reconstructed using the TEX86 paleothermometer for a sequence of early Aptian organic-rich sediments deposited during the first Cretaceous OAE (OAE1a) at Shatsky Rise in the tropical Pacific. SSTs range from ~30 to ~36 °C and include two prominent cooling episodes of ~4 °C. The cooler temperatures reflect significant temperature instability in the tropics likely triggered by changes in carbon cycling induced by enhanced burial of organic matter. SST instability recorded during the early Aptian in the Pacific is comparable to that reported for the late Albian-early Cenomanian in the Atlantic, suggesting that such climate perturbations may have recurred during the Cretaceous with concomitant consequences for biota and the marine environment.

Sea surface temperture reconstruction for ODP Site 189-1172

Relative to the present day, meridional temperature gradients in the Early Eocene age (~56-53 Myr ago) were unusually low, with slightly warmer equatorial regions (Pearson et al., 2007, doi:10.1130/G23175A.1 ) but with much warmer subtropical Arctic (Sluijs et al., 2008, doi:10.1029/2007PA001495) and mid-latitude (Sluijs et al., 2007, doi:10.1038/nature06400) climates. By the end of the Eocene epoch (~34 Myr ago), the first major Antarctic ice sheets had appeared (Zachos et al., 1992, doi:10.1130/0091-7613(1992)020<0569:EOISEO>2.3.CO;2; Barker et al., 2007, doi:10.1016/j.dsr2.2007.07.027), suggesting that major cooling had taken place. Yet the global transition into this icehouse climate remains poorly constrained, as only a few temperature records are available portraying the Cenozoic climatic evolution of the high southern latitudes. Here we present a uniquely continuous and chronostratigraphically well-calibrated TEX86 record of sea surface temperature (SST) from an ocean sediment core in the East Tasman Plateau (palaeolatitude ~65° S). We show that southwest Pacific SSTs rose above present-day tropical values (to ~34° C) during the Early Eocene age (~53 Myr ago) and had gradually decreased to about 21° C by the early Late Eocene age (~36 Myr ago). Our results imply that there was almost no latitudinal SST gradient between subequatorial and subpolar regions during the Early Eocene age (55-50 Myr ago). Thereafter, the latitudinal gradient markedly increased. In theory, if Eocene cooling was largely driven by a decrease in atmospheric greenhouse gas concentration Zachos et al. (2008, doi:10.1038/nature06588), additional processes are required to explain the relative stability of tropical SSTs given that there was more significant cooling at higher latitudes.

Palynology and sea surface temperature reconstruction for the mid-Cretaceous Oceanic Anoxic Event 2 at ODP Hole 210-1276A

The mid-Cretaceous is thought to be a greenhouse world with significantly higher atmospheric pCO2 and sea-surface temperatures as well as a much flatter latitudinal thermal gradient compared to the present. This time interval was punctuated by the Cenomanian/Turonian Oceanic Anoxic Event (OAE-2, ~ 93.5 Myr ago), an episode of global, massive organic carbon burial that likely resulted in a large and abrupt pCO2 decline. However, the climatic consequences of this pCO2 drop are yet poorly constrained. We determined the first, high-resolution sea-surface temperature (SST) record across OAE-2 from a deep-marine sedimentary sequence at Ocean Drilling Program (ODP) Site 1276 in the mid-latitudinal Newfoundland Basin, NW Atlantic. By employing the organic palaeothermometer TEX86, we found that SSTs across the OAE-2 interval were extremely high, but were punctuated by a remarkably large cooling (5-11 °C), which is synchronous with the 2.5-5.5 °C cooling in SST records from equatorial Atlantic sites, and the "Plenus Cold Event". Because this global cooling event is concurrent with increased organic carbon burial, it likely acted in response to the associated pCO2 drop. Our findings imply a substantial increase in the latitudinal SST gradient in the proto-North Atlantic during this period of global cooling and reduced atmospheric pCO2, suggesting a strong coupling between pCO2 and latitudinal thermal gradients under greenhouse climate conditions.

An Investigation into the Failure to Implement a Universal School Lunch Program during WWII and Postwar Reconstruction in Canada: The Case of Ontario

The purpose of this research is to investigate the various social, political and economic factors that contributed to Canada’s failure to implement a universal school lunch program during the 1940s. Although Canada developed several other social welfare programs in the post-war period, it remains one of the only industrialized nations that does not provide hot meals to children in elementary or secondary schools. Data from the province of Ontario, a major site of postwar reconstruction and policy-making, has been taken up to inform the broader national discourse on school lunches from the 1940s. National, Ontario provincial and City of Toronto archival records were collected and analyzed according to common themes, in order to identify key barriers that constrained government support of a hot meal program. Archival records were identified using key words, and were limited to materials created between 1930-1952. Analysis suggests that sufficient need for a hot meal program had not been established during the 1940s. Despite misleading nutrition messages, rates of malnutrition and nutrient-related disease were at an all-time low, and many Ontario school boards did not appear to have the necessary infrastructure required to supply all pupils with hot meals. The Canadian government had already employed significant resources to improve existing social security programs by coupling them with health education. This strategy reflected a shift in understanding malnutrition as a knowledge-based problem, as opposed to income-based. This understanding was further reinforced through the moralized dissemination of nutrition information, which placed blame on women for improperly raising their children. Ultimately, the strong uptake of nutrition as a public health issue in Ontario may have limited prospective responses to solutions already utilized in the public health domain, and directed favour away from a universal school lunch program for Canada.

Improving the Accuracy of 3-D Reconstruction in Robotic Vision Applications

Thesis (Ph.D.)--University of Washington, 2016-08