Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Exposure dating of Late Glacial and pre-LGM moraines in the Cordon de Doña Rosa, Northern/Central Chile (~ 31°S)

Despite the important role of the Central Andes (15–30° S) for climate reconstruction, knowledge about the Quaternary glaciation is very limited due to the scarcity of organic material for radiocarbon dating. We applied 10Be surface exposure dating (SED) on 22 boulders from moraines in the Cordon de Doña Rosa, Northern/Central Chile (~31° S). The results show that several glacial advances in the southern Central Andes occurred during the Late Glacial between ~14.7±1.5 and 11.6±1.2 ka. A much more extensive glaciation is dated to ~32±3 ka, predating the temperature minimum of the global LGM (Last Glacial Maximum: ~20 ka). Reviewing these results in the paleoclimatic context, we conclude that the Late Glacial advances were most likely caused by an intensification of the tropical circulation and a corresponding increase in summer precipitation. High-latitude temperatures minima, e.g. the Younger Dryas (YD) and the Antarctic Cold Reversal (ACR) may have triggered individual advances, but current systematic exposure age uncertainties limit precise correlations. The absence of LGM moraines indicates that moisture advection was too limited to allow significant glacial advances at ~20 ka. The tropical circulation was less intensive despite the maximum in austral summer insolation. Winter precipitation was apparently also insufficient, although pollen and marine studies indicate a northward shift of the westerlies at that time. The dominant pre-LGM glacial advances in Northern/Central Chile at ~32 ka required lower temperatures and increased precipitation than today. We conclude that the westerlies were more intense and/or shifted equatorward, possibly due to increased snow and ice cover at higher southern latitudes coinciding with a minimum of insolation.

LGM and Late Glacial glacier advances in the Cordillera Real and Cochabamba (Bolivia) deduced from 10Be surface exposure dating

Surface exposure dating (SED) is an innovative tool already being widely applied for moraine dating and for Late Quaternary glacier and climate reconstruction. Here we present exposure ages of 28 boulders from the Cordillera Real and the Cordillera Cochabamba, Bolivia. Our results indicate that the local Last Glacial Maximum (LGM) in the Eastern Cordilleras occurred at ~22–25 ka and was thus synchronous to the global temperature minimum. We were also able to date several Late Glacial moraines to ~11–13 ka, which likely document lower temperatures and increased precipitation ("Coipasa" humid phase). Additionally, we recognize the existence of older Late Glacial moraines re-calculated to ~15 ka from published cosmogenic nuclide data. Those may coincide with the cold Heinrich 1 event in the North Atlantic region and the pronounced "Tauca" humid phase. We conclude that (i) exposure ages in the tropical Andes may have been overestimated so far due to methodological uncertainties, and (ii) although precipitation plays an important role for glacier mass balances in the tropical Andes, it becomes the dominant forcing for glaciation only in the drier and thus more precipitation-sensitive regions farther west and south.

Surface uplift and climate change - the geomorphic evolution of the western escarpment of the Andes of northern Chile between the Miocene and present

