COMPARISON OF MEASUREMENTS OF ABSORBED DOSE TO WATER USING A WATER CALORIMETER AND IONIZATION CHAMBERS FOR CLINICAL RADIOTHERAPY PHOTON AND ELECTRON BEAMS

With the development of the water calorimeter direct measurement of absorbed dose in water becomes possible. This could lead to the establishment of an absorbed dose rather than an exposure related standard for ionization chambers for high energy electrons and photons. In changing to an absorbed dose standard it is necessary to investigate the effect of different parameters, among which are the energy dependence, the air volume, wall thickness and material of the chamber. The effect of these parameters is experimentally studied and presented for several commercially available chambers and one experimental chamber, for photons up to 25 MV and electrons up to 20 MeV, using a water calorimeter as the absorbed dose standard and the most recent formalism to calculate the absorbed dose with ion chambers.^ For electron beams, the dose measured with the calorimeter was 1% lower than the dose calculated with the chambers, independent of beam energy and chamber.^ For photon beams, the absorbed dose measured with the calorimeter was 3.8% higher than the absorbed dose calculated from the chamber readings. Such differences were found to be chamber and energy independent.^ The results for the photons were found to be statistically different from the results with the electron beams. Such difference could not be attributed to a difference in the calorimeter response. ^

Reduction of setup uncertainties and tumor motion in the radiation therapy of lung tumors

Because the goal of radiation therapy is to deliver a lethal dose to the tumor, accurate information on the location of the tumor needs to be known. Margins are placed around the tumor to account for variations in the daily position of the tumor. If tumor motion and patient setup uncertainties can be reduced, margins that account for such uncertainties in tumor location in can be reduced allowing dose escalation, which in turn could potentially improve survival rates. ^ In the first part of this study, we monitor the location of fiducials implanted in the periphery of lung tumors to determine the extent of non-gated and gated fiducial motion, and to quantify patient setup uncertainties. In the second part we determine where the tumor is when different methods of image-guided patient setup and respiratory gating are employed. In the final part we develop, validate, and implement a technique in which patient setup uncertainties are reduced by aligning patients based upon fiducial locations in projection images. ^ Results from the first part indicate that respiratory gating reduces fiducial motion relative to motion during normal respiration and setup uncertainties when the patients were aligned each day using externally placed skin marks are large. The results from the second part indicate that current margins that account for setup uncertainty and tumor motion result in less than 2% of the tumor outside of the planning target volume (PTV) when the patient is aligned using skin marks. In addition, we found that if respiratory gating is going to be used, it is most effective if used in conjunction with image-guided patient setup. From the third part, we successfully developed, validated, and implemented on a patient a technique for aligning a moving target prior to treatment to reduce the uncertainties in tumor location. ^ In conclusion, setup uncertainties and tumor motion are a significant problem when treating tumors located within the thoracic region. Image-guided patient setup in conjunction with treatment delivery using respiratory gating reduces these uncertainties in tumor locations. In doing so, margins around the tumor used to generate the PTV can be reduced, which may allow for dose escalation to the tumor. ^

EVALUATION OF THE EFFECTIVENESS OF ANISOTROPIC ANALYTICAL ALGORITHM IN FLATTENED AND FLATTENING-FILTER-FREE BEAMS FOR HIGH ENERGY LUNG DOSE DELIVERY USING THE RADIOLOGICAL PHYSICS CENTER LUNG PHANTOM

The effectiveness of the Anisotropic Analytical Algorithm (AAA) implemented in the Eclipse treatment planning system (TPS) was evaluated using theRadiologicalPhysicsCenteranthropomorphic lung phantom using both flattened and flattening-filter-free high energy beams. Radiation treatment plans were developed following the Radiation Therapy Oncology Group and theRadiologicalPhysicsCenterguidelines for lung treatment using Stereotactic Radiation Body Therapy. The tumor was covered such that at least 95% of Planning Target Volume (PTV) received 100% of the prescribed dose while ensuring that normal tissue constraints were followed as well. Calculated doses were exported from the Eclipse TPS and compared with the experimental data as measured using thermoluminescence detectors (TLD) and radiochromic films that were placed inside the phantom. The results demonstrate that the AAA superposition-convolution algorithm is able to calculate SBRT treatment plans with all clinically used photon beams in the range from 6 MV to 18 MV. The measured dose distribution showed a good agreement with the calculated distribution using clinically acceptable criteria of ±5% dose or 3mm distance to agreement. These results show that in a heterogeneous environment a 3D pencil beam superposition-convolution algorithms with Monte Carlo pre-calculated scatter kernels, such as AAA, are able to reliably calculate dose, accounting for increased lateral scattering due to the loss of electronic equilibrium in low density medium. The data for high energy plans (15 MV and 18 MV) showed very good tumor coverage in contrast to findings by other investigators for less sophisticated dose calculation algorithms, which demonstrated less than expected tumor doses and generally worse tumor coverage for high energy plans compared to 6MV plans. This demonstrates that the modern superposition-convolution AAA algorithm is a significant improvement over previous algorithms and is able to calculate doses accurately for SBRT treatment plans in the highly heterogeneous environment of the thorax for both lower (≤12 MV) and higher (greater than 12 MV) beam energies.

