ALBERT

Description: We compare differences and similarities in the annual stratospheric HNO3 cycle derived from ground‐based measurements at the South Pole during 1993 and 1995, after correcting an error in earlier published profile retrievals for 1993 which led to under estimation of mixing ratios. The data series presented here provide profiling over the range ∼16–48 km, and cover the fall‐winter‐spring cycle in the behavior of HNO3 in the extreme Antarctic with a large degree of temporal overlap. With the exception of one gap of 20 days, the combined data sets cover a full annual cycle. The record shows an increase in HNO3 above 30 km occurring about 20 days before sunset, which appears to be the result of higher altitude heterogeneous conversion of NOx as photolysis diminishes. Both years show a strong increase in HNO3 beginning about polar sunset, in a layer peaking at about 25 km, as additional NOx is heterogeneously converted to nitric acid. When temperatures drop to the polar stratospheric cloud (PSC) formation range near the end of May, gas phase HNO3 is rapidly reduced in the lower stratosphere, although at least 2–3 weeks of temperatures ≤192 K appear to be required to complete most of the gas‐phase removal at the upper end of the depletion range (22–25 km). Despite a significant difference in residual sulfate loading from the explosion of Mount Pinatubo, there appears to be little gross difference in the timing and effects of PSC formation in removing gas phase HNO3 in these 2 years, though removal may be more rapid in 1995. Incorporation of gas phase HNO3 into PSCs appears to be nearly complete up to ∼25 km by midwinter. We also see a repeat of the formation of gas phase HNO3 in the middle stratosphere in early midwinter of 1995 with about the same timing as in 1993, suggesting that this phenomenon is driven by a repetition of dynamical transport and appropriate temperatures and pressures in the polar night, and not (as has been suggested) by ion‐based heterogeneous chemistry that requires triggering by large relativistic electron fluxes. High‐altitude HNO3 production peaks during a period of ∼20 days, but appears to persist for up to ∼40 days in the 40–45 km range, ceasing well before sunrise. This HNO3 descends rapidly throughout the production period, at a rate in good agreement with theoretically determined midwinter subsidence rates. As noted in earlier studies, later warming of this region above PSC evaporation temperatures does not cause reappearance of large amounts of HNO3, indicating that most PSCs gravitationally sink out of the stratosphere before early spring. We present evidence that smaller PSCs do evaporate to ∼1 to 3.5 ppbv of HNO3 in the lower stratosphere, however, working downward from ∼25 km as temperatures rise during the late winter. There is a delay of ∼15 days after sunrise before photolysis causes significant depletion in the altitude range below ∼30 km, where subsidence has carried virtually all higher‐altitude HNO3 by polar sunrise. Some continued subsidence and photolysis combine to keep mixing ratios less than ∼5 ppbv below 30 km until the final breakdown of the vortex in November brings larger amounts of HNO3 with air from lower latitudes.

Description: Published

Description: 17739-17750

Description: 5A. Ricerche polari e paleoclima

Description: JCR Journal

Description: Seismic hazard modeling is a multidisciplinary science that aims to forecast earthquake occurrence and its resultant ground shaking. Such models consist of a probabilistic framework that quantifies uncertainty across a complex system; typically, this includes at least two model components developed from Earth science: seismic source and ground motion models. Although there is no scientific prescription for the forecast length, the most common probabilistic seismic hazard analyses consider forecasting windows of 30 to 50 years, which are typically an engineering demand for building code purposes. These types of analyses are the topic of this review paper. Although the core methods and assumptions of seismic hazard modeling have largely remained unchanged for more than 50 years, we review the most recent initiatives, which face the difficult task of meeting both the increasingly sophisticated demands of society and keeping pace with advances in scientific understanding. A need for more accurate and spatially precise hazard forecasting must be balanced with increased quantification of uncertainty and new challenges such as moving from time‐independent hazard to forecasts that are time dependent and specific to the time period of interest. Meeting these challenges requires the development of science‐driven models, which integrate all information available, the adoption of proper mathematical frameworks to quantify the different types of uncertainties in the hazard model, and the development of a proper testing phase of the model to quantify its consistency and skill. We review the state of the art of the National Seismic Hazard Modeling and how the most innovative approaches try to address future challenges.

Description: Published

Description: e2019RG000653