The Western Escarpment of the Andes at 18.30°S (Arica area, northern Chile) is a classical example for a transient state in landscape evolution. This part of the Andes is characterized by the presence of >10,000 km2 plains that formed between the Miocene and the present, and >1500 m deeply incised valleys. Although processes in these valleys scale the rates of landscape evolution, determinations of ages of incision, and more importantly, interpretations of possible controls on valley formation have been controversial. This paper uses morphometric data and observations, stratigraphic information, and estimates of sediment yields for the time interval between ca. 7.5 Ma and present to illustrate that the formation of these valleys was driven by two probably unrelated components. The first component is a phase of base-level lowering with magnitudes of∼300–500 m in the Coastal Cordillera. This period of base-level change in the Arica area, that started at ca. 7.5 Ma according to stratigraphic data, caused the trunk streams to dissect headward into the plains. The headward erosion interpretation is based on the presence of well-defined knickzones in stream profiles and the decrease in valley widths from the coast toward these knickzones. The second component is a change in paleoclimate. This interpretation is based on (1) the increase in the size of the largest alluvial boulders (from dm to m scale) with distal sources during the last 7.5 m.y., and (2) the calculated increase in minimum fluvial incision rates of ∼0.2 mm/yr between ca. 7.5 Ma and 3 Ma to ∼0.3 mm/yr subsequently. These trends suggest an increase in effective water discharge for systems sourced in the Western Cordillera (distal source). During the same time, however, valleys with headwaters in the coastal region (local source) lack any evidence of fluvial incision. This implies that the Coastal Cordillera became hyperarid sometime after 7.5 Ma. Furthermore, between 7.5 Ma and present, the sediment yields have been consistently higher in the catchments with distal sources (∼15 m/m.y.) than in the headwaters of rivers with local sources (<7 m/m.y.). The positive correlation between sediment yields and the altitude of the headwaters (distal versus local sources) seems to reflect the effect of orographic precipitation on surface erosion. It appears that base-level change in the coastal region, in combination with an increase in the orographic effect of precipitation, has controlled the topographic evolution of the northern Chilean Andes.

Glacier response to the change in atmospheric circulation in the eastern Mediterranean during the Last Glacial Maximum

In this study, we document glacial deposits and reconstruct the glacial history in the Karagöl valley system in the eastern Uludağ in northwestern Turkey based on 42 cosmogenic 10Be exposure ages from boulders and bedrock. Our results suggest the Last Glacial Maximum (LGM) advance prior to 20.4 ± 1.2 ka and at least three re-advances until 18.6 ± 1.2 ka during the global LGM within Marine Isotope Stage-2. In addition, two older advances of unknown age are geomorphologically well constrained, but not dated due to the absence of suitable boulders. Glaciers advanced again two times during the Lateglacial. The older is exposure dated to not later than 15.9 ± 1.1 ka and the younger is attributed to the Younger Dryas (YD) based on field evidence. The timing of the glaciations in the Karagöl valley correlates well with documented archives in the Anatolian and Mediterranean mountains and the Alps. These glacier fluctuations may be explained by the change in the atmospheric circulation pattern during the different phases of North Atlantic Oscillation (NAO) winter indices.

Late Weichselian local ice dome configuration and chronology in Northwestern Svalbard: early thinning, late retreat

The chronology and configuration of the Svalbard Barents Sea Ice Sheet (SBSIS) during the Late Weichselian (LW) are based on few and geographically scattered data. Thus, the timing and configuration of the SBSIS has been a subject of extensive debate. We present provenance data of erratic boulders and cosmogenic 10Be ages of bedrock and boulders from Northwest Spitsbergen (NWS), Svalbard to determine the thickness, configuration and chronology during the LW. We sampled bedrock and boulders of mountain summits and summit slopes, along with erratic boulders from coastal locations around NWS. We suggest that a local ice dome over central NWS during LW drained radially in all directions. Provenance data from erratic boulders from northern coastal lowland Reinsdyrflya suggest northeastward ice flow through Liefdefjorden. 10Be ages of high-elevation erratic boulders in central NWS (687–836 m above sea level) ranging from 18.3 ± 1.3 ka to 21.7 ± 1.4 ka, indicate that the centre of a local ice dome was at least 300 m thicker than at present. 10Be ages of all high-elevation erratics (>400 m above sea level, central and coastal locations) indicate the onset of ice dome thinning at 25–20 ka. 10Be ages from erratic boulders on Reinsdyrflya ranging from 11.1 ± 0.8 ka to 21.4 ± 1.7 ka, indicate an ice cover over the entire Reinsdyrflya during LW and a complete deglaciation prior to the Holocene, but apparently later than the thinning in the mountains. Lack of moraine deposits, but the preservation of beach terraces, suggest that the ice covering this peninsula possibly was cold-based and that Reinsdyrflya was part of an inter ice-stream area covered by slow-flowing ice, as opposed to the adjacent fjord, which possibly was filled by a fast-flowing ice stream. Despite the early thinning of the ice sheet (25–20 ka) we find a later timing of deglaciation of the fjords and the distal lowlands. Several bedrock samples (10Be) from vertical transects in the central mountains of NWS pre-date the LW, and suggest either ice free or pervasive cold-based ice conditions. Our reconstruction is aligned with the previously suggested hypothesis that a complex multi-dome ice-sheet-configuration occupied Svalbard and the Barents Sea during LW, with numerous drainage basins feeding fast ice streams, separated by slow flowing, possibly cold-based, inter ice-stream areas.

