IDENTIFYING AND MONITORING THE ROLES OF CAVITATION IN HEATING FROM HIGH-INTENSITY FOCUSED ULTRASOUND

For high-intensity focused ultrasound (HIFU) to continue to gain acceptance for cancer treatment it is necessary to understand how the applied ultrasound interacts with gas trapped in the tissue. The presence of bubbles in the target location have been thought to be responsible for shielding the incoming pressure and increasing local heat deposition due to the bubble dynamics. We lack adequate tools for monitoring the cavitation process, due to both limited visualization methods and understanding of the underlying physics. The goal of this project was to elucidate the role of inertial cavitation in HIFU exposures in the hope of applying noise diagnostics to monitor cavitation activity and control HIFU-induced cavitation in a beneficial manner. A number of approaches were taken to understand the relationship between inertial cavitation signals, bubble heating, and bubble shielding in agar-graphite tissue phantoms. Passive cavitation detection (PCD) techniques were employed to detect inertial bubble collapses while the temperature was monitored with an embedded thermocouple. Results indicate that the broadband noise amplitude is correlated to bubble-enhanced heating. Monitoring inertial cavitation at multiple positions throughout the focal region demonstrated that bubble activity increased prefocally as it diminished near the focus. Lowering the HIFU duty cycle had the effect of maintaining a more or less constant cavitation signal, suggesting the shielding effect diminished when the bubbles had a chance to dissolve during the HIFU off-time. Modeling the effect of increasing the ambient temperature showed that bubbles do not collapse as violently at higher temperatures due to increased vapor pressure inside the bubble. Our conclusion is that inertial cavitation heating is less effective at higher temperatures and bubble shielding is involved in shifting energy deposition at the focus. The use of a diagnostic ultrasound imaging system as a PCD array was explored. Filtering out the scattered harmonics from the received RF signals resulted in a spatially- resolved inertial cavitation signal, while the amplitude of the harmonics showed a correlation with temperatures approaching the onset of boiling. The result is a new tool for detecting a broader spectrum of bubble activity and thus enhancing HIFU treatment visualization and feedback.

PHOTOREFRACTIVE CRYSTAL-BASED ACOUSTO-OPTIC IMAGING IN THE NEAR-INFRARED AND ITS APPLICATIONS

Acousto-optic (AO) sensing and imaging (AOI) is a dual-wave modality that combines ultrasound with diffusive light to measure and/or image the optical properties of optically diffusive media, including biological tissues such as breast and brain. The light passing through a focused ultrasound beam undergoes a phase modulation at the ultrasound frequency that is detected using an adaptive interferometer scheme employing a GaAs photorefractive crystal (PRC). The PRC-based AO system operating at 1064 nm is described, along with the underlying theory, validating experiments, characterization, and optimization of this sensing and imaging apparatus. The spatial resolution of AO sensing, which is determined by spatial dimensions of the ultrasound beam or pulse, can be sub-millimeter for megahertz-frequency sound waves.A modified approach for quantifying the optical properties of diffuse media with AO sensing employs the ratio of AO signals generated at two different ultrasound focal pressures. The resulting “pressure contrast signal” (PCS), once calibrated for a particular set of pressure pulses, yields a direct measure of the spatially averaged optical transport attenuation coefficient within the interaction volume between light and sound. This is a significant improvement over current AO sensing methods since it produces a quantitative measure of the optical properties of optically diffuse media without a priori knowledge of the background illumination. It can also be used to generate images based on spatial variations in both optical scattering and absorption. Finally, the AO sensing system is modified to monitor the irreversible optical changes associated with the tissue heating from high intensity focused ultrasound (HIFU) therapy, providing a powerful method for noninvasively sensing the onset and growth of thermal lesions in soft tissues. A single HIFU transducer is used to simultaneously generate tissue damage and pump the AO interaction. Experimental results performed in excised chicken breast demonstrate that AO sensing can identify the onset and growth of lesion formation in real time and, when used as feedback to guide exposure parameters, results in more predictable lesion formation.

SCANNING ACOUSTIC MICROSCOPE FOR CHARACTERIZATION OF ARTERIAL PLAQUE