Description: 6T. Studi di pericolosità sismica e da maremoto

Description: JCR Journal

Description: An edited version of this paper was published by AGU. Copyright (2016) American Geophysical Union.

Description: The accuracy of earthquake locations and their correspondence with subsurface geology depends strongly on the accuracy of the available seismic velocity model. Most modern methods to construct a velocity model for earthquake location are based on the inversion of passive source seismological data. Another approach is the integration of high-resolution geological and geophysical data to construct deterministic velocity models in which earthquake locations can be directly correlated to the geological structures. Such models have to be kinematically consistent with independent seismological data in order to provide precise hypocenter solutions. We present the Altotiberina (AT) seismic model, a three-dimensional velocity model for the Upper Tiber Valley region (Northern Apennines, Italy), constructed by combining 300 km of seismic reflection profiles, 6 deep boreholes (down to 5 km depth), detailed data from geological surveys and direct measurements of P- and S-wave velocities performed in situ and in laboratory. We assess the robustness of the AT seismic model by locating 11,713 earthquakes with a non-linear, global-search inversion method and comparing the probabilistic hypocenter solutions to those calculated in three previously published velocity models, constructed by inverting passive seismological data only. Our results demonstrate that the AT seismic model is able to provide higher-quality hypocenter locations than the previous velocity models. Earthquake locations are consistent with the subsurface geological structures and show a high degree of spatial correlation with specific lithostratigraphic units, suggesting a lithological control on the seismic activity evolution.

Description: Published

Description: 8113-8135

Description: 4T. Sismicità dell'Italia

Description: JCR Journal

Description: [1] We present the first intercomparison between the two most comprehensive records of gas‐phase HNO3 profiles in the Antarctic stratosphere, covering the greater part of 1993 and 1995. We compare measurements by the Stony Brook Ground‐Based Millimeter‐wave Spectrometer (GBMS) at the South Pole with Version 5 HNO3 data from the Microwave Limb Sounder (MLS) aboard the Upper Atmospheric Research Satellite. Trajectory tracing was used to select MLS measurements in the 70°–80°S latitude band that sampled air observed by the GBMS during passage over the Pole. When temperatures were near the HNO3 condensation range, additional screening was performed to select MLS measurements that sampled air parcels within 1.5 K of the temperature they experienced over the Pole. Quantitative comparisons are given at 7 different potential temperature levels spanning the range ∼19–30 km. Agreement between the data sets is quite good between 465 and 655 K (∼20–25 km) during a large fraction of the year. Agreement is best during winter and spring, when seasonally averaged differences are generally within 1 ppbv below ∼25 km. At higher altitudes, and during summer and fall, the agreement becomes worse, and GBMS measurements can exceed MLS values by more than 3 ppbv. We provide evidence that differences occurring in the lower stratosphere during fall are due to lack of colocation between the two data sets during a period of strong poleward gradients in HNO3. Remaining discrepancies between GBMS and MLS V5 HNO3 measurements are thought to be due to instrumental or retrieval biases.

Description: Published

Description: id 4809

Description: 5A. Ricerche polari e paleoclima

Description: JCR Journal

Description: Producing probabilistic subseasonal forecasts of extreme events up to six weeks in advance is crucial for many economic sectors. In agribusiness, this time scale is particularly critical because it allows for mitigation strategies to be adopted for counteracting weather hazards and taking advantage of opportunities. For example, spring frosts are detrimental for many nut trees, resulting in dramatic losses at harvest time. To explore subseasonal forecast quality in boreal spring, identified as one of the most sensitive times of the year by agribusiness end users, we build a multisystem ensemble using four models involved in the Subseasonal to Seasonal Prediction project (S2S). Two-meter temperature forecasts are used to analyze cold spell predictions in the coastal Black Sea region, an area that is a global leader in the production of hazelnuts. When analyzed at the global scale, the multisystem ensemble probabilistic forecasts for near-surface temperature are better than climatological values for several regions, especially the tropics, even many weeks in advance; however, in the coastal Black Sea, skill is low after the second forecast week. When cold spells are predicted instead of near-surface temperatures, skill improves for the region, and the forecasts prove to contain potentially useful information to stakeholders willing to put mitigation plans into effect. Using a cost–loss model approach for the first time in this context, we show that there is added value of having such a forecast system instead of a business-as-usual strategy, not only for predictions released 1–2 weeks ahead of the extreme event, but also at longer lead times.

