Understanding the Influence of Photon Energy on 6MV Non Reference Dosimetry and CT Dosimetry using TLD and OSLD

Measurement of the absorbed dose from ionizing radiation in medical applications is an essential component to providing safe and reproducible patient care. There are a wide variety of tools available for measuring radiation dose; this work focuses on the characterization of two common, solid-state dosimeters in medical applications: thermoluminescent dosimeters (TLD) and optically stimulated luminescent dosimeters (OSLD). There were two main objectives to this work. The first objective was to evaluate the energy dependence of TLD and OSLD for non-reference measurement conditions in a radiotherapy environment. The second objective was to fully characterize the OSLD nanoDot in a CT environment, and to provide validated calibration procedures for CT dose measurement using OSLD. Current protocols for dose measurement using TLD and OSLD generally assume a constant photon energy spectrum within a nominal beam energy regardless of measurement location, tissue composition, or changes in beam parameters. Variations in the energy spectrum of therapeutic photon beams may impact the response of TLD and OSLD and could thereby result in an incorrect measure of dose unless these differences are accounted for. In this work, we used a Monte Carlo based model to simulate variations in the photon energy spectra of a Varian 6MV beam; then evaluated the impact of the perturbations in energy spectra on the response of both TLD and OSLD using Burlin Cavity Theory. Energy response correction factors were determined for a range of conditions and compared to measured correction factors with good agreement. When using OSLD for dose measurement in a diagnostic imaging environment, photon energy spectra are often referenced to a therapy-energy or orthovoltage photon beam – commonly 250kVp, Co-60, or even 6MV, where the spectra are substantially different. Appropriate calibration techniques specifically for the OSLD nanoDot in a CT environment have not been presented in the literature; furthermore the dependence of the energy response of the calibration energy has not been emphasized. The results of this work include detailed calibration procedures for CT dosimetry using OSLD, and a full characterization of this dosimetry system in a low-dose, low-energy setting.

Evaluation of a wellness program at a large private sector petrochemical company in Texas

The study purpose was to analyze the effects Integrated Health Solutions (IHS), an employee wellness program that has been implemented for one year on the corporate campus of a major private sector petrochemical company in Houston, TX, has on employee health. ^ Chronic diseases are the leading causes of morbidity and mortality in the United States and are the most preventable of all health problems. The costs of chronic diseases in the working-age adult population include not only health problems and a decrease in quality of life, but also an increase the cost of health care and costs to businesses and employers, both directly and indirectly. These emerging costs to employers as well as the fact that adults now spend the majority of waking hours at the office have increased the interest in worksite health promotion programs that address many of the behavioral factors that lead to chronic conditions. Therefore, implementing and evaluating programs that are aimed at promoting health and decreasing the prevalence of chronic diseases at worksites is very important. ^ Data came from existing data that were collected by IHS staff during employee biometric screenings at the company in 2010 and 2011. Data from employees who participated in screenings in both 2010 and 2011 were grouped into a cohort by IHS staff. ^ One-tailed t-tests were conducted to determine if there were significant improvements in the biometric measures of body fat percentage, BMI, waist circumference, systolic and diastolic blood pressures, total, HDL, and LDL cholesterol levels, triglycerides, blood glucose levels, and cardiac risk ratios. Sensitivity analysis was conducted to determine if there were differences in program outcomes when stratified by age, gender, job type, and time between screenings. ^ Mean differences for the variables from 2010 to 2011 were small and not always in the desired direction for health improvement indicators. Through conducting t-tests, it was found that there were significant improvements in HDL, cardiac risk ratio, and glucose levels. There were significant increases in cholesterol, LDL, and diastolic blood pressures. For the IHS program, it appears that gender, job type, and time between screenings were possible modifiers of program effectiveness. When program outcome measures were stratified by these factors, results suggest that corporate employees had better outcomes than field employees, males had better outcomes overall than females, and more positive program effects were seen for employees with less time between their two screenings. ^ Recommendations for the program based on the results include ensuring validity of instruments and initial and periodic training of measurement procedures and equipment handling, using normative data or benchmarks to decrease chances for biased estimates of program effectiveness, measuring behaviors as well as biometric and physiologic statuses and changes, and collecting level of engagement data.^

