Sea ice thickness is affected by processes that fall into two categories – thermodynamic and dynamic. Thermodynamic processes serve to increase or decrease thickness through phase change; dynamic processes serve to redistribute that thickness both horizontally and vertically. Sea ice models account separately for both thermodynamic and dynamic processes in order to determine ice thickness and predict how it will respond in a changing climate (; ; ). These models are evaluated against sea ice thickness observations, but current state-of-the-art, basin-wide observations capture only overall ice thickness (; ) and are unable to distinguish between thermodynamic and dynamic processes, which are independently affected through different mechanisms in a changing climate. In order to properly predict how sea ice will respond in a changing climate, these independent processes must be individually evaluated within models, requiring independent observations of each.

A few studies have investigated the effects of thermodynamics and dynamics in sea ice observations on a basin-wide scale (; ). These papers calculate changes to sea ice volume or concentration due to dynamics and define the residual between these dynamic effects and overall changes as the sum of thermodynamic and any other effects. completed this analysis from a volume perspective. Volumetric sea ice advection is calculated at regional resolution using sea ice motion and thickness products, while thermodynamic effects are taken to be the residual between this advection and overall regional volume change from the sea ice thickness product. The authors go on to compare these results to model output and investigate trends in both the model and observed data. The results of the initial partitioning between dynamic and thermodynamic effects show that, on a regional spatial scale and in the mean between 2002–2019, thermodynamic effects are larger than dynamic effects, which vary in sign from region to region. However, results on a sub-regional scale or from the central Arctic region are lacking. investigated the effects of thermodynamics and dynamics on sea ice concentration over the entire annual cycle and in both polar regions. The authors estimated sea ice concentration advection using sea ice motion vectors determined from observed passive microwave brightness temperatures. The residual difference between this advection term and the overall changes in sea ice concentration is taken to be the sum of thermodynamic effects and ridging. The results show that dynamics play a significant role in maintaining the observed sea ice cover. In particular, they found mechanical redistribution to be an important sink for sea ice concentration in the central Arctic.

In addition to these observation-based studies, model studies offer useful context for basin-wide analysis. examined the sea ice volume budget across the models within the World Climate Research Programme (WCRP) Coupled Model Intercomparison Project Phase 6 (CMIP6; ). In the basin-wide multi-model mean taken from 1960–1989, basal thermodynamic growth is found to be the dominant source term within the volume budget, while dynamics serve as a volume sink at -30% of basal thermodynamic growth. investigated trends in thermodynamic volume growth as represented by the National Center for Atmospheric Research (NCAR) Community Earth System Model Large Ensemble (CESM-LE; ) in order to ascertain the impact of reduced fall ice thickness on thermodynamic growth. Indeed, CESM-LE results suggest that, in our current climate, thinner fall ice leads to higher thermodynamic growth rates throughout the winter.

The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC; ) presented an opportunity for partitioning thermodynamic and dynamic growth from a Lagrangian perspective along the Transpolar Drift over a full year from October 2019 to September 2020. The study by analyzed data from airborne electromagnetic (AEM) surveys and an ice mass balance buoy network to characterize the annual cycle of both dynamic and thermodynamic sea ice thickness contributions experienced by the ice pack surrounding the MOSAiC drift station. Thermodynamic growth was modeled using ice mass balance buoy temperature profiles and subtracted from overall ice growth captured by the airborne electromagnetic survey data in order to calculate dynamic sea ice effects as a residual. Overall, the dynamic contribution of 0.1 m out of the 1.1 m growth amounts to 10 %. Offering a potential window into basin-wide partitioning of thermodynamics versus dynamics, compared National Aeronautics and Space Administration (NASA) Ice, Cloud and Land Elevation Satellite (ICESat-2) data collected over the MOSAiC drift station to ice mass balance buoy thicknesses collected during the field experiment. They found that the mode of ICESat-2-derived sea ice thickness over this region represented level ice under the effects of thermodynamics only, while mean and median sea ice thickness included sporadic deformation events which increased sea ice thickness under the effects of sea ice dynamics. Comparing the mean and median observations against the mode observations, the authors conclude that dynamics accounted for 35.6 % of the mean sea ice thickness increase and 42.6 % of the median sea ice thickness increase over a region enclosed by a 50 km radius around the Polarstern research vessel.

