Combining borehole temperature and meteorologic data to constrain past climate change

Geothermal observations from a suite of boreholes in western Utah, USA, combined with meteorologic data at nearby weather stations are used to test the hypothesis that temperatures in the earths subsurface contain an accurate record of recent climate change. The change in air temperature over the last hundred years successfully predicts detailed subsurface temperature profiles to better than ±0.05°C, indicating that ground temperatures tract air temperatures over long periods and that climate change signals are conducted into, and recorded in, the solid earth by the process of heat conduction. We combine borehole temperature data with meteorologic data from the nearest weather station to determine the time averaged difference between surface ground temperature and surface air temperature for borehole-weather station pairs and to infer the long term mean air temperature prior to the observational record. For our western Utah sites the preobservational mean temperature is close to the average surface air temperature for this century suggesting that up to 0.5°C of warming deduced from the last 100 years of weather station data may be attributed to recovery from a cool period at the turn of the century.     

Evidence of climatic warming in the southern Urals region derived from borehole temperatures and meteorological data

Thirty borehole temperature–depth profiles in the central and southern Urals, Russia were scrutinized for evidence of ground surface temperature histories. We explored two inversion schemes: a simple ramp inversion in which solutions are parameterized in terms of an onset time and magnitude of change and a more sophisticated functional space inverse algorithm in which the functional form of the solution is left unspecified. To enhance and potentially identify latitudinal differences in the ground surface temperature signal, we subdivided the data into three groups based on geographic proximity and simultaneously inverted the borehole temperature–depth logs. The simultaneous inversions highlighted 13 temperature–depth logs that could not both fit a common ground surface temperature history and a priori models within reasonable bounds. Our results confirm that this is an effective way to reduce site-specific noise from an ensemble of boreholes. Each inversion scheme gives comparable results indicating locally variable warming on the order of 1°C starting between 1800 and 1900 AD. Similarly surface air temperature records from 12 nearby meteorological stations exhibit locally variable warming also on the order of 1°C of warming during the 20th century. To explore the degree to which borehole temperatures and surface air temperature (SAT) time series are responding to the same signal, we average the SAT data into the same three groups and used these averages as a forcing function at the Earth's surface to generate synthetic transient temperature profiles. Root mean square (RMS) misfits between these synthetic temperature profiles and averaged temperature–depth profiles are low, suggesting that first-order curvature in borehole temperatures and variations in SAT records are correlated.     

Comparison of petrophysical properties of impactites for four meteoritic impact structures

We reanalyzed and compared unique data sets, which we obtained in the frame of combined petrophysical and geothermal investigations within scientific drilling projects on four impact structures: the Puchezh–Katunki impact structure (Vorotilovo borehole, Russia), the Ries impact structure (Noerdlingen‐73 borehole, Germany), the Chicxulub impact structure (ICDP Yaxcopoil‐1 borehole, Mexico), and the Chesapeake impact structure (ICDP‐USGS‐Eyreville borehole, USA). For a joined interpretation, we used the following previously published data: thermal properties, using the optical scanning technique, and porosities, both measured on densely sampled halfcores of the boreholes. For the two ICDP boreholes, we also used our previously published P‐wave velocities measured on a subset of cores. We show that thermal conductivity, thermal anisotropy, porosity, and velocity can be correlated with shock metamorphism (target rocks of the Puchezh–Katunki and Ries impact structures), and confirm the absence of shock metamorphism in the samples taken from megablocks (Chicxulub and Chesapeake impact structure). The physical properties of the lithic impact breccias and suevites are influenced mainly by their impact‐related porosity. Physical properties of lower porosity lithic impact breccias and suevites are also influenced by their chemical composition. These data allow for a distinction between different types of breccias due to differences concerning the texture and chemistry and the different amounts of melt and rock clasts.     

Long wavelength ground surface temperature history from continuous temperature logs in the Transylvanian Basin

Continuous temperature logs to depths between 750 and 1400 m in the Transylvanian Basin, Romania, in many cases show temperature gradient variations with depth which cannot be explained by depth variations in thermal conductivity, topography and ground water flow. The only possible responsible agent seems to be past surface temperature variations. The temperature logs from nine boreholes have been interpreted individually and jointly by least squares inverse modelling with the surface temperature history and background heat flux as unknown parameters. All the temperature profiles are consistent with a temperature rise at the end of the last glaciation. The effects of borehole depth, of a wrong choice of thermal conductivity, and the level of uncorrelated random noise were examined using synthetic examples.     

