<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Pacific‐Arctic Connections: Assessing Flow Through Bering Strait in Context With Dynamic Ocean Topography and Surface Stress</title></titleStmt>
			<publicationStmt>
				<publisher>JGR Oceans</publisher>
				<date>08/01/2024</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10610155</idno>
					<idno type="doi">10.1029/2024JC021132</idno>
					<title level='j'>Journal of Geophysical Research: Oceans</title>
<idno>2169-9275</idno>
<biblScope unit="volume">129</biblScope>
<biblScope unit="issue">8</biblScope>					

					<author>Annika Margevich</author><author>Mary‐Louise Timmermans</author><author>Seth Danielson</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[Bering Strait is the only ocean gateway connecting the Pacific and Arctic oceans. The ∼1Sv northward flow of Pacific water through the strait to the Arctic Ocean has been increasing by ∼0.01Sv/yr since 1990. Monthly dynamic ocean topography (DOT), wind, and sea‐ice data at Bering Strait are analyzed in context with the long‐term record of flow through the strait to investigate local drivers. Ocean transport is found to be proportional to the across‐strait slope in DOT, suggesting some component of the flow is in geostrophic balance. Along‐strait ocean surface stresses, which modulate the across‐strait DOT slope via Ekman transport, are analyzed in the presence of a seasonally varying ice cover. It is shown that northward interior ocean flow under sea ice in winter results in southward surface stresses, and westward Ekman transport that slows the geostrophic component of the northward ocean flow. As the number of open water days local to Bering Strait increase each year, we find no trend in the annual mean surface stress, that is, the loss of sea ice is not leading to increased northward wind stress input that would enhance northward ocean flow. This analysis is consistent with the theory that changes in both the atmosphere and ocean non‐local to Bering Strait are likely driving the increased transport from the Pacific into the Arctic via Bering Strait.</p>]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="1.">Introduction</head><p>Bering Strait lies between Alaska and Russia and connects the Northeastern Pacific Ocean to the Arctic Ocean. It is the only source of Pacific inflow to the Arctic Ocean and, therefore, is particularly influential to the downstream Chukchi Shelf and Canada Basin. The strait is &gt;50 m deep and &gt;85 km wide, supporting a northward throughflow of approximately 1 Sv = 10 6 m 3 /s <ref type="bibr">(Roach et al., 1995;</ref><ref type="bibr">Woodgate, 2018)</ref>. The Pacific inflow serves as a source of relatively fresh water to the western Arctic Ocean, accounting for approximately one-third of the freshwater input to the Arctic <ref type="bibr">(Aagaard &amp; Carmack, 1989;</ref><ref type="bibr">Serreze et al., 2006)</ref>. Bering Strait inflow also carries approximately 3-6 &#215; 10 20 J/yr (relative to freezing temperature) of heat in the spring and summer from the Pacific to the Arctic <ref type="bibr">(Shaw et al., 2009;</ref><ref type="bibr">Shimada et al., 2006;</ref><ref type="bibr">Woodgate et al., 2010)</ref>. This incoming warm and fresh water influences the larger Arctic system. For example, the freshwater influences the regional stratification, and the incoming heat has been shown to be important in initiating the sea-ice albedo feedback <ref type="bibr">(Lu et al., 2020;</ref><ref type="bibr">Serreze et al., 2016;</ref><ref type="bibr">Woodgate et al., 2010)</ref>. Over the last three decades, volume transport through Bering Strait has been increasing at &gt;0.01 Sv/yr <ref type="bibr">(Woodgate, 2018)</ref>. Here we explore drivers of the Bering Strait inflow such that we might better understand how future climate change will influence the Pacific-Arctic system. The forcing behind the Bering Strait throughflow is often split into two components: a pressure head component and a local-wind component. The pressure head component relates to higher sea surface heights (SSH) in the Pacific than in the Arctic <ref type="bibr">(Aagaard et al., 1981;</ref><ref type="bibr">Coachman &amp; Aagaard, 1966)</ref>. <ref type="bibr">Aagaard et al. (2006)</ref> argue that the pressure head is of steric origin. More recent studies have used a variety of in situ, model, remotely sensed, and reanalysis data to explore how the Pacific-Arctic pressure head operates on different temporal and spatial scales and how the large-scale forcing relates to flow through Bering Strait. <ref type="bibr">Danielson et al. (2014)</ref> propose that Arctic and sub-Arctic atmospheric pressure fields (e.g., the Siberian High and Aleutian Low) both contribute to Bering Strait flow variations through a combination of direct wind forcing and altered SSH gradients that arise from cross-shore Ekman transport in both the Bering Sea and the East Siberian Sea (ESS). Along-shore winds drive coastal Ekman convergences and divergences that contribute to regional SSH gradients and associated alongshore and through-strait flows <ref type="bibr">(Danielson et al., 2014)</ref>. <ref type="bibr">Peralta-Ferriz and Woodgate (2017)</ref> examine the ESS in greater detail, identifying a temporally shifting control of the Bering Strait flow that seasonally relocates between the ESS and the Bering Sea. <ref type="bibr">Nguyen et al. (2020)</ref> use an adjoint model to explore how local and remote atmospheric forcings relate to flow through Bering Strait. They further affirm that on monthly time scales, both the ESS and the Bering Sea are influential in setting variability in flow through Bering Strait. Similarly, <ref type="bibr">Zhang et al. (2020)</ref> use a model to show that winds over the ESS and the Bering Sea influence the SSH gradient between the two basins which ultimately influences the flow through Bering Strait. Another proposed theory for the pressure head posits that it is set up between the North Pacific and the North Atlantic where the lack of deep water formation in the Pacific allows higher sea surface elevations there than in the Atlantic <ref type="bibr">(Cessi, 2020)</ref>. The consensus, though, seems to be that the Bering Sea and the ESS are important centers of action with respect to the large-scale forcings behind the Pacific-Arctic pressure head that partially drives the Bering Strait throughflow.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>RESEARCH ARTICLE</head><p>The second component forcing Bering Strait throughflow is associated with the local winds. This component has been identified through the correlation between local winds and transport through Bering Strait. In the winter when volume transport is weaker, that is, &gt;0.4 Sv, there is a strong correlation between northward transport and the strong southward winds <ref type="bibr">(Woodgate, 2018)</ref>. This suggests that Ekman transport may play a role, where strong winds from the north reduce northward transport by increasing westward Ekman transport across the strait <ref type="bibr">(Woodgate, 2018)</ref>. In the summer, when flow through the strait is &gt;1.3 Sv, there is only a moderate correlation between northward transport and the weak northward surface winds. This seasonality in the correlation between transport and the local winds has led to the assumption that the pressure head is the primary mechanism driving flow through the strait in summer and local wind forcing is the dominant driver for the flow in winter (Woodgate, 2018). While it is likely that local winds influence Bering Strait throughflow in the winter, the presence of sea ice has not yet been considered when exploring these local surface dynamics.</p><p>Bering Strait is sea-ice covered from around late November until June, where the influence of this seasonal ice cover on the local surface dynamics contributing to flow through Bering Strait is largely unknown. The complicated effect of sea-ice cover in the presence of a sustained ocean flow has been explored, for example, in the context of the anticyclonic wind-driven Beaufort Gyre (e.g., <ref type="bibr">Meneghello et al., 2018)</ref>. When the ice pack is less concentrated and more mobile (summer), anticyclonic wind forcing drives convergent Ekman transport, spinning up the gyre. Conversely, in the presence of a concentrated, relatively immobile ice pack (winter), friction between the anticyclonic geostrophic ocean flow and the slower sea ice drives a divergent Ekman transport, relaxing the gyre <ref type="bibr">(Dewey et al., 2018;</ref><ref type="bibr">Meneghello et al., 2018;</ref><ref type="bibr">Proshutinsky et al., 2002)</ref>. Analogous seasonal dynamics may have some influence on northward flow through Bering Strait.</p><p>In this study, we focus on the Bering Strait region and analyze dynamic ocean topography (DOT) in addition to mooring data, wind reanalysis and sea ice data to assess the Bering Strait flow dynamics. We explore how the seasonal presence of sea ice influences the surface stress in Bering Strait and how the loss of sea ice might be affecting ocean transport. In Section 2, the data are described, and the methods for characterizing the across-and along-strait DOT slopes and the surface stresses are explained. Section 3 outlines the equations of motion in the across-and along-strait directions and assesses the across-strait balance using observations. Section 4 discusses how the presence of sea ice relates to the correlation between local winds and northward transport and assesses whether the loss of sea ice in Bering Strait relates to the observed increase in northward transport. Finally, Section 5 summarizes and discusses these results.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.">Data and Methods</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.1.">Dynamic Ocean Topography</head><p>The monthly DOT data are from the Centre for Polar Observation and Modeling (CPOM), which uses satellite altimetry data from the Envisat and CryoSat-2 missions <ref type="bibr">(Armitage et al., 2016</ref><ref type="bibr">(Armitage et al., , 2017) )</ref> (Figure <ref type="figure">1</ref>). Satellite altimetry measures the distance from the satellite to the surface of the ocean. DOT is calculated by subtracting the altimetry reading from the distance between the satellite and an assumed reference geoid (the GOCO03s geoid in this case). In the presence of sea ice, DOT can be estimated by using leads and other open-water areas to determine the distance between the satellite and the ocean surface rather than the ice surface. The CPOM DOT data span 2003-2014 at monthly time intervals. The spatial coverage is 60&#176;N-81.5&#176;N on a 0.75&#176;&#215; 0.25&#176;grid. The "pole hole" is defined as north of 81.5&#176;N to ensure compatibility between the two satellite missions. Average uncertainties associated with the sea-surface height measurements used to calculate DOT are &gt;1.1 cm for the monthly means and &gt;0.9 cm for the annual means <ref type="bibr">(Armitage et al., 2016)</ref>.</p><p>In recent years, two additional DOT data sets have been released, both using the CryoSat-2 satellite mission and covering up to 88&#176;N monthly from 2011 through 2020 <ref type="bibr">(Doglioni et al., 2023;</ref><ref type="bibr">Lin et al., 2023)</ref>. Differences between these data include differing geoid models and grids. <ref type="bibr">Doglioni et al. (2023)</ref> use the Optimal Geoid for Modeling Ocean Circulation (OGMOC) geoid model and is on a 0.75&#176;&#215; 0.25&#176;grid while <ref type="bibr">Lin et al. (2023)</ref> use the GOCO03s geoid model and a polar stereographic projection. Each of these data sets has a different grid through Journal of Geophysical Research: Oceans The across-strait slope, &#951; x r , is the change in DOT between the Russian coast and the Alaskan coast at the narrowest point of Bering Strait. To compute this across-strait slope, we rotate the frame of reference such that &#951; x r is approximately perpendicular to the DOT contours through Bering Strait (Figure <ref type="figure">1b</ref>), which is a 45&#176;clockwise rotation. Variables with a superscript of r indicate the value is calculated in this rotated frame. It follows that the along-strait slope, &#951; y r , is perpendicular to &#951; x r . We rotate the wind, ice, and water velocities to be in the same frame of reference. From the DOT field, we calculate a monthly time series of &#951; x r using a linear regression over the three grid boxes in the across-strait direction at every monthly interval (Figure <ref type="figure">2a</ref>) from which we can determine the seasonality of the across-strait slope (Figure <ref type="figure">2b</ref>) and the monthly anomalies (Figure <ref type="figure">2c</ref>). To compute &#951; y r , we define a transect perpendicular to the across-strait slope. We identify the latitude along the transect at which the DOT change is steepest and take this to be the center of the along-strait DOT slope. We then compute &#951; y r from a linear regression over the points within &gt;1&#176;latitude north and south of the center of the slope over the along-strait transect.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.2.">Mooring Data</head><p>Ocean velocity data through Bering Strait are from year-round bottom tethered ocean moorings, specifically the A3 mooring located just north of Bering Strait (Figure <ref type="figure">1b</ref>), which is one of several long-term moorings in the Bering Strait region <ref type="bibr">(Woodgate, 2018;</ref><ref type="bibr">Woodgate et al., 2015;</ref><ref type="bibr">Woodgate &amp; Peralta-Ferriz, 2021)</ref>. Because the Diomede Islands split Bering Strait into a western and eastern channel, data from the A3 mooring is taken to approximate water properties and transport through Bering Strait as a whole. For this reason, it is referred to as the "climate site" <ref type="bibr">(Woodgate, 2018)</ref>. From 1990 through 2006, velocities in the lower layer at &gt;45 m depth were measured by an Aanderaa single point current meter (RCM) every hour, where the water depth at A3 is &gt;57 m  <ref type="bibr">(Woodgate, 2018)</ref>. Beginning in 2007, velocities throughout the water column are measured by an acoustic Doppler current profiler (ADCP) every 30 min. To maintain continuity between the RCM and ADCP data, we use velocity data from the second lowest-depth bin of the ADCP data as bottom velocity and interpolate it onto an hourly grid. As a note, the lowest-depth bin of the ADCP data is neglected due to recording erroneously low velocity values <ref type="bibr">(Woodgate, 2018)</ref>. We de-tide the hourly velocities using the t tide function in Matlab <ref type="bibr">(Pawlowicz et al., 2002)</ref> and take daily averages of the detided hourly data. Because we use a rotated frame of reference for the across-and along-strait slopes, we rotate the bottom velocities clockwise by 45&#176;for consistency (Figure <ref type="figure">3a</ref>); terms along-strait and across-strait refer to the coordinate axes drawn in Figure <ref type="figure">1b</ref>. Using the rotated bottom velocity and assuming flow is barotropic and horizontally homogeneous, we estimate volume transport from the A3 mooring by multiplying the rotated velocity by a cross-section area of 4.25 km 2 <ref type="bibr">(Woodgate, 2018)</ref>.</p><p>Sea ice motion data are from under-ice-tracking ADCPs mounted on the year-round bottom tethered ocean moorings at sites A3 and A4 (Figure <ref type="figure">1b</ref>). The A4 mooring has under-ice-tracking ADCP data from mid-2002 to the present, while the A3 mooring has ADCP data spanning mid-2007 to the present. Because the A3 site is considered the climate site, we use the under-ice-tracking ADCP data from A3 for 2007-2021. We then fill any gaps in the A3 time series over 2003-2021 with data from the under-ice-tracking ADCP on the A4 mooring, as there is relatively good agreement in sea ice motion between the two moorings (Figure <ref type="figure">S1</ref> in Supporting Information S1). Similar to the bottom velocity data, we de-tide hourly sea ice velocities using the t tide function in Matlab <ref type="bibr">(Pawlowicz et al., 2002)</ref> and take daily averages of the detided hourly data. We apply a threshold of 0.7 m/s for sea ice speeds, as values exceeding this threshold tend to be associated with short-lived spurious fluctuations. Additionally, we only consider under-ice-tracking ADCP data at times when the sea ice concentration is greater than 0.15 <ref type="bibr">(Cavalieri et al., 1991)</ref>. Smoothed daily time series for the resulting sea ice motion data are shown in Figure <ref type="figure">3c</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.3.">Wind and Sea-Ice Concentration Data</head><p>The 10 m wind field data are from the ERA5 Reanalysis hourly fields with a resolution of 0.25&#176;&#215; 0.25&#176; <ref type="bibr">(Hersbach et al., 2023)</ref>. Daily sea ice concentration data are obtained from the National Oceanic and Atmospheric Administration (NOAA) and National Snow and Ice Data Center (NSIDC) Climate Data Record of Passive Microwave Sea Ice Concentration, Version 4 <ref type="bibr">(Meier et al., 2021)</ref>. We regrid the sea ice concentration data onto the same 0.25&#176;&#215; 0.25&#176;grid as the 10 m ERA5 winds. We then compute the mean 10 m winds and sea ice concentration over 65&#176;N-66.75&#176;N and 189.5&#176;E-193.5&#176;E (dashed box in Figure <ref type="figure">1b</ref>). Smoothed daily time series for these variables are shown in Figure <ref type="figure">3b</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.4.">Quantifying Surface-Ocean Stresses</head><p>Bering Strait is ice-covered between approximately November and June each year. To account for the effect of seasonal sea ice cover in Bering Strait on the local surface stress, we characterize surface stress as</p><p>where SIC is the ice concentration local to Bering Strait; &#961; is the density of water, taken as &#961; = 1,027.5 kg/m 3 ; C ice is the ice-ocean drag coefficient, taken as 0.0055 <ref type="bibr">(Meneghello et al., 2017)</ref>; &#961; a is air density, taken as &#961; a = 1.25 kg/ m 3 ; C air is the air-ocean drag coefficient, taken as 0.00125 <ref type="bibr">(Meneghello et al., 2017)</ref>; u air is the wind velocity; and u rel is the relative velocity between the ice and the ocean, computed as u rel = u ice -(u geo + u ek ) <ref type="bibr">(Meneghello et al., 2017;</ref><ref type="bibr">Yang, 2006</ref><ref type="bibr">Yang, , 2009))</ref>. The characterization of relative velocity follows analyses of the Beaufort Gyre where the ocean velocity is separated into an Ekman and a geostrophic component <ref type="bibr">(Meneghello et al., 2017;</ref><ref type="bibr">Yang, 2006</ref><ref type="bibr">Yang, , 2009))</ref>. For Bering Strait, there are additional components to the flow, such as tides and other barotropic flows, in addition to the balanced (geostrophic) and Ekman flows. Here we refer to u geo as the "interior flow" to encompass all components of the ocean flow not associated with friction at the boundaries. Note that the DOT product includes ocean geostrophic velocity components at the same temporal and spatial resolutions <ref type="bibr">(Armitage et al., 2016</ref><ref type="bibr">(Armitage et al., , 2017))</ref>; however, due to the likely underestimation of these flows (inferred in Section 3.1), we do not use them in our calculations. To understand how this quantification of surface stress is sensitive to different interior flows, we examine three cases: 0 m/s, 0.1 m/s interior flow to the north, and interior flow equal to the bottom velocities (Figure <ref type="figure">3d</ref>), the latter of which is a reasonable assumption given that flow through Bering Strait is approximately barotropic below 10-15 m depth <ref type="bibr">(Roach et al., 1995)</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.">Dynamics of Flow Through Bering Strait</head><p>The monthly across-strait slope calculated from the DOT field indicates a correlation with the monthly alongstrait bottom velocity from 2003 to 2014 (R = 0.53, p = 1.34e-11; Figure <ref type="figure">2</ref>). This suggests that the volume transport through the strait may be related in part to changes in the DOT field across the strait. We explore this further by analyzing the dynamics that govern transport with the steady state, depth-integrated, horizontal components of the momentum balance (e.g., <ref type="bibr">Coachman &amp; Aagaard, 1966)</ref>:</p><p>where f is the Coriolis parameter ( f = 1.46 &#215; 10 -4 s -1 at the latitude of Bering Strait); U r = + 0 -H u r dz; V r = + 0 -H v r dz, and u r and v r are the across-and along-strait components of the ocean velocity, respectively; g is Journal of Geophysical Research: Oceans 10.1029/2024JC021132 gravitational acceleration; H is depth (approximately 50 m in Bering Strait). Lastly, &#964; sfc r is the rotated surface stress; &#964; bot r is the rotated bottom stress, given by &#964; r bot = &#961;C bot |u r bot |(u r bot ) where C bot is the ocean-bottom drag coefficient, taken as 0.0023 <ref type="bibr">(Couto et al., 2020)</ref>, and u r bot is the rotated bottom ocean velocity. In the next section, we will simplify Equation 2 and analyze how the simplified balance holds given observations local to Bering Strait.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.1.">Geostrophic Balance in the Across-Strait Direction</head><p>The across-strait balance is predominantly geostrophic, that is, the dominant balance is between the first two terms of Equation 2 (Figure <ref type="figure">4a</ref>). Using the DOT data and the mooring data, we assess the geostrophic balance in the across-strait direction by quantifying how the across-strait slope compares to observed monthly bottom velocities (rotated, such that they are in the along-strait direction) through the strait. Over 2003-2014, there is a Journal of Geophysical Research: Oceans 10.1029/2024JC021132</p><p>moderate correlation (R = 0.53, p = 1.3e-11) between the across-strait slope and along-strait bottom velocities (Figure <ref type="figure">4b</ref>). This moderate correlation in part reflects the similar seasonal cycles of the across-strait slope to that of the bottom velocities (Figure <ref type="figure">2b</ref>) where the across-strait slope is maximized in summer when bottom velocities are faster and the across-strait slope is lower in winter when the bottom velocities are slower. While ocean transport and &#951; x r have strong seasonalities, the relation between the across-strait slope and the along-strait velocities does not vary based on sea ice cover, suggesting that the geostrophically balanced component of flow in the strait does not vary seasonally.</p><p>The correlation between the across-strait slope and the bottom velocity at Bering Strait suggests that some component of flow is in geostrophic balance. However, an estimate of the geostrophic flow using the across-strait slope (i.e., the second term of Equation <ref type="formula">2</ref>) is only about 40% of the magnitude of the observed depth-integrated along-strait flow (Figure <ref type="figure">4b</ref>). This might be due to the coarse resolution of the DOT product in the vicinity of Bering Strait. There are three grid boxes across the strait, meaning there are likely unresolved variations in the DOT field across Bering Strait. Another reason for the underestimation is that there are additional components to the bottom velocity through Bering Strait beyond the geostrophic component, such as those associated with timevarying, barotropic flow, strait geometry, and friction, that we are unable to characterize and remove before assessing the dynamical balance.</p><p>The along-strait balance (Equation <ref type="formula">3</ref>) has also been used to explore the dynamics of flow through Bering Strait (e.g., <ref type="bibr">Couto et al., 2020;</ref><ref type="bibr">Woodgate, 2018)</ref>. We completed a similar analysis for the along-strait balance, examining both the geostrophic balance and a balance between the Coriolis term and the stress term. We did not find a strong relation in either case. Therefore, we focus on the across-strait slope and processes that may influence transport through modifying &#951; r</p><p>x . There are several mechanisms that could influence the magnitude of &#951; r x such as local surface stress, buoyancy gradients, and large scale pressure differences between the Pacific and Arctic Oceans. For the remainder of this paper, we will focus on understanding the influence of local sea ice conditions on &#951; r x and transport through Bering Strait via the local surface stress.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.2.">Ekman Transport in the Across-Strait Direction</head><p>We concluded in the previous section that some component of flow through Bering Strait is balanced by the across-strait slope, and it is well known that a locally applied along-strait surface stress can generate an acrossstrait slope via Ekman transport. Our goal is to explore how the local along-strait surface stress relates to the across-strait slope, and how this relation might vary based on local sea ice conditions. In the absence of sea ice, the wind stress corresponds to the across-strait slope (Figure <ref type="figure">5a</ref>) where the across-strait slope increases as the alongstrait winds increase to the north.</p><p>We find a statistically significant northward trend in local winds (99% confidence) over 2003-2021 (Figure <ref type="figure">3b</ref>). Following the scaling introduced by Woodgate (2018, their Section 5.2), we estimate that a northward wind trend of 0.45 &#177; 0.12 m/s/decade (Figure <ref type="figure">3b</ref>) would correspond to a bottom velocity trend of 2.8 &#215; 10 -4 &#177; 2.1 &#215; 10 -5 m/ s/decade, which is two orders of magnitude smaller than the observed trend in bottom velocity. This suggests that the observed increase in local winds to the north does not account for the observed northward increase in bottom velocity. Note also that local winds are generally to the south while transport is to the north (Figure <ref type="figure">3</ref>), so the wind effect is to modulate the northward ocean flow.</p><p>When estimating the trend in bottom velocity associated with increased northward winds, we characterized the surface stress as wind stress on open ocean. However, this is not the case in Bering Strait year-round as sea ice is present from roughly November to June. When sea ice is present, the surface stress can strongly depend on the strength of the interior flow, which can be of the same order of magnitude as the ice motion (Figures <ref type="figure">5b</ref>, <ref type="figure">5c</ref>, and <ref type="figure">5d</ref>); in the absence of sea ice, the ocean flow can be neglected relative to the much faster wind speeds in the calculation of surface stress. For sea ice concentrations close to one, in the case of no interior flow, surface stresses are generally northward (Figure <ref type="figure">5b</ref> and dotted black line in Figure <ref type="figure">5d</ref>). Since the winds over Bering Strait are generally to the south, even while sea ice motion is to the north, the northward ocean flow must be playing a role in driving sea ice to the north and any calculation of surface stress must assume non-zero ocean flow. As the interior ocean flow increases to the north (Figures 5c and 5d Using the same scaling from Woodgate (2018, their Section 5.2) as before, we find a statistically significant northward trend in surface stress assuming an interior flow of 0.1 m/s (90% confidence) over 2003-2021 (Figure <ref type="figure">3d</ref>). We estimate the northward trend in surface stress (assuming an interior flow of 0.1 m/s) of 0.008 &#177; 0.005 N/m 2 /decade (Figure <ref type="figure">3b</ref>) would correspond to a bottom velocity trend of 7.1 &#215; 10 -3 &#177; 4.5 &#215; 10 -3 N/ m 2 /decade, which is approximately 35% of the observed trend in bottom velocity. For surface stress calculated assuming the interior flow is equal to the bottom velocity, the surface stress trend accounts for &gt; 27% of the observed increase in bottom velocity. This suggests that changes in sea ice may contribute to the observed changes in bottom velocity, which we explore further in the next section.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="4.">Influence of Sea Ice on Transport via Surface Stress</head><p>Long-term northward trends in surface stress suggest that changes in surface stress may be contributing, in part, to the observed increase in northward transport through Bering Strait. We explore how the sensitivity of surfaceocean stress to sea-ice conditions might translate to influencing transport through Bering Strait in context with the loss of regional sea ice over recent decades. To do this, we first analyze the variability in along-strait winds, sea ice motion, bottom velocity, and across-strait Ekman transport, M x,r ek = &#964; y,r sfc / f , over the course of a seasonal cycle. We focus on 2006 as a representative year to show how the relation between variables differs under ice cover versus open water on daily to weekly timescales.</p><p>When sea ice is absent, we find that northward increases in wind lead to increases in northward flow through the strait via the inferred increased Ekman transport to the east (Figure <ref type="figure">6</ref>). This is consistent with the previous discussion of winds influencing variability in northward transport during days of open water. When sea ice is present in Bering Strait, we find that large increases in wind (to the north or south) are associated with increases in Journal of Geophysical Research: Oceans 10.1029/2024JC021132 sea ice motion, inferred Ekman transport, and bottom velocity in the respective directions. Additionally, the magnitude of the inferred Ekman transport in the presence of sea ice varies according to the assumed interior flow. When we assume no interior flow, the inferred Ekman transport under sea ice is stronger to the east than the cases of 0.1 m/s ocean flow to the north or taking ocean flow into account using the measured bottom velocity. Additionally, all cases illustrate variability that is consistent with a relationship to volume transport (i.e., inferred eastward Ekman transport and northward bottom velocities occur simultaneously). The greater inferred westward Ekman transport under sea ice than open water conditions in all cases could be one reason we observe weaker transport in the winter. We also find that all cases of interior flow show larger variability in inferred Ekman transport between sea ice-free and sea ice-covered conditions. However, we observe comparable variations in bottom velocity between when sea ice is present versus absent (Figure <ref type="figure">6</ref>). This could indicate that remote effects of local wind-forcing changes have a dominant influence on the flow. Given that surface stress is sensitive to the presence of sea ice, it follows that the loss of sea ice might lead to changes in transport through increased interaction with the winds via increased open water days.</p><p>The number of open water days in the Bering Strait region have been increasing over the last several decades <ref type="bibr">(Serreze et al., 2016)</ref>. We explore the possibility of whether the increase in the number of open water days has an effect on the transport through associated changes in surface stress imparted to the ocean. We use a SIC threshold of 15% over the region specified in Section 2.3 to define the day of sea ice retreat and advance locally each year <ref type="bibr">(Cavalieri et al., 1991)</ref>, where days of open water are those between sea ice retreat and sea ice advance. However, we find no trend in the annual mean along-strait surface stress assuming an interior flow of 0.1 m/s over icecovered days (Figure <ref type="figure">7a</ref>) or over open water days (Figure <ref type="figure">7b</ref>). We do find a decrease in the cumulative southward surface stress at the strait during the ice-covered period, due to the trend of fewer days with ice cover each Journal of Geophysical Research: Oceans 10.1029/2024JC021132 year (Figure <ref type="figure">7c</ref>). This is compensated by an increase in the southward surface stress at the strait during the ice-free period as a result of the increasing trend in open water days per year (Figure <ref type="figure">7d</ref>). In other words, the observed increase in Bering Strait throughflow is not associated with a reduction in southward surface stress due to the increasing number of open water days at the strait. A similar result is found for the case of assuming the interior flow is the bottom velocity, though the trends have higher p-values (Figure <ref type="figure">S3</ref> in Supporting Information S1).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="5.">Summary and Discussion</head><p>Bering Strait is the only source of Pacific inflow for the Arctic Ocean, and the volume of Pacific water entering the Arctic has been increasing at 0.01 Sv/yr over the last few decades <ref type="bibr">(Woodgate, 2018)</ref>. We find that some component of flow through Bering Strait is balanced with the across-strait DOT gradient year-round. This is consistent with the findings of <ref type="bibr">Cherniawsky et al. (2005)</ref> who examined only ice-free periods. In the framework of geostrophic balance, the across-strait slope balances some component of the along-strait velocity through Bering Strait, and several processes can contribute to an across-strait slope.</p><p>We focus on how local winds and sea ice conditions contribute to the along-strait surface stress and generate an across-strait slope via inferred Ekman transport. We find that as the along-strait winds increase to the north, the across-strait slope steepens due to increased Ekman transport to the east, and northward transport through the strait increases. Of importance is that local winds in Bering Strait are generally to the south such that increases in volume transport correspond to a weakening of winds to the south as opposed to strengthening of winds to the north. One added complexity in this system is that Bering Strait is ice covered from around November to June, when ocean surface stress is modulated by the presence of sea ice. We find that quantifying the surface stress local to Bering Strait under ice cover is challenging given that the surface stress characterization is sensitive to the interior ocean flow. When we assume no interior flow, the sea ice motion and surface stress are northward, while the winds are predominantly to the south. This implies that the northward ocean flow, at least in part, drives the Journal of Geophysical Research: Oceans 10.1029/2024JC021132 sea ice and, therefore, should be considered when calculating surface stress. For interior flows increasing to the north, the surface stress increases to the south (i.e., westward Ekman transport increases) and northward ocean transport decreases. Similar to how sea ice in the Beaufort Gyre contributes to gyre spin down through divergent Ekman transport, sea ice in Bering Strait weakens northward flow through enhanced westward Ekman transport. In both systems, the presence of sea ice contributes to the reduced strength of each during the winter. Furthermore, this suggests that the strong correlation between local winds and the observed bottom velocities in Bering Strait might be misleading in the presence of sea ice as both the ocean flow and local winds influence the sea ice motion and surface stress. It is possible that the signal in local winds correlating to bottom velocity is indicative of the relationship between changes in the large-scale atmospheric circulation and the pressure head component driving flow through Bering Strait <ref type="bibr">(Danielson et al., 2014;</ref><ref type="bibr">Nguyen et al., 2020;</ref><ref type="bibr">Zhang et al., 2020)</ref>.</p><p>The sensitivity of surface stress to the presence of sea ice suggests that with the increasing number of open water days in Bering Strait, there may be changes in surface stress that would influence transport. However, we do not find a trend in the annual mean surface stress. It would appear that changes in the duration of ice cover local to Bering Strait are not leading to changes in the stress input that affect transport through the strait via changes in Ekman transport. It is possible that these processes are contributing to changes in remote regions associated with the pressure head, thereby affecting the across-strait slope and the volume transport through Bering Strait.</p></div><note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_0"><p>21699291, 2024, 8, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2024JC021132 by University Of Alaska Fairbanks, Wiley Online Library on [19/05/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License</p></note>
		</body>
		</text>
</TEI>
