Steric sea-level variations inferred from combined Topex/Poseidon altimetry and GRACE gravimetry

In this study, we propose to estimate the steric sea-level variations over a < 2-year period (April 2002 through December 2003) by combining global mean sea level (GMSL) based on Topex/ Poseidon (T/P) altimetry with time-variable geoid averaged over the oceans, as observed by the GRACE (Gravity Recovery and Climate Experiment) satellite. In effect, altimetry-derived GMSL changes results from two contributions: Steric (thermal plus salinity) effects due to sea water density change and ocean mass change due to water exchange with atmosphere and continents. On the other hand, GRACE data over the oceans provide the ocean mass change component only. The paper first discusses the corrections to apply to the GRACE data. Then the steric contribution to the GMSL is estimated using GRACE and T/P data. Comparison with available thermal expansion based on in situ hydrographic data is performed.

Effects of ocean acidification on growth and physiology of Ulva lactuca (Chlorophyta) in a rockpool-scenario, link to supplementary data

Rising atmospheric CO2-concentrations will have severe consequences for a variety of biological processes. We investigated the responses of the green alga Ulva lactuca (Linnaeus) to rising CO2-concentrations in a rockpool scenario. U. lactuca was cultured under aeraton with air containing either preindustrial pCO2 (280µatm) or for the end of the 21st century predicted (700µatm) pCO2 for 31 days. We addressed the following question: Will elevated CO2-concentrations affect photosynthesis (net photosynthesis, rETR(max), Fv/Fm, pigment composition) and growth of U. lactuca in rockpools with limited water exchange? Two phases of the experiment were distinguished: In the initial phase (day 1-4) the Seawater Carbonate System (SWCS) of the culture medium could be adjusted to the selected atmospheric pCO2 condition by continuous aeration with target pCO2 values. In the second phase (day 4-31) the SWCS was largely determined by the metabolism of the growing U. lactuca biomass. In the initial phase, Fv/Fm and rETR(max) were only slightly elevated at high CO2-concentrations whereas growth was significantly enhanced. After 31 days the Chl a content of the thalli was significantly lower under future conditions and the photosynthesis of thalli grown under preindustrial conditions was not dependent on external carbonic anhydrase. Biomass increased significantly at high CO2-concentrations. At low CO2-concentrations most adult thalli disintegrated between day 14 and 21, whereas at high CO2-concentrations most thalli remained integer until day 31. Thallus disintegration at low CO2-concentrations was mirrored in a drastic decline in seawater DIC and HCO3-. Accordingly, the SWCS differed significantly between the treatments. Our results indicated a slight enhancement of photosynthetic performance and significantly elevated growth of U. lactuca at future CO2-concentrations. The accelerated thallus disintegration at high CO2-concentrations under conditions of limited water exchange indicates additional CO2 effects on the life cycle of U. lactuca when living in rockpools.

(Table 2) Stable isotopes in foraminifers at DSDP Leg 78

Analyses of stable isotopes of monospecific planktonic foraminifers (G. quadrilobatus group) and monogeneric benthic foraminifers (Cibicidoides spp.) from late Neogene Atlantic Site 502 and Pacific Site 503 were conducted in order to determine the paleoceanographic changes resulting from the late Neogene uplift of the Panama Isthmus and from climatic cooling. In general, results at each site are similar to those from previous studies for the late Miocene and late Pliocene time interval, documenting the late Miocene (6 Ma) shift in carbon isotopes and the inferred growth of permanent Northern Hemisphere continental ice sheets beginning about 3.2 Ma. Comparison of Atlantic-Pacific planktonic-benthic isotope data for four stratigraphic intervals (~6-8, ~5-6, ~3-5, and ~2-3 Ma) suggests that increasing isolation of Atlantic and Pacific low-latitude waters may be related to the emergence of the Panama Isthmus. The contrast between Atlantic and Pacific benthic foraminiferal d13C increased in two steps from 0.60 per mil to 1 per mil (the modern contrast) at about 6 Ma and 3 Ma. The first increase (0.15 per mil) may represent the end of previously limited deep-water communication between the Atlantic and Pacific at the present location of Panama. The second increase (0.25 per mil) may be due to increased production of North Atlantic Deep Water. This probably reflects the development of modern deep-sea circulation. The d18O of planktonic foraminifers begins to increase in Atlantic Site 502 at 4.2 Ma and may reflect the increasing salinity of the North Atlantic Ocean arising from diminishing surface-water exchange across Panama. This increase is clearly shown by contrasting the d18O of Atlantic and Pacific planktonic foraminifers, as well as the d18O of planktonic and benthic foraminifers at Site 502. This inferred increase in surface-water salinity begins at the time of increasing provinciality of Atlantic and Pacific planktonic foraminifers.

