Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.

Late Neogene carbonate sedimentation and stable isotope record of ODP Site 121-758

The evolution of oceanic and climatic conditions the northeast Indian Ocean during the last 7 m.y. is revealed in the sediments from Site 758. We present detailed and continuous records of d18O and d13C from planktonic foraminifers, weight percent calcium carbonate, weight percent coarse fraction, magnetic susceptibility, and geomagnetic reversals. Sample spacing of the records ranges from 3 to 10 cm and is equivalent to an average time interval of 2000 to 6000 yr. Despite the fact that core recovery ranged between 100% and 105%, recovery gaps as large as 2.7 m occurred at nearly every break between advanced hydraulic piston cores. Approximately 12% of the late Neogene sequence was not recovered in each of the two holes drilled at Site 758. To circumvent the discontinuity introduced by the gaps, a composite depth section was constructed from multiple cores taken from offset holes at Site 758. The resulting composite depth section extends continuously from 0 to 116 mbsf, from the Holocene to the upper Miocene. A detailed chronostratigraphy is based on geomagnetic reversals which extend from the Brunhes Chron to Chron 6, and on d18O stages 1 through 105, which span from 0 to 2.5 Ma. The d18O record is dominated by a ~40-k.y. cycle in the late Pliocene and early Pleistocene, and is followed by a change to a ~100-k.y. cycle in the late Pleistocene. The mid-Pleistocene transition between these two modes of variability occurs between d18O stages 25 and 22 (between 860 and 800 Ka). Thirteen major volcanic ash horizons from the Indonesian arc are observed throughout the sedimentary section and are dated by their relative position within the geomagnetic reversals and the d18O chronostratigraphy. Since 5 Ma, there has been a long-term decline in weight percent CaCO3 and CaCO3 mass accumulation rates, and an associated rise in non-CaCO3 mass accumulation rates. We attribute these changes to a decrease in CaCO3 productivity and an increase in terrigenous sedimentation through enhanced riverine input. Such input may be linked to rapid tectonic uplift of the Himalayas and the Tibetan Plateau via mechanisms such as the intensification of the monsoonal rains, increased fluvial erosion, and regional glaciation. The long-term increase in percent coarse fraction since 5 Ma suggests a gradual increase in CaCO3 preservation. Higher frequency fluctuations in CaCO3 preservation are superimposed on the long-term trend and are related to climate fluctuations. The abrupt drop (-50%) in CaCO3 accumulation at 3.4 Ma signals a dramatic decrease in CaCO3 production that occurred over much of the Indian Ocean.

Reconstruction of the provenance and distribution of terrigenous sediments off NE South America over the last 30ka

We investigate the redistribution of terrigenous materials in the northeastern (NE) South American continental margin during slowdown events of the Atlantic Meridional Overturning Circulation (AMOC). The compilation of stratigraphic data from 108 marine sediment cores collected across the western tropical Atlantic shows an extreme rise in sedimentation rates off the Parnaíba River mouth (about 2°S) during Heinrich Stadial 1 (HS1, 18-15 ka). Sediment core GeoB16206-1, raised offshore the Parnaíba River mouth, documents relatively constant 143Nd/144Nd values (expressed as epsilonNd(0)) throughout the last 30 ka. Whereas the homogeneous epsilonNd(0) data support the input of fluvial sediments by the Parnaíba River from the same source area directly onshore, the increases in Fe/Ca, Al/Si and Rb/Sr during HS1 indicate a marked intensification of fluvial erosion in the Parnaíba River drainage basin. In contrast, the epsilonNd(0) values from sediment core GeoB16224-1 collected off French Guiana (about 7°N) suggest Amazon-sourced materials within the last 30 ka. We attribute the extremely high volume of terrigenous sediments deposited offshore the Parnaíba River mouth during HS1 to (i) an enhanced precipitation in the catchment region and (ii) a reduced North Brazil Current, which are both associated with a weakened AMOC.

Palynological and geochemical data of sediment core GeoB3910-2 located offshore Northeast Brazil