Other studies used estimates of sea ice drift vectors to relate dynamics to sea ice thickness growth, again from a Lagrangian perspective. calculated shear and divergence terms averaged over a region north of the Canadian Arctic Archipelago using estimated ice motion vectors during winter from 2011 through 2015. These terms and a constant thermodynamic growth term were linearly regressed to overall sea ice thickness change from European Space Agency (ESA) CryoSat-2. This analysis showed that divergence and shear led to 42 % to 56 % of overall thickness change averaged across the region in question during those winters, with the remaining change due to thermodynamic effects. The sea ice deformation effects of a 2015 winter storm in the Transpolar Drift north of Svalbard were examined by , who analyzed AEM measurements of thickness before and after the storm during the Norwegian Young Sea ICE (N-ICE2015) expedition. By tracking individual features in the measured sea ice distribution, they were able to relate divergence and shear to changes in sea ice deformation. In multiplying the effects of this single storm by the climatological average of 10 to 20 storms per winter, the authors predict a 5 % to 10 % volume increases due to deformation in the region. The study by also took advantage of AEM thickness measurements and satellite synthetic aperture radar observations of an unusually large polynya north of Greenland in 2018 to determine a relationship between deformation and thickness changes. Over the 65 000 km2 polynya and over the 1 month of analysis, deformation of ice was found to account for an average of 50 % of the thickness increase and in some cases as much as 90 % of the thickness increase.

analyzed RADARSAT Geophysical Processor System (RGPS) sea ice motion-vector-derived Eulerian estimates of deformation from 1996 to 2000 over a much larger area, though independent of any sea ice thickness measurements. They report that seasonal ice experiences more deformation than multi-year ice, possibly due to its decreased thickness and strength. They also report a decrease in sea ice divergence as the winter growth season progresses, potentially via the same mechanism of increasing thickness and strength. Even without a link to sea ice thickness, the findings of this study allow for the extrapolation of more localized, short-term thickness effect results to a larger spatial scale.

Existing studies have yielded either regional and long-term or localized and short-term partitioning of thermodynamic and dynamic growth, but a large-scale and longer-term dataset at sub-regional resolution is lacking – especially from a Eulerian perspective for easier comparison to model outputs. In this study, we fill this knowledge gap by presenting a monthly, Arctic-basin-wide, and 25 km resolution Eulerian estimation of thermodynamic sea ice thickness growth and dynamic, advection, and deformation effects on wintertime sea ice thickness. A difficulty inherent to large-scale partitioning of thermodynamic and dynamic effect is the large-scale characterization of basal thermodynamic growth. The Stefan's Law Integrated Conducted Energy (SLICE) retrieval methodology allows for daily and basin-wide retrieval of wintertime thermodynamic sea ice growth rate using passive microwave brightness temperatures . Here, the retrieved thermodynamic growth rate is used in conjunction with overall sea ice thickness changes from the radar altimeter aboard the ESA CryoSat-2 satellite to estimate dynamic sea ice effects during the CryoSat-2 era beginning in 2010. With overall sea ice thickness growth provided by the Alfred Wegener Institute (AWI) CryoSat-2 and ESA Soil Moisture and Ocean Salinity (SMOS) combination sea ice thickness product (CS2SMOS; ) and thermodynamic growth provided by SLICE, Arctic-basin-wide sea ice thickness changes due to dynamics are calculated as the residual difference between the two. The effects of advection are also estimated using a sea ice motion vector dataset allowing for the calculation of deformation thickness as a residual of overall dynamic thickness effects and advection thickness effects.

This analysis requires a sea ice thickness dataset, a sea ice motion estimation dataset, and a thermodynamic sea ice growth retrieval.

Sea ice thickness is affected by thermodynamic processes and dynamic processes. Thermodynamic processes serve to change sea ice thickness through molecular phase change, and dynamic processes serve to change local sea ice thickness through the mechanical processes of advection and deformation (ridging or lead formation). An Eulerian governing equation for sea ice thickness sums thermodynamic and dynamic processes: 2∂H∂t=f(t,H,x)-∇⋅(uH), where H is plane slab sea ice thickness; t is time; f is a function of time, thickness, and position x describing thermodynamic sea ice thickness increase; and u is the ice motion vector. The second term on the right-hand side represents dynamic sea ice thickness processes.

We aimed to partition weekly basin-wide observations of overall changes to the sea ice thickness field, ∂H∂t, into its components thermodynamic growth, f(t,H,x), and dynamic effects, -∇⋅(uH). The result is basin-wide estimates of thermodynamic and dynamic process effects on sea ice thickness. At each weekly time step, ∂H∂t was calculated in centered-difference fashion as half the difference between the CS2SMOS sea ice thickness field from the previous time step and that from the following time step. We estimated dynamics, -∇⋅(uH), as the residual difference between total growth, ∂H∂t, and expected thermodynamic growth, f(t,H,x), retrieved using the SLICE methodology with the CS2SMOS sea ice thickness from the previous time step as the initial thickness. Due to the 95 % or greater sea ice concentration constraint on SLICE, only grid cells that meet this condition at a given time step per an AMSR-E/AMSR2 passive microwave sea ice concentration product (; ) are included in this analysis. The results shown are time averages over various time periods. Grid cells that do not contain 95 % sea ice concentration for over 40 % of the time steps within the averaging time period are discarded and not considered. Figure a shows the portion of total time steps that each grid cell shows over 95 % sea ice concentration. Also shown in Fig.  is a map of regions used in Sect. . Regions are defined similar to those shown in to aid comparison.

