National Institute of Polar Research

Abstract. We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km3) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.

We present a new common stratigraphic timescale for the North Greenland Ice Core Project (NGRIP) and GRIP ice cores. The timescale covers the period 7.9–14.8 kyr before present and includes the Bølling, Allerød, Younger Dryas, and early Holocene periods. We use a combination of new and previously published data, the most prominent being new high‐resolution Continuous Flow Analysis (CFA) impurity records from the NGRIP ice core. Several investigators have identified and counted annual layers using a multiparameter approach, and the maximum counting error is estimated to be up to 2% in the Holocene part and about 3% for the older parts. These counting error estimates reflect the number of annual layers that were hard to interpret, but not a possible bias in the set of rules used for annual layer identification. As the GRIP and NGRIP ice cores are not optimal for annual layer counting in the middle and late Holocene, the timescale is tied to a prominent volcanic event inside the 8.2 kyr cold event, recently dated in the DYE‐3 ice core to 8236 years before A. D. 2000 (b2k) with a maximum counting error of 47 years. The new timescale dates the Younger Dryas‐Preboreal transition to 11,703 b2k, which is 100–150 years older than according to the present GRIP and NGRIP timescales. The age of the transition matches the GISP2 timescale within a few years, but viewed over the entire 7.9–14.8 kyr section, there are significant differences between the new timescale and the GISP2 timescale. The transition from the glacial into the Bølling interstadial is dated to 14,692 b2k. The presented timescale is a part of a new Greenland ice core chronology common to the DYE‐3, GRIP, and NGRIP ice cores, named the Greenland Ice Core Chronology 2005 (GICC05). The annual layer thicknesses are observed to be log‐normally distributed with good approximation, and compared to the early Holocene, the mean accumulation rates in the Younger Dryas and Bølling periods are found to be 47 ± 2% and 88 ± 2%, respectively.

We present new analytical data of major and trace elements for the geological MPI‐DING glasses KL2‐G, ML3B‐G, StHs6/80‐G, GOR128‐G, GOR132‐G, BM90/21‐G, T1‐G, and ATHO‐G. Different analytical methods were used to obtain a large spectrum of major and trace element data, in particular, EPMA, SIMS, LA‐ICPMS, and isotope dilution by TIMS and ICPMS. Altogether, more than 60 qualified geochemical laboratories worldwide contributed to the analyses, allowing us to present new reference and information values and their uncertainties (at 95% confidence level) for up to 74 elements. We complied with the recommendations for the certification of geological reference materials by the International Association of Geoanalysts (IAG). The reference values were derived from the results of 16 independent techniques, including definitive (isotope dilution) and comparative bulk (e.g., INAA, ICPMS, SSMS) and microanalytical (e.g., LA‐ICPMS, SIMS, EPMA) methods. Agreement between two or more independent methods and the use of definitive methods provided traceability to the fullest extent possible. We also present new and recently published data for the isotopic compositions of H, B, Li, O, Ca, Sr, Nd, Hf, and Pb. The results were mainly obtained by high‐precision bulk techniques, such as TIMS and MC‐ICPMS. In addition, LA‐ICPMS and SIMS isotope data of B, Li, and Pb are presented.

The last two abrupt warmings at the onset of our present warm interglacial period, interrupted by the Younger Dryas cooling event, were investigated at high temporal resolution from the North Greenland Ice Core Project ice core. The deuterium excess, a proxy of Greenland precipitation moisture source, switched mode within 1 to 3 years over these transitions and initiated a more gradual change (over 50 years) of the Greenland air temperature, as recorded by stable water isotopes. The onsets of both abrupt Greenland warmings were slightly preceded by decreasing Greenland dust deposition, reflecting the wetting of Asian deserts. A northern shift of the Intertropical Convergence Zone could be the trigger of these abrupt shifts of Northern Hemisphere atmospheric circulation, resulting in changes of 2 to 4 kelvin in Greenland moisture source temperature from one year to the next.