Description: Published

Description: 237–254

Description: 4A. Oceanografia e clima

Description: JCR Journal

Description: Mineral dust plays an important role in the atmospheric radiation budget as well as in the ocean carbon cycle through fertilization and by ballasting of settling organic matter. However, observational records of open‐ocean dust deposition are sparse. Here, we present the spatial and temporal evolution of Saharan dust deposition over 2 years from marine sediment traps across the North Atlantic, directly below the core of the Saharan dust plume, with highest dust fluxes observed in summer. We combined the observed deposition fluxes with model simulations and satellite observations and argue that dust deposition in the Atlantic is predominantly controlled by summer rains. The dominant depositional pathway changes from wet deposition in summer to dry deposition in winter. Wet deposition has previously been suggested to increase the release of dust‐derived nutrients and their bioavailability, which may be a key contributor to surface‐ocean productivity in remote and oligotrophic parts of the oceans.

Description: Streams and rivers are important components of the carbon cycle as they transport and transform dissolved organic matter (DOM). Using high‐resolution Fourier‐transform ion cyclotron resonance mass spectrometry, we studied the spatial distribution of DOM at the molecular level at more than 100 sites across a stream network during summer and winter baseflow. We developed a model approximating the time DOM spent in the fluvial network, a key constraint on the biogeochemical processing of DOM. Discharge‐weighted travel time explained the compositional changes of DOM, which differed markedly in summer and winter. We attribute these seasonal differences to variation in source material, putatively reflecting the dynamics of freshly produced DOM in summer and DOM with an imprint of leaf litter in winter. Hydrological mixing was an important driver of the spatial dynamics of DOM. From the convergence rate of DOM compound intensities to the network‐wide average, we inferred the spatial distribution of sources within the catchment. Finally, we estimated network‐wide apparent mass transfer coefficients (vf app) of individual DOM compounds, which describe the vertical velocity at which DOM compounds are removed by biotic and abiotic processes. We identified the oxidative state of carbon as an important factor explaining vf app, which we consequently attribute to biological uptake of thermodynamically favorable DOM compounds. This work contributes to our understanding of the spatial processes, temporal constraints, and chemical properties of DOM that regulate the transformation and diagenesis of DOM at the fluvial network scale.

Description: The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component contributing to the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The ice stream geometry and surface velocities in the onset region of the NEGIS are not yet sufficiently well reproduced by ice sheet models. We present an assessment of the basal conditions of the onset region in a systematic analysis of airborne ultra-wideband radar data. Our data yield a new detailed model of ice-thickness distribution and basal topography in the upstream part of the ice stream. We observe a change from a smooth to a rougher bed where the ice stream widens from 10 to 60 km, and a distinct roughness anisotropy, indicating a preferred orientation of subglacial structures. The observation of off-nadir reflections that are symmetrical to the bed reflection in the radargrams suggests that these structures are elongated subglacial landforms, which in turn indicate potential streamlining of the bed. Together with basal water routing pathways, our observations hint to two different zones in this part of the NEGIS: an accelerating and smooth upstream region, which is collecting water, with reduced basal traction, and in the further downstream part, where the ice stream is slowing down and is widening, with a distribution of basal water towards the shear margins. Our findings support the hypothesis that the NEGIS is strongly interconnected to the subglacial water system in its onset region, but also to the subglacial substrate and morphology.

Description: Changes in ocean gateway configuration can induce basin‐scale rearrangements in ocean current characteristics. However, there is large uncertainty in the relative timing of the Oligocene/Miocene subsidence histories of the Greenland‐Scotland Ridge (GSR) and the Fram Strait (FS). By using a climate model, we investigate the temperature and salinity changes in response to the subsidence of these two key ocean gateways during early to middle Miocene. For a singular subsidence of the GSR, we detect warming and a salinity increase in the Nordic Seas and the Arctic Ocean. As convection sites shift to the north of Iceland, North Atlantic Deep Water (NADW) is formed at cooler temperatures. The associated deep ocean cooling and upwelling of deep waters to the Southern Ocean surface can cause a cooling in the southern high latitudes. These characteristic responses to the GSR deepening are independent of the FS being shallow or deep. An isolated subsidence of the FS gateway for a deep GSR shows less pronounced warming and salinity increase in the Nordic Seas. Arctic temperatures remain unaltered, but a stronger salinity increase is detected, which further increases the density of NADW. The increase in salinity enhances the contribution of NADW to the abyssal ocean at the expense of the colder southern source water component. These relative changes largely counteract each other and cause a negligible warming in the upwelling regions of the Southern Ocean.

