Arctic Sea Ice : Decadal Simulations and Future Scenarios Using BESM-OA

Important international reports and a significant number of scientific publications have reported on the abrupt decline of Arctic sea ice and its impact on the Global Climate System. In this paper, we evaluated the ability of the newly implemented Brazilian Earth System Model (BESM-OA) to represent Arctic sea ice and sensitivity to CO2 forcing, using decadal simulations (1980-2012) and future scenarios (2006-2100). We validated our results with satellite observations and compared them to Coupled Model Intercomparison Project, Phase 5 (CMIP5) for the same numerical experiment. BESM results for the seasonal cycle are consistent with CMIP5 models and observations. However, almost all models tend to overestimate sea ice extent in March compared to observations. The correct evaluation of minimum record of sea ice, in terms of time, spatial and area remains a limitation in Coupled Global Climate Models. Looking to spatial patterns, we found a systematic model error in September sea ice cover between the Beaufort Sea and East Siberia for most models. Future scenarios show a decrease in sea ice extent in response to an increase in radiative forcing for all models. From the year 2045 onwards, all models show a dramatic shrinking in sea ice and ice free conditions at the end of the melting season. The projected future sea ice loss is explained by the combined effects of the amplified warming in northern hemisphere high latitudes and feedbacks processes.


Introduction
Sea ice is an important and complex component of the global climate system acting both as an indicator as well as an amplifier of climate change [1]- [3].Notz and Marotzke [4] and Doescher et al. [5], indicate that sea ice cover is a more robust indicator of climate change than temperature trends alone, because sea ice changes depend on integrated changes in atmospheric and ocean variables with non linear impacts, on various temporal and spatial scales under global climate forcing.
Over the last 30 years, abrupt changes in sea ice have become evident in the Arctic, especially in the summer months of 2007 and 2012 when Sea Ice Extent (SIE) reached a minimum record extent of 4.2 × 10 6 km 2 and 3.4 × 10 6 km 2 , respectively.Satellite data have shown that the sea-ice loss has happened faster than forecasted and is unprecedented in the past 1.5 millennia [1] [6]- [8].
Sea ice age and sea ice thickness also have decreased rapidly resulting in a sea ice more sensitive to dynamic and thermodynamic forcing [1] [7].There is an agreement among scientists about the direct relationship between the shrinking of the Arctic sea ice and global warming.According to Holland and Bitz [3] and Curry et al. [9], the range of simulated polar warming is from 1.5 to 4.5 times the global mean warming and is widely related to the sea ice-albedo feedback mechanism.Most of the climate models agree that the global air temperature will continue to rise, particularly in northern high latitudes and the Arctic will become ice free in the summer in approximately 30 years, as a response to an increase in atmospheric greenhouse gas concentrations [1] [5].The impacts of melting sea ice in recent and future decades have not yet been fully understood and accurately quantified.Nonetheless, recent studies suggest that sea ice loss is linked to cold winter extremes in the northern continents, hot summer extremes over mid-latitudes continents, as well as wet summers and flooding in Eurasia [10]- [13].
Besides the important role of sea ice in the climate system, knowing the dynamics and geographical sea ice cover is also essential for human activities such as navigation, oil exploration and fishery [14] [15].According to Cochran et al. [16] and Meier et al. [15], changes in Arctic threaten the infrastructure, health and safety of the Arctic indigenous people as well as present a significant risk to local marine biodiversity According to Whiteman et al. [17], sea ice changes will affect all nations, not just those in the world's far north, and all should be concerned about changes that are happening in the Arctic region.In that sense, Global Climate Models, even with inherent uncertainties and limitations are powerful tools for better understanding the changes in sea ice as well as providing future scenarios to guide decision markers, governments and local communities among others.
The recent development of the Brazilian Earth System Model (BESM) is an effort of several institutions and researchers lead by the Brazilian National Institute for Space Research (INPE) to build a multidisciplinary research framework with the intent to understand the causes of global climate change, its effects and its impacts on society.The BESM model, also aims to contribute to Program for Climate Model Diagnostics and Intercomparison (PCMDI) with short-term and long-term simulations, as well as to provide futures scenarios of climate change [18].Based on several studies and reports [5] [14] [17] [19] and understanding the importance of sea ice in the global climate system and the global economy, BESM simulations are expected to contribute with, among other variables, sea ice short and long-term simulations.BESM simulations can also be useful for future studies on ocean-atmosphere-sea ice coupling processes and impacts of sea ice loss around the world.
In this paper, we evaluated decadal simulations (1980-2012) and future scenarios (2006-2100) of SIE as simulated by two versions of the BESM and by other Coupled General Circulation Models participating in the Coupled Model Intercomparison Project, Phase 5 (CMIP5).Our goal is to evaluate the first results on the ability of BESM to represent past and future sea ice changes and sensitivity of the sea ice to the radiative forcing, using the Taylor protocol [20] [21].The paper is structured as follows: first, we present the data sources in Section 2.Then, in Section 3, we examine the seasonal cycle, the spatial pattern and the minimum records of the Arctic sea ice, comparing the BESM decadal simulations to satellite observations and other CMIP5 models.In Section 4, we investigate the future scenarios for two different versions of BESM and the CMIP5 models, using two different scenarios, the Representative Concentration Pathway RCP4.5 and RCP8.5.We discuss the results and also indicate possible causes to explain the differences between the sea ice variation using BESM versions 2.3 and 2.5.Finally, in Section 5 we present our conclusions and lay out our recommendations for future work.

