1. Introduction

Since 1970, the global average temperatures have increased by 0.2˚C decade⁻¹, because of the emissions of greenhouse gases [1]. Also, recent studies show that the climate change in the Arctic regions, has been at least twice as fast as the global average [2]. The decline in Arctic sea ice is the most dramatic indicator of the Arctic warming. Dynamical and physical changes such as thermal and radiative forcing in both Arctic atmosphere and ocean impact on the sea ice and the atmospheric changes act as a driver for the sea ice changes [3]. Feedbacks due to the reduction of ice concentration, surface albedo and ice thickness are primarily responsible for the amplification of the temperature signal [4] [5] [6]. The rate at which changes are occurring in the Arctic makes it a place that is quickly changing and where the environment and its people are likely to be strongly impacted.

Wind plays a key role in the Arctic and sub-Arctic atmosphere and the water cycle. For instance, the wind changes in seashore areas are very important because of their impact on hydrodynamical aspects [7] [8] [9]. Any changes in weather patterns due to wind shifts in the Arctic would be important, both in regional and global climatic changes. Understanding the coastal wind changes is one of the key factors to expect the future environmental changes linked to greenhouse gas changes [10]. The recent studies related to the atmospheric circulation anomalies and trends during winter have shown that summer sea ice decreases over the east part of Europe and north of Alaska are associated with the recent strong anti-cyclonic circulation [11]. Also, the study of the composite mean surface wind stress and sea ice anomaly vectors during summer revealed that the yearly September sea-ice area decreasing is due to the anti-cyclonic circulation over the Arctic during the ice-loss years whereas the cyclonic circulation is linked to ice-gain years [12].

This study can provide a better understanding and denote an assessment for recent local wind trends in Hudson Bay, located in the Canadian Arctic and Sub-Arctic regions (Figure 1). Wind circulation is one of the major atmospheric elements which has important potential impacts on local weather and climate change in Hudson Bay which as an example may cause the increase of the frequency of blizzards. However, the effects of climate change on Hudson Bay winds are not well understood and need more research. The seasonal ice reduction in Hudson Bay affects the surface albedo and near-surface meteorological and physical structures such as air temperature, humidity and cloudiness.

Hudson Bay Climate

Hudson Bay (Figure 1) is the largest inland sea in the world extending between

51 and 63 degrees north latitude and has an area of total 637,000 square kilometres. It is bounded by Nunavut territory on the north and west, Manitoba and Ontario on the south and Quebec on the east. Hudson Bay is connected with the Atlantic Ocean via Hudson Strait from the northeast and via Foxe Chanel to the north with the Arctic Ocean. It is located in the Canadian Arctic and Sub-Arctic regions. Its climate depends largely on the water surface, and it is anomalously cold in comparison with the other regions with the same latitudes, linked to the influence of the seasonal variation of its ice cover even during some parts of summer [13] [14] [15]. The freezing occurs in late October and November and in January and February the Bay is completely covered with the ice pack [16] [17]. Ice melting begins in May and the ice rapidly decreases by June, when cloudiness and fog increase. From mid-August to late October, it is free of ice [16] [18] [19]. The southwest part of Hudson Bay is the region having ice cover extending into the summer due to ice advection by winds and counter-clockwise gyre ocean currents [20] [21].

2. Methodology

This study investigates climate change and the increase in wind speeds over Hudson Bay and its relation to sea ice decline. This study uses NARR (The North American Regional Reanalysis Model) data over the Hudson Bay region. The NARR model [22] uses a very high-resolution NCEP Eta Model (the grid resolution is approximately 0.3 × 0.3 degrees or 32 km at the lowest latitude) together with the Regional Data Assimilation System (RDAS) which significantly, assimilates the precipitation along with other variables. The improvements in the model and assimilation have resulted in a dataset with substantial improvements in the accuracy of temperature, winds and precipitation compared to the NCEP-DOE Global Reanalysis 2. The outputs include 8 times daily data at 29 vertical levels. The NCEP monthly Dataset [23] is available from 1948 and NARR data is available from 1979 for 8-times daily, for every 3 Universal Time Coordinate.