Evaluation of Polymer Gel Dosimeters for Measurements of Dose and LET in Proton Beams

This project assessed the effectiveness of polymer gel dosimeters as tools for measuring the dose deposited by and LET of a proton beam. A total of three BANG® dosimeter formulations were evaluated: BANG®-3-Pro-2 BANGkits™ for dose measurement and two BANG®-3 variants, the LET-Baseline and LET-Meter dosimeters, for LET measurement. All dosimeters were read out using an OCT scanner. The basic characteristics of the BANGkits™ were assessed in a series of photon and electron irradiations. The dose-response relationship was found to be sigmoidal with a threshold for response of approximately 15 cGy. The active region of the dosimeter, the volume in which dosimeter response is not inhibited by oxygen, was found to make up roughly one fourth of the total dosimeter volume. Delivering a dose across multiple fractions was found to yield a greater response than delivering the same dose in a single irradiation. The dosimeter was found to accurately measure a dose distribution produced by overlapping photon fields, yielding gamma pass rates of 95.4% and 93.1% from two planar gamma analyses. Proton irradiations were performed for measurements of proton dose and LET. Initial irradiations performed through the side of a dosimeter led to OCT artifacts. Gamma pass rates of 85.7% and 89.9% were observed in two planar gamma analyses. In irradiations performed through the base of a dosimeter, gel response was found to increase with height in the dosimeter, even in areas of constant dose. After a correction was applied, gamma pass rates of 94.6% and 99.3% were observed in two planar gamma analyses. Absolute dose measurements were substantially higher (33%-100%) than the delivered doses for proton irradiations. Issues encountered while calibrating the LET-Meter gel restricted analysis of the LET measurement data to the SOBP of a proton beam. LET-Meter overresponse was found to increase linearly with track-average LET across the LET range that could be investigated (1.5 keV/micron – 3.5 keV/micron).

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

Comparison of Tumor Shrinkage and Cumulative Dose Distribution for Lung Cancers

Background: The physical characteristic of protons is that they deliver most of their radiation dose to the target volume and deliver no dose to the normal tissue distal to the tumor. Previously, numerous studies have shown unique advantages of proton therapy over intensity-modulated radiation therapy (IMRT) in conforming dose to the tumor and sparing dose to the surrounding normal tissues and the critical structures in many clinical sites. However, proton therapy is known to be more sensitive to treatment uncertainties such as inter- and intra-fractional variations in patient anatomy. To date, no study has clearly demonstrated the effectiveness of proton therapy compared with the conventional IMRT under the consideration of both respiratory motion and tumor shrinkage in non-small cell lung cancer (NSCLC) patients. Purpose: This thesis investigated two questions for establishing a clinically relevant comparison of the two different modalities (IMRT and proton therapy). The first question was whether or not there are any differences in tumor shrinkage between patients randomized to IMRT versus passively scattered proton therapy (PSPT). Tumor shrinkage is considered a standard measure of radiation therapy response that has been widely used to gauge a short-term progression of radiation therapy. The second question was whether or not there are any differences between the planned dose and 5D dose under the influence of inter- and intra-fractional variations in the patient anatomy for both modalities. Methods: A total of 45 patients (25 IMRT patients and 20 PSPT patients) were used to quantify the tumor shrinkage in terms of the change of the primary gross tumor volume (GTVp). All patients were randomized to receive either IMRT or PSPT for NSCLC. Treatment planning goals were identical for both groups. All patients received 5 to 8 weekly repeated 4-dimensional computed tomography (4DCT) scans during the course of radiation treatments. The original GTVp contours were propagated to T50 of weekly 4DCT images using deformable image registration and their absolute volumes were measured. Statistical analysis was performed to compare the distribution of tumor shrinkage between the two population groups. In order to investigate the difference between the planned dose and the 5D dose with consideration of both breathing motion and anatomical change, we re-calculated new dose distributions at every phase of the breathing cycle for all available weekly 4DCT data sets which resulted 50 to 80 individual dose calculations for each of the 7 patients presented in this thesis. The newly calculated dose distributions were then deformed and accumulated to T50 of the planning 4DCT for comparison with the planned dose distribution. Results: At the end of the treatment, both IMRT and PSPT groups showed mean tumor volume reductions of 23.6% ( 19.2%) and 20.9% ( 17.0 %) respectively. Moreover, the mean difference in tumor shrinkage between two groups is 3% along with the corresponding 95% confidence interval, [-8%, 14%]. The rate of tumor shrinkage was highly correlated with the initial tumor volume size. For the planning dose and 5D dose comparison study, all 7 patients showed a mean difference of 1 % in terms of target coverage for both IMRT and PSPT treatment plans. Conclusions: The results of the tumor shrinkage investigation showed no statistically significant difference in tumor shrinkage between the IMRT and PSPT patients, and the tumor shrinkage between the two modalities is similar based on the 95% confidence interval. From the pilot study of comparing the planned dose with the 5D dose, we found the difference to be only 1%. Overall impression of the two modalities in terms of treatment response as measured by the tumor shrinkage and 5D dose under the influence of anatomical change that were designed under the same protocol (i.e. randomized trial) showed similar result.