Warming permafrost in European mountains

Here we present the first systematic measurements of European mountain permafrost temperatures from a latitudinal transect of six boreholes extending from the Alps, through Scandinavia to Svalbard. Boreholes were drilled in bedrock to depths of at least 100 m between May 1998 and September 2000. Geothermal profiles provide evidence for regional-scale secular warming, since all are nonlinear, with near-surface warm-side temperature deviations from the deeper thermal gradient. Topographic effects lead to variability between Alpine sites. First approximation estimates, based on curvature within the borehole thermal profiles, indicate a maximum ground surface warming of +1 °C in Svalbard, considered to relate to thermal changes in the last 100 years. In addition, a 15-year time series of thermal data from the 58-m-deep Murtèl–Corvatsch permafrost borehole in Switzerland, drilled in creeping frozen ice-rich rock debris, shows an overall warming trend, but with high-amplitude interannual fluctuations that reflect early winter snow cover more strongly than air temperatures. Thus interpretation of the deeper borehole thermal histories must clearly take account of the potential effects of changing snow cover in addition to atmospheric temperatures.     

Subsurface temperature—depth profiles, anomalies due to climatic ground surface temperature changes or groundwater flow effects

Climatic temperature changes at the ground surface propagate downward to the subsurface creating transient disturbances to the temperature—depth (T(z)) profile. Due to the poor thermal diffusivity of rocks the disturbances are preserved long times in the bedrock, and in a conductive regime it is possible to reveal the ground surface temperature (GST) history from borehole temperature data with inversion techniques. Geothermal temperature measurements thus provide a source of palaeoclimatic information which so far has not been utilized extensively. Inversion of GST history is, however, not straightforward and any disturbing effects should be excluded before the data can be utilized in inversion. Groundwater flow is of special importance in this respect because it is a common phenomenon in bedrock and convection often produces temperature—depth profiles resembling those affected by palaeoclimatic GST changes. In interpreting temperature—depth (T(z)) logs it is therefore not always clear whether the recorded vertical gradient variations should be attributed to the effects of palaeoclimatic ground surface temperature (GST) changes or to groundwater circulation. Using several synthetic T(z) profiles and applying general least squares inversion techniques we simulate a situation of “misinterpreting” the curvature of the T(z) profile in terms of palaeoclimatic GST changes, although it is actually produced by convective heat transfer due to groundwater flow. For comparison the opposite case is also studied, namely, genuine palaeoclimatic effects are misinterpreted as being due to disturbances caused by groundwater flow. A homogeneous half-space model is used to model T(z) profiles disturbed conductively by GST changes during the time interval 10–10000 yr B.P. and a one-dimensional porous layer model is applied for convective heat transfer calculations. The results indicate that a given T(z) profile can be attributed to either of these effects with reasonable parameter values. In addition to the synthetic T(z) profiles, a case history from a 958 m deep drill hole at Lavia, southwestern Finland, is presented. Special care is needed in analyzing T(z) data. A knowledge of geothermal data, such as temperature, thermal conductivity and diffusivity is not necessarily adequate for determining which of the phenomena (or whether a combination of them) provides the most probable interpretation of a T(z) profile. Additional information on the hydrogeological properties of the drilled strata is essential.     

Potential causes of differences between ground and surface air temperature warming across different ecozones in Alberta, Canada

Analysis and modelling of temperature anomalies from 25 selected deep wells in Alberta show that the differences between GST (ground surface temperature) warming for the northern Boreal Forest ecozone and the combined Prairie Grassland ecozone and Aspen Parkland transition region to the south occur during the latter half of this century. This corresponds with recent changes in surface albedo resulting from permanent land development in the northern areas and also to increases in natural forest fires in the past 20 years. Differences between GST and SAT (surface air temperature) warming are much higher in the Boreal Forest ecozone than in the Prairie Grassland ecozone and Aspen Parkland transition region. Various hypotheses which could account for the existing differences between the GST and SAT warming in the different ecozones of Alberta, and western Canada in general, are tested. Analysis of existing data on soil temperature, hydrological piezometric surfaces, snowfall and moisture patterns, and land clearing and forest fires, indicate that large areas of Alberta, characterised by anomalous GST warming, have experienced widespread changes to the surface landscape in this century. It is postulated that this has resulted in a lower surface albedo with a subsequent increase in the absorption of solar energy. Heat flow modelling shows that, after climatic SAT warming, permanent clearing of the land is the most effective and likely cause of the observed changes in the GST warming. The greater GST warming in the Boreal Forest ecozone in the latter half of this century is related to landscape change due to land development and increasing forest fire activity. It appears to account for a portion of the observed SAT warming in this region through a positive feedback loop with the overlying air. The anthropogenic effect on regional climatic warming through 20th century land clearing and landscape alteration requires further study. In future, more accurate quantification of these various forcings will be necessary in order to distinguish between, and to detect, the variety of natural and anthropogenic influences and on climate.     