CHARACTERIZATION OF THE COUNT RATE PERFORMANCE AND EVALUATION OF THE EFFECTS OF HIGH COUNT RATES ON MODERN GAMMA CAMERAS

CHARACTERIZATION OF THE COUNT RATE PERFORMANCE AND EVALUATION OF THE EFFECTS OF HIGH COUNT RATES ON MODERN GAMMA CAMERAS Michael Stephen Silosky, B.S. Supervisory Professor: S. Cheenu Kappadath, Ph.D. Evaluation of count rate performance (CRP) is an integral component of gamma camera quality assurance and measurement of system dead time (τ) is important for quantitative SPECT. The CRP of three modern gamma cameras was characterized using established methods (Decay and Dual Source) under a variety of experimental conditions. For the Decay method, input count rate was plotted against observed count rate and fit to the paralyzable detector model (PDM) to estimate τ (Rates method). A novel expression for observed counts as a function of measurement time interval was derived and the observed counts were fit to this expression to estimate τ (Counts method). Correlation and Bland-Altman analysis were performed to assess agreement in estimates of τ between methods. The dependencies of τ on energy window definition and incident energy spectrum were characterized. The Dual Source method was also used to estimate τ and its agreement with the Decay method under identical conditions and the effects of total activity and the ratio of source activities were investigated. Additionally, the effects of count rate on several performance metrics were evaluated. The CRP curves for each system agreed with the PDM at low count rates but deviated substantially at high count rates. Estimates of τ for the paralyzable portion of the CRP curves using the Rates and Counts methods were highly correlated (r=0.999) but with a small (~6%) difference. No significant difference was observed between the highly correlated estimates of τ using the Decay or Dual Source methods under identical experimental conditions (r=0.996). Estimates of τ increased as a power-law function with decreasing ratio of counts in the photopeak to the total counts and linearly with decreasing spectral effective energy. Dual Source method estimates of τ varied as a quadratic with the ratio of the single source to combined source activities and linearly with total activity used across a large range. Image uniformity, spatial resolution, and energy resolution degraded linearly with count rate and image distorting effects were observed. Guidelines for CRP testing and a possible method for the correction of count rate losses for clinical images have been proposed.

AN EVALUATION OF THE CONSISTENCY OF IMRT PATIENT-SPECIFIC QA TECHNIQUES

To ensure the integrity of an intensity modulated radiation therapy (IMRT) treatment, each plan must be validated through a measurement-based quality assurance (QA) procedure, known as patient specific IMRT QA. Many methods of measurement and analysis have evolved for this QA. There is not a standard among clinical institutions, and many devices and action levels are used. Since the acceptance criteria determines if the dosimetric tools’ output passes the patient plan, it is important to see how these parameters influence the performance of the QA device. While analyzing the results of IMRT QA, it is important to understand the variability in the measurements. Due to the different form factors of the many QA methods, this reproducibility can be device dependent. These questions of patient-specific IMRT QA reproducibility and performance were investigated across five dosimeter systems: a helical diode array, radiographic film, ion chamber, diode array (AP field-by-field, AP composite, and rotational composite), and an in-house designed multiple ion chamber phantom. The reproducibility was gauged for each device by comparing the coefficients of variation (CV) across six patient plans. The performance of each device was determined by comparing each one’s ability to accurately label a plan as acceptable or unacceptable compared to a gold standard. All methods demonstrated a CV of less than 4%. Film proved to have the highest variability in QA measurement, likely due to the high level of user involvement in the readout and analysis. This is further shown by how the setup contributed more variation than the readout and analysis for all of the methods, except film. When evaluated for ability to correctly label acceptable and unacceptable plans, two distinct performance groups emerged with the helical diode array, AP composite diode array, film, and ion chamber in the better group; and the rotational composite and AP field-by-field diode array in the poorer group. Additionally, optimal threshold cutoffs were determined for each of the dosimetry systems. These findings, combined with practical considerations for factors such as labor and cost, can aid a clinic in its choice of an effective and safe patient-specific IMRT QA implementation.