High resolution palynological and geochemical data of sediment core GeoB 3910-2 (located offshore Northeast Brazil) spanning the period between 19 600 and 14 500 calibrated year bp (19.6-14.5 ka) show a land-cover change in the catchment area of local rivers in two steps related to changes in precipitation associated with Heinrich Event 1 (H1 stadial). At the end of the last glacial maximum, the landscape in semi-arid Northeast Brazil was dominated by a very dry type of caatinga vegetation, mainly composed of grasslands with some herbs and shrubs. After 18 ka, considerably more humid conditions are suggested by changes in the vegetation and by Corg and C/N data indicative of fluvial erosion. The caatinga became wetter and along lakes and rivers, sedges and gallery forest expanded. The most humid period was recorded between 16.5 and 15 ka, when humid gallery (and floodplain) forest and even small patches of mountainous Atlantic rain forest occurred together with dry forest, the latter being considered as a rather lush type of caatinga vegetation. During this humid phase erosion decreased as less lithogenic material and more organic terrestrial material were deposited on the continental slope of northern Brazil. After 15 ka arid conditions returned. During the humid second phase of the H1 stadial, a rich variety of landscapes existed in Northeast Brazil and during the drier periods small pockets of forest could probably survive in favorable spots, which would have increased the resilience of the forest to climate change.

Grain size distribution in surface sediments samples of the continental shelf and slope of southeast Greenland

During "Meteor" Cruise 6/1966 in the northwest Atlantic a systematic survey of the bottom topography of the southeast Greenland continental margin was undertaken. Eighty-seven profiles transverse to the shelf edge at distances of 3-4 nautical miles and two longitudinal profiles parallel to the coast were carried out with the ELAC Narrow Beam Echo-Sounder giving a reliable record of even steep slopes. On the basis of the echo soundings the topography and morphology of the continental shelf and slope are evaluated. A detailed bathymetric chart and a serial profile chart were designed as working material for the morphological research. These maps along with the original echograms are morphometrically evaluated. The analysis of the sea bottom features is the basis of a subsequent morphogenetical interpretation, verified and extended by means of interpretation of magnetic data and sediment analysis (grain size, roundness, lithology). The results of the research are expressed in a geomorphological map. The primary findings can be summarized as follows: 1) The southeast Greenland shelf by its bottom topography can be clearly designated as a glacially formed area. The glacial features of the shelf can be classified into two zones nearly parallel to the coast: glacial erosion forms on the inner shelf and glacial accumulation forms on the outer shelf. The inner shelf is characterized by the rugged and hummocky topography of ice scoured plains with clear west/east slope asymmetry. On the outer shelf three types of glacial accumulation forms can be recognized: ice margin deposits with clearly expressed terminal moraines, glacial till plains and glaciomarine outwash fans. Both zones of the shelf can be subdivided into two levels of relief. The ice scoured plains, with average depths of 240 meters (m), are dissected to a maximum depth of 1060 m (Gyldenloves Trough) by trough valleys, which are the prolongations of the Greenland fjords. The banks of the outer shelf, with an average depth of 180 m, surround glacial basins with a maximum depth of 670 meters. 2) The sediments of the continental shelf can be classified as glacial due to their grain size distribution and the degree of roundness of the gravel particles. The ice margin deposits on the outer shelf can be recognized by their high percentage of gravels. On the inner shelf a rock surface is suggested, intermittently covered by glacial deposits. In the shelf troughs fine-grained sediments occur mixed with gravels. 3) Topography and sediments show that the southeast Greenland shelf was covered by an ice sheet resting on the sea floor during the Pleistocene ice-age. The large end moraines along the shelf edge probably indicate the maximum extent of the Wurm shelf ice resting on the sea floor. The breakthroughs of the end moraines in front of the glacial basins suggest that the shelf ice has floated further seaward over the increasing depths. 4) Petrographically the shelf sediments consist of gneisses, granites and basalts. While gneisses and granites occire on the nearby coast, basalt is not known to exist here. Either this material has been drifted by icebergs from the basalt province to the north or exists on the southeast Greenland shelf itself. The last interpretation is supported bythe high portion of basalt contained in the sediment samples taken and the strong magnetic anomalies probably caused by basaltic intrusions. 5) A magnetic profile allows the recognition of two magnetically differing areas which approximately coincide with the glacial erosion and accumulation zones. The inner shelf shows a strong and variable magnetic field because the glacially eroded basement forms the sea floor. The outer shelf is characterized by a weak and homogenous magnetic field, as the magnetized basement lies at greater depthy, buried by a thick cover of glacial sediments. The strong magnetic anomalies of the inner shelf are probably caused by dike swarms, similar to those observed further to the north in the Kangerdlugssuaq Fjord region. This interpretation is supported by the high basalt content of the sediment samples and the rough topography of the ice scoured plains which correlates in general with the magnetic fluctuations. The dike structures of the basement have been differentially eroded by the shelf ice. 6) The continental slope, extending from the shelf break at 313 m to a depth of 1270 m with an average slope of 11°, is characterized by delta-shaped projections in front of the shelf basins, by marginal plateaus, ridges and hills, by canyons and slumping features. The projections could be identified as glaciomarine sediment fans. This conclusion is supported by the strong decrease of magnetic field intensity. The deep sea hills and ridges with their greater magnetic intensities have to be regarded as basement outcrops projecting through the glaciomarine sediment cover. The upper continental rise, sloping seaward at about 2°, is composed of wide sediment fans and slump material. A marginal depression on the continental rise running parallel to the shelf edge has been identified. In this depression bottom currents capable of erosion have been recorded. South of Cape Farvel the depression extends to the accumulation zone of the "Eirik" sedimentary ridge. 7) By means of a study of the recent marine processes, postglacial modification of the ice-formed relief can be postulated. The retention effect of the fjord troughs and the high velocity of the East Greenland stream prevents the glacial features from being buried by sediments. Bottom currents capable of active erosion have only been found in the marginal depression on the continental rise. In addition, at the time of the lowest glacio-eustatic sea level, the shelf bottom was not situated in the zone of wave erosion. Only on the continental slope and rise bottom currents, sediment slumps and turbidity currents have led to significant recent modifications. Considering these results, the geomorphological development of the southeast Greenland continental terrace can be suggested as follows: 1. initial formation of a "peneplain", 2. fluvial incision, 3. submergence, and finally 4. glacial modification.