A comparison of five different analyses in the interpretation of five borehole temperature data sets

To assess the effectiveness of various methods in inferring ground surface temperature history (GSTH) from perturbed borehole temperature data, five synthetic and real data sets were prepared and disseminated to interested groups for analysis. This report summarizes the results obtained by five different groups using different inverse methods. The five inverse methods yield similar results inasmuch as the salient features of the GSTH are concerned. However, the estimated GSTHs differ in details of the timing and the magnitude of changes. This is attributed to the differences in parameterization schemes and to the different a priori constraints imposed on the GSTH. Better agreement among the inverse methods may be obtainable when information other than that contained in the borehole data is incorporated.     

Climate change of the last millennium inferred from borehole temperatures: regional patterns of climatic changes in the Czech Republic — Part III

Accurate temperature–depth profiles may help to assess the temperature variations associated with the climate changes in the past. Ninety-eight ground surface temperature histories inverted from the temperature–depth borehole logs drilled on the territory of the Czech Republic [Bodri, L., ermák, V., 1995. Climate changes of the last millennium inferred from borehole temperatures: results from the Czech Republic — Part I. Global Planet. Change 11, pp. 111–125; Bodri, L., ermák, V., 1997. Climate changes of the last two millennia inferred from borehole temperatures: results from the Czech Republic — Part II. Global Planet. Change 14, pp. 163–173.] are used to reconstruct the regional patterns of the respective climate change. The climate was mapped for the following periods: 1100–1300 A.D. (Little Climatic Optimum), 1400–1500 A.D., 1600–1700 A.D. (main phase of the Little Ice Age), and for the most recent climate trend after year 1960. Comparison of the obtained maps with the meteorological observations and proxy climatic reconstructions confirmed good applicability of the “geothermal” paleoclimatic reconstructions for the regional studies.     

Geochemical studies of the SUBO 18 (Enkingen) drill core and other impact breccias from the Ries crater,Germany

Suevite and melt breccia compositions in the boreholes Enkingen and Polsingen are compared with compositions of suevites from other Ries boreholes and surface locations and discussed in terms of implications for impact breccia genesis. No significant differences in average chemical compositions for the various drill cores or surface samples are noted. Compositions of suevite and melt breccia from southern and northeastern sectors of the Ries crater do not significantly differ. This is in stark contrast to the published variations between within‐crater and out‐of‐crater suevites from northern and southern sectors of the Bosumtwi impact structure, Ghana. Locally occurring alteration overprint on drill cores—especially strong on the carbonate‐impregnated suevite specimens of the Enkingen borehole—does affect the average compositions. Overall, the composition of the analyzed impact breccias from Ries are characterized by very little macroscopically or microscopically recognized sediment‐clast component; the clast populations of suevite and impact melt breccia are dominated consistently by granitic and intermediate granitoid components. The Polsingen breccia is significantly enriched in a dioritic clast component. Overall, chemical compositions are of intermediate composition as well, with dioritic‐granodioritic silica contents, and relatively small contributions from mafic target components. Selected suevite samples from the Enkingen core have elevated Ni, Co, Cr, and Ir contents compared with previously analyzed suevites from the Ries crater, which suggest a small meteoritic component. Platinum‐group element (PGE) concentrations for some of the enriched samples indicate somewhat elevated concentrations and near‐chondritic ratios of the most immobile PGE, consistent with an extraterrestrial contribution of 0.1–0.2% chondrite‐equivalent.     

Lunar heat flow and regolith structure inferred from interferometric observations at a wavelength of 49.3 cm