Antarctic and Southern Ocean (ASO) marine ecosystems have been changing for at least the last 30 years, including in response to increasing ocean temperatures and changes in the extent and seasonality of sea ice; the magnitude and direction of these changes differ between regions around Antarctica that could see populations of the same species changing differently in different regions. This article reviews current and expected changes in ASO physical habitats in response to climate change. It then reviews how these changes may impact the autecology of marine biota of this polar region: microbes, zooplankton, salps, Antarctic krill, fish, cephalopods, marine mammals, seabirds, and benthos. The general prognosis for ASO marine habitats is for an overall warming and freshening, strengthening of westerly winds, with a potential pole-ward movement of those winds and the frontal systems, and an increase in ocean eddy activity. Many habitat parameters will have regionally specific changes, particularly relating to sea ice characteristics and seasonal dynamics. Lower trophic levels are expected to move south as the ocean conditions in which they are currently found move pole-ward. For Antarctic krill and finfish, the latitudinal breadth of their range will depend on their tolerance of warming oceans and changes to productivity. Ocean acidification is a concern not only for calcifying organisms but also for crustaceans such as Antarctic krill; it is also likely to be the most important change in benthic habitats over the coming century. For marine mammals and birds, the expected changes primarily relate to their flexibility in moving to alternative locations for food and the energetic cost of longer or more complex foraging trips for those that are bound to breeding colonies. Few species are sufficiently well studied to make comprehensive species-specific vulnerability assessments possible. Priorities for future work are discussed.

With the advent of the Heliophysics/Geospace System Observatory (H/GSO), a complement of multi-spacecraft missions and ground-based observatories to study the space environment, data retrieval, analysis, and visualization of space physics data can be daunting. The Space Physics Environment Data Analysis System (SPEDAS), a grass-roots software development platform (www.spedas.org), is now officially supported by NASA Heliophysics as part of its data environment infrastructure. It serves more than a dozen space missions and ground observatories and can integrate the full complement of past and upcoming space physics missions with minimal resources, following clear, simple, and well-proven guidelines. Free, modular and configurable to the needs of individual missions, it works in both command-line (ideal for experienced users) and Graphical User Interface (GUI) mode (reducing the learning curve for first-time users). Both options have "crib-sheets," user-command sequences in ASCII format that can facilitate record-and-repeat actions, especially for complex operations and plotting. Crib-sheets enhance scientific interactions, as users can move rapidly and accurately from exchanges of technical information on data processing to efficient discussions regarding data interpretation and science. SPEDAS can readily query and ingest all International Solar Terrestrial Physics (ISTP)-compatible products from the Space Physics Data Facility (SPDF), enabling access to a vast collection of historic and current mission data. The planned incorporation of Heliophysics Application Programmer's Interface (HAPI) standards will facilitate data ingestion from distributed datasets that adhere to these standards. Although SPEDAS is currently Interactive Data Language (IDL)-based (and interfaces to Java-based tools such as Autoplot), efforts are under-way to expand it further to work with python (first as an interface tool and potentially even receiving an under-the-hood replacement). We review the SPEDAS development history, goals, and current implementation. We explain its "modes of use" with examples geared for users and outline its technical implementation and requirements with software developers in mind. We also describe SPEDAS personnel and software management, interfaces with other organizations, resources and support structure available to the community, and future development plans. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s11214-018-0576-4) contains supplementary material, which is available to authorized users.