Data Sources
This study uses short-term simulations (decadal hindcasts) and long-term simulations (future scenarios) of 11 state-of-the-art General Circulation Models (GCMs) and Earth System Models (ESMs), seen in Table 1.The numerical experiment design follows the CMIP5 protocol, for decadal data and future projections based on the Taylor protocol [20] [21].CMIP is an international effort of the scientific community to provide simulations of many different climate models in order to better understand past and future climate changes as well to provide a scientific data set for the Intergovernmental Panel on Climate Change (IPCC).
The BESM ensemble members of the decadal simulations were integrated for 10 years, each with initial conditions (IC) on 1 -10 December of the years 1960, 1965, 1970, 1975, 1980, 1985, 1990, 1995, 2000 and 2005.Three of these ensembles (1960, 1990 and 2005) were extended for an extra 20 years for each of the 10 members, completing 30 years long integrations each.These simulations used atmospheric CO2 concentrations derived from in situ air samples collected at Mauna Loa Observatory, Hawaii [18].The Atmospheric model initial conditions for each ensemble member used the National Centers for Environmental Prediction (NCEP-NCAR) reanalysis fields for the 0000 UTC of each day from 1 to 10 December of the chosen years.The ocean initial states were chosen from the same dates from a spinup run of MOM4p1 that used prescribed atmospheric fields of momentum, solar radiation, air temperature, and freshwater described in Nobre et al. [18].
The future scenarios are defined by the Representative Concentration Pathways (RCPs) and each RCP defines a specific emissions trajectory and subsequent radiative forcing.The radiative forcing values in the year 2100 relative to pre-industrial values are 4.5 W•m −2 and 8.5 W•m −2 for RCP4.5 and RCP8.5 respectively, which include the period from 2006 to 2100.The CO 2 concentration in the year 2100 for each RCP is approximately 600 ppm and 1300 ppm for RCP4.5 and RCP8.5, respectively.
We compared BESM results with CMIP5 models simulations using the same numerical experiment setup.Still, the models differ in spatial resolution, physical component and parameterizations.For decadal simulations we chose to work with time series from 1980 to 2012 due to the availability of satellite observations for comparison.The SSM/I (Special Sensor Microwave Imager) satellite observations obtained from the National Snow and Ice Data Center (NSIDC) were used to validate the numerical simulations.For all simulations we calculated the SIE, defined as the area where the sea ice concentration is greater than 15% in a grid.

