Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle ...failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.
Abstract Ambient seismic recordings taken at broad locations across Ross Ice Shelf and a dense array near West Antarctic Ice Sheet (WAIS) Divide, Antarctica, show pervasive temporally variable ...resonance peaks associated with trapped seismic waves in near-surface firn layers. These resonance peaks feature splitting on the horizontal components, here interpreted as frequency-dependent anisotropy in the firn and underlying ice due to several overlapping mechanisms driven by ice flow. Frequency peak splitting magnitudes and fast/slow axes were systematically estimated at single stations using a novel algorithm and compared with good agreement with active source anisotropy measurements at WAIS Divide determined via active sources recorded on a 1 km circular array. The approach was further applied to the broad Ross Ice Shelf (RIS) array, where anisotropy axes were directly compared with visible surface features and ice shelf flow lines. The near-surface firn, depicted by anisotropy above 30 Hz, was shown to exhibit a novel plastic stretching mechanism of anisotropy, whereby the fast direction in snow aligns with accelerating ice shelf flow.
This paper presents the study of modern pollen analogs from the Balearic Islands. While similar studies have been largely applied to mainland areas, research focused on modern vegetation dynamics on ...Mediterranean islands remains very rare. In this research, we combine vegetation surveys, pollen analysis and multivariate statistics to understand landscape composition. The main objectives of are: (1) to examine pollen-vegetation relationships in relation to environmental and land-use variables; (2) to understand modern pollen representation in a mosaic landscape structure; and (3) to propose pollen indicators that characterize the primary vegetation types from the Balearic Islands to better interpret past pollen records in Mediterranean island environments. Pollen results and Redundancy Analysis (RDA) distinguish three major groups: (a) Holm oak and box formations; (b) maquis and garrigues; and (c) anthropogenic and open habitats. Landscape form, mean decadal rainfall, mean decadal temperature, fire activity, trampling, slope percentage, wet/flooded soil, saline soil, distance to agropastoral cells, gHM index, domestic herbivory presence, agropastoral use, and soil type are the major variables explaining modern pollen assemblage variation in our research. Poaceae undiff., Plantago sp., Apiaceae undiff., Cerealia-t, and Cichorieae are highly correlated to human activities but should be interpreted cautiously when occurring in low values. Quercus ilex-t, Hypericum, and Buxus are correlated to humid locations while Pistacia, Pinus, Juniperus-t, and Olea to high mean decadal temperatures. Our study indicates how pollen analysis and multivariate analysis are powerful tools for characterizing the mosaic landscape, with special focus on the main vegetation types of the Balearic Islands.
This open access volume is the first comprehensive assessment of the Hindu Kush Himalaya (HKH) region. It comprises important scientific research on the social, economic, and environmental pillars of ...sustainable mountain development and will serve as a basis for evidence-based decision-making to safeguard the environment and advance people’s well-being. The compiled content is based on the collective knowledge of over 300 leading researchers, experts and policymakers, brought together by the Hindu Kush Himalayan Monitoring and Assessment Programme (HIMAP) under the coordination of the International Centre for Integrated Mountain Development (ICIMOD). This assessment was conducted between 2013 and 2017 as the first of a series of monitoring and assessment reports, under the guidance of the HIMAP Steering Committee: Eklabya Sharma (ICIMOD), Atiq Raman (Bangladesh), Yuba Raj Khatiwada (Nepal), Linxiu Zhang (China), Surendra Pratap Singh (India), Tandong Yao (China) and David Molden (ICIMOD and Chair of the HIMAP SC). This First HKH Assessment Report consists of 16 chapters, which comprehensively assess the current state of knowledge of the HKH region, increase the understanding of various drivers of change and their impacts, address critical data gaps and develop a set of evidence-based and actionable policy solutions and recommendations. These are linked to nine mountain priorities for the mountains and people of the HKH consistent with the Sustainable Development Goals. This book is a must-read for policy makers, academics and students interested in this important region and an essentially important resource for contributors to global assessments such as the IPCC reports. ; Constitutes the first comprehensive assessment of the Hindu Kush Himalaya region, providing an authoritative overview of the region Assembles the collective knowledge of over 300 leading researchers, practitioners, experts, and policymakers Combines the current state of knowledge of the Hindu Kush Himalaya region in one volume Offers Open Access to a set of practically oriented policy recommendations
Coseismic landslides pose immediate and prolonged hazards to mountainous communities, and provide a rare opportunity to study the effect of large earthquakes on erosion and sediment budgets. By ...mapping landslides using high-resolution satellite imagery, we find that the 25 April 2015 Mw7.8 Gorkha earthquake and aftershock sequence produced at least 25,000 landslides throughout the steep Himalayan Mountains in central Nepal. Despite early reports claiming lower than expected landslide activity, our results show that the total number, area, and volume of landslides associated with the Gorkha event are consistent with expectations, when compared to prior landslide-triggering earthquakes around the world. The extent of landsliding mimics the extent of fault rupture along the east-west trace of the Main Himalayan Thrust and increases eastward following the progression of rupture. In this event, maximum modeled Peak Ground Acceleration (PGA) and the steepest topographic slopes of the High Himalaya are not spatially coincident, so it is not surprising that landslide density correlates neither with PGA nor steepest slopes on their own. Instead, we find that the highest landslide density is located at the confluence of steep slopes, high mean annual precipitation, and proximity to the deepest part of the fault rupture from which 0.5–2Hz seismic energy originated. We suggest that landslide density was determined by a combination of earthquake source characteristics, slope distributions, and the influence of precipitation on rock strength via weathering and changes in vegetation cover. Determining the relative contribution of each factor will require further modeling and better constrained seismic parameters, both of which are likely to be developed in the coming few years as post-event studies evolve. Landslide mobility, in terms of the ratio of runout distance to fall height, is comparable to small volume landslides in other settings, and landslide volume-runout scaling is consistent with compilations of data on larger slope failures. In general, the size ratios of landslide source area to full landslide area are smaller than global averages, and hillslope length seems to largely control runout distance, which we propose reflects a topographic control on landslide mobility in this setting. We find that landslide size dictates runout distance and that more than half of the landslide debris was deposited in direct connection with stream channels. Connectivity, which is defined as the spatial proximity of landslides to fluvial channels, is greatest for larger landslides in the high-relief part of the High Himalaya. Although these failures are less abundant than those at lower elevations, they may have a disproportionate impact on sediment dynamics and cascading hazards, such as landslide reactivation by monsoon rainfall and landslide dams that lead to outburst floods. The overall high fluvial connectivity of coseismic landsliding in the Gorkha event suggests coupling between the earthquake cycle and sediment/geochemical budgets of fluvial systems in the Himalaya.
The timing of sea ice retreat and advance in Arctic coastal waters variessubstantially from year to year. Various activities, ranging from marinetransport to the use of sea ice as a platform for ...industrial activity orwinter travel, are affected by variations in the timing of breakup andfreeze-up, resulting in a need for indicators to document the regional andtemporal variations in coastal areas. The primary objective of this study isto use locally based metrics to construct indicators of breakup andfreeze-up in the Arctic and subarctic coastal environment. The indicatorsdeveloped here are based on daily sea ice concentrations derived fromsatellite passive-microwave measurements. The “day of year” indicators aredesigned to optimize value for users while building on past studiescharacterizing breakup and freeze-up dates in the open pack ice. Relativeto indicators for broader adjacent seas, the coastal indicators generallyshow later breakup at sites known to have landfast ice. The coastalindicators also show earlier freeze-up at some sites in comparison withfreeze-up for broader offshore regions, likely tied to earlier freezing ofshallow-water regions and areas affected by freshwater input from nearbystreams and rivers. A factor analysis performed to synthesize the localindicator variations shows that the local breakup and freeze-up indicatorshave greater spatial variability than corresponding metrics based onregional ice coverage. However, the trends towards earlier breakup andlater freeze-up are unmistakable over the post-1979 period in thesynthesized metrics of coastal breakup and freeze-up and the correspondingregional ice coverage. The findings imply that locally defined indicatorscan serve as key links between pan-Arctic or global indicators such assea ice extent or volume and local uses of sea ice, with the potential toinform community-scale adaptation and response.
Accurate modeling of cryospheric surface albedo is essential for ourunderstanding of climate change as snow and ice surfaces regulate the globalradiative budget and sea-level through their albedo and ...massbalance. Although significant progress has been made using physicalprinciples to represent the dynamic albedo of snow, models of glacier icealbedo tend to be heavily parameterized and not explicitly connected withphysical properties that govern albedo, such as the number and size of airbubbles, specific surface area (SSA), presence of abiotic and biotic lightabsorbing constituents (LACs), and characteristics of any overlyingsnow. Here, we introduce SNICAR-ADv4, an extension of the multi-layertwo-stream delta-Eddington radiative transfer model with theadding–doubling solver that has been previously applied to represent snowand sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnelreflectance and transmittance between overlying snow and higher-densityglacier ice, scattering by air bubbles of varying sizes, and numerous typesof LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadbandalbedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03for bare and bubble-free ice under direct light. Our results indicate thatrepresenting ice with a density of 650 kg m-3 as snow with norefractive Fresnel layer, as done previously, generally overestimates theBBA by an average of 0.058. However, because mostnaturally occurring ice surfaces are roughened “white ice”, we recommendmodeling a thin snow layer over bare ice simulations. We find optimalagreement with measurements by representing cryospheric media with densitiesless than 650 kg m-3 as snow and larger-density media as bubbly icewith a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedoimpacts from LACs with changing ice SSA, with peak impact per unit mass ofLACs near SSAs of 0.1–0.01 m2 kg-1. For bare, bubble-free ice, LACsactually increase the albedo. SNICAR-ADv4 represents smooth transitionsbetween snow, firn, and ice surfaces and accurately reproduces measuredspectral albedos of a variety of glacier surfaces. This work paves the wayfor adapting SNICAR-ADv4 to be used in land ice model components of Earthsystem models.
Determining the rates, mechanisms, and geographic variability of relative sea-level (RSL) change following the Last Glacial Maximum (LGM) provides insight into the sensitivity of ice sheets to ...climate change, the response of the solid Earth and gravity field to ice-mass redistribution, and constrains statistical and physical models used to project future sea-level rise. To do so in a scientifically robust way requires standardized datasets that enable broad spatial comparisons that minimize bias. As part of a larger goal to develop a unified, spatially-comprehensive post-LGM global RSL database, in this special issue we provide a standardized global synthesis of regional RSL data that resulted from the first ‘Geographic variability of HOLocene relative SEA level (HOLSEA)’ meetings in Mt Hood, Oregon (2016) and St Lucia, South Africa (2017). The HOLSEA meetings brought together sea-level researchers to agree upon a consistent protocol to standardize, interpret, and incorporate realistic uncertainties of RSL data. This special issue provides RSL data from ten geographical regions including new databases from Atlantic Europe and the Russian Arctic and revised/expanded databases from Atlantic Canada, the British Isles, the Netherlands, the western Mediterranean, the Adriatic, Israel, Peninsular Malaysia, Southeast Asia, and the Indian Ocean. In total, the database derived from this special issue includes 5634 (5290 validated) index (n = 3202) and limiting points (n = 2088) that span from ∼20,000 years ago to present. Progress in improving the standardization of sea-level databases has also been accompanied by advancements in statistical and analytical methods used to infer spatial patterns and rates of RSL change from geological data that have a spatially and temporally sparse distribution and geochronological and elevational uncertainties. This special issue marks the inception of a unified, spatially-comprehensive post-LGM global RSL database.
Display omitted
The 2016 National Land Cover Database (NLCD) product suite (available on www.mrlc.gov), includes Landsat-based, 30 m resolution products over the conterminous (CONUS) United States ...(U.S.) for land cover, urban imperviousness, and tree, shrub, herbaceous and bare ground fractional percentages. The release of NLCD 2016 provides important new information on land change patterns across CONUS from 2001 to 2016. For land cover, seven epochs were concurrently generated for years 2001, 2004, 2006, 2008, 2011, 2013, and 2016. Products reveal that land cover change is significant across most land cover classes and time periods. The land cover product was validated using existing reference data from the legacy NLCD 2011 accuracy assessment, applied to the 2011 epoch of the NLCD 2016 product line. The legacy and new NLCD 2011 overall accuracies were 82% and 83%, respectively, (standard error (SE) was 0.5%), demonstrating a small but significant increase in overall accuracy. Between 2001 and 2016, the CONUS landscape experienced significant change, with almost 8% of the landscape having experienced a land cover change at least once during this period. Nearly 50% of that change involves forest, driven by change agents of harvest, fire, disease and pests that resulted in an overall forest decline, including increasing fragmentation and loss of interior forest. Agricultural change represented 15.9% of the change, with total agricultural spatial extent showing only a slight increase of 4778 km2, however there was a substantial decline (7.94%) in pasture/hay during this time, transitioning mostly to cultivated crop. Water and wetland change comprised 15.2% of change and represent highly dynamic land cover classes from epoch to epoch, heavily influenced by precipitation. Grass and shrub change comprise 14.5% of the total change, with most change resulting from fire. Developed change was the most persistent and permanent land change increase adding almost 29,000 km2 over 15 years (5.6% of total CONUS change), with southern states exhibiting expansion much faster than most of the northern states. Temporal rates of developed change increased in 2001–2006 at twice the rate of 2011–2016, reflecting a slowdown in CONUS economic activity. Future NLCD plans include increasing monitoring frequency, reducing latency time between satellite imaging and product delivery, improving accuracy and expanding the variety of products available in an integrated database.