Ice and snow measurements during Nathaniel B. Palmer cruise NBP0709 (SIMBA experiment)

The Sea Ice Mass Balance in the Antarctic (SIMBA) experiment was conducted from the RVIB N.B. Palmer in September and October 2007 in the Bellingshausen Sea in an area recently experiencing considerable changes in both climate and sea ice cover. Snow and ice properties were observed at 3 short-term stations and a 27-day drift station (Ice Station Belgica, ISB) during the winter-spring transition. Repeat measurements were performed on sea ice and snow cover at 5 ISB sites, each having different physical characteristics, with mean ice (snow) thicknesses varying from 0.6 m (0.1 m) to 2.3 m (0.7 m). Ice cores retrieved every five days from 2 sites and measured for physical, biological, and chemical properties. Three ice mass-balance buoys (IMBs) provided continuous records of snow and ice thickness and temperature. Meteorological conditions changed from warm fronts with high winds and precipitation followed by cold and calm periods through four cycles during ISB. The snow cover regulated temperature flux and controlled the physical regime in which sea ice morphology changed. Level thin ice areas had little snow accumulation and experienced greater thermal fluctuations resulting in brine salinity and volume changes, and winter maximum thermodynamic growth of ~0.6 m in this region. Flooding and snow-ice formation occurred during cold spells in ice and snow of intermediate thickness. In contrast, little snow-ice formed in flooded areas with thicker ice and snow cover, instead nearly isothermal, highly permeable ice persisted. In spring, short-lived cold air episodes did not effectively penetrate the sea ice nor overcome the effect of ocean heat flux, thus favoring net ice thinning from bottom melt over ice thickening from snow-ice growth, in all cases. These warm ice conditions were consistent with regional remote sensing observations of earlier ice breakup and a shorter sea ice season, more recently observed in the Bellingshausen Sea.

Thorium and radiocarbon measurements on sediment cores from the Northseast Atlantic

Radiocarbon and 230Thexcess data from six NE Atlantic box cores are considered. The cores form a transect from the Porcupine Abyssal Plain over the East Thulean Rise to the southern end of Feni Drift. The chronology for the cores is established from bulk sediment carbonate radiocarbon data and reveals that sections exhibiting constant accumulation rates can be identified in all the cores, with rates of 3.0-3.5 cm/kyr on the plain through the Holocene and late Holocene rates of 4.3-6.6 cm/kyr elsewhere. Five out of the six cores show accumulations of more 230Thexcess than is produced in the overlying water column, with the greatest inventories (up to 225% of production) in the cores from the rise and drift. A size fraction comparison between two cores from the plain and rise reveals that the higher overall accumulation rates and 230Thexcess inventories in the off-plain cores are due to an increased fine (<5 µm) component fraction, whereas the flux of coarser material is similar to that received on the plain. This suggests that the higher fluxes of materials observed are physically (rather than biogeochemically) driven and also that drift formation has been continuously active in the late Holocene. Sections of all the cores where regular accumulation is defined by the radiocarbon data are modeled first by a linear radiocarbon age/depth model and second by a constant rain (230Thexcess)0 model prorated for the observed core inventories. These modeling approaches yield historical mass accumulation rate estimates which are generally in reasonable agreement (±30%), but the differences observed appear to be well organized in time rather than random.