Glacier advances in northeastern Turkey before and during the global Last Glacial Maximum

Abstract Our study in the Başyayla Valley in northeastern Anatolia showed evidence of four glacier advances that built terminal and lateral moraines. Surface exposure dating of boulders on these moraines showed that the Maximum Ice Extent (MIE) was asynchronous with the global Last Glacial Maximum (LGM; 22.1 ± 4.3 thousand years; ka). The local {MIE} took place at least 57.0 ± 3.5 ka ago. The extent of the Başyayla Glacier during this advance is not known exactly because the boulders are only preserved on a lateral moraine. The next advance was prior to 41.5 ± 2.5 ka, and it descended down the valley to approximately 2320 m above sea level (m a.s.l.), with a glacier length of 5.3 km. During the early global LGM, the Başyayla Glacier extended for a distance of 4.9 km down to approx. 2430 m a.s.l. The last recorded advance occurred during the global LGM. This extension was 0.7 km smaller than the local {MIE} and its terminus reached 2490 m a.s.l. only. The exposure ages of boulders in a retreat position at an altitude of approx. 3045 m a.s.l. indicate that the valley has remained ice-free since the Lateglacial period. Therefore, the Lateglacial extent was limited to the cirque system in the uppermost part of the catchment. Furthermore, Holocene glacier oscillations seem to be either absent or restricted to solifluction in the whole catchment and to rock glacier movements in the southern tributary of the Başyayla Valley system.

Timing of retreat of the Reuss Glacier (Switzerland) at the end of the Last Glacial Maximum

We used cosmogenic 10Be and 36Cl to establish the timing of the onset of deglaciation after the Last Glacial Maximum of the Reuss Glacier, one of the piedmont lobes of the Alpine ice cap that reached the northern Alpine foreland in Switzerland. In this study, we sampled erratic boulders both at the frontal position in the foreland (Lenzburg and Wohlen, canton Aargau) and at the lateral Alpine border position (Seeboden moraine, Rigi, canton Schwyz). The minimum age for the beginning of retreat is 22.2 ± 1.0 ka at the frontal (terminal) position and 20.4 ± 1.0 ka at the lateral position. These ages are directly comparable with exposure ages from the other piedmont lobes in the northern Alpine foreland. Our data from the mountain called Rigi, do not support the hypothesis that boulders located external to the Seeboden moraine were deposited prior to the last glacial cycle. We present a first exposure age from an erratic boulder in a retreat position in the Alpine foreland. The Reuss Glacier was approximately 12 km behind the maximal extent no later than at 18.6 ± 0.9 ka.

Minor inheritance inhibits the calibration of the 10Be production rate from the AD 1717 Val Ferret rock avalanche, European Alps

To calibrate the in situ 10Be production rate, we collected surface samples from nine large granitic boulders within the deposits of a rock avalanche that occurred in AD 1717 in the upper Ferret Valley, Mont Blanc Massif, Italy. The 10Be concentrations were extremely low and successfully measured within 10% analytical uncertainty or less. The concentrations vary from 4829 ± 448 to 5917 ± 476 at g−1. Using the historical age exposure time, we calculated the local and sea level-high latitude (i.e. ≥60°) cosmogenic 10Be spallogenic production rates. Depending on the scaling schemes, these vary between 4.60 ± 0.38 and 5.26 ± 0.43 at g−1 a−1. Although they correlate well with global values, our production rates are clearly higher than those from more recent calibration sites. We conclude that our 10Be production rate is a mean and an upper bound for production rates in the Massif region over the past 300 years. This rate is probably influenced by inheritance and will yield inaccurate (e.g. too young) exposure ages when applied to surface-exposure studies in the area. Other independently dated rock-avalanche deposits in the region that are approximately 103 years old could be considered as possible calibration sites.