Geochemical composition of the Trincheira Complex in the southwestern Amazon craton, Brazil

We document the first-known Mesoproterozoic ophiolite from the southwestern part of the Amazon craton, corresponding to the Trincheira Complex of Calymmian age, and propose a tectonic model that explains many previously enigmatic features of the Precambrian history of this key craton, and discuss its role in the reconstruction of the Columbia supercontinent. The complex comprises extrusive rocks (fine-grained amphibolites derived from massive and pillowed basalts), mafic-ultramafic intrusive rocks, chert, banded iron formation (BIFs), pelites, psammitic and a smaller proportion of calc-silicate rocks. This sequence was deformed, metasomatized and metamorphosed during the development of the Alto Guaporé Belt, a Mesoproterozoic accretionary orogen. The rocks were deformed by a single tectonic event, which included isoclinal folding and metamorphism of the granulite-amphibolite facies. Layered magmatic structures were preserved in areas of low strain, including amygdaloidal and cumulate structures. Metamorphism was pervasive and reached temperatures of 780-853°C in mafic granulites and 680-720°C in amphibolites under an overall pressure of 6.8 kbar. The geochemical composition of the extrusive and intrusive rocks indicates that all noncumulus mafic-ultramafic rocks are tholeiitic basalts. The mafic-ultramafic rocks display moderately to strongly fractionation of light rare earth elements (LREE), near-flat heavy rare earth elements (HREE) patterns and moderate to strong negative high field strength elements (HFSE) anomalies (especially Nb), a geochemical signature typical of subduction zones. The lowest units of mafic granulites and porphyroblastic amphibolites in the Trincheira ophiolite are similar to the modern mid-ocean ridge basalt (MORB), although they locally display small Ta, Ti and Nb negative anomalies, indicating a small subduction influence. This behavior changes to an island arc tholeiites (IAT) signature in the upper units of fine-grained amphibolites and amphibole rich-amphibolites, characterized by progressive depletion in the incompatible elements and more pronounced negative Ta and Nb anomalies, as well as common Ti and Zr negative anomalies. Tectono-magmatic variation diagrams and chondrite-normalized REE and primitive mantle normalized patterns suggest a back-arc to intra-oceanic island arc tectonic regime for the eruption of these rocks. Therefore, the Trincheira ophiolite appears to have originated in an intraoceanic supra-subduction setting composed of an arc-back-arc system. Accordingly, the Trincheira Complex is a record of oceanic crust relics obducted during the collision of the Amazon craton and the Paraguá block during the Middle Mesoproterozoic. Thus, the recognition of the Trincheira ophiolite and suture significantly changes views on the evolution of the southern margin of the Amazon craton, and how it can influence the global tectonics and the reconstruction of the continents.

XRF data, magnetic susceptibility, greyscale values and revised depths of sediments from ODP Leg 154 sites

The distinctly cyclic sediments recovered during ODP Leg 154 played an important role in constructing the astronomical time scale and associated astro(bio)chronology for the Miocene, and in deciphering ocean-climate history. The accuracy of the timescale critically depends on the reliability of the shipboard splice used for the tuning and on the tuning itself. New high-resolution colour- and magnetic susceptibility core scanning data supplemented with limited XRF-data allow improvement of the stratigraphy. The revised composite record results in an improved astronomical age model for ODP Site 926 between 5 and 14.4 Ma. The new age model is confirmed by results of complex amplitude demodulation of the precession and obliquity related cycle patterns. Different values for tidal dissipation are applied to improve the fit between the sedimentary cycle patterns and the astronomical solution. Due to the improved stratigraphy and tuning, supported by the results of amplitude demodulation, the revised time scale yields more reliable age estimates for planktic foraminiferal and calcareous nannofossil events. The results of this study highlight the importance of stratigraphy for timescale construction.