Abstract. Data from both laboratory studies and atmospheric measurements are used to develop an empirical parameterization for the immersion freezing activity of natural mineral dust particles. Measurements made with the Colorado State University (CSU) continuous flow diffusion chamber (CFDC) when processing mineral dust aerosols at a nominal 105% relative humidity with respect to water (RHw) are taken as a measure of the immersion freezing nucleation activity of particles. Ice active frozen fractions vs. temperature for dusts representative of Saharan and Asian desert sources were consistent with similar measurements in atmospheric dust plumes for a limited set of comparisons available. The parameterization developed follows the form of one suggested previously for atmospheric particles of non-specific composition in quantifying ice nucleating particle concentrations as functions of temperature and the total number concentration of particles larger than 0.5 μm diameter. Such an approach does not explicitly account for surface area and time dependencies for ice nucleation, but sufficiently encapsulates the activation properties for potential use in regional and global modeling simulations, and possible application in developing remote sensing retrievals for ice nucleating particles. A calibration factor is introduced to account for the apparent underestimate (by approximately 3, on average) of the immersion freezing fraction of mineral dust particles for CSU CFDC data processed at an RHw of 105% vs. maximum fractions active at higher RHw. Instrumental factors that affect activation behavior vs. RHw in CFDC instruments remain to be fully explored in future studies. Nevertheless, the use of this calibration factor is supported by comparison to ice activation data obtained for the same aerosols from Aerosol Interactions and Dynamics of the Atmosphere (AIDA) expansion chamber cloud parcel experiments. Further comparison of the new parameterization, including calibration correction, to predictions of the immersion freezing surface active site density parameterization for mineral dust particles, developed separately from AIDA experimental data alone, shows excellent agreement for data collected in a descent through a Saharan aerosol layer. These studies support the utility of laboratory measurements to obtain atmospherically relevant data on the ice nucleation properties of dust and other particle types, and suggest the suitability of considering all mineral dust as a single type of ice nucleating particle as a useful first-order approximation in numerical modeling investigations.

The low energy particle (LEP) instrument onboard GEOTAIL is designed to make comprehensive observations of plasma and energetic electrons and ions with fine temporal resolution in the terrestrial magnetosphere (mainly magnetotail) and in the interplanetary medium. It consists of three units of sensors (LEP-EA, LEP-SW and LEP-MS) and a common electronics (LEP-E). The Energy-per-charge Analyzers (EA) measure three-dimensional velocity distributions of electrons (with EA-e) and ions (with EA-i), simultaneously and separately, over the. energy-per-charge range of several eV/q to 43 keV/q. Emphasis in the EA design is laid on the large geometrical factor to measure tenuous plasma in the magnetotail with sufficient counting statistics in the high-time-resolution measurement. On the other hand, the Solar Wind ion analyzer (SW) has smaller geometrical factor, but fine angular and energy resolutions, to measure energy-per-charge spectra of the solar wind ions. In both EA and SW sensors, the complete three-dimensional velocity distributions can only be obtained in a period of four spins, while the velocity moments up to the third order are calculated onboard every spin period (nominally, 3 sec). The energetic-ion Mass Spectrometer (MS) can provide three-dimensional determinations of the ion composition. In this paper, we describe the instrumentation and present some examples of the inflight measurements.

Reproducible climate reconstructions of the Common Era (1 CE to present) are key to placing industrial-era warming into the context of natural climatic variability. Here we present a community-sourced database of temperature-sensitive proxy records from the PAGES2k initiative. The database gathers 692 records from 648 locations, including all continental regions and major ocean basins. The records are from trees, ice, sediment, corals, speleothems, documentary evidence, and other archives. They range in length from 50 to 2000 years, with a median of 547 years, while temporal resolution ranges from biweekly to centennial. Nearly half of the proxy time series are significantly correlated with HadCRUT4.2 surface temperature over the period 1850-2014. Global temperature composites show a remarkable degree of coherence between high- and low-resolution archives, with broadly similar patterns across archive types, terrestrial versus marine locations, and screening criteria. The database is suited to investigations of global and regional temperature variability over the Common Era, and is shared in the Linked Paleo Data (LiPD) format, including serializations in Matlab, R and Python.