36Cl production rate from K-spallation in the European Alps (Chironico landslide, Switzerland)

The abundant production of in situ cosmogenic 36Cl from potassium renders 36Cl measurements in K-rich rocks or minerals, such as K-feldspars, potentially useful for precisely dating rock surfaces, either in single-nuclide or in multi-nuclide studies, for example combined with 10Be measurements in quartz. However, significant discrepancies in experimentally calibrated 36Cl production rates from spallation of potassium (36PK-sp), referenced to sea-level/high-latitude (SLHL), limit the accuracy of 36Cl dating from K-rich lithologies. We present a new 36Cl calibration using K-feldspars, in which K-spallation is the most dominant 36Cl production pathway (>92% of total 36Cl), thus minimizing uncertainties from the complex multi-pathway 36Cl production systematics. The samples are derived from boulders of an ∼13.4 ka-old landslide in the Swiss Alps (∼820 m, 46.43°N, 8.85°E). We obtain a local 36PK-sp of 306 ± 16 atoms 36Cl (g K)−1 a−1 and an SLHL 36PK-sp of 145.5 ± 7.7 atoms 36Cl (g K)−1 a−1, when scaled with a standard scaling protocol (‘Lm’). Applying this SLHL 36PK-sp to determine 36Cl exposure ages of K-feldspars from 10Be-dated moraine boulders yields excellent agreement, confirming the validity of the new SLHL 36PK-sp for surface exposure studies, involving 36Cl in K-feldspars, in the Alps.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.

Search for satellites near comet 67P/Churyumov-Gerasimenko using Rosetta/OSIRIS images

Context. The European Space Agency Rosetta mission reached and started escorting its main target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, at the beginning of August 2014. Within the context of solar system small bodies, satellite searches from approaching spacecraft were extensively used in the past to study the nature of the visited bodies and their collisional environment. Aims. During the approaching phase to the comet in July 2014, the OSIRIS instrument onboard Rosetta performed a campaign aimed at detecting objects in the vicinity of the comet nucleus and at measuring these objects' possible bound orbits. In addition to the scientific purpose, the search also focused on spacecraft security to avoid hazardous material in the comet's environment. Methods. Images in the red spectral domain were acquired with the OSIRIS Narrow Angle Camera, when the spacecraft was at a distance between 5785 km and 5463 km to the comet, following an observational strategy tailored to maximize the scientific outcome. From the acquired images, sources were extracted and displayed to search for plausible displacements of all sources from image to image. After stars were identified, the remaining sources were thoroughly analyzed. To place constraints on the expected displacements of a potential satellite, we performed Monte Carlo simulations on the apparent motion of potential satellites within the Hill sphere. Results. We found no unambiguous detections of objects larger than similar to 6 m within similar to 20 km and larger than similar to 1 m between similar to 20 km and similar to 110 km from the nucleus, using images with an exposure time of 0.14 s and 1.36 s, respectively. Our conclusions are consistent with independent works on dust grains in the comet coma and on boulders counting on the nucleus surface. Moreover, our analysis shows that the comet outburst detected at the end of April 2014 was not strong enough to eject large objects and to place them into a stable orbit around the nucleus. Our findings underline that it is highly unlikely that large objects survive for a long time around cometary nuclei.

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Aims. We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material (tensile, shear, and compressive strengths). Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies (overhangs, collapsed structures, boulders, cliffs, and Philae's footprint) and mechanical considerations. Results. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 degrees) are covered by a fine material and contain a few large (>10 m) and isolated boulders; ii) intermediate-slope terrains (20-45 degrees) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from <1m to 10m for the majority of them; and iii) high-slope terrains (45-90 degrees) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.