Reefs shift from net accretion to net erosion along a natural environmental gradient

Coral reefs persist in an accretion-erosion balance and ocean acidification resulting from anthropogenic CO2 emissions threatens to shift this balance in favor of net reef erosion. Corals and calcifying algae, largely responsible for reef accretion, are vulnerable to environmental changes associated with ocean acidification, but the direct effects of lower pH on reef erosion has received less attention, particularly in the context of known drivers of bioerosion and natural variability. This study examines the balance between reef accretion and erosion along a well-characterized natural environmental gradient in Kane'ohe Bay, Hawai'i using experimental blocks of coral skeleton. Comparing before and after micro-computed tomography (µCT) scans to quantify net accretion and erosion, we show that, at the small spatial scale of this study (tens of meters), pH was a better predictor of the accretion-erosion balance than environmental drivers suggested by prior studies, including resource availability, temperature, distance from shore, or depth. In addition, this study highlights the fine-scale variation of pH in coastal systems and the importance of microhabitat variation for reef accretion and erosion processes. We demonstrate significant changes in both the mean and variance of pH on the order of meters, providing a local perspective on global increases in pCO2. Our findings suggest that increases in reef erosion, combined with expected decreases in calcification, will accelerate the shift of coral reefs to an erosion-dominated system in a high-CO2 world. This shift will make reefs increasingly susceptible to storm damage and sea-level rise, threatening the maintenance of the ecosystem services that coral reefs provide.

Geometric parameters of shoreface profiles of Laptev Sea erosion coast (Table 1)

Field investigations of the Laptev Sea shoreface morphology were carried out (1) off erosional shores composed of unconsolidated sediments, (2) off the modern delta shores of the Lena River, and (3) off rocky shores. It was found that profiles off erosional shores had a concave shape. This shape is not well described by commonly applied power functions, a feature, which is in disagreement with the generally accepted concept of the equilibrium shape of shoreface profiles. The position of the lower shoreface boundary is determined by the elevation of the coastal lowland inundated during the last transgression (at -5 to -10 m) and may easily be recognized by a sharp, an order of magnitude decrease in the mean inclination of the sea floor. The mean shoreface inclination depends on sediment grain-size and ranges from 0.0022 to 0.033. The concave shape of the shoreface did not change substantially during the last 20-30 years, which indicates that shoreline retreat did not slow down and hence suggests continued intensive coastal erosion in the 21st century. The underwater part of the Lena River delta extends up to 35 km offshore. Its upper part is formed by a shallow and up to 18-km wide bench, which reaches depths of 2-3 m along the outer edge. The evolution of the delta was irregular. Whereas some parts of the delta are advancing rapidly (58 m/year), other parts are eroding. Comparison of measured profiles with older bathymetric data gave an opportunity to evaluate the changes of the underwater delta over past decades. Bathymetric surveys of the seabed around the delta can thus contribute towards a quantification of the sediment budget of the river-sea system. In addition, some sections of the Laptev Sea coast are composed of bedrock that has a comparatively low resistance to wave erosion. These sections may supply a considerable amount of sediment, especially if the cliffs are high. This source must therefore also be taken into account when assessing the contribution of shore erosion to the Laptev Sea sediment budget.

Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

