The Greenland ice sheet (GrIS) stores the largest amount of freshwater in the Northern Hemisphere and has been recently losing mass at an increasing rate. An eddy‐permitting ocean general circulation model is forced with realistic estimates of freshwater flux from the GrIS. Two approaches are used to track the meltwater and its trajectory in the ocean. We show that freshwater from western and eastern GrIS have markedly different fates, on a decadal time scale. Freshwater from west Greenland predominantly accumulates in Baffin Bay before being exported south down the Labrador shelf. Meanwhile, GrIS freshwater entering the interior of the Labrador Sea, where deep convection occurs, comes predominantly (∼80%) from east Greenland. Therefore, hosing experiments, which generally assume a uniform freshwater flux spatially, will not capture the true hydrographic response and regional impacts. In addition, narrow boundary currents are important for freshwater transport and distribution, requiring simulations with eddy‐resolving resolution.

No abstract

No abstract

Abstract. Sea ice thickness evolution within the Canadian Arctic Archipelago (CAA) is of great interest to science, as well as local communities and their economy. In this study, based on the NEMO numerical framework including the LIM2 sea ice module, simulations at both 1/4 and 1/12 • horizontal resolution were conducted from 2002 to 2016. The model captures well the general spatial distribution of ice thickness in the CAA region, with very thick sea ice (∼ 4 m and thicker) in the northern CAA, thick sea ice (2.5 to 3 m) in the west-central Parry Channel and M'Clintock Channel, and thin (< 2 m) ice (in winter months) on the east side of CAA (e.g., eastern Parry Channel, Baffin Island coast) and in the channels in southern areas. Even though the configurations still have resolution limitations in resolving the exact observation sites, simulated ice thickness compares reasonably (seasonal cycle and amplitudes) with weekly Environment and Climate Change Canada (ECCC) New Ice Thickness Program data at first-year landfast ice sites except at the northern sites with high concentration of old ice. At 1/4 to 1/12 • scale, model resolution does not play a significant role in the sea ice simulation except to improve local dynamics because of better coastline representation. Sea ice growth is decomposed into thermodynamic and dynamic (including all non-thermodynamic processes in the model) contributions to study the ice thickness evolution. Relatively smaller thermodynamic contribution to ice growth between December and the following April is found in the thick and very thick ice regions, with larger contributions in the thin ice-covered region. No significant trend in winter maximum ice volume is found in the northern CAA and Baffin Bay while a decline (r 2 ≈ 0.6, p < 0.01) is simulated in Parry Channel region. The two main contributors (thermodynamic growth and lateral transport) have high interannual variabilities which largely balance each other, so that maximum ice volume can vary interannually by ±12 % in the northern CAA, ±15 % in Parry Channel, and ±9 % in Baffin Bay. Further quantitative evaluation is required.

A coupled 1-D sea ice-ocean physical-biogeochemical model was developed to investigate the processes governing ice algal and phytoplankton blooms in the seasonally ice-covered Arctic Ocean. The 1-D column is representative of one grid cell in 3-D model applications and provides a tool for parameterization development. The model was applied to Resolute Passage in the Canadian Arctic Archipelago and assessed with observations from a field campaign during spring of 2010. The factors considered to limit the growth of simulated ice algae and phytoplankton were light, nutrients, and in the case of ice algae, ice melt. In addition to the standard simulation, several model experiments were conducted to determine the sensitivity of the simulated ice algal bloom to parameterizations of light, mortality, and pre-bloom biomass. Model results indicated that: (1) ice algae limit subsequent pelagic productivity in the upper 10 m by depleting nutrients to limiting levels; (2) light availability and pre-bloom biomass determine the onset timing of the ice algal bloom; (3) the maximum biomass is relatively insensitive to the pre-bloom biomass, but is limited by nutrient availability; (4) a combination of linear and quadratic parameterizations of mortality rate is required to adequately simulate the evolution of the ice algal bloom; and (5) a sinking rate for large detritus greater than a threshold of ∼10 m d -1 effectively strips the surface waters of the limiting nutrient (silicate) after the ice algal bloom, supporting the development of a deep chlorophyll maximum.

Greenland ice sheet meltwater runoff has been increasing in recent decades, especially in the southwest and the northeast. To determine the impact of this accelerating meltwater flux on Baffin Bay, we examine eight numerical experiments using an ocean‐sea ice model: Nucleus for European Modelling of the Ocean. Enhanced runoff causes shoreward increasing sea surface height and strengthens the stratification in Baffin Bay. The changes in sea surface height reduces the southward transport through the Canadian Arctic Archipelago and strengthens the gyre circulation within Baffin Bay. The latter leads to further freshening of surface waters as it produces a larger northward surface freshwater transport across Davis Strait. Increasing the meltwater runoff leads to a warming and shallowing of the west Greenland Irminger water on the northwest Greenland shelf. These warmer waters can now more easily enter fjords on the Greenland coast and thus provide additional heat to accelerate the melting of marine‐terminating glaciers.

Labrador Sea Water (LSW) is one of the main contributors to the lower limb of the Atlantic Meridional Overturning Circulation. In this study, we explore the sensitivity of LSW formation to model resolution, Greenland melt, absence of high‐frequency atmospheric phenomena, and changes in precipitation. We use five numerical model simulations at both (1/4)° and (1/12)° resolutions. A kinematic subduction approach is used to obtain the LSW formation rate over the period 2004 to 2016. The control simulation, with (1/4)° resolution, showed a mean annual production rate of 1.9 Sv (1 Sv = 106 m3/s) in the density range of 27.68–27.80 kg/m3 for the period 2004–2016. Deep convection events that occurred during 2008, 2012, and 2014–2016 were captured. We found that with (1/4)° resolution the LSW formation rate is 19% larger compared with its counterpart at (1/12)° resolution. The presence of Greenland melt and an increase in the precipitation impact the denser LSW layer replenishment but do not decrease the overall LSW formation rate nor the maximum convection depth. A dramatic response was found when filtering the atmospheric forcing, which induced a decrease of 44% in heat loss over the Labrador Sea, strong enough to halt the deep convection and decrease the LSW formation rate by 89%. Even if our experiment was extreme, a decrease in the storms crossing the Labrador Sea with a consequent reduction in the winter heat loss might be a bigger threat to deep convection and LSW formation in the future than the expected increases in the freshwater input.

A set of numerical simulations (with horizontal resolutions of 1/4° and 1/12°) is conducted to study the Pacific Water pathway in the Arctic Ocean and the freshwater content in Beaufort Gyre. Passive tracer tags the Pacific Water entering through Bering Strait into the Arctic Ocean and further reveals its circulation routes and spatial distribution. Both the 1/4° and 1/12° simulations show Pacific Water mainly follows the Transpolar Drift over the integration period of 2002–2016, with a limited amount being able to flow eastward along the Alaskan coast to enter the Canadian Arctic Archipelago. However, the circulation pattern of Pacific Water within the Beaufort Gyre is quite different with a stronger and tighter anticyclonic circulation in the 1/12° simulation corresponding to the difference in freshwater content. The 1/12° simulation successfully reproduces the overall recent increasing trend in the freshwater content in the Beaufort Gyre, while the 1/4° simulation fails to maintain the high freshwater content state after 2007. Budget analysis suggests that this difference in Beaufort Gyre freshwater storage is mainly caused by lateral advection. The lateral freshwater flux is decomposed into two components due to the slow‐varying circulation and mesoscale eddies. The difference in the capability to resolve eddies in the two simulations causes the difference in the temporal evolution of both components of the lateral flux.