The dynamic term within Eq. () can be further decomposed to form 3∇⋅(uH)=(∇H)⋅u+H(∇⋅u), where the first term on the right represents changes to local sea ice thickness due to advection, i.e., the movement of ice transporting ice of a new thickness into a grid cell, and the second term on the right represents sea ice thickness changes due to deformation processes caused by divergence of the ice motion vector field. The deformation effect does not include advection and therefore can be considered Lagrangian dynamics – i.e., the dynamic effect as observed by a Lagrangian drifter. Our estimates of dynamic effects were partitioned into advection effects and residual effects using the sea ice motion vector product. At every weekly time step, the first term on the right-hand side of Eq. () is calculated as the mean of the three time steps centered on the current time step (in order to maintain temporal resolution with the centered-difference scheme used to calculated overall thickness change) using the CS2SMOS ice thickness field and sea ice motion vectors and is taken to be advection effect. The residual difference of the overall dynamic effect and this advection effect includes the effects of ice deformation and any other effects that are not accounted for in SLICE or the calculation of advection. These additional effects include lateral growth, snow ice formation, and any frazil or new ice growth that occurs above 95 % sea ice concentration and is not captured within SLICE. The uncertainty in this residual is a summation of the uncertainties in CS2SMOS, SLICE, and the advection calculation. This approach was taken over calculating the deformation effect using the motion vectors and advection effect as a residual because motion vector divergence was found to be significantly more noisy than the motion vector fields themselves. As motion vectors are not available within the Canadian Arctic Archipelago, advection effects and the residual effects are not available in this region.

In summary, the governing equation expressed in Eq. () is conceptually reconstructed as follows: 4∂(CS2SMOS)∂t=SLICE+dynamiceffect, and the dynamic effect is further decomposed into advection and residual effects with Eq. () conceptually reconstructed as follows: 5dynamic effect=(∇CS2SMOS)⋅ice motion vector+deformation effect.

The results are presented monthly rather than weekly to reduce noise. Weekly integrations of the terms within Eqs. () and () were summed to monthly temporal resolution. The dataset covers November through April, due to SLICE assumptions described in Sect.  and availability of CS2SMOS data, spanning the 10 winters beginning in the years 2010 and 2012–2020. The 2010 data begin on 15 November rather than 1 November along with the availability of CryoSat-2 data, and the winter beginning in 2011 is not included due to a gap between the availability of passive microwave data from the earlier AMSR-E and later AMSR2. At each weekly time step, i, the following steps were completed at each point on the 25 km EASE-Grid 2.0:

6∂(CS2SMOS)∂t=12⋅(CS2SMOSi+1-CS2SMOSi-1)

Figure  shows the wintertime mean total sea ice thickness growth, dynamic effects, and thermodynamic growth for Arctic sea ice thickness across the entire 10-year analysis period. Total growth and thermodynamic effects are always positive, though with varying magnitudes, while the magnitude and sign of dynamic effect vary across the Northern Hemisphere sea ice. As expected, thermodynamic thickness growth is highest in the seasonal ice zones, in some areas greater than 0.3 m per month, and inversely proportional to the climatological sea ice thickness, leading to less than 0.1 m per month of thermodynamic growth in much of the central Arctic. Dynamic effects decrease sea ice thickness over 63 % of the area exhibiting ice during the study period and increase sea ice thickness over the remaining area. An increase in sea ice thickness due to dynamics occurs off the Siberian coast in the Chukchi Sea, where the Beaufort Gyre tends to deposit advected ice, and similarly north of the Canadian Arctic Archipelago, where both the Beaufort Gyre and Transpolar Drift tend to deposit advected ice. The highest positive dynamic effects of greater than 0.1 m per month occur just north of the central Canadian Arctic Archipelago and in the Chukchi and East Greenland seas. A decrease in sea ice thickness due to dynamics, in some areas less than -0.2 m per month, occurs in the coastal regions of the Kara and Laptev seas, from where the Transpolar Drift tends to remove ice, and similarly in the coastal regions of the Beaufort Sea due to a similar effect of the Beaufort Gyre.

Figure  shows the wintertime mean advection effect and residual effects across the 10-year analysis period. Negative advection effect dominates the Arctic sea ice, covering 70 % of the study area. The exception primarily occurs where the Beaufort Gyre advects thick ice from north of the Canadian Arctic Archipelago to the Beaufort Sea. Here, advection effects greater than 0.2 m per month can be found. The most significant negative advection effects, less than -0.2 m per month, occur in the coastal Laptev Sea. Residual effects are negative over 48 % of the study area. The largest residual effects of greater than 0.1 m per month occur where the ice motion tends to deposit and form ice ridges, north of the central Canadian Arctic Archipelago in the Chukchi Sea and north of Greenland in the Transpolar Drift. The Barents and Kara seas are dominated by strong negative deformation effects, some areas with less than -0.3 m per month. Coupled with the positive advection effects in the Beaufort Sea are negative deformation effects in this region.