Description: The Antarctic Circumpolar Current (ACC) is the world's largest current system connecting all three major basins of the global ocean. Our knowledge of glacial‐interglacial changes in ACC dynamics in the southeast Pacific is not well constrained and presently only based on reconstructions covering the last glacial cycle. Here we use a combination of mean sortable silt grain size of the terrigenous sediment fraction (10–63 μm, "Sortable Silt") and X‐ray fluorescence scanner‐derived Zr/Rb ratios as flow strength proxies to examine ACC variations at the Pacific entrance to the Drake Passage (DP) in the vicinity of the Subantarctic Front. Our results indicate that at the DP entrance, ACC strength varied by ~6–16% on glacial‐interglacial time scales, yielding higher current speeds during interglacial times and reduced current speeds during glacials. We provide evidence that previous observations of a reduction in DP throughflow during the last glacial period are part of a consistent pattern extending for at least the last 1.3 Ma. The orbital‐scale cyclicity follows well‐known global climate changes from prevailing ca. 41‐kyr cycles in the early part of the record (1.3 Ma to 850 ka; marine isotope stage 21) across the mid‐Pleistocene transition into the middle and late Pleistocene 100‐kyr world. A comparison to a bottom water flow record from the deep western boundary current off New Zealand (Ocean Drilling Program Site 1123) reveals anti‐phased changes between the two sites. The enhanced supply of deep water along the DP and into the Atlantic Ocean during interglacials corresponds to a weakened flow of the SW Pacific deep western boundary current.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 33(4), (2020): 1535-1545, doi:10.1175/JCLI-D-19-0547.1.

Description: In a transient warming scenario, the North Atlantic is influenced by a complex pattern of surface buoyancy flux changes that ultimately weaken the Atlantic meridional overturning circulation (AMOC). Here we study the AMOC response in the CMIP5 experiment, using the near-geostrophic balance of the AMOC on interannual time scales to identify the role of temperature and salinity changes in altering the circulation. The thermal wind relationship is used to quantify changes in the zonal density gradients that control the strength of the flow. At 40°N, where the overturning cell is at its strongest, weakening of the AMOC is largely driven by warming between 1000- and 2000-m depth along the western margin. Despite significant subpolar surface freshening, salinity changes are small in the deep branch of the circulation. This is likely due to the influence of anomalously salty water in the subpolar intermediate layers, which is carried northward from the subtropics in the upper limb of the AMOC. In the upper 1000 m at 40°N, salty anomalies due to increased evaporation largely cancel the buoyancy increase due to warming. Therefore, in CMIP5, temperature dynamics are responsible for AMOC weakening, while freshwater forcing instead acts to strengthen the circulation in the net. These results indicate that past modeling studies of AMOC weakening, which rely on freshwater hosing in the subpolar gyre, may not be directly applicable to a more complex warming scenario.

Description: We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. We also thank John Marshall for helpful discussions on the driving mechanisms of the AMOC, and three anonymous reviewers whose comments greatly improved the manuscript. This work was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Award 80NSSC17K0372, and by National Science Foundation Award OCE-1433132.

Description: 2020-07-20

Description: Methane emissions from northern high latitude wetlands constitute a major uncertainty in the atmospheric methane (CH4) budget during the Holocene. To reconstruct northern wetland methane emissions, we used an empirical model based on syntheses of observations of peat initiation from more than 3600 radiocarbon-dated basal peat ages, plant-macrofossil-derived peatland type from more than 250 peat cores from sites across the northern high latitudes, and observed CH4 emissions averaged from modern-day wetland types in order to explore the effects of wetland expansion and changes in wetland type. Peatland basal ages and plant macrofossil records showed the widespread formation of fens in major northern wetland complexes before 8000 BP. After 8000 BP, new fen formation continued, but widespread peatland succession (to bogs) and permafrost aggradation also occurred. Reconstructed CH4 emissions from peatlands increased rapidly between 10,600 BP and 6900 BP due to fen formation and expansion, then stabilized after 5000 BP at 42 ± 25 Tg CH4 y-1, as high methane-emitting fens transitioned to lower methane-emitting bogs and permafrost peatlands. Permafrost formation in northern peatlands after 1000 BP decreased CH4 emissions by 20% to 34 ± 21 Tg y-1 by the present day. Warming temperatures, changes in peatland hydrology, and permafrost thaw will likely change the magnitude of northern peatland emissions in the future.