Ordinary and duplicate analysis from major and trace elements of ODP Holes 146-888B and 168-1027B

Oceanic sediments deposited at high rate close to continents are dominated by terrigenous material. Aside from dilution by biogenic components, their chemical compositions reflect those of nearby continental masses. This study focuses on oceanic sediments coming from the juvenile Canadian Cordillera and highlights systematic differences between detritus deriving from juvenile crust and detritus from old and mature crust. We report major and trace element concentrations for 68 sediments from the northernmost part of the Cascade forearc, drilled at ODP Sites 888 and 1027. The calculated weighted averages for each site can then be used in the future to quantify the contribution of subducted sediments to Cascades volcanism. The two sites have similar compositions but Site 888, located closer to the continent, has higher sandy turbidite contents and displays higher bulk SiO2/Al2O3 with lower bulk Nb/Zr, attributed to the presence of zircons in the coarse sands. Comparison with published data for other oceanic sedimentary piles demonstrates the existence of systematic differences between modern sediments deriving from juvenile terranes (juvenile sediments) and modern sediments derived from mature continental areas (cratonic sediments). The most striking systematic difference is for Th/Nb, Th/U, Nb/U and Th/Rb ratios: juvenile sediments have much lower ratios than cratonic sediments. The small enrichment of Th over Nb in cratonic sediments may be explained by intracrustal magmatic and metamorphic differentiation processes. In contrast, their elevated Th/U and Nb/U ratios (average values of 6.87 and 7.95, respectively) in comparison to juvenile sediments (Th/U ~ 3.09, Nb/U ~ 5.15) suggest extensive U and Rb losses on old cratons. Uranium and Rb losses are attributed to long-term leaching by rain and river water during exposure of the continental crust at the surface. Over geological times, the weathering effects create a slow but systematic increase of Th/U with exposure time.

Geochemical composition of gabbroic rocks from ODP Hole 179-1105A, southwest Indian Ridge

Oxide-free olivine gabbro and gabbro, and oxide olivine gabbro and gabbro make up the bulk of the gabbroic suite recovered from Ocean Drilling Program (ODP) Leg 179 Hole 1105A, which lies 1.2 km away from Hole 735B on the eastern transverse ridge of the Atlantis II Fracture Zone, Southwest Indian Ridge. The rocks recovered during Leg 179 show striking similarities to rocks recovered from the uppermost 500 m of Hole 735B during ODP Leg 118. The rocks of the Atlantis platform were likely unroofed as part of the footwall block of a large detachment fault on the inside corner of the intersection of the Southwest Indian Ridge and the Atlantis II Transform at ~11.5 Ma. We analyzed the lithologic, geochemical, and structural stratigraphy of the section. Downhole lithologic variation allowed division of the core into 141 lithologic intervals and 4 main units subdivided on the basis of predominance of oxide gabbroic vs. oxide-free gabbroic rocks. Detailed analyses of whole-rock chemistry, mineral chemistry, microstructure, and modes of 147 samples are presented and clearly show that the gabbroic rocks are of cumulate origin. These studies also indicate that geochemistry results correlate well with downhole magnetic susceptibility and Formation MicroScanner (FMS) resistivity measurements and images. FMS images show rocks with a well-layered structure and significant numbers of mappable layer contacts or compositional contrasts. Downhole cryptic mineral and whole-rock chemical variations depict both "normal" and inverse fine-scale variations on a scale of 10 m to <2 m with significant compositional variation over a short distance within the 143-m section sampled. A Mg# shift in whole-rock or Fo contents of olivine of as much as 20-30 units over a few meters of section is not atypical of the extreme variation in downhole plots. The products of the earliest stages of basaltic differentiation are not represented by any cumulates, as the maximum Fo content was Fo78. Similarly, the extent of fractionation represented by the gabbroic rocks and scarce granophyres in the section is much greater than that represented in the Atlantis II basalts. The abundance of oxide gabbros is similar to that in Hole 735B, Unit IV, which is tentatively correlated as a similar unit or facies with the oxide gabbroic units of Hole 1105A. Oxide phases are generally present in the most fractionated gabbroic rocks and lacking in more primitive gabbroic rocks, and there is a definite progression of oxide abundance as, for example, the Mg# of clinopyroxene falls below 73-75. Coprecipitation of oxide at such early Mg#s cannot be modeled by perfect fractional crystallization. In situ boundary layer fractionation may offer a more plausible explanation for the complex juxtaposition of oxide- and nonoxide-bearing more primitive gabbroic rocks. The geochemical signal may, in part, be disrupted by the presence of mylonitic shear zones, which strike east-west and dip both to the south and north, but predominantly to the south away from the northern rift valley where they formed. Downhole deformation textures indicate increasing average strain and crystal-plastic deformation in units that contain oxides. Oxide-rich zones may represent zones of rheologic weakness in the cumulate section along which mylonitic and foliated gabbroic shear zones nucleate in the solid state at high temperature, or the oxide may be a symptom of former melt-rich zones and hypersolidus flow, as predicted during study of Hole 735B.