BESM-OA Model
In this work, we used two versions of the BESM Coupled Ocean Atmosphere (BESM-OA) model: BESM-OA V2.3 for decadal and RCP simulations and BESM-OA V2.5 for RCP simulations only.The main differences between these two versions are the microphysics scheme proposed by Ferrier et al. [31] and a new surface layer scheme based on Jimenez and Dudhia [32] described by Capistrano et al. [33] [34].
Both BESM versions used in this research are composed of the INPE/CPTEC atmospheric general circulation model (AGCM) coupled to NOAA/GFDL's Modular Ocean Model version 4p1 (MOM4p1) oceanic general circulation model (OGCM) via GFDL's Flexible Modular System [18] [35] [36].The INPE/CPTEC AGCM has a spectral horizontal resolution truncated at triangular wave number 62, giving an equivalent grid size of 1.8758 degrees of latitude and longitude and 28 sigma levels unevenly spaced in the vertical (i.e., T062L28).The exchanges of heat, moisture and momentum between the surface and atmosphere in INPE/CPTEC AGCM over the ocean and continents are computed differently by various physical processes that define the surface fluxes.
The ocean model MOM4p1 [35] from GFDL, includes the Sea Ice Simulator (SIS), described in Winton [37].The SIS is a dynamical model with three vertical layers (two ice and one snow), and five ice thickness categories.The elastic-viscous-plastic technique of Hunke and Dukowicz [38] is used to calculate ice internal stresses, and the thermodynamics is a modified Semtner's three-layer scheme [39].SIS calculates the concentration, thickness, temperature, brine content, and snow cover of an arbitrary number of sea ice thickness categories (including open water) as well the motion of the complete pack.Additionally, the model is responsible for calculating ice/ocean fluxes and communicating fluxes between the ocean and atmosphere models globally.
The MOM4p1 horizontal grid resolution is set to 1˚ in the longitudinal direction, and in the latitudinal direction the grid spacing is 1/4˚ in the tropical region (10˚S -10˚N), decreasing uniformly to 1˚ at 45˚ and to 2˚ at 90˚ in both hemispheres.For the vertical axis, 50 levels are adopted with a 10 m resolution in the upper 220 m, increasing gradually to about 370 m of grid spacing in deeper layers.We used FMS to coupling MOM4p1 and CPTEC/AGCM.Thus, wind stress fields are computed, using Monin-Obukhov scheme within MOM4p1, from the field 10 meters above the above the ocean surface.Adjustments were done to the Monin-Obukhov boundary layer scheme, whose parameters were tuned according to the wind fields output by the CPTEC AGCM.The AGCM receives the following two fields from the coupler: sea surface temperature (SST) and ocean albedo from ocean and sea ice models at an hourly rate (coupling time step).Adjustments were also made to ocean shortwave penetration parameters due to the CTPEC AGCM supply of visible and infrared short wave radiation.The coupling variables supplied by the AGCM are as follows: freshwater (liquid and solid precipitation), specific humidity, heat, vertical diffusion of velocity components, momentum fluxes, and surface pressure.
The microphysics of Ferrier et al. [31] used in BESM-OA V2.5, replaced the Large Scale Precipitation scheme used in BESM-OA V2.3.This new microphysics scheme computes changes in water vapor, cloud water, rain, cloud ice and precipitation ice.Also, BESM-OA V2.5 uses a new surface layer scheme based on Jimenez and Dudhia [32] and described by Capistrano et al. [33] [34].In this scheme, the surface and the first AGCM level values are used to assess wind, air temperature and humidity at 10 m.The changes introduced lead to a more consistent surface layer formulation that resulted in a near-surface wind, air temperature and humidity more consistent with observations than previous BESM version.This occurs mainly over the ocean, where those variables are important to compute the heat fluxes at ocean-atmosphere interface.

Seasonal Cycle
Seasonal melt-freeze transitions are important to continuously monitor sea ice over the Arctic.Sea ice formation, growth and decay are closely related to air temperature, ocean heat content, albedo and heat fluxes and hence can vary strongly from month to month [5] [40].
Thus, we present in this section the Arctic seasonal cycle of sea ice, in order to better understand the differences between the models studied, with a focus on the performance of the BESM-OA V2.All of the models were able to represent the seasonal cycle of the Arctic sea ice.Large oscillations between Based on our results, Sortberg et al. [42] and Karlsson and Svensson [43], we suggest the following scheme to explain the differences between winter and summer model's performance in representing SIE.First, the presence of sea ice affects strongly the sea ice albedo, which has a key influence on the energy budget and is directly linked to the cloud-albedo effect and cloud-radiation effect.Clouds are linked with the energy budget by reflecting shortwave radiation back to space, trapping Longwave (LW) radiation and radiating it back to the surface, providing one of more the strongest feedbacks in the climate system [44].Second, in wintertime, the amount of solar radiation is low or non-existent and the ability of the clouds to reemitted LW to the surface presents a positive cloud radiative effect on the surface energy budget.On the other hand, during the seasons with solar radiation, the positive greenhouse effect is competing with a negative cloud albedo effect, because the clouds decrease the amount of incident solar radiation at the surface.Finally, recent publications using CMIP3 and CMIP5 models [42] [43] suggest that models generally have the tendency to underestimate the amount of LW radiation reemitted back to the surface in winter.As a consequence of these processes, the models tend to overestimate SIE in wintertime.Additionally, the annual amplitude of sea ice cover depends inversely on the model's sea ice albedo [42] [43].BESM-OA V2.3 results agree with this scheme, as both downward and upward LW radiation at the surface are underestimated in winter.The ensemble mean is lower than the mean of the observations by approximately 30 W•m −2 (not shown).Another notable example is related to the MPI-ESM-LR model's performance, which presents a high sea ice albedo and low annual amplitude of sea ice.According to Wild et al. [45], the bias in LW radiation depends on the climate conditions and is not geographically uniform, with higher (smaller) bias in cold and dry climates (warm and humid climates) with low (high) downward LW radiation emission.
Stroeve et al. [1], Knutti and Sedlacek [46] and Li et al. [44], assessed the evolution between CMIP3 and CMIP5 and showed an improvement in the Arctic sea ice prediction and radiation in CMIP5.Nevertheless, a better representation of sea ice depends also on improvements in the representation of the Arctic sea-ice albedo, clouds, cloud-radiation effects and feedback processes.However, when looking at separate months, the correlation drops, as consequence in the same month time series, only the interannual variability is being evaluated.For the month of March, we suggest the previously described scheme (radiation effect) to explain the low correlation between observation and models.To understand the low correlation in September, we suggest a relation with sea ice thickness.According to Shu et al. [47] and Stroeve [1], the sea ice thickness simulated in the CMIP5 models is too thin, resulting in enhanced sea ice melt and an underestimation fo SIE in summertime, as shown in Figure 1.BESM-OA V2.3 agrees with Shu et al. [47], showing an underestimate in sea ice thickness, notably greater in the marginal sea ice zones (not shown).