X-ray fluorescence scanning of discrete samples on the Schwalbenberg II loess-paleosol sequence, raw data, magnetic susceptibility, organic carbon and U-ratio as grainsize index

The Schwalbenberg II loess-paleosol sequence (LPS) denotes a key site for Marine Isotope Stage (MIS 3) in Western Europe owing to eight succeeding cambisols, which primarily constitute the Ahrgau Subformation. Therefore, this LPS qualifies as a test candidate for the potential of temporal high-resolution geochemical data obtained X-ray fluorescence (XRF) scanning of discrete samplesproviding a fast and non-destructive tool for determining the element composition. The geochemical data is first contextualized to existing proxy data such as magnetic susceptibility (MS) and organic carbon (Corg) and then aggregated to element log ratios characteristic for weathering intensity [LOG (Ca/Sr), LOG (Rb/Sr), LOG (Ba/Sr), LOG (Rb/K)] and dust provenance [LOG (Ti/Zr), LOG (Ti/Al), LOG (Si/Al)]. Generally, an interpretation of rock magnetic particles is challenged in western Europe, where not only magnetic enhancement but also depletion plays a role. Our data indicates leaching and top-soil erosion induced MS depletion at the Schwalbenberg II LPS. Besides weathering, LOG (Ca/Sr) is susceptible for secondary calcification. Thus, also LOG (Rb/Sr) and LOG (Ba/Sr) are shown to be influenced by calcification dynamics. Consequently, LOG (Rb/K) seems to be the most suitable weathering index identifying the Sinzig Soils S1 and S2 as the most pronounced paleosols for this site. Sinzig Soil S3 is enclosed by gelic gleysols and in contrast to S1 and S2 only initially weathered pointing to colder climate conditions. Also the Remagen Soils are characterized by subtle to moderate positive excursions in the weathering indices. Comparing the Schwalbenberg II LPS with the nearby Eifel Lake Sediment Archive (ELSA) and other more distant German, Austrian and Czech LPS while discussing time and climate as limiting factors for pedogenesis, we suggest that the lithologically determined paleosols are in-situ soil formations. The provenance indices document a Zr-enrichment at the transition from the Ahrgau to the Hesbaye Subformation. This is explained by a conceptual model incorporating multiple sediment recycling and sorting effects in eolian and fluvial domains.

Clay mineralogy of the Cenozoic Pagodroma Group, East Antarctica

The Cenozoic Pagodroma Group in the northern Prince Charles Mountains, East Antarctica, is a glaciomarine succession of fjordal character, comprising four uplifted formations of different ages. The composition of the <2 µm fraction of sediments of the Pagodroma Group was analysed in order to help identify source areas, past weathering conditions and glacial regimes. Both clay and non-clay minerals have been quantified. The assemblage of the upper Oligocene to lower Miocene Mount Johnston Formation is characterised by the dominance of illite and intermediate concentrations of chlorite. Similar to that assemblage is the clay mineral suite of the middle Miocene Fisher Bench Formation, where illite and chlorite together account for 95% of the clay minerals. The middle to upper Miocene Battye Glacier Formation is the only formation with significant and persistent smectite concentrations, although illite is still dominant. The kaolinite concentration is also high and is even higher than that of chlorite. The clay fraction of the upper Pliocene to lower Pleistocene Bardin Bluffs Formation is characterised by maximum kaolinite concentrations and relatively low illite and chlorite concentrations. The bulk of the clay fraction in each formation can be explained by the physical weathering and erosion of a nearby source under glacial conditions. In the case of Mount Johnston Formation and Fisher Bench Formation this source may be situated in the metavolcanic and gneissic rocks of Fisher Massif. The sediments of the Bardin Bluffs Formation indicate a local source within the Amery Oasis, where Proterozoic granitoid rocks and gneisses, and Permo-Triassic fluvial rocks of the Amery Group are exposed. These results suggest a strong local imprint on the glacial sediments as northwards flowing ice eroded the bedrock in these areas. The origin of the clay fraction of the Battye Glacier Formation is a matter of debate. The smectite and kaolinite content most easily can be explained by erosion of sources largely hidden beneath the ice upstream. Less likely, these clay minerals reflect climatic conditions that were much warmer and wetter than today, facilitating chemical weathering.

Extensive management of field margins enhances their potential to mitigate soil erosion