Organic geochemistry of lower Cretaceous sediments at DSDP Hole 77-535

At Site 535, the four lithologic units of Cretaceous age are controlled by two types of sedimentologic facies: (1) the massive light-colored limestones or marly limestones in which the total organic carbon (TOC) content is low and the organic matter more or less oxidized and (2) laminated dark facies in which the TOC content is higher and associated with a well-preserved organic matter of Type II origin. Very little typical Type III organic matter occurs in the whole series from late Berriasian to Aptian and Cenomanian. Fluctuations from oxidizing to reducing environments of deposition are proposed to account for the variations in properties of the Type II organic matter between the different facies. Dark laminated layers are good but immature potential source rocks: petroleum potential is often higher than 2 kg HC/t of rock.

Organic carbon and total nitrogen stocks in soils of the Lena River Delta

The Lena River Delta, which is the largest delta in the Arctic, extends over an area of 32 000 km**2 and likely holds more than half of the entire soil organic carbon (SOC) mass stored in the seven major deltas in the northern permafrost regions. The geomorphic units of the Lena River Delta which were formed by true deltaic sedimentation processes are a Holocene river terrace and the active floodplains. Their mean SOC stocks for the upper 1 m of soils were estimated at 29 kg/m**2 ± 10 kg/m**2 and at 14 kg/m**2 ± 7 kg/m**2, respectively. For the depth of 1 m, the total SOC pool of the Holocene river terrace was estimated at 121 Tg ± 43 Tg, and the SOC pool of the active floodplains was estimated at 120 Tg ± 66 Tg. The mass of SOC stored within the observed seasonally thawed active layer was estimated at about 127 Tg assuming an average maximum active layer depth of 50 cm. The SOC mass which is stored in the perennially frozen ground at the increment 50-100 cm soil depth, which is currently excluded from intense biogeochemical exchange with the atmosphere, was estimated at 113 Tg. The mean nitrogen (N) stocks for the upper 1 m of soils were estimated at 1.2 kg/m**2 ± 0.4 kg/m**2 for the Holocene river terrace and at 0.9 kg/m**2 ± 0.4 kg/m**2 for the active floodplain levels, respectively. For the depth of 1 m, the total N pool of the river terrace was estimated at 4.8 Tg ± 1.5 Tg, and the total N pool of the floodplains was estimated at 7.7 Tg ± 3.6 Tg. Considering the projections for deepening of the seasonally thawed active layer up to 120 cm in the Lena River Delta region within the 21st century, these large carbon and nitrogen stocks could become increasingly available for decomposition and mineralization processes.

Geochemistry and mineralogy of ODP Leg 107 basalts