Spatial Pattern
Several studies have compared the observed SIE variation using climate models and CMIP data sets in a seasonal cycle or time series approach [1] [6] [41].This type of analysis is important to know the model's ability to predict SIE.However, when considering only SIE, the information related to spatial pattern is lost.Analyzing spatial patterns avoids overconfidence in the predictions and excludes compensation of errors of opposite sign in different regions [48].Cavalieri and Parkinson [49] also show the importance of evaluating the Arctic Ocean by regions.The authors, using satellite data set to analyze sea ice variability and trends from 1979 to 2010, found that trends for nine distinct regions in the Arctic are not homogeneous and indicated the complex nature of the Arctic climate system by regions.
Figure 3 shows September average Sea Ice Cover (SIC) observations over the study area.The spatial difference between modeled and observed SIC in September average is shown in Figure 4. September was chosen because commonly is when sea ice reaches its annual minimum over the Arctic.
Despite obvious inter-model differences depicted in Figure 4, there is a reasonable agreement between all the models.Most models tend to well represent SIC in the central Arctic, whereas the opposite occurs in marginal ice zones.There is a general tendency to underestimate SIC in areas such as the Beaufort Sea and the East Siberian Sea (except for MPI-ESM-LR and HadCM3) suggesting a systematic model error in this region.However, NCAR-CCSM4 overestimates SIC in both Laptev and Kara Sea (Figure 4).This may reflect the NCAR model's overestimation observed for September and shown in Figure 1.
The SIC in the region between Canada and Greenland is well represented by BESM-OA V2.3, GFDL-CM2.1 and MPI-ESM-LR models, while SIC between East Siberia and the Barents Seas is only well simulated by   4), it is clear that the amplitude of its annual cycle is smaller than both the others models and satellite observations (Figure 1).This reveals a certain deficiency in representing physical processes between ocean-atmosphere-sea ice, although the good representation during summer.
The detailed analysis of simulated SIC by regions using Climate Models is justifiable due to both economic and scientific reasons.Economically, as a result of sea ice loss maritime transports may gain two new routes with the opening of the "Northwest Passage" in Northern Canada and Greenland, and the "Northeast Passage" between Northern Russia and Norway [14].This is considered a hot topic because these passages could lead to fast and cheaper ship transport between Europe and North America.Scientifically, this is relevant because the importance of properly account for the dynamical/thermodynamical processes taking place in shaping SIC over the Arctic region.
It is instructive to compare Figure 1 with Figure 4, analyzing only SIE in September for GFDL-CM2.1 and CanCM4 models (Figure 1).This can induce overconfidence in how well the models agree (SIE in both models are approximately 5.6 × 10 6 km 2 ).However, when we look at the spatial patterns in Figure 4, we find quite different SIE distributions.CanCM4 model shows a large area of high negative values (especially between Greenland and Canada), whereas GFDL-CM2.1 shows a small area of high negative values only in parts of Beaufort Sea and East Siberian Sea.Thus, even some climate models showing a good performance in simulating SIE during summertime do not necessarily simulate a realistic spatial sea ice distribution.