Figure  shows the uncertainty in the mean effects shown in Figs.  and . Per Eq. (), uncertainties for mean values across time are reduced through the averaging operation. When reduced in this manner, uncertainty in the thermodynamic effect is small. The dynamic effect uncertainty, which includes the effects of uncertainty in both thermodynamic and CS2SMOS overall thickness change, is greater than thermodynamic growth uncertainty and highest in the lower latitudes. The uncertainty in CS2SMOS is the larger source of uncertainty in the dynamics term because the differencing of thickness between two separate time steps means overall uncertainty from both time steps are included. This is the case even when the covariance between time steps is included. Advection uncertainty is affected by a similar mechanism, as it includes the effects of spatial derivatives in the CS2SMOS field in addition to the motion vector uncertainty. Finally, residual effects have the highest uncertainty, as they include the uncertainty of overall change, thermodynamic growth, and advection summed in quadrature. In some areas, uncertainty has a similar magnitude to the effects themselves.

To investigate how the budgetary terms interact on a regional scale and facilitate comparison with and , Fig.  depicts mean total monthly sea ice volume changes due to each effect over the regions shown in Fig.  and summed across the entire Arctic. Also included in this figure is the mean monthly volume contribution from areas with less than or equal to the 95 % sea ice concentration threshold required for inclusion in the other terms to illustrate the effect of this condition. Each budgetary term was summed across each region and each month before a mean across these months was taken. The ≤95% sea ice concentration contribution was calculated by a centered-difference scheme similar to Eq. () but with the factor of its sea ice concentration included. The uncertainty shown on Fig.  is a sum of uncertainties from all grid cells within each region for a given effect, which are then binned and summed by month using Eq. () and averaged over time using Eq. () to yield a monthly uncertainty for each effect and region. Applying Eq. () to the uncertainty summation across each region yielded unreasonably low uncertainty, so a basic sum was chosen.

Across the Arctic, thermodynamic growth accounts for 1936±45 km3 per month of volume growth. The dynamic effect across the Arctic is a volume sink of -583±359 km3 per month, which is -30% of thermodynamic growth. Advection accounts for -126±229 km3 per month of this dynamic effect, while residual effects are -457±397 km3 per month. Volume changes in areas that do not have more than 95 % ice concentration and therefore are not included in the other terms are 247±52 km3 per month or 13 % that of thermodynamic growth. In all regions, thermodynamic growth is the strongest effect. The central Arctic has the highest thermodynamic volume growth at 320±12 km3 per month and is the only region with a positive dynamic effect at 56±101 km3 per month. The regions with the strongest dynamic effects are the Laptev Sea, the Kara Sea, and Baffin Bay, all regions with significant polynyas, with -135±24, -142±22, and -162±19 km3 per month, respectively. The remaining regions all have negative dynamic effects, but they are nearer to zero. The central Arctic has the strongest advection and residual effects among the regions, with -106±77 km3 per month advection effect and 163±127 km3 per month residual effect, and it is the only region with a positive residual effect. The only regions with a positive advection effect are the Beaufort and East Greenland seas and Baffin Bay, with 46±34, 53±39, and 3±13 km3 per month, respectively. The remaining regions have negative advection effects that are closer to zero. The strongest negative residual effect is in Baffin Bay, with -165±20 km3 per month. In the Barents and East Greenland seas, volume contributions in grid cells with less than or equal to 95 % sea ice concentration are not insignificant. Other regions have similar contributions from these grid cells below the sea ice concentration threshold, but they make up a smaller portion of the total volume changes.

Figure  shows mean dynamics and mean residual effects in relative terms as ratios to thermodynamic growth. To calculate these metrics, the ratio of dynamics and residual to thermodynamics was calculated at each time step, and then a time mean of these ratios was taken across all time steps in this study. These plots look similar to those in Fig. , as the scaling quantity, thermodynamic growth, is fairly uniform across most of the Arctic. Nevertheless, the relative importance of dynamics and other residual effects to thermodynamic growth is an important result. Much of the Arctic shows a slightly negative impact of total dynamics relative to thermodynamic growth. The areas north of the Canadian Arctic Archipelago and Greenland show the highest relative impact of dynamics, with some grid cells showing dynamics with over twice the impact relative to thermodynamics. The coastal regions of the Kara and Laptev seas show significant negative impacts of dynamics, with magnitudes nearly equal to thermodynamics. The greatest relative importance of residual effects is also found in the areas north of the Canadian Arctic Archipelago and Greenland, though skewing more towards the central Arctic. In these regions, residual effects can be twice that of thermodynamics. The largest positive relative impact of residual effects is found between Svalbard and Greenland, where the Transpolar Drift causes ridging in thick ice that is not experiencing large thermodynamic growth.

