A Brief History
Over the centuries, sea ice observations have come in many forms, designed for different needs and covering different scales. The earliest observations began thousands of years ago to meet the needs of Arctic Indigenous peoples. By using sea ice as a platform for subsistence hunting, traveling over ice for long distances, and sharing experiences from one generation to the next, their knowledge of ice behavior, stability, and characteristics was formed. A rich and immeasurable understanding of sea ice continues to grow to this day, as sea ice is still a cultural livelihood for Arctic Indigenous communities (Krupnik et al., 2010; Huntington et al., 2017).
During the age of exploration, mariners navigated through the outer reaches of the Arctic ice pack, drafting its periphery onto coarse maps of lands and seas. The rise of routine, quantitative sea ice observations came with the explosion of whaling operations and expeditions. Iterative maps, reworked time and time again, began to formulate a clearer picture of the marginal ice zone and the seasonality of Arctic sea ice. The observational record expanded immensely during this time, so much so that these early observations of sea ice coverage were combined with modern-day satellite data to extend the record back to 1750 (Divine and Dick, 2006) and 1850 (Walsh et al., 2016). In recent decades, paleoclimate proxies from marine sediments, land ice cores, and coastal records have been used to reconstruct the sea ice record back even further to thousands of years (Polyak et al., 2010; Abram et al., 2013), revealing new insight on the natural variability of Arctic sea ice.
From the late nineteenth century and through the twentieth century, the Fram Expedition (1893–1896), drifting ice stations, autonomous drifting buoys, airborne campaigns, and submarine surveys bolstered the foundation of sea ice observations. A multinational, coordinated endeavor called the International Polar Year (IPY) began in 1882 with the aim of collecting geophysical observations of the polar regions year-round (Barr and Lüdecke, 2010). IPY was the driving force for many of these field activities. Collectively, they produced the first quantitative records of sea ice circulation (Thorndike and Colony, 1982) and sea ice properties (Nansen, 1902), including ridge size and distribution (Romanov, 1995), melt ponds (Nazintsev, 1964), ice draft (Gossett, 1996), and snow depth and density (Warren et al., 1999). Carefully pieced together, these data sets provided the first pan-Arctic, multidecadal, year-round time series; some have since served as baselines from which long-term changes in Arctic sea ice properties have been measured (e.g., Webster et al., 2014; Kwok, 2018).
In the mid-twentieth century, the age of remote sensing profoundly advanced our observational capabilities. Polar-orbiting satellites collect near-continuous data over vast spatial scales, readily filling in gaps in the pan-Arctic picture. These satellites also relay observations from buoys on drifting sea ice in real time. In the late 1970s, passive microwave remote sensing initiated the longest, most consistent observational record of sea ice on the pan-Arctic scale. This iconic record spans more than 43 years, monitoring sea ice day and night, in cloudy and clear skies, and has proven instrumental in unveiling the acute sensitivity of Arctic sea ice to global warming (Fox-Kemper et al., 2021). The passive microwave record continues to serve as a valuable metric against which to test the ability of climate models to accurately simulate the Earth system (Notz and SIMIP Community, 2020).
In recent decades, technology has honed our remote-sensing capabilities. Airborne campaigns have become testbeds of lidar altimetry (Ketchum, 1971), synthetic aperture radar (Holmes et al., 1984), and electromagnetic soundings (Kovacs et al., 1987). Sea ice properties previously unresolvable from air and space now include ridges (Fredensborg Hansen et al., 2021), melt ponds (Wright and Polashenski, 2018; Farrell et al., 2020), freeboard and thickness (Ricker et al., 2017; Kwok et al., 2019), sea ice age (Tschudi et al., 2020), snow depth (Rostosky et al., 2018; Kwok et al., 2020), and leads (Reiser et al., 2020). Conjointly, the temporal resolution of satellite data is ever refining, enabling the collection of process-oriented information such as divergence, convergence, shear, and vorticity of the ice pack along with melt pond evolution and more. Radio triangulation was used early on to locate drifting buoys to within ~25 km. With the advent of satellites, Doppler positioning of buoys increased location accuracy to ~300 m, then ~2–3 m with the Global Positioning System, and centimeter-scale with the newer geodetic-quality global navigation satellite system (GNSS) buoys. The Iridium constellation of satellites now allows continuous temporal coverage by drifting buoys.
This brings us to today, a time of plentiful Arctic sea ice observations. In the sections that follow, we demonstrate current observational capabilities across local, synoptic, and pan-Arctic scales using a wide range of in situ and remote-sensing tools. We showcase sea ice conditions in 2018, with an emphasis on the Beaufort and Chukchi Seas, where significant summertime ice loss has occurred in recent decades (Frey et al., 2015). We present complementary data sets collected from field and airborne campaigns, buoy deployments, and Lagrangian tracking by satellites. The results reveal that, through coordinated efforts across communities, a richer understanding of Arctic sea ice system behavior is readily within reach.
The Issue of Scale
Our capability to observe Arctic sea ice is better than ever before, but many knowledge gaps in the Arctic sea ice system remain—and filling in these gaps is a balancing act. There are pressing needs to observe Arctic sea ice across a multitude of spatial and temporal scales, from improving the development of retrieval algorithms and data assimilation in weather and sea ice forecasts to monitoring climate-scale changes. Fundamental to all measurements are their observational scales, the times and spaces at which the measurements can be made (Blöschl and Sivapalan, 1995). The observational scale can be a point measurement, such as a snow depth value obtained with a ruler, or a multi-kilometer footprint of a satellite sensor, as with passive microwave retrievals of sea ice concentration. Although the boundaries between scales are becoming increasingly blurred (Figure 1), the scales at which observations are made often differ from the scales at which processes occur (Blöschl and Sivapalan, 1995). For satellite retrievals, the spatial resolution of gridded products differs from that of a satellite’s footprint. Furthermore, what is measured within a satellite footprint may not statistically represent the average of the variable being derived. Collectively, these issues pose considerable challenges in the appropriate interpretation of in situ and remotely sensed observations.
FIGURE 1. Examples of sea ice processes that occur over multiple spatial and temporal scales. The scales at which measurements are made often do not match the scales at which phenomena occur. The green background shading indicates scales over which in situ measurements are made, while the purple shading represents scales at which remote sensing retrievals are possible. Adapted from Blöschl and Sivapalan (1995) for the Arctic sea ice environment. > High res figure
|
As a case example, consider the relationship between snow depth and sea ice thickness in the Arctic in mid-spring. On point measurement scales (<0.5 m footprint), snow depth and the thickness of smooth sea ice beneath it are often negatively correlated, with locally thin sea ice having locally thick snow, and vice-versa. This inverse relationship arises due to snow’s high insulating capacity, which inhibits heat flux and sea ice growth (Sturm et al., 2002) and the greater susceptibility of wind-blown snow to be deposited in surface depressions, such as on top of refrozen melt ponds (see Figure 11d in Perovich et al., 2003).
In contrast, airborne retrievals of snow depth and sea ice thickness (Kurtz et al., 2015) show the opposite result. The airborne retrievals, as 40 m averages, yield a positive correlation between snow depth and sea ice thickness, meaning that on thicker ice, there is deeper snow, and vice versa. Although the two sets of observations seemingly contradict one another, both are correct for the spatial scales over which they are measured. Their differences underscore the importance of using the process scale to guide the interpretation of different observational scales (Blöschl and Sivapalan, 1995). Point measurements are effective in capturing local interactions and isolating mechanisms, making them central to gaining a deep understanding of physical processes and improving their representation in models; however, their broad-scale application across regions is limited.
Larger-scale measurements, such as airborne and satellite retrievals, are an integration of numerous processes, with large-scale processes having a stronger influence than small-scale processes. This can be readily seen in the pan-Arctic distribution of snow on Arctic sea ice (Webster et al., 2018). The deepest snow is found in regions with higher snowfall rates, older sea ice, and rougher sea ice. Unlike point measurements, larger-scale observations are better suited for monitoring pan-Arctic distributions and assessing the ability of Earth system models to simulate the cumulative effect of different processes on variables.
A Sample in Time
Temporal sampling of Arctic sea ice has greatly improved over the decades. Historically, most in situ measurements were made during spring-summer when the polar day permitted airborne surveys and landings on sea ice, and the summer retreat of ice allowed the passage of ships. From these early observations, weather patterns were quantified and related to patterns in sea ice drift around the Beaufort Gyre and along the Transpolar Drift Stream (Figure 2a). The discovery of the persistent Beaufort High in sea level pressure resulted from these observations.
FIGURE 2. Maps of buoy drifts overlaid on satellite images of increasing resolution for June 2018. The yellow boxes indicate the areas zoomed-in from one figure panel to the next. (a) Blue dots show buoy positions for April 1–June 30, and the green and red lines show the drifts of buoy Clusters 1 and 2 during the same period, overlaid on a 1 km resolution MODIS true-color scene from June 22. The thick, light blue lines and arrows are schematics depicting the circulation of the Beaufort Gyre and the Transpolar Drift Stream. (b) The Cluster 2 drift overlaid on a 10 m resolution Sentinel-2 scene from June 24. (c) The Cluster 2 drift overlaid onto a ~0.4 m resolution WorldView-3 scene from June 28. (d) A WorldView-3 scene from June 28, where, upon closer inspection, individual buoys can be identified. ©2018 DigitalGlobe NextView License. > High res figure
|
Deployment of autonomous buoys (Polashenski et al., 2011; Lei et al., 2016; Wang et al., 2016; Perovich et al., 2021) and moored observatories (Spreen et al., 2020) enabled year-round observations of surface air pressure and temperature, upper ocean temperature and salinity, and other geophysical parameters. The temporal resolution of buoy observations was previously limited to the frequency of satellites passing overhead, resulting in sub-daily sampling on a timescale of synoptic atmospheric and oceanic processes. Now, with continuous coverage of the North Pole by Iridium, observations are considerably more frequent and can reveal short-lived processes. For example, hourly observations can capture the fracturing, ridging, and rafting of sea ice (Figures 2 and 3d), as well as the inertial oscillations in sea ice motion (Figure 2b). Processes that quickly thicken the ice relative to slow thermodynamic growth can readily be monitored by analyzing the areal geometric changes between buoy arrays over nested spatial scales (e.g., the Distributed Network in Nicolaus et al., 2022). Such high-frequency observations provide new opportunities for studying atmosphere-ice-ocean interactions: the impact of tides on sea ice surface roughness, the role of tides and wind-driven circulation in the evolution and distribution of sea ice thickness, and the influence of net divergence or convergence on geographic differences in summer melt processes, to name a few.
FIGURE 3. The evolution of the total areal surface fraction for (a) Cluster 1 and (b) Cluster 2. (c) The melt pond fraction, given as the ponded ice fraction. (d) A subset of pond evolution from WorldView images of Cluster 1. Panels in (d) show pond formation on multiyear (top row) and first-year (bottom row) ice. Each panel is 400 m x 400 m. Dates are presented as “month-day.” ©2018 DigitalGlobe NextView License. > High res figure
|
The amount of data transmitted by buoys has also increased dramatically, from a few bytes to hundreds of bytes. In the near future, the Iridium Certus system will enable the transmission of thousands of bytes. Together with higher data transmission, current electronics and batteries permit the collection of hourly observations. Further advances in low-power electronics and significant improvements in the power density of batteries have spurred initiatives within the World Meteorological Organization and Intergovernmental Oceanographic Commission to deploy buoys that transmit observations every 10 minutes.
Spatial Heterogeneity
The heterogeneity of Arctic sea ice is at its most vibrant display during the melt season. From May to September, the vast expanse of the snow-covered ice pack transforms into a mosaic of floes composed of bare and ponded sea ice. Melt ponds decrease the surface albedo and increase the amount of absorbed and transmitted sunlight (Perovich et al., 2002a). Features of the sea ice cover, like melt ponds, are highly variable in their spatial distribution, and they often change with time-varying processes. The scales over which these features change dictate the type of sensors that can be used to detect them and the frequency at which they can be sampled. For example, fully resolving individual melt ponds requires observing sea ice at the meter scale or finer. To capture pond drainage events, daily temporal resolution is needed at a minimum.
To explore spatial heterogeneity in sea ice observations, consider the evolution of melt ponds at the locations of two drifting buoy clusters in summer 2018 (Figures 2 and 3). Clusters 1 and 2 were tracked with high-resolution optical satellite (WorldView) imagery as they drifted across the Beaufort and Chukchi Seas. The images were analyzed using the open source sea ice processing algorithm (Wright and Polashenski, 2018) to detect the areal fraction of each image that falls into one of four surface categories: (1) snow and bare sea ice, (2) gray (thin or slushy) ice, (3) melt ponds or submerged ice, and (4) ocean.
Figure 3d shows a bird’s-eye view of the surface conditions. These subsets reveal how truly heterogeneous the Arctic sea ice cover can be on meter-to-kilometer scales. The smooth sea ice on the lower portion is representative of undeformed first-year sea ice, while the rougher, more variable surface in the upper portion is characteristic of multiyear sea ice. Both ice types coexist on the same floe, but the progression of melt ponding differs between the two. This arises from dissimilar surface roughness and ice permeability, which affects the areal coverage and persistence of melt ponds. Yet, even across floes of smooth first-year sea ice, there can be contrasting discrepancies in melt pond coverage, with widespread ponding on one floe and complete absence of ponds on another.
The different timing and distribution of melt ponds lead to locally variable rates of melt. If melt ponds form early in one location, that surface will undergo greater melt earlier in the season than unponded ice due to positive albedo feedback. If an area has larger melt ponds, its melt rate will be higher than that in areas with smaller melt ponds. At Cluster 1 (Figure 3d), melt pond formation began earlier on the multiyear ice portion. Despite this earlier formation, the melt ponding on the smooth, first-year ice was more extensive by late June in 2018, and the first-year ice soon deteriorated into a skeletal-like structure, full of holes where ponds once were. Shortly thereafter, the first-year portion of the floe broke up, notably weeks earlier than the multiyear portion. When considering sea ice mass balance, whether the timing of formation or areal coverage of melt ponds is more important remains an open question, but this may be better understood with future coordinated tracking of drifting buoys by satellites.
Beyond the kilometer scale is the mesoscale (~10 km–1,000 km), a scale at which local heterogeneity starts to diminish and geographic differences and predominant sea ice conditions become increasingly important. Continuing with our melt pond example, we turn to the drifting buoy clusters to examine the scaling behavior of spatial heterogeneity. The two clusters remained ~150 km apart (Figure 2) as they drifted with the winds and ocean currents in May-August in 2018, and hence, they were subjected to similar synoptic conditions: similar snowfall, rainfall, temperatures, and radiative forcing. Accordingly, the average snow depths in their broader areas were alike, 0.19 m vs 0.18 m for Clusters 1 and 2, respectively. Even so, the physical state of the sea ice differed between the two clusters. In mid-April 2018, NASA’s Operation IceBridge airborne mission flew over both sites, measuring considerably thicker ice at Cluster 1 (2.2 ± 1.7 m), while Cluster 2 had relatively thinner sea ice (1.4 ± 1.1 m). Cluster 1’s location was also slightly northward and farther from the ice edge, which, together with having thicker ice, may have contributed to its longer survival during the melt season.
While there are some differences in the evolution of surface conditions between the two clusters at the 15 km × 15 km scale—notably that there was more open water and thinner, smoother sea ice at Cluster 2—the overall evolution of melt ponds was not significantly different between the two. These results, consistent with previous findings, reveal an important relationship between sea ice surface heterogeneity and scale: the coverage of melt ponds is highly variable across small spatial scales (tens of meters), but becomes increasingly consistent across larger scales (tens of kilometers) (Perovich et al., 2002b; Wright and Polashenski, 2018; and others). This phenomenon is known as the “aggregate scale,” or the scale at which the observed property is statistically representative of the larger region, and has been found to be tens of kilometers during the melt season (Perovich et al., 2002b). In essence, the local scale heterogeneity (Figure 3d; ~400 m × 400 m) collectively combines into a single aggregate scale (~15 km × 15 km) whose result reveals itself in the time series of the ponded ice fraction (Figure 3c).
Synthesizing Disparate Observations
Drifting buoys, moorings, and field campaigns typically have higher frequency sampling than air- and spaceborne missions, and the lower sampling frequency of satellite products limits their ability to detect short-lived events. In Figure 4, we combine disparate pieces of information to highlight the value of high frequency sampling and the advantages of co-deploying complementary instruments, a common practice of collaborators in the International Arctic Buoy Programme (IABP; https://iabp.apl.uw.edu/). The data in Figure 4 come from several sources, including compiled data sets found in the Lagrangian ice parcel database (Horvath et al., 2021) and in ice mass balance (IMB; Perovich et al., 2021) and WArming and iRradiance Measurements (WARM; Hill et al., 2018) buoy data.
FIGURE 4. A compilation of satellite, model, and buoy data from Cluster 2 during 2018. (a) ARTIST Sea Ice (ASI) concentrations. (b) Air temperatures from the seasonal ice mass balance (SIMB), WArming and iRradiance Measurements (WARM) buoy, and derived by satellite (Atmospheric Infrared Sounder). Storm events are marked by cyan asterisks and satellite-derived dates of early onset (yellow) and continuous (red) melt by vertical bars. (c) Snow depth from the SIMB and Lagrangian snow-evolution model (SnowModel-LG). (d) Under-ice photos from the WARM buoy. (e) In-ice and ocean temperatures with ice drafts from the SIMB and simulated sea ice thickness from Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS). (f) Ocean temperature from the WARM buoy. (g) Ocean salinity from the WARM buoy. (h) Photosynthetically active radiation (PAR) from the WARM buoy. > High res figure
|
In Figure 4b, consider the effects of temporal resolution on observed temperature. The buoy data exhibit a distinct diurnal cycling of temperatures, while the satellite-derived temperatures show a more slowly varying behavior and no diurnal cycling whatsoever. Different temporal sampling and different measurement techniques can produce substantially different results for the same phenomenon. The date of melt onset is an especially insightful example of the large differences that can arise between in situ and satellite-derived measurements (Figure 4b). In mid-May, in situ temperatures exceeded 0°C and surface melt was observed. Meanwhile, the satellite retrievals show a melt onset date three weeks later (Figure 3d). Deriving air temperature by satellite requires several assumptions, including those about emissivity, a property that is strongly influenced by the surface state and atmospheric water content (Ulaby et al., 1981; Jackson et al., 2006). Thus, discrepancies between in situ and satellite data sets can arise if assumptions in satellite retrievals do not truly reflect the ever-changing surface and atmospheric conditions. A three- to four-week discrepancy in the timing of melt onset changes the interpretation of cause and effect when it comes to interactions among the atmosphere, sea ice, and ocean, especially regarding primary productivity and melt processes.
And timing is everything for sea ice algae. Ice algae are a fundamental component of the Arctic marine food web (Kohlbach et al., 2016; Wiedmann et al., 2020) and are highly sensitive to the thickness of snow and sea ice due to their adaptation to low irradiance levels (Leu et al., 2015). At Cluster 2, the timing and duration of snow melt (Figure 4c) coincided with a notable ice algal bloom (Victoria Hill, Old Dominion University, pers. comm., 2021), resulting in less photosynthetically active radiation (PAR) in the upper ocean (Figure 4h). As melt pond formation progressed (Figure 3), irradiance levels returned to near-normal levels. The ice algae bloom gradually diminished (Figure 4d), but it had left its mark: the long duration of low ocean PAR levels substantially disrupted phytoplankton productivity in the upper water column for the season (Victoria Hill, Old Dominion University, pers. comm., 2021).
Auxiliary information, such as ocean temperature and salinity (Figure 4f,g), provides added insight on potential nutrient loading and also on environmental conditions leading to abrupt sea ice changes. This case study illustrates how coordinated deployments enable single discipline and cross-discipline analyses, which can greatly aid system science investigations.
Scalability and its Impacts
One of the most powerful traits of Arctic sea ice is its high albedo, which results in reflection of ~40%–80% of solar radiation back into space, compared to only ~7% reflected by the open ocean (Perovich et al., 2011). Here, we estimate solar heat input into the ice-free and ice-covered Arctic Ocean to examine how small-scale heterogeneity manifests across spatial scales (Figure 5). We estimate solar heat input using a daily mean downwelling solar flux from ERA5 reanalysis (Hersbach et al., 2018), prescribed surface albedo as in Perovich et al. (2011), stages of melt and freeze onset dates derived from satellite data (Markus et al., 2009), and sea ice concentrations also derived from satellite data. This exercise was performed using three sea ice concentration retrievals, each covering a different spatial scale: WorldView (WV) ~15 km scene with ~0.4 m resolution (Figure 3), ARTIST Sea Ice (ASI) 6.25 km grid (Spreen et al., 2008), and Climate Data Record (CDR) 25 km grid (Meier et al., 2021).
FIGURE 5. (left) Total annual solar heat input into the ice and ocean in 2018 from Climate Data Record (CDR) sea ice concentrations. Tracks of the drifting buoy clusters are drawn in cyan, and the region of the Beaufort and Chukchi Seas is outlined by the light gray box. (right) Equivalent amounts of sea ice melt based on the mean value of total annual solar heat input, an ice density of 900 kg m3, latent heat of fusion of 0.335 MJ km–1, and different sea ice concentration retrievals. Going from small to large scales, the estimates are derived from the buoy clusters (Buoys WV; Buoys CDR; Buoys ASI), the Beaufort and Chukchi Sea region (B&C CDR; B&C ASI), and the Arctic, averaged over 67°N–90°N (Arctic CDR; Arctic ASI). The inlaid bar chart displays the equivalent amounts of ice melt (meters) at the end of 2018. In reality, a small amount of heat would be lost through turbulent fluxes and other processes, and less ice melt would occur. B&C = Beaufort and Chukchi Sea region. ASI = ARTIST Sea Ice. CDR = Climate Data Record. WV = WorldView. > High res figure
|
First, consider the cumulative heat input on the smallest scale (6.25 km to 25 km) at the buoy clusters, where coincident pixels from the three sea ice concentration retrievals were evaluated (Figure 5, right: black, orange, and cyan lines). The estimates from the buoy clusters use the same solar flux and the same prescribed albedo for sea ice and open water; the differences arise solely from the sea ice concentrations, which are within ~1%–2% of one another throughout the year except for July, when ASI and CDR show 25% less sea ice coverage than the WV retrieval. The 25% difference in July sea ice coverage results in ~70–90 MJ m–2 more heat going into the sea ice and open ocean. This is equivalent to ~0.2–0.3 m of additional sea ice melt in a region that had, on average, 1.7 m thick ice prior to melt onset.
There are several possibilities for the discrepancies in the solar heat input estimates. For one, the passive microwave retrievals (ASI, CDR) may erroneously detect melt ponds as open water. In the range of microwave electromagnetic radiation, melt ponds are sufficiently deep to mask the emissivity of the underlying sea ice and produce an emissivity similar to that of open ocean (Ulaby et al., 1981). Thus, the passive microwave retrievals may underestimate sea ice concentrations and overestimate solar absorption when melt ponds are present. Previous work (Rösel et al., 2012) suggests a low bias in passive microwave sea ice concentrations by as much as ~40% due to melt ponds. Approaches blending higher-resolution optical imagery with passive microwave may be a fruitful way forward in improving sea ice retrievals during the melt season, as the use of optical imagery alone is limited by clouds in summer.
Discrepancies in the estimated solar heat input could also stem from a sampling issue. Here, we examined grid cells of the satellite products at each buoy cluster, with the centers askew from one another. Moreover, the clusters are located in a markedly dynamic environment; the marginal ice zone is susceptible to large changes in ice coverage over relatively short periods, which can result in retrieval differences if satellites collect data at different times. To explore the differences in spatial resolution, we performed the same exercise using the ASI and CDR products on regional and pan-Arctic scales. The higher resolution product (ASI) yields more solar heat input (30 MJ m–2) for both the Beaufort and Chukchi Sea regions and Arctic-wide. This equates to ~0.1 m of additional sea ice melt, which is distributed fairly uniformly across the Arctic, across ice concentrations, and between multiyear and first-year ice regions.
While there are significant differences between passive microwave products and their input data (Ivanova et al., 2015; Meier and Stewart, 2019), the differences in the solar heat input estimates shown here are still noteworthy. Collectively, they point to the need to continually revisit and refine remotely sensed sea ice observations as technology advances and our ability to resolve surface processes improves. Much like our understanding of the Arctic sea ice system, the techniques used to document sea ice conditions are ever-evolving, and there is added value in modernizing past observations to retain and reestablish a baseline from which long-term changes can be measured.
Auspicious Directions
Arctic sea ice loss is expected to continue over the next decades, with most model projections showing the first ice-free (extent under 1 million square kilometers) Arctic summer to occur by 2050 (Notz and SIMIP Community, 2020). As sea ice loss continues and technologies advance, autonomous systems and remote sensing will play increasingly larger roles in sea ice observations. As shown in the drifting cluster examples, harmonization between different observational scales and techniques augments the science that can be achieved with observational investments. Instrument arrays with complementary sensors are of particular interest since they enable cross-disciplinary science. They can capture the concomitant changes in atmospheric, sea ice, biological, and oceanic conditions, allowing for better understanding of processes on local scales. Coordinated airborne measurements and Lagrangian tracking of drifting arrays by satellite are at the forefront of sea ice observations. These types of observations provide the context necessary for relating local-scale heterogeneity and temporal variability to aggregate-scale properties and evolution. In particular, such fine-scale observations enable the scalability of in situ measurements for relating to coarse-resolution satellite products and model output.
Although future field operations may encounter greater risks with a thinner, less stable Arctic sea ice cover, in situ observation will be equally, if not more, valuable. As noted in Gerland et al. (2019), there are only a few clusters of ice and ocean buoys across Arctic sea ice at any given time, and there are recurring gaps in the observing network during winter when buoys drift into the Transpolar Drift Stream and away from the Eurasian coast. These gaps increase the uncertainties in our analyses of weather, sea ice, and climate change (Inoue et al., 2009; Inoue, 2021), and more efforts are needed to reseed the Arctic Observing Network during winter. Similarly, routine observations of sea ice conditions from ships initiated by the “Ice Watch” Arctic Shipborne Sea Ice Standardization Tool (ASSIST) Data Network (https://icewatch.met.no/) would benefit from the increasing number of ships operating in the Arctic.
New remote-sensing technologies, such as dual-sensor altimetry, will also require considerable ground validation programs to account for the anticipated environmental changes to come. The shift to seasonal, more saline sea ice, more frequent freeze-thaw cycling, and more frequent rainfall (Fox-Kemper et al., 2021) will lead to a more complex vertical substrate of snow and sea ice, posing greater challenges for remote-sensing retrievals in the future. To sufficiently quantify uncertainties and biases for scientific reliability, ground-truthing observations should optimally cover a wide range of snow and sea ice conditions. Ideally, such ground-truthing activities can be done in coordination with other observational efforts, such as drifting array deployments and planned field programs, including ecosystem studies, to foster collaboration and inclusivity across the broader community.
Lastly, models can optimize observations by helping guide observational priorities. Sensitivity experiments and observational assessments are effective in pinpointing sources of deficiencies in models, such as specific processes and physics. Targeting observations for such processes can help advance incomplete representations of heterogeneity, variability, and atmosphere-ice-ocean-ecosystem interactions in models. Model inter-comparisons can elucidate regions and seasons where intra-model spread is large, and subsequently direct where and when observing systems are best placed. By bridging observational and modeling efforts, we can elevate the scientific, operational, and community return on investments, and better anticipate the cascading effects of the changing Arctic sea ice system.
Acknowledgments
We are especially grateful for the assistance of John Sonntag and Kyle Krabill from NASA’s Operation IceBridge team in coordinating overflights of field sites; James S. Hak from the United States Geological Survey and Mike Cloutier and Paul Morin from the University of Minnesota Polar Geospatial Center in acquiring co-located WorldView imagery; Cameron Planck from Cryosphere Innovation for his assistance in providing IMB data (Perovich et al., 2021); Don Perovich for his assistance in processing the IMB sonic range-finder data; and Sean Horvath for providing output of cyclone events, AIRS, SnowModel-LG, and PIOMAS from the Lagrangian ice parcel database (https://tc.copernicus.org/preprints/tc-2021-297/). WARM buoy data was obtained from the NSF Arctic Data Center. WorldView satellite imagery was provided with DigitalGlobe NextView License 2018 through the University of Minnesota Polar Geospatial Center. MODIS reflectance data (MYD09GA; Vermote, 2015) and Copernicus Sentinel-2 (ESA) reflectance data (10.5066/F76W992G) were acquired through the US Geological Survey. We thank two anonymous reviewers and Stephen Warren for their constructive reviews, which greatly improved the manuscript.
Funding
IR is funded by contributors to the US Interagency Arctic Buoy Program, which includes support from the US Coast Guard, Department of Energy, North Slope Borough, NASA (80NSSC21K1148), NOAA (NA15OAR4310154, and NA20OAR4320271), NSF (OPP-1951762), Office of Naval Research (N00014-20-1-2207), and USNIC. MW was funded by NASA’s Interdisciplinary Research in Earth Science program (80NSSC21K0264) and New Investigator Program in Earth Science (80NSSC20K0658). NW was funded by NASA’s Cryospheric Sciences Program (80HQTR21T0063).