Minimum of Sea Ice Extent
Changes in ice extent due to the seasonal cycle are so large that they tend to obscure any signal related to interanual variability.To remove the strong seasonal cycle, we again specifically focus on September since it shows the minimum annual of SIE.According to Doescher et al. [5], the ability to identify real changes in the Arctic Climate System increases when we focus on individual seasons.In this context, Figure 5 [50], such abrupt sea ice loss resulted from a complex interplay between the thermodynamics and dynamics of sea ice, ocean and atmosphere and successful prediction requires careful initialization with ocean and sea ice conditions.
Figure 6 illustrates the spatial distribution of average September SIE (left) and minima of September SIE value found between 1980 and 2010 (right) for all the models evaluated in this work.This figure aimed to show the model's performance to represent the spatial pattern in episodes of low SIE, regardless of year.
Looking at the spatial patterns of the SIE climatological mean and minimum record, it is clear that the climate models are able to reproduce the seasonal cycle of the SIE (Figure 1 and Figure 2) better than they represent the minimum records.Only GFDL-CM2.1 model presents a good spatial agreement of minimum records with observation.This could be explained by two main reasons.First, it was the only one that matched the spatial pattern of the observed minimum, which may lead to a better agreement with observed meteo-oceanographic patterns (not studied here).Second, it may be related to a better representation of the sea ice and feedbacks processes in the parameterizations of the GFDL-CM2.1 model.Two other models (BESM-OA V2.3 and CanCM4) also show a reasonable spatial agreement with observations, although not as well as GFDL-CM2.1.These two models presented an underestimation of SIE, but presented a very good representation of the SIE in  the central Arctic region (Figure 6).The minimum record for the BESM-OA V2.3 shows a deficiency near Greenland and at the north of Canada.We understand that this happens because of the overestimation of Sea Surface Temperature (SST) in that region by the BESM-OA V2.3 model (not shown).Although BESM-OA V2.3 and CanCM4 were capable to capture the correct signature of the SIE minimum record with a decrease in SIE followed by an increase in the following year (Figure 5), the correct estimation of minimum SIE, in time, spatial, area and processes signatures remains a challenge for the modeling community.
Due to the sea ice retreat in recent Septembers months, ice cover in the following spring tends to be thinner, thus vulnerable to melting in summer.According to Doscher et al. [5], each record of low SIE is followed by a partial recovery.Additionally, Tietsche et al. [48] suggest that the minimum record of SIE during a single September is reversible, as the albedo sea ice mechanism is compensated by large scale recovery mechanisms.According to Vihma [51] the sea ice loss increases the heat flux from the ocean to the atmosphere in early winter and autumn.As result of this, a local increase of air temperature, humidity and cloud cover is expected thus reducing the stability of the atmospheric boundary layer.
Hunke et al. [52] evaluated the retrospective and new directions of sea ice models.The authors indicated some deficiencies in the dynamics (e.g.transport processes, dynamic coupling and mechanical redistribution) and thermodynamics (e.g.feedback processes and melt ponds) and suggested that improvements in the sea ice prediction dependent on improvements in the descriptions of the physical processes and characteristics, as well as, extending the models for Earth System Model simulations including biogeochemistry.According to Flocoo et al. [53] and Roeckner et al. [54], one of the processes, poorly represented in sea ice models, is the formation and evolution of melt ponds.Melt ponds affect the heat and mass balances of SIC, mainly by reducing albedo by up to 20%.Consequentially, a reduction of the sea ice volume can reach 40%, leading to further sea ice melt.At the end of the melting season, melt ponds cover up to 50% of the sea ice surface.A better representation of the melt pond scheme will improve the sea ice simulation and is essential for accurate future sea ice projections.

Future Projection of Arctic Sea Ice
The long-term evolution of SIE in the northern hemisphere as simulated by BESM and CMIP5 models, using RCP4.5 and RCP8.5 is shown in Figure 7.The simulations clearly show a decrease in SIE up to 2100, for all simulations and both RCPs.Arctic SIE decline with the increase of the radiative forcing in all models.The BESM-OA V2.3 control experiment (gray lines in Figure 7) reinforce that ice-free conditions only happen when external forcing from anthropogenic sources are include in climate model simulations.These results are in agreement with Stroeve et al. [1].
For September, at the beginning of the series (2006 to present-day), the HadGEM2 model SIE values are close to satellite observations.During March the best representation of the observed data was obtained by MIROC5 and BESM-OA V2.5 models.It is possible to observe a higher inter-annual variability in September than in March for all models, as well as for the early period's satellite observations.The changes in inter-annual variability are important for sea ice prediction and frequency and for assessing the frequency of occurrence of extreme SIE anomalies.
It is noteworthy that the models comparatively show different tendencies for the months of maximum and minimum SIE.For the month of March, the MPI-ESM-LR model presented the lowest values compared to the other models used here, whereas for the month of September the lowest values encountered are those of the BESM-OA V2.5 model.In general, when compared to the other CMIP5 models, the BESM-OA V2.3 model tends to overestimate SIE in March and September, for both RCP simulations.
From the year 2040 onward, all models show a dramatic shrinking in SIE in the RCP 8.5 scenario.This indicates a high sensitivity of sea ice cover in response to an increase in the atmospheric carbon dioxide.The GFDL-CM3 model clearly shows this abrupt decrease in SIE with the RCP8.5 scenario when compared to the RCP4.5 one.In this case, the decline is so strong that at the end of the 21st century the SIE maximum (in March) is similar to the minimum SIE (in September) found in the beginning of the 21st century.If the GFDL-CM3 model is reasonably correct, it means that the Arctic can be ice-free also during the coldest season of the year just after 2100.
For September, ice-free conditions (defined as less than 0.5 × 10 6 km 2 ) are obtained from 2020 in Can-ESM2 model, BESM-OA V2.5 and HAGEM2-2S with the RCP8.5 scenario.According to Chylek et al. [23] the addition of the land-vegetation model and terrestrial oceanic interactive carbon cycle to the coupled atmosphere-ocean in the Can-ESM2 model improved the simulations, although increased the overstimation of atmospheric warming after 1970.That explains the minimum values found here for CAN-ESM2 sea ice projections.Also focusing on the RCP8.5 scenario, most of the models show ice free situations, or episodes, after 2045 for the month of September.The exceptions are BESM-OA V2.3 and NCAR-CCSM4.These two models are a bit more conservative than the others, pointing out to ice free situations to starting after 2060.
It is expected that ice-free conditions will have strong effects on the global climate system though changes in both ocean and atmospheric circulations.It is known that sea ice loss amplifies the effects of radiative forcing by the albedo-sea ice feedback mechanism and cloud effects.It also, affects the meridional and inter-hemispheric temperature gradients that can affect mid-latitude circulation.However the quantification of these effects remains unclear requiring improvements in the global climate models.