Abstract. In spring 2006, the European Arctic was abnormally warm, setting new historical temperature records. During this warm period, smoke from agricultural fires in Eastern Europe intruded into the European Arctic and caused the most severe air pollution episodes ever recorded there. This paper confirms that biomass burning (BB) was indeed the source of the observed air pollution, studies the transport of the smoke into the Arctic, and presents an overview of the observations taken during the episode. Fire detections from the MODIS instruments aboard the Aqua and Terra satellites were used to estimate the BB emissions. The FLEXPART particle dispersion model was used to show that the smoke was transported to Spitsbergen and Iceland, which was confirmed by MODIS retrievals of the aerosol optical depth (AOD) and AIRS retrievals of carbon monoxide (CO) total columns. Concentrations of halocarbons, carbon dioxide and CO, as well as levoglucosan and potassium, measured at Zeppelin mountain near Ny Ålesund, were used to further corroborate the BB source of the smoke at Spitsbergen. The ozone (O3) and CO concentrations were the highest ever observed at the Zeppelin station, and gaseous elemental mercury was also elevated. A new O3 record was also set at a station on Iceland. The smoke was strongly absorbing – black carbon concentrations were the highest ever recorded at Zeppelin – and strongly perturbed the radiation transmission in the atmosphere: aerosol optical depths were the highest ever measured at Ny Ålesund. We furthermore discuss the aerosol chemical composition, obtained from filter samples, as well as the aerosol size distribution during the smoke event. Photographs show that the snow at a glacier on Spitsbergen became discolored during the episode and, thus, the snow albedo was reduced. Samples of this polluted snow contained strongly elevated levels of potassium, sulphate, nitrate and ammonium ions, thus relating the discoloration to the deposition of the smoke aerosols. This paper shows that, to date, BB has been underestimated as a source of aerosol and air pollution for the Arctic, relative to emissions from fossil fuel combustion. Given its significant impact on air quality over large spatial scales and on radiative processes, the practice of agricultural waste burning should be banned in the future.

Abstract A database of surface Antarctic snow isotopic composition is constructed using available measurements, with an estimate of data quality and local variability. Although more than 1000 locations are documented, the spatial coverage remains uneven with a majority of sites located in specific areas of East Antarctica. The database is used to analyze the spatial variations in snow isotopic composition with respect to geographical characteristics (elevation, distance to the coast) and climatic features (temperature, accumulation) and with a focus on deuterium excess. The capacity of theoretical isotopic, regional, and general circulation atmospheric models (including “isotopic” models) to reproduce the observed features and assess the role of moisture advection in spatial deuterium excess fluctuations is analyzed.

Abstract. The EPICA (European Project for Ice Coring in Antarctica) Dome C drilling in East Antarctica has now been completed to a depth of 3260 m, at only a few meters above bedrock. Here we present the new EDC3 chronology, which is based on the use of 1) a snow accumulation and mechanical flow model, and 2) a set of independent age markers along the core. These are obtained by pattern matching of recorded parameters to either absolutely dated paleoclimatic records, or to insolation variations. We show that this new time scale is in excellent agreement with the Dome Fuji and Vostok ice core time scales back to 100 kyr within 1 kyr. Discrepancies larger than 3 kyr arise during MIS 5.4, 5.5 and 6, which points to anomalies in either snow accumulation or mechanical flow during these time periods. We estimate that EDC3 gives accurate event durations within 20% (2σ) back to MIS11 and accurate absolute ages with a maximum uncertainty of 6 kyr back to 800 kyr.

Animal-attached remote sensing, or bio-logging, refers to the deployment of autonomous recording tags on free-living animals, so that multiple variables can be monitored at rates of many times per second, thereby generating millions of data points over periods ranging from hours to years. Rapid advances in technology are allowing scientists to use data-recording units to acquire huge, quantitative datasets of behavior from animals moving freely in their natural environment. In other words, scientists can examine wild animals in the field, behaving normally, with the same rigor that is normally used in the laboratory. The flexibility of such recording systems means that bio-logging science operates at the interface of several biological disciplines, looking at a wide array of aquatic, airborne, and terrestrial species, monitoring not only the physical characteristics of the environment, but also the animal's reactions to it. This approach is critically important in an era when global change threatens the survival of species and where habitat loss is leading to widespread extinctions.