Unstable arterial plaque is likely the key component of atherosclerosis, a disease which is responsible for two-thirds of heart attacks and strokes, leading to approximately 1 million deaths in the United States. Ultrasound imaging is able to detect plaque but as of yet is not able to distinguish unstable plaque from stable plaque. In this work a scanning acoustic microscope (SAM) was implemented and validated as tool to measure the acoustic properties of a sample. The goal for the SAM is to be able to provide quantitative measurements of the acoustic properties of different plaque types, to understand the physical basis by which plaque may be identified acoustically. The SAM consists of a spherically focused transducer which operates in pulse-echo mode and is scanned in a 2D raster pattern over a sample. A plane wave analysis is presented which allows the impedance, attenuation and phase velocity of a sample to be de- termined from measurements of the echoes from the front and back of the sample. After the measurements, the attenuation and phase velocity were analysed to ensure that they were consistent with causality. The backscatter coefficient of the samples was obtained using the technique outlined by Chen et al [8]. The transducer used here was able to determine acoustic properties from 10-40 MHz. The results for the impedance, attenuation and phase velocity were validated for high and low-density polyethylene against published results. The plane wave approximation was validated by measuring the properties throughout the focal region and throughout a range of incidence angles from the transducer. The SAM was used to characterize a set of recipes for tissue-mimicking phantoms which demonstrate indepen- dent control over the impedance, attenuation, phase velocity and backscatter coefficient. An initial feasibility study on a human artery was performed.

QUANTITATIVE THREE DIMENSIONAL ELASTICITY IMAGING

Neoplastic tissue is typically highly vascularized, contains abnormal concentrations of extracellular proteins (e.g. collagen, proteoglycans) and has a high interstitial fluid pres- sure compared to most normal tissues. These changes result in an overall stiffening typical of most solid tumors. Elasticity Imaging (EI) is a technique which uses imaging systems to measure relative tissue deformation and thus noninvasively infer its mechanical stiffness. Stiffness is recovered from measured deformation by using an appropriate mathematical model and solving an inverse problem. The integration of EI with existing imaging modal- ities can improve their diagnostic and research capabilities. The aim of this work is to develop and evaluate techniques to image and quantify the mechanical properties of soft tissues in three dimensions (3D). To that end, this thesis presents and validates a method by which three dimensional ultrasound images can be used to image and quantify the shear modulus distribution of tissue mimicking phantoms. This work is presented to motivate and justify the use of this elasticity imaging technique in a clinical breast cancer screening study. The imaging methodologies discussed are intended to improve the specificity of mammography practices in general. During the development of these techniques, several issues concerning the accuracy and uniqueness of the result were elucidated. Two new algorithms for 3D EI are designed and characterized in this thesis. The first provides three dimensional motion estimates from ultrasound images of the deforming ma- terial. The novel features include finite element interpolation of the displacement field, inclusion of prior information and the ability to enforce physical constraints. The roles of regularization, mesh resolution and an incompressibility constraint on the accuracy of the measured deformation is quantified. The estimated signal to noise ratio of the measured displacement fields are approximately 1800, 21 and 41 for the axial, lateral and eleva- tional components, respectively. The second algorithm recovers the shear elastic modulus distribution of the deforming material by efficiently solving the three dimensional inverse problem as an optimization problem. This method utilizes finite element interpolations, the adjoint method to evaluate the gradient and a quasi-Newton BFGS method for optimiza- tion. Its novel features include the use of the adjoint method and TVD regularization with piece-wise constant interpolation. A source of non-uniqueness in this inverse problem is identified theoretically, demonstrated computationally, explained physically and overcome practically. Both algorithms were test on ultrasound data of independently characterized tissue mimicking phantoms. The recovered elastic modulus was in all cases within 35% of the reference elastic contrast. Finally, the preliminary application of these techniques to tomosynthesis images showed the feasiblity of imaging an elastic inclusion.

EVALUATION OF HARMONIC MOTION ELASTOGRAPHY AND ACOUSTO-­OPTIC IMAGING FOR MONITORING LESION FORMATION BY HIGH INTENSITY FOCUSED ULTRASOUND