Surface Anomalies Temperatures
In this section, we compare Surface Anomalies of Temperatures (SAT) for BESM-OA V2.3 and BESM-OA V2.5 to explain the differences between those versions in SIE presented in Figure 8.
Figure 8 shows SAT and Total Cloud Cover for BESM-OA 2.3 and BESM-OA 2.5 using future scenarios, relative to the period from 2006 to 2100.A marked warming in the northern high-latitudes is observed in both BESM versions, being notably stronger in the RCP8.5 simulation.This warming called Polar Amplification (PA) occurs due to the increase in the atmospheric greenhouse gas concentration, and is accompanied by an expressive reduction in SIE in both simulations (Figure 8).The relationship among air temperature rises and sea ice loss is evident and underpinned statistically [5].PA is associated with several feedback processes as the ice-albedo feedback, temperature, water vapor and clouds.Most of studies indicate that the ice-albedo feedback is the main contributor for PA [2] [3] [9].However Pithan and Maurtsen [55], using CMIP5 simulations found that, the major contributor to PA comes from air temperature feedbacks (as the surface warms, more energy is radiated into space in lower latitudes compared with the Arctic region).
When comparing the warming between the two BESM versions, we observe that higher amplitude values are observed in BESM-OA V2.5, particularly in high latitudes between 75˚N to 90˚N.As a result of these warming discrepancies the SIE is lower in BESM-OA V2.5 than in BESM-OA V2.3 (Figure 7).
The microphysics of Ferrier et al. [31] and the new surface scheme based on Jimenez and Dudhia [32] used in BESM-OA V2.5 produced an improved in the representation of precipitation, wind, air temperature, humidity, and energy balance at the top of the atmosphere (not shown).A better representation of these variables exerts strong influences in coupling ocean-atmosphere-sea ice simulations and teleconnections with higher latitudes.The microphysics adopted in BESM-OA V2.5 produced a decrease in the total cloud cover in the Arctic region (Figure 8).This allows the ocean to absorb more heat from incident shortwave radiation and then contribute to a greater melting of the sea ice.The decrease in total cloud cover and consequent strong increase in SAT is consistent with the SIE reduction showed in Figure 7.
According to Jiang et al. [56], Clouds (in both ice and liquid forms) are important modulators of the climate system and are involved in several feedback processes that affect the global atmospheric circulation and the energy budget.Improving the cloud microphysics in Coupled Climate Models result in an advance in climate prediction and reduce the uncertainties in future projections.As recently pointed out by Eyring et al. [57], the understanding of the role of the clouds in the general atmospheric circulation, climate sensitivity and assessing the response of the cryosphere system to a warming climate, are among the greatest challenges for CMIP6.

Conclusion
In this work, we evaluated the decadal simulation  and assessed the future climate projections (2006-2100) generated by BESM-OA and CMIP5 models.BESM-OA V2.3 results for the seasonal cycle are consistent with satellites observations and the other CMIP5 models, however almost all models tend to overestimate SIE in March in relation to observations.Based on our results and [42] [43] [45], we suggest that the winter Arctic SIE bias is related to a LW radiation bias in climate models.Spatial patterns of climatological averages at the end of the melting season presented a deficiency in capturing the correct signature of the minimum SIE record, as well a model systematic model error between Beaufort Sea and East Siberia (Figure 4 and Figure 6).Future scenarios show an abrupt shrinking of the sea ice and ice-free summer conditions from the year 2045 and onwards, for both RCPs projections.This is a result of the internal climate response to the changing in radiative forcing throughout the years.Polar Amplification and feedback processes explain the rapid Arctic sea ice loss, despite the uncertainties and limitations of Global Climate Models.The sea ice responses are different in CMIP5 models due to differences in the ocean, atmosphere, sea ice conditions, as well the coupling between the components in each model.Future progress in sea ice modeling is essential and requires advances in the parameterizations of climate feedback processes.The climate in the Arctic region will change even more and will induce complex changes in the global climate, thereafter will induce changes in Arctic climate over again.In synthesis, we can say that the Arctic region and its climate are way more complex than forecasted.
3 model.First, to understand the ability of BESM-OA V2.3 to simulate the seasonal cycle in relation to observation and other CMIP5 models, we present in Figure 1 the seasonal cycle of climatological average of SIE from CMIP5 model and observed values for the period from 1980 to 2010.