Seawater Nd and Pb isotopic record in central Pacific and Caribbean sediments

The shoaling and final closure of the Central American Seaway (CAS) resulted in a major change of the global ocean circulation and has been suggested as an essential driver for strengthening of Atlantic Meridional Overturning Circulation (AMOC). The exact timing of CAS closure is key to interpreting its importance. Here we present a reconstruction of deep and intermediate water Nd and Pb isotope compositions obtained from fossil fish teeth and the authigenic coatings of planktonic foraminifera in the eastern equatorial Pacific (Ocean Drilling Program (ODP) Site 1241) and the Caribbean (ODP Sites 998, 999, and 1000) covering the final stages of CAS closure between 5.6 and 2.2 Ma. The data for the Pacific site indicate no significant Atlantic/Caribbean influence over this entire period. The Caribbean sites show a continuous trend to less radiogenic Nd isotope compositions during the Pliocene, consistent with an enhancement of Upper North Atlantic Deep Water (UNADW) inflow and a strengthening of the AMOC. Superimposed onto this long-term trend, shorter-term changes of intermediate Caribbean Nd isotope signatures approached more UNADW-like values during intervals when published reconstructions of seawater salinity suggested complete closure of the CAS. The data imply that significant deep water exchange with the Pacific essentially stopped by 7 Ma and that shallow exchange, which still occurred at least periodically until approximately 2.5 Ma, may have been linked to the strength of the AMOC but did not have any direct effect on the intermediate and deep Caribbean Nd isotope signatures through mixing with Pacific waters.

(Table 1) Content of sand fraction components in sediments of the interglacial MIS 7 optimum from Core LV28-42-4

Distribution of diatoms, radiolarians, planktonic and benthic foraminifers, and sediment components in fraction >0.125 mm was analyzed in a core obtained from the central Sea of Okhotsk within frameworks of the Russian-German KOMEX Project. The core section characterizes the period 190-350 ka, which corresponds to marine-isotopic stages (MIS) 7 to 10. During glacial MIS 10 and MIS 8, the basin accumulated terrigenous material lacking microfossils or containing them in low abundance, which reflects, along with their composition, heavy sea-ice conditions, suppressed bioproductivity, and bottom environment aggressive toward calcium carbonate. Interglacial MIS 9 was characterized by elevated bioproductivity with accumulation of diatomaceous ooze during the climatic optimum (328 to 320 ka). Water exchange with the Pacific was maximal from 328 to 324 ka ago. Environment became moderate and close to the present-day one at the end of the optimum exhibiting possible existence of a dichothermal layer with substantial amounts of surface Pacific water still flowing into the basin. Similar to interglacial MIS 5e and MIS 1, ''old'' Pacific water determined near-bottom environment in the central Sea of Okhotsk during that period, although influx of terrigenous material was higher, probably reflecting more humid climate of the region. Slight warming marked the terminal MIS 8 (approximately 260 ka ago). Paleoceanographic situation during the interglacial MIS 7 was highly variable: from warm-water to almost glacial. The main climatic optimum of MIS 7 occurred within 220-210 ka, when subsurface stratification increased and the dichothermal layer developed. Bottom environment during the studied time interval, except for the optimum of interglacial MIS 9, resembled those characteristic of glacial periods: actively formed ''young'' Okhotsk water displaced ''old'' Pacific deep water.

Synthesis of gadolinium-labeled shell-crosslinked nanoparticles for magnetic resonance imaging applications