During ODP Leg 107, the basement of the Tyrrhenian Sea was drilled at Site 650, located in the Marsili basin, and at Sites 651 and 655, both located in the Vavilov basin. In addition, a lava flow was drilled at Site 654 on the Sardinia rifted margin. Mineral and whole rock major and trace element chemistry, including rare earth element (REE) and Sr and Nd isotopic ratios, were determined in samples of these rocks. Site 654 lava was sampled within uppermost Pliocene postrift sediments. This lava is a basaltic andesite of intraplate affinity, and is analogous to some Plio-Pleistocene tholeiitic lavas from Sardinia. Site 650 basalts, drilled beneath 1.7-1.9-Ma-old basal sediment, are strongly altered and vesicular suggesting a rapid subsidence of the Marsili basin. Based on incompatible trace elements, these basalts show calc-alkaline affinity like some products of the Marsili Seamount and the Eolian arc. The basement of the two sites drilled within Vavilov basin shows contrasting petrologies. Site 655, located along the Gortani ridge in the western part of the basin, drilled a 116-m-thick sequence of basalt flows beneath 3.4-3.6-Ma-old basal sediments. These basalts are chemically relatively homogeneous and show affinity to transitional MORB. Four units consisting of slightly differentiated basaltic lavas, have been identified. Site 655 basalts are geochemically similar to the high Ti lavas from DSDP Leg 42, Site 373 (Vavilov Basin). The basement at Site 651, overlain by 40 m of metalliferous dolostone covered by fossiliferous sediments with an age of 2 Ma, consists of two basalt units separated by a dolerite-albitite intrusive body; serpentinized harzburgites were drilled for 30 m at the base of the hole. The two basalt units of Site 651 are distinct petrochemically, though both show incompatible elements affinity with high-K calc-alkaline/calc-alkaline magmas from Eolian arc. The cpx chemistry and high K/Na ratio of the lower unit lavas suggest a weak alkaline tendency of potassic lineage. Leg 107 basement rock data, together with data from DSDP Site 373 and from dredged samples, indicate that the deepest basins of the central Tyrrhenian Sea are underlain by a complex back-arc basin crust produced by magmas with incompatible element affinities to transitional MORB (Site 655 and DSDP Site 373), and to calc-alkaline and high-K calc-alkaline converging plate margin basalts (Sites 650 and 651). This petrogenetic complexity is in accordance with the back-arc setting of the Vavilov and Marsili basins. Other back-arc basin basalts, particularly those from ensialic basins such as the Bransfield Strait (Antarctica), show a comparable petrogenetic complexity (cf., Sounders and Tarney, 1984).

The petrology of the lower series volcanics of ODP Hole 104-642E

Between 1086.6 and 1229.4 m below seafloor at Site 642 on the Outer Vøring Plateau, a series of intermediate volcanic extrusive flow units and volcaniclastic sediments was sampled. A mixed sequence of dacitic subaerial flows, andesitic basalts, intermediate volcaniclastics, subordinate mid-ocean ridge basalt, (MORB) lithologies, and intrusives was recovered, in sharp contrast to the more uniform tholeiitic T-type MORB units of the overlying upper series. This lower series of volcanics is composed of three chemically distinct groups, (B, A2, A1), rather than the two previously identified. Flows of the dacitic group (B) have trace-element and initial Sr isotope signatures which indicate that their source magma derived from the partial melting of a component of continental material in a magma chamber at a relatively high level in the crust. The relative proportions of crustal components in this complex melt are not known precisely. The most basic group (A2) probably represents a mixture of this material with MORB-type tholeiitic melt. A third group (A1), of which there was only one representative flow recovered, is chemically intermediate between the two groups above, and may suggest a repetition of, or a transition phase in, the mixing processes.

Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Abb. 3 Change of altitude in the EGIG-West-East-profile

The members of the two International Glaciological Greenland Expeditions (EGIG) in 1959 and 1968 determined the altitudes in the West-East-profile by levelling. The comparison of these two measurements gives evidence of the change of altitudes and of the surface-waves. The ways of measurement and of analysis are described here. The graphic representation of the result is discussed here in connexion with the possible causes.