Figure 1 .
Figure 1.Climatology of SIE (1980 to 2010) in the northern hemisphere simulated by BESM-OA V2.3, CMIP5 models and observations.summerand winter are evident, with sea ice growing from autumn and winter reaching a peak in March, and then declining throughout spring and summer as the melting season progresses.However, most models overestimate SIE values in winter (except the MPI-ESM-LR model), and underestimate in summer (except HadCM3 and NCAR-CCSM4).BESM-OA V2.3 ensemble agrees quite well with observations and satisfactorily represents the seasonality of sea ice, although the model's sea ice decays more rapidly than observed in summer and autumn.The observational data (BESM-OA V2.3) shows that Arctic SIE varies between approximately 15 × 10 6 km 2 (18 × 10 6 km 2 ) at winter maximum and 6 × 10 6 km 2 (6 × 10 6 km 2 ) at summer minimum.The difference between the model's performance for winter and summer are in agreement with[1] [26] [41][42].It is clear that the BESM-OA V2.3 model, even with an overestimation during winter presents a very good agreement in summer when SIE reaches critical values.Based on our results, Sortberg et al.[42] and Karlsson and Svensson[43], we suggest the following scheme to explain the differences between winter and summer model's performance in representing SIE.First, the presence of sea ice affects strongly the sea ice albedo, which has a key influence on the energy budget and is directly linked to the cloud-albedo effect and cloud-radiation effect.Clouds are linked with the energy budget by reflecting shortwave radiation back to space, trapping Longwave (LW) radiation and radiating it back to the surface, providing one of more the strongest feedbacks in the climate system[44].Second, in wintertime, the amount of solar radiation is low or non-existent and the ability of the clouds to reemitted LW to the surface presents a positive cloud radiative effect on the surface energy budget.On the other hand, during the seasons with solar radiation, the positive greenhouse effect is competing with a negative cloud albedo effect, because the clouds decrease the amount of incident solar radiation at the surface.Finally, recent publications using CMIP3 and CMIP5 models[42] [43] suggest that models generally have the tendency to underestimate the amount of LW radiation reemitted back to the surface in winter.As a consequence of these processes, the models tend to overestimate SIE in wintertime.Additionally, the annual amplitude of sea ice cover depends inversely on the model's sea ice albedo[42] [43].BESM-OA V2.3 results agree with this scheme, as both downward and upward LW radiation at the surface are underestimated in winter.The ensemble mean is lower than the mean of the observations by approximately 30 W•m −2 (not shown).Another notable example is related to the MPI-ESM-LR

Figure 2
shows a Taylor diagram for September, March and annual climatology of SIE.This diagram is a useful tool to compare observed and simulated data in terms of correlation coefficient, RMS and standard deviation.A shorter distance between model and REF (observed) in a Taylor Diagram indicates a better model's performance.For the annual values (black) all the six models have a correlation with the observations higher than 0.96, while for March and September the correlation coefficient presents low values, as expected.For all models (except MPI-ESM-LR) the correlation in March is smaller than 0.6.The annual cycle of SIE is quite well represented because the seasonal cycle of SIE was well represented by all the models as shown in Figure 1.

Figure 2 .
Figure 2. Taylor diagram of September, March and climatological annual cycle of SIE for the period 1980-2012.The x-axis and y-axis is the standard deviation normalized.The correlation coefficient between observations and each model is given by the azimuthal position.The centered RMS difference between simulated and observed is proportional to their distance one from another.