An ethogram is a catalogue of discrete behaviors typically employed by a species. Traditionally animal behavior has been recorded by observing study individuals directly. However, this approach is difficult, often impossible, in the case of behaviors which occur in remote areas and/or at great depth or altitude. The recent development of increasingly sophisticated, animal-borne data loggers, has started to overcome this problem. Accelerometers are particularly useful in this respect because they can record the dynamic motion of a body in e.g. flight, walking, or swimming. However, classifying behavior using body acceleration characteristics typically requires prior knowledge of the behavior of free-ranging animals. Here, we demonstrate an automated procedure to categorize behavior from body acceleration, together with the release of a user-friendly computer application, "Ethographer". We evaluated its performance using longitudinal acceleration data collected from a foot-propelled diving seabird, the European shag, Phalacrocorax aristotelis. The time series data were converted into a spectrum by continuous wavelet transformation. Then, each second of the spectrum was categorized into one of 20 behavior groups by unsupervised cluster analysis, using k-means methods. The typical behaviors extracted were characterized by the periodicities of body acceleration. Each categorized behavior was assumed to correspond to when the bird was on land, in flight, on the sea surface, diving and so on. The behaviors classified by the procedures accorded well with those independently defined from depth profiles. Because our approach is performed by unsupervised computation of the data, it has the potential to detect previously unknown types of behavior and unknown sequences of some behaviors.

With the Arctic rapidly changing, the needs to observe, understand, and model the changes are essential. To support these needs, an annual cycle of observations of atmospheric properties, processes, and interactions were made while drifting with the sea ice across the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition from October 2019 to September 2020. An international team designed and implemented the comprehensive program to document and characterize all aspects of the Arctic atmospheric system in unprecedented detail, using a variety of approaches, and across multiple scales. These measurements were coordinated with other observational teams to explore cross-cutting and coupled interactions with the Arctic Ocean, sea ice, and ecosystem through a variety of physical and biogeochemical processes. This overview outlines the breadth and complexity of the atmospheric research program, which was organized into 4 subgroups: atmospheric state, clouds and precipitation, gases and aerosols, and energy budgets. Atmospheric variability over the annual cycle revealed important influences from a persistent large-scale winter circulation pattern, leading to some storms with pressure and winds that were outside the interquartile range of past conditions suggested by long-term reanalysis. Similarly, the MOSAiC location was warmer and wetter in summer than the reanalysis climatology, in part due to its close proximity to the sea ice edge. The comprehensiveness of the observational program for characterizing and analyzing atmospheric phenomena is demonstrated via a winter case study examining air mass transitions and a summer case study examining vertical atmospheric evolution. Overall, the MOSAiC atmospheric program successfully met its objectives and was the most comprehensive atmospheric measurement program to date conducted over the Arctic sea ice. The obtained data will support a broad range of coupled-system scientific research and provide an important foundation for advancing multiscale modeling capabilities in the Arctic.

Snow on Antarctic sea ice plays a complex and highly variable role in air‐sea‐ice interaction processes and the Earth's climate system. Using data collected mostly during the past 10 years, this paper reviews the following topics: snow thickness and snow type and their geographical and seasonal variations; snow grain size, density, and salinity; frequency of occurrence of slush; thermal conductivity, snow surface temperature, and temperature gradients within snow; and the effect of snow thickness on albedo. Major findings include large regional and seasonal differences in snow properties and thicknesses; the consequences of thicker snow and thinner ice in the Antarctic relative to the Arctic (e.g., the importance of flooding and snow‐ice formation); the potential impact of increasing snowfall resulting from global climate change; lower observed values of snow thermal conductivity than those typically used in models; periodic large‐scale melt in winter; and the contrast in summer melt processes between the Arctic and the Antarctic. Both climate modeling and remote sensing would benefit by taking account of the differences between the two polar regions.

Invasive alien species are among the primary causes of biodiversity change globally, with the risks thereof broadly understood for most regions of the world. They are similarly thought to be among the most significant conservation threats to Antarctica, especially as climate change proceeds in the region. However, no comprehensive, continent-wide evaluation of the risks to Antarctica posed by such species has been undertaken. Here we do so by sampling, identifying, and mapping the vascular plant propagules carried by all categories of visitors to Antarctica during the International Polar Year's first season (2007-2008) and assessing propagule establishment likelihood based on their identity and origins and on spatial variation in Antarctica's climate. For an evaluation of the situation in 2100, we use modeled climates based on the Intergovernmental Panel on Climate Change's Special Report on Emissions Scenarios Scenario A1B [Nakićenović N, Swart R, eds (2000) Special Report on Emissions Scenarios: A Special Report of Working Group III of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, UK)]. Visitors carrying seeds average 9.5 seeds per person, although as vectors, scientists carry greater propagule loads than tourists. Annual tourist numbers (∼33,054) are higher than those of scientists (∼7,085), thus tempering these differences in propagule load. Alien species establishment is currently most likely for the Western Antarctic Peninsula. Recent founder populations of several alien species in this area corroborate these findings. With climate change, risks will grow in the Antarctic Peninsula, Ross Sea, and East Antarctic coastal regions. Our evidence-based assessment demonstrates which parts of Antarctica are at growing risk from alien species that may become invasive and provides the means to mitigate this threat now and into the future as the continent's climate changes.