Figure  shows the monthly mean overall dynamic effect, thermodynamic growth, advection effect, and residual effects across the 10 winters of data. The thermodynamic growth field, inversely proportional to sea ice thickness, remains consistently low across the central Arctic through most of the winter and higher in the perennial ice zones, though growth decreases as thickness increases with time in these areas. The residual and advection fields sum to an overall dynamic effect that increases with time in winter, dominated by the residual effect, potentially as ridging effects increase with overall ice thickness. The monthly residual effect fields depict a negative residual effect in the westward leg of the Beaufort Gyre that peaks in December, decreases in January, and is nearly absent in February and March. February and March do, however, depict residual effect maxima north of Svalbard and eastern Siberia. Positive advection of ice thickness by the Beaufort Gyre similarly peaks in the early winter, though the broad pattern of the advection field remains consistent throughout the winter.

A detailed investigation of interannual variability upon mean thermodynamic, dynamic, advection, and deformation effects is outside the scope of this work. However, Fig.  shows yearly averages for each of these effects during the study period.

The climatology of ice motion during the CryoSat-2 era as plotted in Fig.  suggests that the patterns of dynamic effect determined here are sound. The Beaufort Gyre and Transpolar Drift both transport ice towards the Canadian Arctic Archipelago, where high positive dynamic effects are found. The high residual effects in this region are likely explained by ridging. Where the Transpolar Drift originates, in the coastal region of the Laptev Sea, strong negative dynamic effects dominate. In these regions, lead formation likely explains the strongly negative residual effects. Between these regions and the Canadian Arctic Archipelago, a couplet of negative advection effect and positive residual effect characterizes a region where the Transpolar Drift tends to move thinner ice towards thicker ice, all the while experiencing ridging and other effects that increase thickness. Where the motion vectors show the Beaufort Gyre transporting ice westward from north of the Canadian Arctic Archipelago, a similar but reversed couplet of high positive advection effect and strong negative residual effect is found. The ice is transported from a region of climatologically thicker ice north of the Canadian Arctic Archipelago to a region of climatologically thinner ice in the Beaufort Sea, leading to a positive advection effect. In this same region, the Beaufort Gyre flow pattern is diverging and accelerating westward. This divergence leads to lead formation and a reduction in mean overall thickness, as reflected by negative residual effects in this region. This matches previous work suggesting high divergence and lead formation in this region (; ; ). As this ice is further advected westward around the Beaufort Gyre, the advection and residual effects return to near zero before reaching areas north of eastern Siberia, where ridging may explain positive dynamic and residual effects. The location and magnitude of these patterns vary year to year with sea ice flow patterns and likely atmospheric conditions, as demonstrated in Fig. .

Our results add spatial detail to previous studies and agree well with previous studies when summed to similar scales. When summed across the Arctic, our results show that dynamic effects are a volume sink with a magnitude that is 30 % of thermodynamic growth. This similar to the CMIP6 multi-model mean of 30 % as reported by . Volume changes that occurred in grid cells containing sea ice concentration beneath the 95 % threshold for inclusion in our results account for 13 % of total volume growth, potentially containing much of the frazil growth making up 19 % of total growth in the CMIP6 models . When our results are summed to the regional scale, they agree well with the previous regional-scale results from , though their results go back to 2002, whereas ours begin in 2018. For each of the six regions included in , mean thermodynamic growth from our results is within 50 km3 per month. Of these six regions, our results also agree on which have the strongest dynamic effects (Laptev and Kara seas), while the rest of the regions show dynamic effects closer to zero – though sometimes with differing signs, perhaps due to the differing time frames.

reported on the spatial and seasonal characteristics of Arctic sea ice deformation in the ice motion vector fields using high-resolution RGPS data from 1997–2000. Though the years in question do not overlap with the analysis period shown here, their results offer context for understanding dynamic sea ice effects. Their analysis showed divergence in the Beaufort Sea and convergence north of eastern Siberia, a pattern reflected here by the residual effects in these same regions. They show that the fraction of deformed ice in these regions decreases over the course of the growth season, a phenomenon also shown by a lessening of the negative deformation effects from November through March in Fig. . These points support our supposition that residual effects are comprised mostly of the effects of deformed ice on sea ice thickness, whether it be through ridging (positive residual effects) or lead formation (negative residual effects). The ice motion vectors used here are not suited for vector calculus calculations of deformation terms, but a future comparison between concurrent observations of vector deformation fields and dynamics would be fruitful for model improvement. We set a 95 % sea ice concentration threshold, at or below which a grid cell is not considered for partitioning of thermodynamic and dynamic changes. In Fig. , we report sea ice volume changes that occurred in grid cells with below this threshold to be 266 km3 per month, or 13 % that of thermodynamic growth. Volume changes due to changes in sea ice concentration between time steps will be at most 5 % but would be included in the residual effects. It should be noted that new ice formation in leads that have opened up due to divergence in the flow field will not be included in the thermodynamic effect; rather, the balance of the new ice thickness and the thickness of the remaining ice in the grid cell will be quantified as a negative residual effect. This type of effect has the potential to be greater than 5 %, as the leads may grow to arbitrary size and not be apparent in the passive microwave sea ice concentration product if new ice fills the lead prior to the next passive microwave brightness temperature observation. It is this initial closing of the lead with new ice that will not be captured as thermodynamic growth here.