We discuss observations of the Moon at a wavelength of 49.3 cm made with the Owens Valley Radio Observatory Interferometer. These observations have been fit to models in order to estimate the lunar dielectric constant, the equatorial subsurface temperature, the latitude dependence of the subsurface temperature, and the subsurface temperature gradient. The models are most consistent with a dielectric constant of 2.52 ± 0.01 (formal errors), an equatorial subsurface temperature of 249_?5⁺⁸K, and a change in the subsurface temperature with latitude (ψ), which is proportional to cos^0.38ψ. Since the temperature of the Moon has been measured by the Apollo Lunar Heat Flow Experiment, we have been able to use our determination of the equatorial temperature to estimate the error in the flux density calibration scale at 49.3cm (608 MHz). This results in a correction factor of 1.03 ± 0.04, which must be applied to the flux density scale. This factor is much different from 1.21 ± 0.09 estimated by Muhleman et al. (1973) from the brightness temperature of Venus and apparently indicates that the observed decrease in the brightness temperature of Venus at long wavelengths is a real effect.The estimates of the temperature gradient, which are based on the measurement of limb darkening, are small and negative (temperature decreases with depth) and may be insignificantly different from zero since they are only as large as their formal errors. We estimate that a temperature gradient in excess of 0.6K/m at 10m depth would have been observed. Thus, a temperature gradient like that measured in situ at the Apollo 15 and 17 landing sites in the upper 2m of the regolith is not typical of the entire lunar frontside at the 10m depths where the 49.3 cm wavelength emission originates. This result may indicate that the mean lunar heat flow is lower than that measured at the Apollo landing sites, that the thermal conductivity is greater at 10m depth than it is at 2m depth, or that the radio opacity is greater at 10m depth than at 2m depth. The negative estimates of the temperature gradient indicate that the Moon appeared limb bright and might be explained by scattering of the emission from boulders or an interface with solid rock. The presence of solid rock at 10m depths will probably cause heat flows like those measured by Apollo to be unobservable by our interferometric method at long wavelengths, since it will cause both the thermal conductivity and radio opacity of the regolith to increase. Thus, our data may be most consistent with a change in the physical properties of the regolith to those of solid rock or a mixture of rock and soil at depths of 7 to 16m. Our results show that future radio measurements for heat flow determinations must utilize wavelengths considerably shorter than 50 cm (25 cm or less) to avoid the rock regions below the regolith.     

Ground temperature histories in eastern and central Canada from geothermal measurements: evidence of climatic change

Inverse and direct methods have been used to analyze a large number of borehole temperature logs in order to infer past climatic changes. Results indicate a warming of 1–2°C in eastern and central Canada during the past 150 years. A period of cooling between 500 and 200 years before present, corresponding to the time of the “Little Ice Age”, has also been identified in the same areas. A regional ground temperature history is estimated for eastern and central Canada from the simultaneous inversion of several temperature logs. The inferred temperature changes appear correlated with the concentration of atmospheric carbon dioxide as reported from a Greenland ice core, and agree with existing meteorological and dendrochronological records for the area.     

Borehole temperatures, climate change and the pre-observational surface air temperature mean: allowance for hydraulic conditions

Joint analysis of surface air temperature series recorded at weather stations together with the inversion of the temperature-depth profiles logged in the near-by boreholes enables an estimate of the conditions existing prior to the beginning of the meteorological observation, the so-called pre-observational mean (POM) temperature.Such analysis is based on the presumption of pure diffusive conditions in the underground. However, in real cases a certain subsurface fluid movement cannot be excluded and the measured temperature logs may contain an advective component. The paper addresses the correction for the hydraulic conditions, which may have perturbed the climate signal penetrating from the surface into the underground. The method accounts for vertical conductive and vertical advective heat transport in a 1-D horizontally layered stratum and provides a simultaneous evaluation of the POM-temperature together with the estimate of the Darcy fluid velocity. The correction strategy is illustrated on a synthetic example and its use is demonstrated on temperature logs measured in four closely spaced boreholes drilled near Tachlovice (located about 15 km SW of Prague, Czech Republic). The results revealed that in a case of moderately advectively affected subsurface conditions (fluid velocities about 10⁻⁹ m/s), the difference between POM-values assessed for a pure conductive approach and for combined vertical conductive/advective approach may amount up to 0.3–0.5 K, the value comparable with the amount usually ascribed to the 20th century climate warming.     