Stable isotopes of water are important climatic tracers used to understand atmospheric moisture cycling and to reconstruct paleoclimate. The combined use of hydrogen and oxygen isotopes in water provides an additional parameter, deuterium excess ( d ), which might reflect ocean surface conditions in moisture source regions for precipitation. The d records from polar ice cores covering glacial‐interglacial cycles were used to reconstruct ocean surface temperatures at the moisture source, enabling elimination of source effects from the conventional isotope thermometer. However, observations of the essential relationship between d in vapor and ocean surface conditions are very limited. To date, theoretical values predicted using simple and atmospheric general circulation models (GCM) have not been validated against the data. Here, we show the isotope ratios of atmospheric water vapor near the ocean surface in middle and high latitudes of the Southern Ocean. Our observations show that d negatively correlates with relative humidity ( h ) above the ocean and correlates with sea surface temperature (SST). Despite the fact that the GCMs would underestimate the absolute value of observed d , the observations and simulation results are consistent for slopes between d versus h and d versus SST, suggesting that d is a reliable index to h and SST over the ocean surface.

Hayabusa2 at the asteroid Ryugu Asteroids fall to Earth in the form of meteorites, but these provide little information about their origins. The Japanese mission Hayabusa2 is designed to collect samples directly from the surface of an asteroid and return them to Earth for laboratory analysis. Three papers in this issue describe the Hayabusa2 team's study of the near-Earth carbonaceous asteroid 162173 Ryugu, at which the spacecraft arrived in June 2018 (see the Perspective by Wurm). Watanabe et al. measured the asteroid's mass, shape, and density, showing that it is a “rubble pile” of loose rocks, formed into a spinning-top shape during a prior period of rapid spin. They also identified suitable landing sites for sample collection. Kitazato et al. used near-infrared spectroscopy to find ubiquitous hydrated minerals on the surface and compared Ryugu with known types of carbonaceous meteorite. Sugita et al. describe Ryugu's geological features and surface colors and combined results from all three papers to constrain the asteroid's formation process. Ryugu probably formed by reaccumulation of rubble ejected by impact from a larger asteroid. These results provide necessary context to understand the samples collected by Hayabusa2, which are expected to arrive on Earth in December 2020. Science , this issue p. 268 , p. 272 , p. 252 ; see also p. 230

Abstract Cloud observations over the past decade from six Arctic atmospheric observatories are investigated to derive estimates of cloud occurrence fraction, vertical distribution, persistence in time, diurnal cycle, and boundary statistics. Each observatory has some combination of cloud lidar, radar, ceilometer, and/or interferometer for identifying and characterizing clouds. By optimally combining measurements from these instruments, it is found that annual cloud occurrence fractions are 58%–83% at the Arctic observatories. There is a clear annual cycle wherein clouds are least frequent in the winter and most frequent in the late summer and autumn. Only in Eureka, Nunavut, Canada, is the annual cycle shifted such that the annual minimum is in the spring with the maximum in the winter. Intersite monthly variability is typically within 10%–15% of the all-site average. Interannual variability at specific sites is less than 13% for any given month and, typically, is less than 3% for annual total cloud fractions. Low-level clouds are most persistent at the observatories. The median cloud persistence for all observatories is 3–5 h; however, 5% of cloud systems at far western Arctic sites are observed to occur for longer than 100 consecutive hours. Weak diurnal variability in cloudiness is observed at some sites, with a daily minimum in cloud occurrence near solar noon for those seasons for which the sun is above the horizon for at least part of the day.