Shell-crosslinked knedel-like nanoparticles (SCKs; knedel is a Polish term for dumplings) were derivatized with gadolinium Shell chelates and studied as robust magnetic-resonance-imaging-active structures with hydrodynamic diameters of 40 +/- 3 nm. SCKs possessing an amphiphilic core-shell morphology were produced from the aqueous assembly of diblock copolymers of poly(acrylic acid) (PAA) and poly(methyl acrylate) (PMA), PAA(52)-b-PMA(128), and subsequent covalent crosslinking by amidation upon reaction with 2,2'-(ethylenedioxy)bis(ethylamine) throughout the shell layer. The properties of these materials, including non-toxicity towards mammalian cells, non-immunogenicity within mice, and capability for polyvalent targeting, make them ideal candidates for utilization within biological systems. The synthesis of SCKs derivatized with Gd-III and designed for potential use as a unique nanometer-scale contrast agent for MRI applications is described herein. Utilization of an amino-functionalized diethylenetriaminepentaacetic acid-Gd analogue allowed for direct covalent conjugation throughout the hydrophilic shell layer of the SCKs and served to increase the rotational correlation lifetime of the Gd. In addition, the highly hydrated nature of the shell layer in which the Gd was located allowed for rapid water exchange; thus, the resulting material demonstrated large ionic relaxivities (39 s(-1) mM(-1)) in an applied magnetic field of 0.47 T at 40 degrees C and, as a result of the large loading capacity of the material, also demonstrated high molecular relaxivities (20 000 s(-1) mM(-1)).

Spatial patterns of aerobic and anaerobic mineralization rates and oxygen penetration dynamics in coral reef sediments

Oxygen consumption rates (OCR), aerobic mineralization and sulfate reduction rates (SRR) were studied in the permeable carbonate reef sediments of Heron Reef, Australia. We selected 4 stations with different hydrodynamic regimes for this study. In situ oxygen penetration into the sediments was measured with an autonomous microsensor profiler. Areal OCR were quantified from the measured oxygen penetration depth and volumetric OCR. Oxygen penetration and dynamics (median penetration depths at the 4 stations ranged between 0.3 and 2.2 cm), OCR (median 57 to 196 mmol C m(-2) d(-1)), aerobic mineralization (median 24 to 176 mmol C m(-2) d(-1)) and SRR (median 9 to 42 mmol C m(-2) d(-1)) were highly variable between sites. The supply of oxygen by pore water advection was a major cause for high mineralization rates by stimulating aerobic mineralization at all sites. However, estimated bottom water filtration rates could not explain the differences in volumetric OCR and SRR between the 4 stations. This suggests that local mineralization rates are additionally controlled by factors other than current driven pore water advection, e.g. by the distribution of the benthic fauna or by local differences in labile organic carbon supply from sources such as benthic photosynthesis. Carbon mineralization rates were among the highest reported for coral reef sediments, stressing the role of these sediments in the functioning of the reef ecosystem.

Top down control in a relatively pristine seagrass ecosystem

The loss of large-bodied herbivores and/or top predators has been associated with large-scale changes in terrestrial, freshwater, and marine ecosystems around the world. Understanding the consequences of these declines has been hampered by a lack of studies in relatively pristine systems. To fill this gap, I investigated the dynamics of the relatively pristine seagrass ecosystem of Shark Bay, Australia. I began by examining the seagrass species distributions, stoichiometry, and patterns of nutrient limitation across the whole of Shark Bay. Large areas were N-limited, P-limited, or limited by factors other than nutrients. Phosphorus-limitation was centered in areas of restricted water exchange with the ocean. Nutrient content of seagrasses varied seasonally, but the strength of seasonal responses were species-specific. Using a cafeteria-style experiment, I found that fast-growing seagrass species, which had higher nutrient content experienced higher rates of herbivory than slow-growing species that are dominant in the bay but have low nutrient content. Although removal rates correlated well with nutrient content at a broad scale, within fast-growing species removal rates were not closely tied to N or P content. Using a combination of stable isotope analysis and animal borne video, I found that green turtles (Chelonia mydas)—one of the most abundant large-bodied herbivores in Shark Bay—appear to assimilate little energy from seagrasses at the population level. There was, however, evidence of individual specialization in turtle diets with some individuals foraging largely on seagrasses and others feeding primarily on macroalgae and gelatinous macroplankton. Finally, I used exclusion cages, to examine whether predation-sensitive habitat shifts by megagrazers (green turtles, dugongs) transmitted a behavior-mediated trophic cascade (BMTC) between sharks and seagrasses. In general, data were consistent with predictions of a behavior-mediated trophic cascade. Megaherbivore impacts on seagrasses were large only in the microhabitat where megaherbivores congregate to reduce predation risk. My study highlights the importance of large herbivores in structuring seagrass communities and, more generally, suggests that roving top predators likely are important in structuring communities—and possibly ecosystems—through non-consumptive pathways.