Composite surface temperature history from simultaneous inversion of borehole temperatures in western Canadian plains

Functional space inversions (FSI) of precise temperature logs from wells located in low conductivity clastic sediments of the western Canadian Sedimentary Basin show evidence of extensive, recent ground surface temperature (GST) warming. Simultaneous inversion of the data, as well as averaging of the individual site reconstructions, indicate that this high magnitude of GST warming exceeds over two times that of globally averaged GST's [Science 282 (1998) 279] and is significantly higher than that of surface temperature histories based on instrumental records and tree ring reconstruction in northern and western Canada [Holocene 7 (1997) 375; Science 278 (1997) 1251; Clim. Res. 12 (1999) 39].     

The effect of target properties on crater morphology: Comparison of central peak craters on the Moon and Ganymede

Abstract— We examine the morphology of central peak craters on the Moon and Ganymede in order to investigate differences in the near‐surface properties of these bodies. We have extracted topographic profiles across craters on Ganymede using Galileo images, and use these data to compile scaling trends. Comparisons between lunar and Ganymede craters show that crater depth, wall slope and amount of central uplift are all affected by material properties. We observe no major differences between similar‐sized craters in the dark and bright terrain of Ganymede, suggesting that dark terrain does not contain enough silicate material to significantly increase the strength of the surface ice. Below crater diameters of ?12 km, central peak craters on Ganymede and simple craters on the Moon have similar rim heights, indicating comparable amounts of rim collapse. This suggests that the formation of central peaks at smaller crater diameters on Ganymede than the Moon is dominated by enhanced central floor uplift rather than rim collapse. Crater wall slope trends are similar on the Moon and Ganymede, indicating that there is a similar trend in material weakening with increasing crater size, and possibly that the mechanism of weakening during impact is analogous in icy and rocky targets. We have run a suite of numerical models to simulate the formation of central peak craters on Ganymede and the Moon. Our modeling shows that the same styles of strength model can be applied to ice and rock, and that the strength model parameters do not differ significantly between materials.     

Variations in ground surface temperature histories in the Thompson Belt, Manitoba, Canada: environment and climate changes

Several temperature–depth profiles recorded at Pipe Mine, 32 km southwest of Thompson, Manitoba, in central Canada, exhibit a marked departure from the equilibrium gradient. These profiles could be interpreted as indicating strong warming (up to 4.5 K) of the ground surface during the last 200 years. All the temperature profiles at Pipe Mine show perturbations stronger than at the others sites in the Thompson Nickel Belt. Temperature profiles recorded near the town of Thompson show a moderate warming (≈1–2 K) trend, while temperature profiles at Soab, 45 km southwest of Pipe Mine, indicate very moderate cooling (≈0.5 K). There was little human activity in this part of Manitoba before the development of the mining camp of Thompson in the late 1950s. Our study shows the variability of ground surface temperature histories at a very local scale (i.e. <1 km) with much stronger signals at some of the Pipe Mine drill holes than at others. These holes are located within 500 m of the highway and a power line built after 1955, at ≈3 km from the now abandoned open pit mine. The ground surface temperature history (GSTH) obtained by the inversion of Pipe Mine temperature profiles suggests that a recent (50 years) and strong (≈1–2 K) ground surface warming is superimposed on a 1–2 K warming trend that started 200 years ago, without any indication of a cold (little ice ages) episode before. The recent warming (40 years) at Pipe Mine is only a local effect and is likely to be related to the presence of the highway. Before 1960, the ground surface temperature history for Pipe is similar to other sites in the Thompson region. Ground surface temperature histories from other profiles within and near the city of Thompson seem less affected by environmental perturbations and their trends are parallel to that of the meteorological records in the Canadian Prairies.     

Estimation of plasma parameters in an accretion column near the surface of accreting white dwarfs from their flux variability