The mean relative dynamic and residual effects over thermodynamic growth shown in Fig.  are useful for understanding the relative importance of these processes and are useful for an eventual comparison across time periods. The highest relative impact of dynamics is found north of the Canadian Arctic Archipelago, where the Beaufort Gyre and Transpolar Drift both deposit ice and thermodynamic growth is limited due to the high thickness of the ice. The relative residual plot tells a different story. With advection removed, this plot shows that a drifting observer would experience strong negative residual deformation effects relative to thermodynamic growth in the Beaufort Sea and increasingly strong positive residual deformation effects relative to thermodynamics in the Transpolar Drift and eastward leg of the Beaufort Gyre. In the areas north of the Canadian Arctic Archipelago, both dynamics and residual effects are equal to and greater than thermodynamic effects, consistent with , who report that 42 %–56 % of mean thickness change is due to deformation. In some areas, the effects of dynamics dominate thermodynamics by a factor of 3. These results confirm the findings of and that dynamics play an important role in shaping the climatological sea ice thickness patterns in the Arctic.

Within the methodology used here, the dynamic effects term represents Eulerian dynamics – i.e., a spatially stationary observer of sea ice thickness would observe changes due to thermodynamics and changes captured by the dynamics term, which includes advection. A Lagrangian observer, who is advecting as described by sea ice motion vector, would not experience changes due to this advection. The Lagrangian observer would only experience changes due to the residual in the framework here. In our interpretation, this residual effect is dominated by deformation effects. In this way, our Eulerian deformation term can be considered Lagrangian dynamics. Two Lagrangian studies of dynamics observed along the MOSAiC drift track offer useful context for validating our Eulerian results (; ). The comparison is necessarily between our monthly Eulerian deformation term – i.e., Lagrangian dynamics – from closest in space and time to the MOSAiC drift track and dynamics, as described within those studies.

A 10 % dynamic sea ice thickness growth relative to total growth along the MOSAiC drift was reported by , lower than the 48 % reported here and the 42.6 % reported by . A likely primary cause of this discrepancy is related to temporal resolution. AEM sea ice thickness distributions were analyzed by across the 50 km buoy network at the beginning of the growth season and at the end of the growth season, estimating a cumulative thermodynamic growth during the season using buoy thicknesses. Through a phenomenon acknowledged by the authors, this does not account for dynamics affecting thermodynamic growth throughout the growth season. If dynamics were to increase thickness, as our results and those reported by show did indeed occur along the MOSAiC drift track, thermodynamic growth rate would decrease, as shown in Eq. (). Without accounting for this effect, thermodynamic growth is overestimated and dynamically driven growth is underestimated, due to the latter being calculated as a residual. The higher temporal resolutions in this work and greatly improve – though do not eliminate – this issue. Indeed, report a cumulative thermodynamic growth of 1 m, whereas we report total thermodynamic growth of 0.72 m. This is significant, especially given measured growth between 14 October and 17 April, whereas our analysis period is from 1 November to 1 April (this time discrepancy alone may cause differences as well).

Another potential reason for discrepancy between our results and those of is the higher overall thickness growth measured by CryoSat-2 relative to the AEM measurements. We report a total mean growth along the MOSAiC drift track of 1.39 m relative to 1.1 m from the AEM survey. Indeed, CS2SMOS shows a mean sea ice thickness of 2.52 m on 1 April 2020 at the MOSAiC location relative to 2.2 m from the AEM surveys. This would manifest as an increase in dynamic effect in our analysis, as dynamics are calculated as a residual when thermodynamics are subtracted from total growth. That there are differences here is not surprising. The satellite measurements are gridded and taken from the nearest grid cell to MOSAiC, while the AEM surveys are centered on the MOSAiC buoy array. On the other hand, though the AEM has higher spatial resolution, the coverage over the 50 km buoy network is not complete. The satellite samples a larger area, although not centered exactly on the MOSAiC buoy array. That our results agree better with is not surprising, given that both studies have used satellite sea ice thickness as the primary dataset. Given that our study aims at a temporally and spatially larger scale while these studies are more focused on a single drift track, we can expect differences in results while using these more localized studies to provide context for our larger-scale study.