Description: Infrastructure and anthropogenic impacts are expanding across the Arctic. A consistent record of human impact is required in order to quantify the changes and to assess climate change impacts on the communities. We derived a first panarctic satellite-based record of expanding infrastructure and anthropogenic impacts along all permafrost affected coasts (100 km buffer) within the H2020 project Nunataryuk based on Sentinel-1/2 satellite imagery. C-band synthetic aperture radar and multi-spectral information is combined through a machine learning framework. Depending on region, we identified up to 50% more information (human presence) than in OpenStreetMap. The combination with satellite records on vegetation change (specifically NDVI from Landsat since 2000) allowed quantification of recent expansion of infrastructure. Most of the expanded human presence occurred in Russia related predominantly to oil/gas industry. The majority of areas with human presence will be subject to thaw by mid-21st century based on ground temperature trends derived from the ESA CCI+ Permafrost time series (1997-2019). Of specific concern in this context are also settlements located directly at permafrost affected coasts. An efficient erosion rate monitoring scheme needs to be developed and combined with settlement records in order to assess the risk for local communities and infrastructure. Relevant progress in the framework of the ESA EO4PAC project will be discussed.

Description: Retrogressive thaw slumps (RTS) are typical landforms indicating processes of rapid thawing and degrading permafrost. Their abundance is increasing in many regions and quantifying their dynamics is of high importance for assessing geomorphic, hydrologic, and biogeochemical impacts of climate change in the Arctic. Here we present a deep-learning (DL) based semantic segmentation framework to detect RTS, using high-resolution multi-spectral PlanetScope, topographic (ArcticDEM elevation and slope), and medium-resolution multi-temporal Landsat Trend data. We created a highly automated processing pipeline, which is designed to allow reproducible results and to be flexible for multiple input data types. The processing workflow is based on the pytorch deep-learning framework and includes a variety of different segmentation architectures (UNet, UNet++, DeepLabV3), backbones and includes common data transformation techniques such as augmentation or normalization. We tested (training, validation) our DL based model in six different regions of 100 to 300 km² size across Canada, and Siberia. We performed a regional cross-validation (5 regions training, 1 region validation) to test the spatial robustness and transferability of the algorithm. Furthermore, we tested different architectures, backbones and loss-functions to identify the best performing and most robust parameter sets. For training the models we created a database of manually digitized and validated RTS polygons. The resulting model performance varied strongly between different regions with maximum Intersection over Union (IoU) scores between 0.15 and 0.58. The strong regional variation emphasizes the need for sufficiently large training data, which is representative of the diversity of RTS types. However, the creation of good training data proved to be challenging due to the fuzzy definition and delineation of RTS. We are further continuing to improve the usability and the functionality to add further datasets and classes. We will show first results from the upscaling beyond small test areas towards large spatial clusters of extensive RTS presence e.g. Peel Plateau in NW Canada.

Description: With the Earth’s climate rapidly warming, the Arctic represents one of the most vulnerable regions to environmental change. These northern high latitude regions experience intensified fire seasons and especially tundra fires are projected to become more frequent and severe. Fires in permafrost regions have extensive impacts, including the initiation of thermokarst (rapid thaw of ice-rich ground), as they combust the upper organic soil layers which provide insulation to the permafrost below. Rapid permafrost thaw is, thus, often observable in fire scars in the first years post-disturbance. In polygonal ice-wedge landscapes, this becomes most prevalent through melting ice wedges and degrading troughs. The further these ice wedges degrade, the more troughs will likely connect and build an extensive hydrological network with changing patterns and degrees of connectivity that influences hydrology and runoff. While subsiding troughs over melting ice wedges may host new ponds, an increasing connectivity may also subsequently lead to more drainage of ponds, which in turn can limit further thaw and help stabilize the landscape. To quantify the changes in such dynamic landscapes over large regions, highly automated methods are needed that allow extracting information on the geomorphic state and changes over time of ice-wedge trough networks from remote sensing data. We developed a computer vision algorithm to automatically derive ice-wedge polygonal networks and the current microtopography of the degrading troughs from very high resolution, airborne laserscanning-based digital terrain models. We represent the networks as graphs (a concept from the computer sciences to describe complex networks) and apply methods from graph theory to describe and quantify hydrological network characteristics of the changing landscape. In fire scars, we especially observe rapidly growing networks and fast micromorphological change in those degrading troughs. In our study, we provide a space-for-time substitution comparing fire scars throughout the Alaskan tundra of up to 70 years since the fire disturbance, to show how this type of disturbed landscape evolves over time.