Geophysical measurements on sediment cores from North Pond, Mid Atlantic Ridge

Die Porosität ist einer der wichtigsten gesteinsphysikalischen Parameter von Sedimenten. Daher werden in dieser Arbeit die Abhängigkeiten zwischen dem zentralen gesteinsphysikalischen Parameter Porosität und den Parametern Wärmeleitfähigkeit und spezifischer elektrischer Widerstand an Sedimenten von North Pond, eines Sedimentbeckens an der westlichen Flanke des Mittelatlantischen Rückens bei ca. 23°N, qualitativ und empirisch untersucht. Die Messung des spezifischen Widerstands wurde mittels eines miniaturisierten Wenner-Arrays mit einem Zentimeter Elektrodenabstand noch an Bord der FS Maria S. Merian auf der Reise MSM11/1 vom 10. Februar bis 13. März 2009 an Sedimentkernhälften durchgeführt, ebenso wie die Messung der Wärmeleitfähigkeit mit einer Nadelsonde. Die Messung der Gamma-Dichte an einem Multi Sensor Core Logger erfolgte im Rahmen dieser Arbeit und wurde wie die Messung der Porosität an Proben im Marum in Bremen durchgeführt. Die zugrundeliegende und zu untersuchende Beziehung zwischen spezifischem Widerstand und Porosität ist Archie's Law, die untersuchte Beziehung zwischen Wärmeleitfähigkeit und Porosität ist das Modell des geometrischen Mittels der Sedimentkomponenten. Die physikalische Beziehung zwischen der Porosität und den beiden Parametern spezifischer Widerstand und Wärmeleitfähigkeit konnte an den Sedimenten von North Pond bestätigt werden.

Tomographic measurements of O. umbonatus and N. truempyi during the PETM and ETM-2 at ODP Sites 208-1262 and 208-1263

Predicting the impact of ongoing anthropogenic CO2 emissions on calcifying marine organisms is complex, owing to the synergy between direct changes (acidification) and indirect changes through climate change (e.g., warming, changes in ocean circulation, and deoxygenation). Laboratory experiments, particularly on longer-lived organisms, tend to be too short to reveal the potential of organisms to acclimatize, adapt, or evolve and usually do not incorporate multiple stressors. We studied two examples of rapid carbon release in the geological record, Eocene Thermal Maximum 2 (~53.2 Ma) and the Paleocene Eocene Thermal Maximum (PETM, ~55.5 Ma), the best analogs over the last 65 Ma for future ocean acidification related to high atmospheric CO2 levels. We use benthic foraminifers, which suffered severe extinction during the PETM, as a model group. Using synchrotron radiation X-ray tomographic microscopy, we reconstruct the calcification response of survivor species and find, contrary to expectations, that calcification significantly increased during the PETM. In contrast, there was no significant response to the smaller Eocene Thermal Maximum 2, which was associated with a minor change in diversity only. These observations suggest that there is a response threshold for extinction and calcification response, while highlighting the utility of the geological record in helping constrain the sensitivity of biotic response to environmental change.

Field measurements on Atka Bay landfast sea ice in 2012

Basal melt of ice shelves may lead to an accumulation of disc-shaped ice platelets underneath nearby sea ice, to form a sub-ice platelet layer. Here we present the seasonal cycle of sea ice attached to the Ekström Ice Shelf, Antarctica, and the underlying platelet layer in 2012. Ice platelets emerged from the cavity and interacted with the fast-ice cover of Atka Bay as early as June. Episodic accumulations throughout winter and spring led to an average platelet-layer thickness of 4 m by December 2012, with local maxima of up to 10 m. The additional buoyancy partly prevented surface flooding and snow-ice formation, despite a thick snow cover. Subsequent thinning of the platelet layer from December onwards was associated with an inflow of warm surface water. The combination of model studies with observed fast-ice thickness revealed an average ice-volume fraction in the platelet layer of 0.25 +/- 0.1. We found that nearly half of the combined solid sea-ice and ice-platelet volume in this area is generated by heat transfer to the ocean rather than to the atmosphere. The total ice-platelet volume underlying Atka Bay fast ice was equivalent to more than one-fifth of the annual basal melt volume under the Ekström Ice Shelf.