The calculations and analyses carried out here are all performed using satellite data on a 25 km EASE-Grid 2.0. Whereas sea ice processes can occur on much smaller scales, the results on a satellite scale are useful for deciphering patterns and trends on an Arctic-basin-wide basis. The CS2SMOS sea ice thickness dataset represents mean sea ice thickness within each grid cell. In actuality, thickness over the grid cell is defined by a distribution rather than a single value. However, without having observed local thickness distributions available at each time step, we have omitted thickness distributions and applied the SLICE retrieval using the mean thickness provided by CS2SMOS. It is likely that implementing thickness distributions would augment our results. Given the non-linear and inverse relationship between thickness and thermodynamic growth rate present in Eq. (), a distribution favoring ice which is thinner than the mean thickness over that which is thicker than the mean within a grid cell would increase thermodynamic growth and decrease dynamic effects (and vice versa). Though the snow–ice interface temperature is not expected to vary as greatly as sea ice thickness across a 25 km grid cell, passive microwave snow–ice interface temperature retrieval also represents the mean across each grid cell. The dynamic, thermodynamic, advection, and residual effects then necessarily represent mean effects over the grid cell area. While few areas within a grid cell will have experienced exactly the effects described by our results, the cell will have experienced these effects on the mean.

Moving through the methodology described in Sect.  and uncertainty calculations shown in Sect. , it is apparent that the uncertainties associated with the various input products stack in such a way that uncertainty increases moving from thermodynamic growth with the lowest uncertainty through the dynamic effect, advection effect, and residual effect with the highest uncertainty. The uncertainty in thermodynamic growth is a summation in quadrature of SLICE uncertainty and that of the input thickness from CS2SMOS. As the dynamic effect is a residual between thermodynamic growth and overall sea ice thickness change from CS2SMOS, uncertainty in this term is a summation in quadrature of uncertainties in these terms. The first term under the radical is uncertainty in overall thickness change rate from CS2SMOS, as calculated in Eq. (). The uncertainty in the advection calculation is an application of Eqs. () to (), with uncertainties in the spatial derivatives of CS2SMOS appearing similar to that of the time derivative in Eq. (). Residual effects, calculated as the difference between the dynamic effect and advection effect, are a summation in quadrature of the uncertainties in these terms, meaning the effects of uncertainty in SLICE, CS2SMOS, and CS2SMOS temporal and spatial derivatives and ice motion are all included.

When compared against buoy data and using the buoy thickness as the a priori initial thickness, report SLICE to have a thermodynamic growth mean bias of 4×10-4 m d-1 and a standard deviation bias of 2.2×10-3 m d-1. The assumption of 2 W m-2 of basal flux from liquid water to solid sea ice leads to additional uncertainty from SLICE. Assuming a density of 917 kg m-2 and a latent heat of fusion of 3.32×105 J kg-1, each 1 W m-2 of basal sensible heat flux from the liquid seawater to solid sea ice is equivalent to a sea ice thermodynamic growth rate of 2.84×10-4 m d-1. If the assumed 2 W m-2 of basal sensible heat flux were removed, sea ice growth would increase by 5.67×10-4 m d-1, and an increase from 2 to 10 W m-2 would decrease thermodynamic sea ice thickness growth by 2.27×10-3 m d-1. The SLICE thermodynamic growth retrieval also does not account for lateral melt and freeze processes or any new or frazil ice growth that occurs above 95 % sea ice concentration.

Because used buoy thickness as the initial a priori thickness when comparing against buoys, the uncertainty associated with initial input thickness was not accounted for. To account for this, we have added the second term under the radical in Eq. (), which accounts for CS2SMOS uncertainty as per Eq. (). The AWI CS2SMOS product merges sea ice thickness retrievals from CryoSat-2 and SMOS into a product that contains reduced uncertainties relative to each instrument's products independently . CryoSat-2 uncertainties are highest over thin ice, while SMOS uncertainties are highest over thick ice, creating the opportunity for synergy. CryoSat-2 uncertainties are made up of observational uncertainties or noise and systemic uncertainties or bias . Observational uncertainties are reduced through spatial averaging on the grid and optimal estimation methodology used to create the CS2SMOS product. While systemic uncertainties affect estimates of absolute thickness, the differencing of thickness between time steps removes them from the estimations of various thickness effects calculated in this work. SMOS uncertainties are caused by uncertainties in the input parameters to the energy budget used to estimate sea ice thickness and are especially high over multi-year ice (MYI), results from which are removed from the optimal interpolation. The AWI CS2SMOS product provides an uncertainty value that estimates observational uncertainties for each individual estimate at each time step and grid cell, which allows our calculations of thermodynamic growth, dynamic effect, advection effect, and deformation effect to also have associated uncertainties at each time step and grid cell.

