1. Introduction
The Arctic sea ice extent has been decreasing under global warming. According to IPCC (2013), the annual mean Arctic sea ice extent decreased over the period 1979–2012 with a rate that was very likely in the range 0.45 to 0.51 million km2 decade−1, and very likely in the range 0.73 to 1.07 million km2 decade−1 for the summer sea ice minimum (perennial sea ice). The IPCC (2013) concluded with high confidence that the average decrease in decadal mean extent of Arctic sea ice has been most rapid in summer, and that the spatial extent has decreased in every season, and in every successive decade since 1979. Satellite observation records show substantial losses in Arctic sea ice, particularly at the time of the minimum extent, which occurs in September at the end of the annual melt season (Hartmann et al. 2013). It is very likely that the Arctic sea ice cover will continue to shrink and thin during the twenty-first century as global mean surface temperature rises (IPCC 2013). Decreased sea ice extent opens the ocean surface up for wind-wave generation, exposing a greater coastal/offshore area to impacts of ocean surface waves (e.g., coastal erosion, shoreline sediment transportation), among other opportunities (e.g., shipping) and impacts [see Stephenson et al. (2011) for examples]. Since open water distances are the controlling variable for wave heights in the Beaufort–Chukchi–Bering (BCB) Seas, larger waves and more common swells are to be expected under scenarios for reduced seasonal sea ice cover, and the increased wave activity could be one feedback mechanism that drives the Arctic system toward an ice-free summer (Thomson and Rogers 2014). All these indicate that scientific information about and assessment of the climate, trends, and variability of the Arctic Ocean surface waves are important for climate change projections and climate prediction on the global and regional scales, and for adaptation and migration of climate change in the Arctic region, among many other applications. However, such information scarcely exists. There is an emerging need to study surface waves in the emerging Arctic Ocean.
However, because of a lack of wave data, there exist few studies on long-term historical changes in ocean surface waves in the Arctic, especially in the BCB Seas. Recently, Francis et al. (2011) analyzed cross-validated satellite altimeter radar observations for the period 1993–2011 and showed a strong trend of increasing surface wave heights in the Arctic Ocean, with the highest trends of up to 0.03–0.04 m yr−1 near the northern Alaskan coast (usually located 100–200 km offshore).
As described in Swail et al. (2007), Environment Canada produced the first comprehensive high-resolution wind and wave reanalysis for this region, namely the MSC (Meteorological Service of Canada, Environment Canada) Beaufort Wind and Wave Reanalysis, initially for the period 1985–2005. This wind and wave reanalysis has been extended to span the period 1970–2013 and is being updated to the present. Since buoy wave observations are limited to a number of buoy sites, and satellite wave observations to the recent couple of decades, wave reanalysis data have been used to analyze historical trends in waves (e.g., Semedo et al. 2011; Wang and Swail 2001, 2002). This study aims to use the MSC Beaufort Wind and Wave Reanalysis data to characterize the climate and trends of wind and waves in the BCB Seas.
The remainder of the paper is arranged as follows. The data and methodology used in this study are described in sections 2 and 3, respectively. Historical changes in wind speed and waves (significant wave height and mean wave period) are presented in sections 4 and 5, respectively. A summary in section 6 completes this study.
2. Data
The wind and wave data analyzed in this study, including surface wind speed (
As detailed in Swail et al. (2007), the wind fields were analyzed using the Interactive Objective Kinematic Analysis (IOKA) system of Oceanweather Inc. (OWI) (Cox et al. 1995). In the wind analysis, all available measured wind data were brought into a background field and assimilated at 6-hourly intervals with the guidance and input of an experienced kinematic analyst working through a graphical interactive tool called the Wind WorkStation. The coastal stations measurements of wind speed were transformed to effective overwater exposure at 10-m elevation, using the station-dependent overwater/overland transformation ratios derived from station surface roughness and anemometer height. The background wind field is taken from the 10-m surface wind-gridded products of the NCEP–NCAR reanalysis (NRA; Kalnay et al. 1996), which covers the period 1948 to the present. The global QuikSCAT wind database (http://manati.star.nesdis.noaa.gov/products/QuikSCAT.php) was used to identify and correct systematic errors in the NRA winds within the domain of the ocean response models.
The wave fields were produced by using the above analyzed wind fields to drive the Oceanweather third-generation (OWI-3G) wave model, with boundary spectra of the basin model being supplied from the GROW (Global Reanalysis of Ocean Waves) hindcast (Swail et al. 2007). The physics of the OWI-3G model are described in Khandekar et al. (1994). The General Bathymetric Chart of the Oceans (GEBCO) provided water depths for most of the basinwide hindcast; depths for the fine grid as well as overlapping regions of the coarse grid were supplied from the Canadian Hydrographic Service (CHS) archive. The two datasets were found to be consistent, so little smoothing was required at their boundary. Grid boxes with greater than 50% ice concentration are considered as land with no wave generation and/or propagation. This simple treatment of ice was implemented in the model because this wave reanalysis was produced prior to much of the most recent developments in the wave–ice interactions. In the coarse BCB domain the ice edge was allowed to change on a weekly basis to represent changing ice conditions during transition periods. The ice concentration data sources include the Goddard Space Flight Center (GFSC) daily data for the BCB domain for November 1970 to December 2000, the Defense Meteorological Satellite Program (DMSP) daily data for the BCB domain for January 2001 to the present (2013), and the Canadian Ice Service (CIS) weekly high-resolution ice concentration dataset for Canadian waters for the period since January 1970.
Swail et al. (2007) validated the MSC Beaufort Wave Reanalysis with the wave measurements from the 12 MEDS (Marine Environmental Data Service) buoys in the Canadian Beaufort for the period 1981–86. They concluded that overall the wave reanalysis reflects the measured day-to-day and storm conditions, although both under- and overestimation are seen in the reanalysis. For 32 storm peaks, the reanalysis shows on average an underestimation of 22 cm in wave height and of 1.05 s in wave period (Swail et al. 2007). For events of greater than 2-m
3. Methodology
As shown in Fig. 1, during the period 1970–2013, the open water exists in the BCB domain only in July to September (with a few grid boxes of open water in October; not shown), and the area of wave data (open water) increases considerably from July to September (larger increases from July to August than from August to September). Thus, we carried out our analysis for July, August, and September, separately.
(a) The mean area of wave data (open water) in the indicated months over the indicated periods. (b) The standardized time series of the area of wave data in the Beaufort–Chukchi–Bering (BCB) Seas domain for July, August, and September. The three time series were standardized using the same mean and standard deviation, which were derived from pooling the three time series together. Thus, the differences between these time series are proportional to the differences in the area of open water between the three months.
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
Considering that the open water area in the BCB domain increased more notably between the last two decades than between the earlier decades (Fig. 1a) and has remarkable interannual variations (Fig. 1b), we characterize trends in three ways. First, we show the 1970–91 and 1992–2013 climatological mean fields of each variable analyzed for each of the three months, and the changes between these two 22-yr periods (Figs. 2–5). Second, we derive the regional mean time series of each variable for each of the three months and estimate trends in these time series (Fig. 6). Note that in each calendar month the sampling region for the regional mean time series increases over time and changes from year to year (Fig. 1b). Third, we estimate trends at each wave grid point for the periods 1971–2013 and 1992–2013 to show the spatial patterns of trend (Fig. 7).
The climatological (1970–91 and 1992–2013) means of monthly mean and maximum wind speed (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
Maps of the climatological mean wind directions (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
The climatological (1970–91 and 1992–2013) means of significant wave height (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
As in Fig. 4, but for the climatological (1970–91 and 1992–2013) means of mean wave period (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
The regional mean time series of wind speed (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
Maps of 1971–2013 and 1992–2013 linear trends in the indicated monthly mean wind speed (
Citation: Journal of Climate 28, 19; 10.1175/JCLI-D-15-0190.1
We analyzed monthly mean and maximum values of wind speed (
Considering that wind speeds, wave heights, and mean wave periods are nonnegative nonnormal data, we used the trend analysis method of Wang and Swail (2001), which is a Mann–Kendall estimator and test (Kendall 1970; Mann 1945) that accounts for the effect of lag-1 autocorrelation on trend estimates and has been found to perform best in comparison with other trend calculation methods (Hartmann et al. 2013; Zhang and Zwiers 2004). The Mann–Kendall estimator is based on ranks (nonparametric) and is thus less vulnerable to gross errors and the associated confidence interval is less sensitive to nonnormality of the parent distribution than the commonly used least squares estimator. The trend estimates in Figs. 6 and 7 (upper panels) are for the 43-yr period 1971–2013 (instead of the 44-yr period 1970–2013 because of accounting for lag-1 autocorrelation). The results are discussed in the next section.
4. Trends in surface wind speed
Figure 2 shows the climatological mean fields of monthly mean and maximum wind speeds for the two 22-yr periods, 1970–91 and 1992–2013, and the changes between these two periods. During July–September, the climatological highest wind speeds over the domain are seen in the Bering Sea (just west of Alaska); the climatological mean wind speeds increase from July to September (Fig. 2).
There were significant changes in the climates of wind speeds from the first to the second period, as shown in Fig. 2. The climatological mean of monthly mean wind speeds increased significantly from the first to the second period over the Bering Sea to north of Alaska in July, over the central part of the domain in August and September, but decreased in the region off the Canadian coasts in August (Fig. 2). The climate of monthly maximum wind speeds show similar changes, except in the region off the Canadian coasts where significant increases are seen only in July maximum wind speeds and significant decreases are seen only in August mean wind speeds (Fig. 2).
As shown in Fig. 3, the climatological mean wind direction also varies notably from July to September, and there are also notable changes in the mean wind direction between the two 22-yr periods, especially in the eastern part of the domain in August and in the northern part in September (Fig. 3). In the area west of the Canadian coast, the July and August mean wind direction has rotated clockwise, with the climatological anticyclonic center seen in this region in August being displaced from near the Alaskan north coast (71.5°N, 152°W) northeastward to around 74°N, 142°W. Another anticyclonic center can be seen in the first period in the region north of the western Siberia, around 75.5°N, 155°E in September, which shows northeastward displacement to north of 76°N in the second period.
The regional mean time series of the four monthly statistics (mean, 90th percentile, 99th percentile, and maximum) of wind speeds, along with the estimates of their 1971–2013 trends, are shown in Figs. 6a–f. Averaged over the entire BCB domain (Figs. 6a–c), all the four statistics of wind speeds show an increasing trend in each of the three months, but the increasing trends are statistically significant (at 5% level) only in September. When expressed in percentage of the corresponding 1970–99 climatological mean, the increases range from 0.1% to 0.3% yr−1. In July and August, the monthly maximum wind speed appears to have greater percentage increases than does the monthly mean, but the opposite is seen in September (Figs. 6a–c, numbers in parentheses). Averaged over the area of waves (open water) (Figs. 6d–f), the wind speed trends are statistically not significant for all statistics in all the three months.
There are also notable spatial variations of trends in wind speed
5. Trends in ocean surface waves
In this section, we present trends in four monthly statistics (mean, 90th percentile, 99th percentile, and maximum) of significant wave heights (
For the two 22-yr periods, 1970–91 and 1992–2013, the climatological mean fields of monthly mean and maximum
It is worth noting that because the wave grid is truncated at 76°N and thus does not account for wave growth north of this region, it is possible that the climatological values of
As shown in Fig. 4, the area of waves increased notably from the first to the second period, with larger increases in August and September than in July. In July, the area of waves is limited to the southern Bering Sea in the 1971–92 period; it extends northward to about 72°N and emerges in the Canadian Beaufort Sea since early 1990s (Fig. 4a). In August and September, the area of waves exists in the Canadian Beaufort Sea throughout the reanalysis period (1970–2013) and extends northwestward in the Bering Sea (Figs. 4b,c).
During July–September, the climatological highest
There were significant changes in the
The regional mean time series of the four monthly
The regional mean monthly mean wave period (
Changes in the local wind speeds alone cannot explain the significant trends in waves (
As shown in Figs. 7d–i, there are notable spatial variations of trends in
The 1992–2013 trend estimates for September mean
Since wave energy flux is proportional to the product of mean wave period and squared
6. Summary and discussion
We have analyzed the historical trends/changes in surface wind speed and ocean surface waves in the Beaufort–Chukchi–Bering (BCB) Seas using the MSC Beaufort Wind and Wave Reanalysis for the period 1970–2013. Considering the changing open water area, we characterize trends in terms of changes in climatological mean fields between two subperiods (1970–91 and 1992–2013), as well as trends in the regional mean time series of each variable for July, August, and September, separately. The 1971–2013 and 1992–2013 trends at each wave grid point were also estimated to show the spatial patterns of trend.
The results show that the climatological means of monthly mean
This study also shows that the regional mean of monthly mean
The results also show that changes in the local wind speeds alone cannot explain the significant changes in waves (
The trend of lengthening mean wave period and increasing wave height imply a trend of increasing wave energy flux, providing a mechanism to break up sea ice and accelerate ice retreat (Thomson and Rogers 2014). This stresses the importance of studying waves in the emerging Arctic Ocean, including characterizing historical climate and trends, as well as projections of possible future wave conditions. Dynamical projections of future scenarios of wave climate in the emerging Arctic Ocean, using the surface winds simulated by phase 5 of the Coupled Model Intercomparison Project (CMIP5; Taylor et al. 2012), are currently ongoing and will be presented in a separate study.
There are other wave parameters, such as wave steepness and the fraction of swells, that are worth analyzing, especially for the changing Arctic Ocean. Currently, the Environment Canada’s Beaufort Wave Reanalysis is being updated to cover the recent years. We plan to assess changes in wave steepness and in the fraction of swells in this region using an updated version of the wave reanalysis dataset.
Acknowledgments
The authors are grateful to Drs. Greg Flato and Mercè Casas-Prat for their helpful internal review of an earlier version of this manuscript. The wind and wave data (monthly statistics) analyzed in this study will be available via request made to the corresponding author: xiaolan.wang@ec.gc.ca.
REFERENCES
Cox, A. T., J. A. Greenwood, V. J. Cardone, and V. R. Swail, 1995: An interactive objective kinematic analysis system. Proc. Fourth Int. Workshop on Wave Hindcasting and Forecasting, Banff, AB, Canada, Atmospheric Environment Service, 109–118.
Francis, O. P., G. G. Planteleev, and V. E. Atkinson, 2011: Ocean wave conditions in the Chukchi Sea from satellite and in situ observations. Geophys. Res. Lett., 38, L24610, doi:10.1029/2011GL049839.
Hartmann, D. L., and Coauthors, 2013: Observations: Atmosphere and surface. Climate Change 2013: The Physical Science Basis, T. F. Stocker et al., Eds., Cambridge University Press, 159–254.
IPCC, 2013: Summary for policymakers. Climate Change 2013: The Physical Science Basis. T. F. Stocker et al., Eds., Cambridge University Press, 3–29.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471, doi:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
Kendall, M. G., 1970: Rank Correlation Methods. 4th ed. Charles Griffin, 202 pp.
Khandekar, M. L., R. Lalbeharry, and V. J. Cardone, 1994: The performance of the Canadian spectral ocean wave model (CSOWM) during the Grand Banks ERS-1 SAR wave spectra validation experiment. Atmos.–Ocean, 32, 31–60, doi:10.1080/07055900.1994.9649489.
Kohout, A. L., M. J. M. Williams, S. M. Dean, and M. H. Meylan, 2014: Storm-induced sea-ice breakup and the implications for ice extent. Nature, 509, 604–607, doi:10.1038/nature13262.
Mann, H. B., 1945: Non-parametric tests against trend. Econometrica, 13, 245–259, doi:10.2307/1907187.
Semedo, A., K. Suselj, A. Rutgersson, and A. Sterl, 2011: A global view on the wind sea and swell climate and variability from ERA-40. J. Climate, 24, 1461–1479, doi:10.1175/2010JCLI3718.1.
Stephenson, S. R., L. C. Smith, and J. A. Agnew, 2011: Divergent long-term trajectories of human access to the Arctic. Nat. Climate Change, 1, 156–160, doi:10.1038/nclimate1120.
Swail, V. R., V. J. Cardone, B. Callahan, M. Ferguson, D. J. Gummer, and A. T. Cox, 2007: The MSC Beaufort Wind and Wave Reanalysis. Proc. 10th Int. Workshop on Wave Hindcasting and Forecasting and Coastal Hazard Symp., North Shore, Oahu, Hawaii, WMO/IOC Joint Technical Commission for Oceanography and Marine Meteorology. [Available online at http://www.waveworkshop.org/10thWaves/Papers/The%20MSC%20Beaufort%20Wind%20and%20Wave%20Reanalysis.pdf.]
Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498, doi:10.1175/BAMS-D-11-00094.1.
Thomson, J., and W. E. Rogers, 2014: Swell and sea in the emerging Arctic Ocean. Geophys. Res. Lett., 41, 3136–3140, doi:10.1002/2014GL059983.
Wang, X. L., and V. R. Swail, 2001: Changes of extreme wave heights in Northern Hemisphere oceans and related atmospheric circulation regimes. J. Climate, 14, 2204–2221, doi:10.1175/1520-0442(2001)014<2204:COEWHI>2.0.CO;2.
Wang, X. L., and V. R. Swail, 2002: Trends of Atlantic wave extremes as simulated in a 40-year wave hindcast using kinematically reanalyzed wind fields. J. Climate, 15, 1020–1035, doi:10.1175/1520-0442(2002)015<1020:TOAWEA>2.0.CO;2.
Zhang, X., and F. Zwiers, 2004: Comment on “Applicability of prewhitening to eliminate the influence of serial correlation on the Mann–Kendall test” by Sheng Yue and Chun Yuan Wang. Water Resour. Res., 40, W03805, doi:10.1029/2003WR002073.