Abundance, biomass, production, grazing, and biometric measurements of prokaryotes, nanoflagellates, ciliates, and dinoflagellates in three areas close to the Antarctic Peninsula in spring and summer 1995/96

Two cruises were carried out during the Austral spring-summer (November 1995 - January 1996: FRUELA 95, and January - February 1996: FRUELA 96), sampling in Bellingshausen Sea, western Bransfield Strait and Gerlache Strait. We investigated whether there were any spatial (among locations) or temporal (between cruises) differences in abundance and biomass of microbial heterotrophic and autotrophic assemblages. Changes in the concentration of chlorophyll a, prokaryotes, heterotrophic and phototrophic nanoflagellates abundance and biomass were followed in the above mentioned locations close to the Antarctic Peninsula. Parallel to these measurements we selected seven stations to determine grazing rates on prokaryotes by protists at a depth coincident with the depth of maximum chlorophyll a concentration. Measuring the disappearance of fluorescent minicells over 48 h assessed grazing by the protist community. From prokaryotes grazing rates, we estimated how much prokaryotic carbon was channeled to higher trophic levels (protists), and whether this prokaryotic carbon could maintain protists biomass and growth rates. In general higher values were reported for Gerlache Strait than for the other two areas. Differences between cruises were more evident for the oligotrophic areas in Bellingshausen Sea and Bransfield Strait than in Gerlache Strait (eutrophic area). Higher values for phototrophic (at least for chlorophyll a concentration) and abundance of all heterotrophic microbial populations were recorded in Bellingshausen Sea and Bransfield Strait during late spring - early summer (FRUELA 95) than in mid-summer (FRUELA 96). However, similar results for these variables were observed in Gerlache Strait as in spring-early summer as well as in mid-summer. Also, we found differences in grazing rates on prokaryotes among stations located in the three areas and between cruises. Thus, during late spring-early summer (FRUELA 95), the prokaryotic biomass consumed from the standing stock was higher in Bellingshausen Sea (26%/day) and Gerlache Strait (18-26%/day) than in Bransfield Strait (0.68-14%/day). During mid-summer (FRUELA 96) a different pattern was observed. The station located in Bellingshausen Sea showed higher values of prokaryotic biomass consumed (11%/day) than the one located in Gerlache Strait (2.3%/day). Assuming HNF as the main prokaryotic consumers, we estimated that the prokaryotic carbon consumed by heterotrophic nanoflagellates (HNF) barely covers their carbon requirements for growth. These results suggest that in Antarctic waters, HNF should feed in other carbon sources than prokaryotes.

Biogenic and lithogenic silica, silicic acid, production of biogenic silica, and irradiance measurements during the Ross Sea Bloom Project 1994-1996

The development of the seasonal phytoplankton bloom in the Ross Sea was studied during two cruises. The first, conducted in November-December 1994, investigated the initiation and rapid growth of the bloom, whereas the second (December 1995-January 1996) concentrated on the bloom's maximum biomass period and the subsequent decline in biomass. Central to the understanding of the controls of growth and the summer decline of the bloom is a quantitative assessment of the growth rate of phytoplankton. Growth rates were estimated over two time scales with different methods. The first estimated daily growth rates from isotropic incorporation under simulated in situ conditions, including 14C, 15N and 32Si uptake measurements combined with estimates of standing stocks of particulate organic carbon, nitrogen and biogenic silica. The second method used daily to weekly changes in biomass at selected locations, with net growth rates being estimated from changes in standing stocks of phytoplankton. In addition, growth rates were estimated in large-volume experiments under optimal irradiances. Growth rates showed distinct temporal patterns. Early in the growing season, short-term estimates suggested that growth rates of in situ assemblages were less than maximum (relative to the temperature-limited maximum) and were likely reduced due to low irradiance regimes encountered under the ice. Growth rates increased thereafter and appeared to reach their maximum as biomass approached the seasonal peak, but decreased markedly in late December. Differences between the major taxonomic groups present were also noted, especially from the isotopic tracer experiments. The haplophyte Phaeocystic antarctica was dominant in 1994 throughout the growing season, and it exhibited the greatest growth rates (mean 0.41/day) during spring. Diatom standing stocks were low early in the growing season, and growth rates averaged 0.100/day. In summer diatoms were more abundant, but their growth rates remained much lower (mean of 0.08/day) than the potential maximum. Understanding growth rate controls is essential to the development of predictive models of the carbon cycle and food webs in Antarctic waters.