Soil erosion is a widespread problem in agricultural landscapes, particularly in regions with strong rainfall events. Vegetated field margins can mitigate negative impacts of soil erosion by trapping eroded material. In this data set, we present data of sediment trapped by 12 field margins during the monsoon season of 2013 in an agricultural landscape in the Haean-myun catchment in South Korea. Prior to the beginning of monsoon season, we equipped a total of 12 sites representing three replicates for each of four different types of field margins ("managed flat", "managed steep", "natural flat" and "natural steep") with Astroturf mats with a size of 34 cm x 25 cm (850 cm**2). The mats (n = 15 / site) were installed at three levels: upslope, immediately before the field margin to quantify the sediments that reach it, in the middle of the field margin to quantify the locally trapped sediments, and after the field margin at the downslope edge to quantify the sediments that leave the field margin to the next field or to the stream. Sediment was collected after each rain event until the end of the monsoon season.

Coastal erosion dynamics and storm characteristics of Bykovsky Peninsula

This study investigates the rate of erosion during the 1951-2006 period on the Bykovsky Peninsula, located north-east of the harbour town of Tiksi, north Siberia. Its coastline, which is characterized by the presence of ice-rich sediment (Ice Complex) and the vicinity of the Lena River Delta, retreated at a mean rate of 0.59 m/yr between 1951 and 2006. Total erosion ranged from 434 m of erosion to 92 m of accretion during these 56 years and exhibited large variability (sigma = 45.4). Ninety-seven percent of the rates observed were less than 2 m/yr and 81.6% were less than 1 m/yr. No significant trend in erosion could be recorded despite the study of five temporal subperiods within 1951-2006. Erosion modes and rates actually appear to be strongly dependant on the nature of the backshore material, erosion being stronger along low-lying coastal stretches affected by past or current thermokarst activity. The juxtaposition of wind records monitored at the town of Tiksi and erosion records yielded no significant relationship despite strong record amplitude for both data sets. We explain this poor relationship by the only rough incorporation of sea-ice cover in our storm extraction algorithm, the use of land-based wind records vs. offshore winds, the proximity of the peninsula to the Lena River Delta freshwater and sediment plume and the local topographical constraints on wave development.

Pollen analysis, particle size analysis, end-member modelling and x-ray fluorescence measurements from sediment cores off Cape Ghir, Morocco

Understanding past human-climate-environment interactions is essential for assessing the vulnerability of landscapes and ecosystems to future climate change. This is particularly important in southern Morocco where the current vegetation is impacted by pastoralism, and the region is highly sensitive to climate variability. Here, we present a 2000-year record of vegetation, sedimentation rate, XRF chemical element intensities, and particle size from two decadal-resolved, marine sediment cores, raised from offshore Cape Ghir, southern Morocco. The results show that between 650 and 850 AD the sedimentation rate increased dramatically from 100 cm/1000 years to 300 cm/1000 years, and the Fe/Ca and pollen flux doubled, together indicating higher inputs of terrestrial sediment. Particle size measurements and end-member modelling suggest increased fluvial transport of the sediment. Beginning at 650 AD pollen levels from Cichorioideae species show a sharp rise from 10% to 20%. Pollen from Atemisia and Plantago, also increase from this time. Deciduous oak pollen percentages show a decline, whereas those of evergreen oak barely change. The abrupt increase in terrestrial/fluvial input from 650 to 850 AD occurs, within the age uncertainty, of the arrival of Islam (Islamisation) in Morocco at around 700 AD. Historical evidence suggests Islamisation led to population increase and development of southern Morocco, including expanded pastoralism, deforestation and agriculture. Livestock pressure may have changed the vegetation structure, accounting for the increase in pollen from Cichorioideae, Plantago, and Artemisia, which include many weedy species. Goats in particular may have played a dominant role as agents of erosion, and intense browsing may have led to the decline in deciduous oak; evergreen oak is more likely to survive as it re-sprouts more vigorously after browsing. From 850 AD to present sedimentation rates, Fe/Ca ratios and fluvial discharge remain stable, whereas pollen results suggest continued degradation. Pollen results from the past 150 years suggest expanded cultivation of olives and the native argan tree, and the introduction of Australian eucalyptus trees. The rapidly increasing population in southern Morocco is causing continued pressure to expand pastoralism and agriculture. The history of land degradation presented here suggests that the vegetation in southern Morocco may have been degraded for a longer period than previously thought and may be particularly sensitive to further land use changes. These results should be included in land management strategies for southern Morocco.