We consider the behavior of matter in the accretion column that emerges under accretion in binary systems near the surface of a white dwarf. The plasma heated in a standing shock wave near the white dwarf surface efficiently radiates in the X-ray energy band. We suggest a method for estimating post-shock plasma parameters, such as the density, temperature, and height of the hot zone, from the power spectrum of its X-ray luminosity variability. The method is based on the fact that the flux variability amplitude for the hot region at various Fourier frequencies depends significantly on its cooling time, which is determined by the parameters of the hot zone in the accretion column. This allows the density and temperature of the hot matter to be estimated. We show that the characteristic cooling time can be efficiently determined from the break frequency in the power spectrum of the X-ray flux variability for accreting white dwarfs. The currently available X-ray instruments do not allow such measurements to be made because of an insufficient collecting area, but this will most likely become possible with new-generation large-area X-ray spectrometers.     

Leaf margins and temperature in the North American flora: Recalibrating the paleoclimatic thermometer

For almost a century, leaf margin analysis has been an important method of estimating Cenozoic paleotemperatures. However, the relation between the proportion of toothed versus entire leaves and temperature has so far been calibrated and characterized using relatively small numbers of opportunistically and irregularly collected floras.Here we present the results of the first systematic spatially-distributed analysis of the relation between leaf margins and temperature for North America.We used species range maps of native North American dicot trees to derive synthetic local floras for each 50 km grid cell in the USA and Canada, and compared the percentage of entire-margined species with the mean annual temperature in each cell. Leaf margin type for each species was assigned using illustrations from appropriate regional North America floras.This analysis confirms that high-rainfall areas of the eastern USA and Canada show a strong inverse relationship between toothed margin percentage and mean annual temperature in the range 0–25 °C. The values we estimate are roughly similar to those obtained from floras in east Asia and previous analyses in North America, but the detailed curve is more complex, showing a distinct change in slope above a mean annual temperature of 20 °C. Wetter climate areas of the coastal western USA do not show any significant correlation between leaf margin percentage and temperature. This may relate to the confounding effect of the diverse topography of the West Coast, or it may be caused by the depleted dicotyledonous tree flora of this area.Like previous studies, this paper demonstrates that there indeed is a strong relationship between leaf margin percentage and temperature on a regional scale in eastern North America. However the relationship does not have the same linear behavior as results obtained from previous analyses. We discuss the implications of the results obtained here for paleoclimatic reconstruction in the content of other leaf margin analyses based on North American sites.     

The extent of temporal smearing in surface-temperature histories derived from borehole temperature measurements

The ability of borehole temperature data to resolve past climatic events is investigated using Backus-Gilbert inversion methods. Two experimental approaches are considered: (1) the data consist of a single borehole temperature profile, and (2) the data consist of climatically-induced temperature transients measured within a borehole during a monitoring experiment. The sensitivity of the data's resolving power to the vertical distribution of the measurements, temperature measurement errors, the inclusion of a local meteorological record, and the duration of a monitoring experiment, are investigated. The results can be used to help interpret existing surface temperature histories derived from borehole temperature data and to optimize future experiments for the detection of climatic signals.     

Overview of analogue science activities at the McGill Arctic Research Station,Axel Heiberg Island,Canadian High Arctic

The Canadian High Arctic contains several of the highest fidelity Mars analogue sites in the world. Situated at nearly 80° north, Expedition Fjord on Axel Heiberg Island is located within a polar desert climate, with the surrounding landscape and conditions providing an invaluable opportunity to examine terrestrial processes in a cold, dry environment. Through the Canadian Space Agency's Analogue Research Network program, scientific activities based out of the McGill Arctic Research Station (M.A.R.S.) are extremely broad in scope, representing physical, biological, and technological investigations. Some of the most unique hydrogeologic features under investigation near M.A.R.S. are a series of cold saline springs that maintain liquid-state flow year round regardless of air temperature. Previous studies have examined their geomorphic relation to discharge-related formations, water chemistry, temperature monitoring, discharge rates, and combined flow/thermal modeling. Recent investigations have identified microbial communities and characterized biological activity within the springs and within permafrost sections, having direct relevance to astrobiological analogue research goals. Another main thrust of research activities based at M.A.R.S. pertains to the detection, mapping, and quantification of subsurface ice deposits. A long-term study is presently underway examining polygonal terrain, comparing surficial patterns found in the region with those identified on Mars, and using surface morphology to estimate ice wedge volumes through a combination of aerial photography interpretation and ground-based geophysical techniques. Other technological developments include the use of in situ microscopy for the detection of biomarkers and improved permafrost drilling techniques. This paper presents an overview of previous studies undertaken at M.A.R.S. over the past decades and will describe in detail both present and upcoming work.