Malignant or benign tumors may be ablated with high‐intensity focused ultrasound (HIFU). This technique, known as focused ultrasound surgery (FUS), has been actively investigated for decades, but slow to be implemented and difficult to control due to lack of real‐time feedback during ablation. Two methods of imaging and monitoring HIFU lesions during formation were implemented simultaneously, in order to investigate the efficacy of each and to increase confidence in the detection of the lesion. The first, Acousto‐Optic Imaging (AOI) detects the increasing optical absorption and scattering in the lesion. The intensity of a diffuse optical field in illuminated tissue is mapped at the spatial resolution of an ultrasound focal spot, using the acousto‐optic effect. The second, Harmonic Motion Imaging (HMI), detects the changing stiffness in the lesion. The HIFU beam is modulated to force oscillatory motion in the tissue, and the amplitude of this motion, measured by ultrasound pulse‐echo techniques, is influenced by the stiffness. Experiments were performed on store‐bought chicken breast and freshly slaughtered bovine liver. The AOI results correlated with the onset and relative size of forming lesions much better than prior knowledge of the HIFU power and duration. For HMI, a significant artifact was discovered due to acoustic nonlinearity. The artifact was mitigated by adjusting the phase of the HIFU and imaging pulses. A more detailed model of the HMI process than previously published was made using finite element analysis. The model showed that the amplitude of harmonic motion was primarily affected by increases in acoustic attenuation and stiffness as the lesion formed and the interaction of these effects was complex and often counteracted each other. Further biological variability in tissue properties meant that changes in motion were masked by sample‐to‐sample variation. The HMI experiments predicted lesion formation in only about a quarter of the lesions made. In simultaneous AOI/HMI experiments it appeared that AOI was a more robust method for lesion detection.

Interior pathways of the North Atlantic meridional overturning circulation.

To understand how our global climate will change in response to natural and anthropogenic forcing, it is essential to determine how quickly and by what pathways climate change signals are transported throughout the global ocean, a vast reservoir for heat and carbon dioxide. Labrador Sea Water (LSW), formed by open ocean convection in the subpolar North Atlantic, is a particularly sensitive indicator of climate change on interannual to decadal timescales. Hydrographic observations made anywhere along the western boundary of the North Atlantic reveal a core of LSW at intermediate depths advected southward within the Deep Western Boundary Current (DWBC). These observations have led to the widely held view that the DWBC is the dominant pathway for the export of LSW from its formation site in the northern North Atlantic towards the Equator. Here we show that most of the recently ventilated LSW entering the subtropics follows interior, not DWBC, pathways. The interior pathways are revealed by trajectories of subsurface RAFOS floats released during the period 2003-2005 that recorded once-daily temperature, pressure and acoustically determined position for two years, and by model-simulated 'e-floats' released in the subpolar DWBC. The evidence points to a few specific locations around the Grand Banks where LSW is most often injected into the interior. These results have implications for deep ocean ventilation and suggest that the interior subtropical gyre should not be ignored when considering the Atlantic meridional overturning circulation.

Three Hundred Years in Annapolis: Phase III Archaeological Investigations of the Anne Arundel County Courthouse Site (18AP63), Annapolis, Maryland

During the summer of 1994, Archaeology in Annapolis conducted archaeological investigations of the city block bounded by Franklin, South and Cathedral Streets in the city of Annapolis. This Phase III excavation was conducted as a means to identify subsurface cultural resources in the impact area associated with the proposed construction of the Anne Arundel County Courthouse addition. This impact area included both the upper and lower parking lots used by Courthouse employees. Investigations were conducted in the form of mechanical trenching and hand excavated units. Excavations in the upper lot area yielded significant information concerning the interior area of the block. Known as Bellis Court, this series of rowhouses was constructed in the late nineteenth century and was used as rental properties by African-Americans. The dwellings remained until the middle of the twentieth century when they were demolished in preparation for the construction of a Courthouse addition. Portions of the foundation of a house owned by William H. Bellis in the 1870s were also exposed in this area. Construction of this house was begun by William Nicholson around 1730 and completed by Daniel Dulany in 1732/33. It was demolished in 1896 by James Munroe, a Trustee for Bellis. Excavations in the upper lot also revealed the remains of a late seventeenth/early eighteenth century wood-lined cellar, believed to be part of the earliest known structure on Lot 58. After an initially rapid deposition of fill around 1828, this cellar was gradually covered with soil throughout the remainder of the nineteenth century. The fill deposit in the cellar feature yielded a mixed assemblage of artifacts that included sherds of early materials such as North Devon gravel-tempered earthenware, North Devon sgraffito and Northem Italian slipware, along with creamware, pearlware and whiteware. In the lower parking lot, numerous artifacts were recovered from yard scatter associated with the houses that at one time fronted along Cathedral Street and were occupied by African- Americans. An assemblage of late seventeenth century/early eighteenth century materials and several slag deposits from an early forge were recovered from this second area of study. The materials associated with the forge, including portions of a crucible, provided evidence of some of the earliest industry in Annapolis. Investigations in both the upper and lower parking lots added to the knowledge of the changing landscape within the project area, including a prevalence of open space in early periods, a surprising survival of impermanent structures, and a gradual regrading and filling of the block with houses and interior courts. Excavations at the Anne Arundel County Courthouse proved this to be a multi-component site, rich in cultural resources from Annapolis' Early Settlement Period through its Modern Period (as specified by Maryland's Comprehensive Historic Preservation Plan (Weissman 1986)). This report provides detailed interpretations of the archaeological findings of these Phase III investigations.