Figure 3 .
Figure 3. Arctic study area and September SIC Climatology (1980-2010) from satellite observations (shaded colors).Dark gray and orange lines refer to the 2007 and 2012 minimum events respectively.
shows the SIE time series of September averages from 1982 to 2014.Arctic sea ice has declined sharply during the last three decades, with record low summer ice cover in September 2007 and 2012 as illustrated in Figure 1.Here, we show the time series of SIE and analyze the ability of the models to represent recent changes.Arctic SIE averages from 1980 to 2010 (Figure 5) show a noticeable decrease in Arctic SIE.September SIE simulated by BESM-OA V2.3 (satellite observations) between 2000 and 2010 was 4.2 × 10 6 km 2 (5.7 × 10 6 km 2 ), while between 1980 and 1990 it was 6 × 10 6 km 2 (7.1 × 10 6 km 2 ), showing a reduction of approximately 30% (19.8%) in SIE.March SIE for example (not shown), although at a slower rate in comparison to September, also decreases with time.The minimum satellite record of SIE occurred in September 2012, with 3.6 × 10 6 km 2 against 4.5 × 10 6 km 2 in BESM-OA V2.3.In 2012, satellite observations (GFDL-CM2.1)presented a decrease of 50.4% (42.3%) of SIE in relation to the 1980s decade.Except for the GFDL-CM2.1 model, no other model was able to well represent the observed 2012 minimum.However, BESM-OA V2.3 and CanCM4, generated episodes of low SIE in September with a magnitude and behavior comparable with the low observed in 2012.These episodes of minimum SIE occurred in 2006 and 2002 for the BESM-OA V2.3 and CanCM4 models, respectively.According to Doscher et al. [5] and Holland et al.

Figure 5 .
Figure 5. Arctic sea ice extent time series of September from 1980 to 2014 for CMIP5 models and observational data.

Figure 6 .
Figure 6.Spatial distribution of SIE average (left) and lowest values of September SIE found between 1980 to 2010, (right), for all CMIP5 models evaluated in this work.

Figure 7 .
Figure 7. Time series of modeled Arctic SIE in September and March from 2006 to 2100, using Representative Concentration Pathways RCP4.5 (solid lines) and RCP8.5 (dash lines).Black lines are the satellite observations and gray lines refer to the control run of the BESM-OA V2.3.During the first 30 years of the series, values from both RCPs are very similar in March and September months.For March, the SIE in the first years of the 21st century ranges from 11.8 × 10 6 km 2 (MPI-ESM-LR) to 18.8 × 10 6 km 2 (BESM-OA V2.3) in the RCP45 simulation.For the RCP8.5 simulation they vary between 11.7 × 10 6 km 2 (MPI-ESM-LR) and 18.4 × 10 6 km 2 (BESM-OA V2.3).For September SIE values vary from approximately 2 × 10 6 km 2 (BESM-OA V2.3) to 6.7 × 10 6 km 2 (MIROC 5).Already during these early years, it is possible to observe the discrepancy between the two different BESM configurations.For all RCP simulations, BESM-OA V2.3 and BESM-OA V2.5 show higher (lower) values in SIE during March (September) when compared to other models used in this work.The models reveal strong amplitudes in SIE between different seasons.Both BESM simulations clearly present the higher values in March for all years.However, for September the higher SIE was found in MIROC5 model (similar amplitude was observed in MPI-ESM-LR model).It is possible to observe a higher inter-annual variability in September than in March for all models, as well as for the early period's satellite observations.The changes in inter-annual variability are important for sea ice prediction and frequency and for assessing the frequency of occurrence of extreme SIE anomalies.It is noteworthy that the models comparatively show different tendencies for the months of maximum and minimum SIE.For the month of March, the MPI-ESM-LR model presented the lowest values compared to the other models used here, whereas for the month of September the lowest values encountered are those of the BESM-OA V2.5 model.In general, when compared to the other CMIP5 models, the BESM-OA V2.3 model tends to overestimate SIE in March and September, for both RCP simulations.From the year 2040 onward, all models show a dramatic shrinking in SIE in the RCP 8.5 scenario.This indicates a high sensitivity of sea ice cover in response to an increase in the atmospheric carbon dioxide.The GFDL-CM3 model clearly shows this abrupt decrease in SIE with the RCP8.5 scenario when compared to the RCP4.5 one.In this case, the decline is so strong that at the end of the 21st century the SIE maximum (in March) is similar to the minimum SIE (in September) found in the beginning of the 21st century.If the GFDL-CM3 model is reasonably correct, it means that the Arctic can be ice-free also during the coldest season of the year just after 2100.For September, ice-free conditions (defined as less than 0.5 × 10 6 km 2 ) are obtained from 2020 in Can-ESM2 model, BESM-OA V2.5 and HAGEM2-2S with the RCP8.5 scenario.According to Chylek et al.[23] the addition of the land-vegetation model and terrestrial oceanic interactive carbon cycle to the coupled atmosphere-ocean in the Can-ESM2 model improved the simulations, although increased the overstimation of atmospheric warming after 1970.That explains the minimum values found here for CAN-ESM2 sea ice projections.

Figure 8 .
Figure 8. Surface Anomalies Temperature (SAT) and total cloud cover from January 2006 to December 2100 for BESM-OA V2.3 and BESM-OA V2.For SAT, the red lines represent the average for latitudes between 75˚N to 90˚N.Green and blue lines are for latitudes between 45˚N to 75˚N and 25˚N to 45˚N, respectively.Latitudes between 0˚N to 90˚N are represented by the black lines.