Description: Using our custom visualization tool for multitemporal Landsat satellite imagery we discovered, to our knowledge, an undocumented mega-landslide in far-east Siberia, which occurred in summer 2017 (https://bit.ly/2WYRLM1; 61.55°N; 170.01°E). To create and visualize this unique dataset, we processed temporal trends of multispectral indices of 〉100,000 Landsat images for a period from 2000-2019 using the freely available Google Earth Engine cloud processing platform (https://ingmarnitze.users.earthengine.app/view/hotspottcvisapp). The megaslide has a size of 3.66 km² and using the ArcticDEM data we estimate a volume movement of ~20 Mm³. With this size and volume, the landslide is among the largest globally known in recent decades. The landslide is accompanied by a smaller one (0.31 km², 1 Mm³) about 600 m further east, which already occurred in summer 2015. The large landslide caused the formation of several small lakes by blocking two valleys with debris and within newly formed crevasses near the hilltop, which are still persisting as of August 2021. As this event occurred in a remote valley far from any settlement, no visible damage to infrastructure or human livelihoods was detected. The remoteness has likely led to being not detected, like many similar, albeit a lot smaller, erosion features in the Arctic permafrost region. In this presentation we will show the main properties of the landslide, potential trigger mechanisms in the traditional sense. As this region is located along the fringes of permafrost presence we will discuss its potential connection to the rapidly warming climate in the high latitudes. Further, we will discuss how such a large event remained undetected for several years. We discuss and highlight the value of our landscape change visualization tool based on Landsat trend analysis (see Nitze et al., AGU 2020), which helped us to identify this extreme event. With more and more available data sources, this tool in addition to automated image analysis (e.g. deep-learning) or seismic analysis will help to uncover the hidden processes and dynamics of the Earth’s surface.

Description: Several decades of research have provided insight into patterns of and controls on thermokarst initiation and expansion, yet studies tend to focus on individual types of thermokarst (i.e., thaw lake formation and subsequent drainage) in particular regions. Today, we are left with uneven knowledge about abrupt permafrost thaw both conceptually and regionally. The goal of this presentation is to summarize recent advancements in monitoring thermokarst and its impact on soil, vegetation, and water while also framing a call to action for the next decade of research. Over the next decade, permafrost researchers must align their efforts on several fronts to not only increase our knowledge about changing permafrost but to align this knowledge with key community and policy needs. To support climate change planning and adaptation, northern communities need future thaw vulnerability mapped at scales relevant to their needs, which will require a suite of downscaled and new mapping and remote sensing products. Thermokarst predisposition maps based on circumpolar datasets greatly overestimate the area vulnerable to thermokarst, which can lead to poor planning and climate anxiety. In some situations, existing mapping products may be useful for downscaling with more detailed input data. In other situations, entirely new approaches may be required to support local action. A second key need for community relevant research is the ability to detect and monitor early warning indicators of thermokarst. Such information is needed to support scenario planning and to help mitigate the risks to social, cultural, and physical infrastructure created by permafrost change. We are evaluating the potential for using changes in vegetation, wetting/drying and topography as early warning indicators of thermokarst, all of which can be remotely sensed. Finally, integrating fine-scale disturbances such as thermokarst into large scale models remains a key challenge but critical for supporting sound climate policy. While a diversity of permafrost modeling approaches is necessary, we outline guiding principles that will help enhance model comparisons, assimilation of simulated data across spatiotemporal scales, and the ability for policy decisions to be rapidly informed by emerging science on permafrost change.