The island-scale internal wave climate of Moorea, French Polynesia

Analysis of five-year records of temperatures and currents collected at Moorea reveal strong internal wave activity at predominantly semi-diurnal frequencies impacting reef slopes at depths 30m around the entire island. Temperature changes of 1.5C to 3C are accompanied by surges of upward and onshore flow and vertical shear in onshore currents. Superimposed on annual temperature changes of approximately 3C, internal wave activity is high from Oct-May and markedly lower from Jun-Sep. The offshore pycnocline is broadly distributed with continuous stratification to at least 500m depth, and a subsurface fluorescence maximum above the strong nutricline at approximately 200m. Minimum buoyancy periods range from 4.8 to 6min, with the maximum density gradient occurring at 50 to 60m depth in summer and deepening to approximately 150 to 200m in winter. The bottom slope angle around all of Moorea is super-critical relative to the vertical stratification angle suggesting that energy propagating into shallow water is only a portion of total incident internal wave energy. Vertical gradient Richardson numbers indicate dominance by density stability relative to current shear with relatively limited diapycnal mixing. Coherence and lagged cross-correlation of semi-diurnal temperature variation indicate complex patterns of inter-site arrival of internal waves and no clear coherence or lagged correlation relationships among island sides. Semi-diurnal and high frequency internal wave packets likely arrive on Moorea from a combination of local and distant sources and may have important impacts for nutrient and particle fluxes in deep reef environments. © 2012 American Geophysical Union. All Rights Reserved.

Ocean structure and climate (Eastern North Atlantic): in situ measurement and remote sensing

Structure and climate of the east North Atlantic are appraised within a framework of in situ measurement and altimeter remote sensing from 0 degree - 60 degree N. Long zonal expendable bathythermograph /conductivity-temperature-depth probe sections show repeating internal structure in the North Atlantic Ocean. Drogued buoys and subsurface floats give westward speeds for eddies and wavelike structure. Records from longterm current meter deployments give the periodicity of the repeating structure. Eddy and wave characteristics of period, size or wavelength, westward propagation speed, and mean currents are derived at 20 degree N, 26 degree N, 32.5 degree N, 36 degree N and 48 degree N from in situ measurements in the Atlantic Ocean. It is shown that ocean wave and eddy-like features measured in situ correlate with altimeter structure. Interior ocean wave crests or cold dome-like temperature structures are cyclonic and have negative surface altimeter anomalies; mesoscale internal wave troughs or warm structures are anticyclonic and have positive surface height anomalies. Along the Eastern Boundary, flows and temperature climate are examined in terms of sla and North Atlantic Oscillation (NAO) Index. Longterm changes in ocean climate and circulation are derived from sla data. It is shown that longterm changes from 1992 to 2002 in the North Atlantic Current and the Subtropical Gyre transport determined from sla data correlate with winter NAO Index such that maximum flow conditions occurred in 1995 and 2000. Minimum circulation conditions occurred between 1996-1998. Years of extreme negative winter NAO Index result in enhanced poleward flow along the Eastern Boundary and anomalous winter warming along the West European Continental Slope as was measured in 1990, 1996, 1998 and 2001.

Particulate and dissolved aluminum and titanium in the upper water column of the Atlantic Ocean