We investigate the available historical archived model dataset by analysing mean values to understand the general structure and anomaly patterns and then study shifts and departures from climatological mean values from 1981 to 2010 (normal period) as a base, following World Meteorological Organization (WMO) Policy, which suggests using the latest decade for the thirty year average. Studying the anomalies can describe the climate variability and make it possible to understand the changes and anomaly structures on the regional scale when the high-resolution NARR data is used. Besides, the anomalies give a reference frame which allows more meaningful comparisons between locations and more accurate calculations of trends. Also, this study uses the statistical analysis such as t-tests for temperature at 2 m, wind speed at 10m and surface albedo for mean values during winter anomaly period (1998-2015) with the values during climatology normal periodic time to examine for an existence of a significant difference between anomaly period and climatology normal periodic time for mentioned parameters over Hudson Bay. For this stage, at first the multipledimension NetCDF data set for the study period (December-February 1998 to 2015) from NCEP/NARR website (https://www.esrl.noaa.gov/psd/) have been downloaded. Then, the obtain Hudson Bay shapefile for the coastal boundary ArcGIS from Canadian census website (http://www12.statcan.gc.ca/census-recensement/2011/geo/bound-limit/bound-limit-2011-eng.cfm) have been extracted. In the next step the extract information for latitude and longitude Hudson Bay shapefile also extracted and then the information from both NetCDF data set and Hudson Bay shapefile for the coastal boundary ArcGIS from Canadian census website have been combined. In the last step the data have been exported to the excel files for farther analysis for different variables.

3. The Recent Warm Trend in the Hudson Bay Area

The time series of yearly average1000 mb and 850 mb air temperatures from1948 to 2015 shows the rapid increase from late 1998 in the Hudson Bay area (Latitude: 50˚N-65˚N and Longitude: 265˚E-285˚E) and a significant increase in average yearly datasets for 1000 mb and 850 mb air temperatures in the Hudson Bay region during 2010 (Figure2(a) and Figure2(b)). Figure2(c) shows a time series of yearly mean sea surface temperature over the area of study for the same time as Figure2(a) and Figure2(b). This Figuredisplays a dramatic change from 1998 over the study area. This maximum change corresponds to the lowest value in the North Atlantic Oscillation index (NAO) which is shown in Figure2(d).

Additionally, the results of some models such as NCEP/NCAR R1, NCEP CSFR, NCEP/DOE R2, ERA-Interim, NOAA ERSST, COBE, CRUTEM4, JRA-55 are presented in Figure 3. All models show the increase in 2 m air temperature from late 1998 in the Hudson Bay region. Also, the long-term daily maximum and minimum Hovmöller diagrams (from 1948-2015) show this rapid change from late 1998 (figures have been not presented here).

4. Yearly Historical Ice Coverage in Different Part of Hudson Bay

The yearly historical total accumulated Hudson Bay ice coverage during 1980-2015, using Canadian Ice Service data (http://iceweb1.cis.ec.gc.ca/IceGraph/page1.xhtml) for the weeks of 0101-1225 for all seasons from 1980 to 2015 for north (a), north-west (b), east (c) and central (d) parts of Hudson Bay are examined. This results show the negative yearly trends for all parts of Hudson Bay changing from largest to lowest values in the north ~ −3.25(+/−1.86)% per decade, north-west ~ −2.36(+/−1.2)% per decade, east ~ −1.99(+/−1.56)% per decade and central part ~ −1.57(+/−1.24)% per decade respectively from 1980 to 2015 (Figures are not presented here).

5. Winter NCEP/NARR Physical and Meteorological Data Analysis in Hudson Bay Area

This section presents results for winter time (DJF) using NCEP/NARR reanalysis data for meteorological parameters such as 1000 mb wind vector, 10 m wind speed, mean sea level pressure, zonal and meridional winds at 850 hPa, vertical profiles of mean and anomalies of zonal wind, polar jet stream and vertical velocity (omega) at 700 hPa mean and anomalies for the Hudson Bay region from December-February 1998 to 2015 (study period). Here, anomalies for departures from 1981-2010 (normal period) are analysed.

5.1. Wind Vector Pattern at 1000 mb

Figure 4(a) shows mean DJF wind vector at 1000 mb during study period and Figure 4(b) shows the wind vector anomaly map of departures from the normal period. Areas with cyclonic and anti-cyclonic circulations are shown with C and A. The mean near surface wind circulation over Hudson Bay is cyclonic. The

anomaly map shows the increasing of the anti-cyclonic circulation over Hudson Bay during winter. In the east boundary and north west the contrast causing to increase the horizontal wind shear, between cyclonic and anti-cyclonic anomalies can be seen.

5.2. Mean Sea Level Pressure (MSL)

MSL pressure (Figure 5(a)) shows the existence of two cyclonic trough lines (red dash line) near the western and eastern boundaries of Hudson Bay and the anomaly map (Figure 5(b)) shows that the eastern trough has been deepened (pressure decreasing by 20 Pa in the north east to around 65 Pa in the east part of James Bay) whereas the western trough has been weakening (nearly 40 Pa increasing). The decreases (trough intensification with cyclonic circulation development) for the eastern trough are mostly over the areas with the increase in vertical motion of the Hudson Bay (see later in Section 5.6). This Figurealso shows that horizontal pressure gradient (the change in pressure with distance which is proportional to acceleration in momentum equation) has been significantly increased (the area with the contrast between cyclonic and anti-cyclonic anomalies) over Hudson Bay eastern boundary and some parts in the north west. It is remarkable that 10 m wind analysis also shows the increase of the wind speed (due to momentum equation and the pressure gradient determining the wind velocity) mostly over eastern boundary and in north west area (see later in Section 5.3).

5.3. Wind Speed Anomalies at 10 m

Figure 6 shows wind speed anomalies at 10m during study period departure from normal period. The wind speed has been increased over the eastern

boundary at around 0.2 - 1.4 m∙s⁻¹ and in north west area at around 0.2 - 0.8 m∙s⁻¹. The maximum wind speed increase is around 1.4 m∙s⁻¹ in the eastern part in the Umiujaq area.

5.4. Zonal and Meridional Winds Anomalies at 850 hPa

Figure 7(a) shows the decreasing of zonal wind at level of 850 hPa in all parts of Hudson Bay region from 0.2 m∙s⁻¹ in the north and east whereas the maximum values are in the north-central and southern boundary around 0.5 m∙s⁻¹. The zonal wind decreasing promises to decrease the westerly and rather colder currents. Figure 7(b) shows remarkable increase in the meridional wind mostly in

the eastern (at around 0.3 m∙s⁻¹), south (at around 0.1 m∙s⁻¹) and west (around 0.1 m∙s⁻¹) boundaries which shows the increase of rather warmer south-north currents in the east and south parts of Hudson Bay. The central part does not show any meridional anomaly.

5.5. Vertical Profiles of Mean and Anomalies of Zonal Wind and Polar Jet Stream

Figure8(a) shows vertical profile of zonal wind in Hudson Bay area. The position of jet core in Hudson Bay region is around 250 mb. Figure8(b) shows the anomaly map for the same time which shows the weakening of the zonal polar jet at around 0.8 m∙s⁻¹ in the core mostly in 50˚N - 65˚N. Figure8(c) shows mean North America map projection for the wind vector polar jet position at 200 hPa during the study period. This Figureshows that the core of jet is below Hudson Bay with the maximum of ~45 m∙s⁻¹. Figure8(d) shows the wind vector jet stream mean at same level and Figure8(e) shows the anomaly map departure from normal period for the same time but in map projection over Hudson Bay region. Comparison of the last two later Figures shows that the wind direction has been changed from northwesterly (in mean map) to rather warmer south easterly flows during anomaly period. South easterly flows have been intensified in Hudson Bay with the maximum value ~1.05 m∙s⁻¹ in core over south west of James Bay area during the anomaly period. Therefore, rather colder westerlies have been weakening during this anomaly time.

5.6. Vertical Velocity at 700 hPa

Figure 9(a) shows mean omega at 700 mb during study period and Figure 9(b) shows omega anomaly map for the same time departure from normal period. The maximum increase in vertical motion is in the southern part of Hudson Bay in the James Bay at around −0.015 Pa∙s⁻¹. The area with cyclonic circulation

controls with negative omega anomalies and anti-cyclonic circulations are associated with positive anomalies of omega (vorticity map is not presented here). The recent climate change study [24] in Baffin Island as a part of the Canadian Arctic hasshown that the recent increase in 1000 mb air temperature can be linked to vertical motion increasing (negative omega) and likely local increase in instability.

5.7. 2 m Air Temperature Anomaly Map

Figure10 shows the 2 m air temperature anomaly map. This Figureshows the increase of 2 m air temperature almost in the Hudson Bay region. The increase in 2 m air temperature in the boundary region in the north, east and south boundaries is remarkable with the maximum value at around 3.5 K in the east and around 5.5 K in the northwest area.

5.8. Albedo at the Surface

Figure 11(a) shows the albedo mean during December 1998-February 2015. Figure 11(b) shows the anomaly map for the same time departure from normal period. The mean albedo map shows the albedo values are changing from lower values around 35 percent in some narrow parts in the east to the maximum 75 percent in the north-west in Hudson Bay. Figure 11(b) shows that the albedo values have been decreased mostly in the north-west, eastern, some parts in the southern, east-central and also James Bay region around 5% along with ~1 K increase in 2 m air temperature to the maximum value at around 15% - 18% along with ~4 K increase in 2 m air temperature.

6. T-Test and Linear Regression Analysis during Winter Anomaly Period over Hudson Bay

The t-tests for mean values of temperature at 2 m, wind speed at 10m and surface albedo during winter anomaly period (1998-2015) as the study period with

the normal period show a much smaller p-values rather than the significance level set at 0.05. So consequently, t-test shows a significant difference during anomaly period (1998-2015) rather than the mean climatology (1981-2010) for mentioned parameters. The linear regression analysis of the winter temperature anomaly at 2 m against albedo anomaly shows the meaningful relationship between the decreased albedo and the increasing in temperature anomaly (R² = 0.59, n = 699, two-tail t-test at p-value nearly zero) where the slope fits at (−0.11 * albedo anomaly + 1.087) in Kelvin degree * % (Figure 12(a)). Also, the linear regression analysis of the winter wind speed anomaly at 10m against air temperature anomaly at 2 m shows the meaningful relationship between the increased wind speed anomaly and the increasing in temperature anomaly (R² = 0.44, n = 699, two-tail t-test at p-value nearly zero) and the slope fits at (0.22 * temperature anomaly −0.17) in ms⁻¹ K * % (Figure 12(b)). Figure 12(c) shows the linear regression analysis of the anomaly winter accumulation total evaporation at surface against albedo anomaly which has the meaningful relationship between the decreased albedo and the increasing in winter accumulation total evaporation at surface (R² = 0.97, n = 699, two-tail t-test at p-value nearly zero) and the slope fits at (−0.0062 * albedo anomaly −6 * 10⁻⁵) in kg∙m⁻² * %. This shows the increased ice melt is increasing evaporation in a statistically significant way.

7. Summary and Conclusions

The physical and meteorological high-resolution data set from NARR (The North American Regional Reanalysis Model) and NCEP (National Center for Atmospheric Research) have been investigated to understand the recent wind changes along with the possible relation to other physical and meteorological factors such as surface albedo, sea ice decline, temperature and evaporation in Hudson Bay. Here, at the first step, the long-term daily maximum and minimum

Hovmöller diagrams (figures have been not presented here) and yearly mean surface temperatures and yearly mean sea surface temperature analysis from 1948 to 2015 have been analysed. These results clearly have denoted a rapid warmer time for the Hudson Bay area from late 1998. For the next step, parameters such as surface albedo, 1000 mb wind vectors, MSL pressure, zonal and meridional winds at 850 hPa, polar jet, vertical velocity at 700 hPa and 2 m air temperature, during wintertime from 1998 to 2015, for mean values and anomalies with departure from normal period have been examined. The results indicate that during the anomaly period, the near-surface air temperatures at 2 m have been increased in most parts of Hudson Bay with remarkable increases in the boundary regions in the north, east and south with the maximum values at around 4 K in the east and at around 6 K in the north-west area. Seasonal winter time analysis of the 10 m wind speed anomaly has shown the increases of wind speed mostly over south of Hudson Bay in James Bay area, east and north-west boundaries of Hudson Bay with the greatest values over the eastern part at around 1.4 m∙s⁻¹ in maximum value. The wind speed has been increased over the north-west area at around 0.2 - 0.8 m∙s⁻¹. The wind vector at 1000 mb analysis has been revealed anti-cyclonic anomaly pattern over the most parts of Hudson Bay. Recent studies of the atmospheric circulation anomalies and trends during winter have shown, a link between accelerating of the summer sea ice decreases over the eastern areas of Europe and north of Alaska associated with the recent strong anti-cyclonic circulation [11]. Also, the study of the composite mean surface wind stress and sea ice anomaly vectors during summer with the composites of ice-loss and ice-gain years have revealed the contribution of the yearly September sea ice area changes to the anti-cyclonic circulation over the Arctic during the ice-loss years whereas the cyclonic circulation is linked to the ice-gain years [12]. Also, wind vector anomalies at 1000 mb, have shown the contrast pattern between cyclonic and anti-cyclonic local wind circulations over Hudson Bay area mostly over east, north and north-west. This can be associated with the change in mean sea level pressure patterns and pressure gradient increase. The study results indicate that over the areas with the contrast structure between cyclonic and anti-cyclonic circulation would result in increasing the horizontal wind shear leads to speedup the near-surface wind. Likewise, the wind vector analysing at different pressure levels shows the change in wind direction mostly from northwesterly (zonal wind weakening) to south and easterlies. The weakening of the zonal winds can be linked to the reduction of the pole to mid-latitude temperature gradient causing the enhanced Arctic warming [25] [26]. The zonal wind speed at 850 hpa along with its vertical structure has been weakened versus the meridional wind intensification at 850mb mostly over the south and east during winter.

Upper atmospheric transport including jet stream like the polar jet wind vectors at 200 mb during anomaly time has revealed the changes in magnitude and position. During winter anomaly time the polar jet at 200 mb has been shifted mostly from rather mean colder northwesterly flows to the rather warmer south and easterly anomaly flows over Hudson Bay areas. These results agree with the recent studies related to rapid Arctic warming linked to more common high-amplitude (wavy) jet stream and smaller westerlies structure of weaker zonal upper-level winds, and a more meridional flow direction [27].

Surface albedo has also been decreased mostly in the north-west, eastern, some part in the southern, east-central part and James Bay region around 5% (~1 K increase in 2 m air temperature) to maximum value at around 15% - 18% (~4 K increase in 2 m air temperature) during winter anomaly period.

Additionally, the yearly historical ice coverage from 1980 to 2015 for north, north-west, east and central parts of Hudson Bay has shown a negative yearly trend for all parts of Hudson Bay changing from largest to lowest values in the north (~ −3.25(+/− 1.86)% per decade), north-west (~ −2.36(+/−1.2)% per decade), east (~ −1.99(+/−1.56)% per decade) and central part (~ −1.57(+/−1.24)% per decade) respectively.

The statistical analysis such as t-tests for the temperature at 2 m, wind speed at 10 m and surface albedo for mean values during winter anomaly period (1998-2015) with mean climatology values (1981-2010) has shown an existence of a significant difference between anomaly period and climatology periodic time for mentioned parameters over Hudson Bay. The linear regression analysis of the anomaly winter temperature at 2 m against albedo anomaly and between the anomaly wind speed at 10 m against anomaly 2 m air temperature represents a meaningful relationship between these variables. Also, a linear regression analysis of the anomaly accumulation of total evaporation at surface against albedo anomaly shows a meaningful relationship between the decreased albedo and the increasing in winter accumulation total evaporation at the surface. This shows the increased ice melt can increase the evaporation in a meaningful statistical way. This study has shown many changes in winter atmospheric and surface conditions significantly in the north, north-west and eastern boundaries which have been quite dramatic over the anomaly period over Hudson Bay.

Acknowledgements

Thanks are given to NOAA/ESRL PSD, Physical Science Division, Boulder Colorado web page through http://www.esrl.noaa.gov/psd/, Canadian ice service webpage through http://iceweb1.cis.ec.gc.ca/IceGraph/page1.xhtml and Canadian census website through http://www12.statcan.gc.ca/census-recensement/2011/geo/bound-limit/bound-limit-2011-eng.cfm. I would like to thank Prof. Bello and Prof. Higuchi for suggestion and comments. Great thanks to Prof. Peter Taylor for useful comments and editing the text.