(Table 2) Relative abundance of planktonic foraminifera species in Pliocene - Early Quaternary sediments from ODP Hole 111-677A

Planktonic foraminifera from Pliocene - Early Quaternary sediments of ODP Hole 111-677A were studied in detail. It was shown that the majority of detected zonal taxa are reliable biostratigraphic reference points. Between 30 and 210 m in the core zones of planktonic foraminifera from PL1b to Pt1 (according to the W.A. Berggren scale) were distinguished. Changes of planktonic foraminifera complexes from sediments of Hole 111-677A are closely associated with climate-controlled development of surface water masses of the Eastern Equatorial Pacific during 4.6-0.65 million years ago. Sharp decrease in equatorial-tropical species about 3.4 million years ago correlated with cessation of surface water exchange between tropical regions of the Pacific and Atlantic oceans due to formation of the Central American isthmus. The paleotemperature method of M.S Barash was used for reconstructing surface temperatures. Maximum temperatures were reconstructed in late Early Pliocene (26.4°C) and in Late Pliocene (26.6°C) and minimum ones - in the beginning of Early Pliocene (18.4°C), in the middle of Late Pliocene (19.6°C). Cold events occurred: 4.6-4.3, 2.8-2.5, and 1.7-1.2 million years ago, and warm: 4.3, 4.18-3.4, 2.5-2.3, and 1 million years ago. In general, the middle of Early Pliocene, the middle of late Pliocene and early Pleistocene are characterized by cold-water conditions, and the end of Early and the end of Late Pliocene - by warm-water conditions.

Distribution of planktonic and benthic foraminifers in bottom sediments from Core LV28-42-5, Sea of Okhotsk

Distribution of diatoms and planktonic and benthic foraminifers, as well as correlation of components of sandy grain size fraction were studied in the Quaternary sediment core LV28-42-5 (720 cm long) col¬lected on the southeastern slope (1045 m depth) of the Institute of Oceanology Rise, Sea of Okhotsk. This study allowed to reconstruct principle features of paleoceanographic evolution. In the course of penultimate and last continental glaciations (isotope stages 6 and 4-2) and during the later period of the last interglacial (substages 5.d-5.a) the following conditions were characteristic of this area: low temperatures of surface water, terrigenous sediment accumulation including coarse grained ice-rafted material, minimum bioproductivity and microfossil content in sediments, low sea level, reduced water exchange with the ocean, low position of old deep Pacific water. During the interglacial optimum (substage 5.e), as well as in the last deglaciation and Holocene (stage 1) water temperature and bioproductivity increased, sea level rose, and active surface water exchange between the Sea of Okhotsk and the Pacific Ocean and the Sea of Japan took place. This resulted in intensive inflow of the old deep Pacific water into the Sea of Okhotsk and elevation of its upper boundary by few hundred meters. During the later intervals of these warm periods a dichothermal structure of the upper water layer formed and diatom oozes accumulated.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

Radiocarbon dating, oxygen isotopes, and sea-level reconstruction of sediment core GeoB5844-2 from the northern end of the Red Sea

Based on a radiocarbon and paleomagnetically dated sediment record from the northern Red Sea and the exceptional sensitivity of the regional changes in the oxygen isotope composition of sea water to the sea-level-dependent water exchange with the Indian Ocean, we provide a new global sea-level reconstruction spanning the last glacial period. The sea-level record has been extracted from the temperature-corrected benthic stable oxygen isotopes using coral-based sea-level data as constraints for the sea-level/oxygen isotope relationship. Although, the general features of this millennial-scale sea-level records have strong similarities to the rather symmetric and gradual Southern Hemisphere climate patterns, we observe, in constrast to previous findings, pronounced sea level rises of up to 25 m to generally correspond with Northern Hemisphere warmings as recorded in Greenland ice-core interstadial intervals whereas sea-level lowstands mostly occur during cold phases. Corroborated by CLIMBER-2 model results, the close connection of millennial-scale sea-level changes to Northern Hemisphere temperature variations indicates a primary climatic control on the mass balance of the major Northern Hemisphere ice sheets and does not require a considerable Antarctic contribution.