This study investigates the oceanic behavior of the lithogenic trace elements Al and Ti in the upper 200 m of the Atlantic Ocean. The distribution of both metals in the dissolved and particulate phases was assessed along an E-W transect in the eastern tropical North Atlantic (December 2009) and along a meridional Atlantic transect (April-May 2010). The surface water concentrations of particulate and dissolved Al and Ti reflected the previously observed pattern of atmospheric inputs into the Atlantic Ocean. Subsurface minima at stations with pronounced fluorescence maxima were observed, suggesting a link between biological productivity and the removal of both dissolved and particulate Al and Ti. This may include uptake mechanisms, adsorption and aggregation processes on biogenic particle surfaces and the formation of large, fast sinking biogenic particles, e.g., fecal pellets. Residence times in the upper water column (100 m) of the tropical and subtropical North Atlantic were estimated to range in the order of days to weeks in the particulate phases (Al: 3-22 days, Ti: 4-37 days) and were 0.9-3.8 years for Al and 10-31 years for Ti in the dissolved phases. Longer residence times in both phases in the South Atlantic are consistent with lower biological productivity and decreased removal rates. In the upper water column, Al was predominantly present in the dissolved form, whereas Ti mostly occurred in the particulate form. Largest deviations in the partition coefficients between the particulate and dissolved phases were found in the surface waters, together with excess dissolved Al over Ti compared to the crustal source. This likely reflects elevated dissolution of Al compared to Ti from atmospheric mineral particles.

Phytoplankton carbon fixation chlorophyll-biomass and diagnostic pigments in the Atlantic Ocean

We have made daily measurements of phytoplankton pigments, size-fractionated (<2 and >2-μm) carbon fixation and chlorophyll-a concentration during four Atlantic Meridional Transect (AMT) cruises in 2003–04. Surface rates of carbon fixation ranged from <0.2-mmol C m−3 d−1 in the subtropical gyres to 0.2–0.5-mmol C m−3 d−1 in the tropical equatorial Atlantic. Significant intercruise variability was restricted to the subtropical gyres, with higher chlorophyll-a concentrations and carbon fixation in the subsurface chlorophyll maximum during spring in either hemisphere. In surface waters, although picoplankton (<2-μm) represented the dominant fraction in terms of both carbon fixation (50–70%) and chlorophyll-a (80–90%), nanoplankton (>2-μm) contributions to total carbon fixation (30–50%) were higher than to total chlorophyll-a (10–20%). However, in the subsurface chlorophyll maximum picoplankton dominated both carbon fixation (70–90%) and chlorophyll-a (70–90%). Thus, in surface waters chlorophyll-normalised carbon fixation was 2–3 times higher for nanoplankton and differences in picoplankton and nanoplankton carbon to chlorophyll-a ratios may lead to either higher or similar growth rates. These low chlorophyll-normalised carbon fixation rates for picoplankton may also reflect losses of fixed carbon (cell leakage or respiration), decreases in photosynthetic efficiency, grazing losses during the incubations, or some combination of all these. Comparison of nitrate concentrations in the subsurface chlorophyll maximum with estimates of those required to support the observed rates of carbon fixation (assuming Redfield stoichiometry) indicate that primary production in the chlorophyll maximum may be light rather than nutrient limited.

Seasonal dynamics of the carbonate system in the Western English Channel

We present over 900 carbonate system observations collected over four years (2007–2010) in the Western English Channel (WEC). We determined CO2 partial pressure (pCO2), Total Alkalinity (TA) and Dissolved Inorganic Carbon (DIC) along a series of 40 km transects, including two oceanographic stations (L4 and E1) within a sustained coastal observatory. Our data follow a seasonal pattern of CO2 undersaturation from January to August, followed by supersaturation in September–October and a return to near-equilibrium thereafter. This pattern is explained by the interplay of thermal and biological sinks in winter and spring–summer, respectively, followed by the breakdown of stratification and mixing with deeper, high-CO2 water in autumn. The drawdown of DIC and inorganic N between March and June with a C:N ratio of 8.7–9.5 was consistent with carbon over-consumption during phytoplankton growth. Monthly mean surface pCO2 was strongly correlated with depth integrated chlorophyll a highlighting the importance of subsurface chlorophyll a maxima in controlling C-fluxes in shelf seas. Mixing of seawater with riverine freshwater in near-shore samples caused a reduction in TA and the saturation state of calcite minerals, particularly in winter. Our data show that the L4 and E1 oceanographic stations were small, net sinks for atmospheric CO2 over an annual cycle (−0.52±0.66 mol C m−2 y−1 and −0.62±0.49 mol C m−2 y−1, respectively).

Low colloidal associations of aluminium and titanium in surface waters of the tropical Atlantic