A potential mechanism for error occurs in the relationship between lead frequency and the snow–ice interface retrieval results. Leads and areas of lower sea ice concentrations contain open seawater exposed at the surface. Seawater has significantly lower emissivity in the microwave band than sea ice, therefore reducing passive microwave brightness temperatures in these regions. To the extent that leads or open water cover a grid cell, these lower brightness temperatures would then artificially reduce the retrieved snow–ice interface temperature and cause erroneously large thermodynamic growth rates. Via Eq. (), erroneously high thermodynamic growth without a change to the CS2SMOS estimates leads to erroneously lower dynamic effects. This phenomenon is difficult to spot because negative dynamic effects are expected in regions with high lead frequency. We restrict our analysis to sea ice concentrations of greater than 95 %, as retrieved by established passive microwave methods. As such, the highest possible open water fraction within a grid cell is 5 %. Assuming that emissivity of a satellite field of view is a linear sum of scene type emissivities weighted by area fraction, the effect of 5 % open water by area on satellite-retrieved snow–ice interface temperature can be approximated. Using an approximate open water emissivity at 6.9 GHz of 0.56 and a sea ice emissivity of 0.98, the emissivity of a 95 % sea ice concentration is 0.959. This reduction in emissivity from 0.98 for a 100 % sea ice concentration scene equates to a 5.25 K reduction in brightness temperatures for a 250 K snow–ice interface temperature. Propagating this difference through the retrieval algorithm per leads to a reduction of the retrieved snow–ice interface temperature of 6 K. In a scenario with thin ice and a small temperature gradient across the ice, this difference could be significant.

Sea ice models, including those contained within global climate models, account for sea ice thickness and volume changes through separate thermodynamic and dynamic processes. These processes are affected by different mechanisms in a changing climate, meaning independent observations of each are essential for model validation. In this study, we present a monthly, Arctic-basin-wide, and 25 km resolution Eulerian estimation of thermodynamic, dynamic, advection, and residual effects on wintertime sea ice thickness from 2010–2021. By retrieving thermodynamic sea ice thickness growth via a simple model driven by passive-microwave-based snow–ice interface temperature observations and differencing this growth on a weekly basis from overall sea ice thickness growth from a satellite altimeter–passive microwave combination sea ice thickness product , we show new spatial detail in these effects with a spatial resolution beyond the regional studies available to date. Using a sea ice motion product , we also separated the overall dynamic effect into its Eulerian, independent component effects of advection and residual effects.

When summed to a basin-wide total, our results show that dynamic effects are a sea ice volume sink with a magnitude that is 30 % of thermodynamic growth, similar to the results of a recent model-based study . Regional totals are also in line with previous estimates . However, our sub-regional results show significant local deviations from these basin-wide and regional results. The highest impact of dynamic effect relative to thermodynamic effect is found north of the Canadian Arctic Archipelago, where dynamic effects account for twice and sometimes 3 times the thickness growth of thermodynamics. Similarly, residual effects are highest relative to thermodynamics in these regions as well, with residual effects more than doubling thermodynamics here and slightly farther north, near the North Pole. This is likely due to ridging in these regions.

Thermodynamic growth is lowest in the central Arctic, lower than 0.1 m per month, and highest in the seasonal ice zones, often greater than 0.3 m per month. The highest positive dynamic effects of greater than 0.1 m per month are found north of the Canadian Arctic Archipelago, where the Transpolar Drift and Beaufort Gyre deposit ice. Strong negative dynamic effects of less than -0.2 m per month are found where the Transpolar Drift originates. The residual and advection effect fields are dominated by couplets with opposite signs between the two. The Beaufort Sea is characterized by positive advection effects of 0.1 m per month and negative residual effects of similar magnitude, while most other regions are characterized by negative advection effects, sometimes as low as -0.2 m per month, and positive residual effects, often greater than 0.1 m per month. A seasonal cycle is also shown for all thickness effects, the most prominent feature of which is an increasing positive residual thickness effect and overall dynamic thickness effect as the winter season progresses. A potential mechanism for this is increasing ice thickness resisting lead formation and making more ice volume available for ridging.

Monthly results compare well with a recent study of the Lagrangian dynamic and thermodynamic effects on sea ice thickness along the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) drift track during the winter of 2019–2020. Where our data show Lagrangian dynamics accounting for 48 % of growth in the grid cells nearest the drifting study area during this time period, found similar results of 42.6 % over a similar spatial scale. This lends confidence to our larger spatial- and temporal-scale results.

Next steps for these data include further interrogation of trends and patterns. There may be a relation to atmospheric conditions or patterns like the Arctic Oscillation or trends related to the changing climate. An additional step will be the comparison of these results to those given by sea ice and global climate models.