The distribution of dissolved, soluble and colloidal fractions of Al and Ti was assessed by ultrafiltration studies in the upper water column of the eastern tropical North Atlantic. The dissolved fractions of both metals were found to be dominated by the soluble phase smaller than 10 kDa. The colloidal associations were very low (0.2–3.4%) for Al and not detectable for Ti. These findings are in some contrast to previous estimations for Ti and to the predominant occurrence of both metals as hydrolyzed species in seawater. However, low tendencies to form inorganic colloids can be expected, as in seawater dissolved Al and dissolved Ti are present within their inorganic solubility levels. In addition, association with functional organic groups in the colloidal phase is unlikely for both metals. Vertical distributions of the dissolved fractions showed surface maxima with up to 43 nM of Al and 157 pM of Ti, reflecting their predominant supply from atmospheric sources to the open ocean. In the surface waters, excess dissolved Al over dissolved Ti was present compared to the crustal source, indicating higher solubility and thus elevated inputs of dissolved Al from atmospheric mineral particles. At most stations, subsurface minima of Al and Ti were observed and can be ascribed to scavenging processes and/or biological uptake. The dissolved Al concentrations decreased by 80–90% from the surface maximum to the subsurface minimum. Estimated residence times in the upper 100 m of the water column ranged between 1.6 and 4 years for dissolved Al and between 14 and 17 years for dissolved Ti. The short residence times are in some contrast to the low colloidal associations of Al and Ti and the assumed role of colloids as intermediates in scavenging processes. This suggests that either the removal of both metals occurs predominantly via direct transfer of the hydrolyzed species into the particulate fraction or that the colloidal phase is rapidly turned over in the upper water column.

Methanol acetaldehyde and acetone in the surface waters of the Atlantic Ocean

Oceanic methanol, acetaldehyde, and acetone concentrations were measured during an Atlantic Meridional Transect (AMT) cruise from the UK to Chile (49°N to 39°S) in 2009. Methanol (48–361 nM) and acetone (2–24 nM) varied over the track with enrichment in the oligotrophic Northern Atlantic Gyre. Acetaldehyde showed less variability (3–9 nM) over the full extent of the transect. These oxygenated volatile organic compounds (OVOCs) were also measured subsurface, with methanol and acetaldehyde mostly showing homogeneity throughout the water column. Acetone displayed a reduction below the mixed layer. OVOC concentrations did not consistently correlate with primary production or chlorophyll-a levels in the surface Atlantic Ocean. However, we did find a novel and significant negative relationship between acetone concentration and bacterial leucine incorporation, suggesting that acetone might be removed by marine bacteria as a source of carbon. Microbial turnover of both acetone and acetaldehyde was confirmed. Modeled atmospheric data are used to estimate the likely air-side OVOC concentrations. The direction and magnitude of air-sea fluxes vary for all three OVOCs depending on location. We present evidence that the ocean may exhibit regions of acetaldehyde under-saturation. Extrapolation suggests that the Atlantic Ocean represents an overall source of these OVOCs to the atmosphere at 3, 3, and 1 Tg yr−1 for methanol, acetaldehyde, and acetone, respectively.

Variations in GDGT distributions through the water column in the South East Atlantic Ocean

The TetraEther indeX of 86 carbon atoms (TEX86) temperature proxy is widely used in reconstructions of past sea surface temperature. Most current calibrations are based on surface sediment distributions of the glycerol dialkyl glycerol tetraether lipids (GDGTs) that comprise TEX86 and assume that these GDGTs are exported from the upper mixed layer. However, GDGT export from deeper waters could impact sedimentary GDGT distributions and therefore TEX86 paleothermometry. Here we examine GDGT distributions in suspended particulate matter (SPM) and underlying sediments collected from the Southeast Atlantic Ocean. Our results reveal different GDGT distributions - specifically the ratio between GDGTs bearing 2 vs. 3 cyclopentyl moieties, [2/3] ratios - between surface, subsurface (>50-200 m) and deep water (>200 m) SPM, which suggests the occurrence of in situ (deep) production that is not apparent when considering TEX86. The GDGT distributions in sediments match those of subsurface waters rather than surface waters, suggesting that they have not been preferentially derived from the upper mixed layer; this is consistent with GDGT abundances being highest in shallow subsurface SPM (˜100 to 200 m). It remains unclear what governs the different [2/3] ratios throughout the water column, but it is likely related to a combination of temperature and thaumarchaeotal community structure.