A new high-resolution atmosphere-ocean coupled general circulation model named MIROC4h has been developed, and its performance in a 120-year control experiment (including a 50-year spin-up) under the ...present conditions (the year 1950) is examined. The results of the control experiment by MIROC4h are compared with simulations of preindustrial conditions carried out for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4) using the previous high- and medium-resolution versions of the model, called MIROC3h and MIROC3m, respectively. A major change in MIROC4h is a doubling of the resolution of the atmospheric component to 0.5625°, compared to 1.125° for MIROC3h. The oceanic components of MIROC4h and MIROC3h are eddy-permitting, with a horizontal resolution of 0.28125° (zonal) × 0.1875° (meridional). In MIROC3m, the horizontal resolution is 2.8125° for the atmospheric component and 1.40625° (zonal) × 0.56°–1.4° (meridional) for the ocean component. Compared with MIROC3h and MIROC3m, many improvements have been achieved; for example, errors in the surface air temperature and sea surface temperature are smaller, there is less drift of the ocean water temperature in the subsurface-deep ocean, and the frequency of heavy rain is comparable to observations. The fine horizontal resolution in the atmosphere makes orographic wind and its effects on the ocean more realistic than those of the former models, and the treatment of coastal upwelling motion in the ocean has been improved. Phenomena in the atmosphere and ocean related to the El Niño and southern oscillation are now closer to observations than was obtained by MIROC3h and MIROC3m. The effective climate sensitivity for CO2 doubling is calculated to be about 5.7 K, which is much larger than the value obtained using the IPCC AR4 models, and is mainly due to a decrease in the low-level clouds at low latitudes.
We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, ...spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.
The correlation between transport via the Bering Strait throughflow (BTF) and sea surface salinity (SSS) in the Bering Sea has been examined mainly using an atmosphere–ocean–ice coupled climate model ...that has an eddy-permitting ocean component. The SSS anomaly in the northwestern Bering Sea is high from winter to spring when the BTF transport anomaly is large in the cold season. Similar features can be seen in an observation dataset and two kinds of ocean data assimilation product. BTF transport is strongly correlated with sea surface height (SSH) in the northeastern Bering Sea, the southwestern Chukchi Sea and the East Siberian Sea. The SSH along the Russian coast in the Arctic Ocean is uncorrelated with the SSH in the Bering Sea, the meaning being that the Arctic SSH affects the BTF and the SSS independently of the SSH in the Bering Sea. The low SSH along the Siberian coast is correlated with easterly wind anomalies over the Laptev Sea and north of the New Siberian Islands. The relationship between the low Siberian-coast SSH and the high SSS in the Bering Sea, however, is not confirmed in 10 years of satellite-derived SSH. Mixed-layer salt budget analysis has revealed that the high SSS anomalies are mainly caused by the increases of horizontal advective salt convergence north of 62.5°N, and by the decreases of sea-ice melting south of 62.5°N, through the strengthening of the near-surface northeastward currents. In the warm season, these two factors fade and the salinization disappears.
Decadal climate predictability is examined in hindcast experiments by a multi-model ensemble using three versions of the coupled atmosphere-ocean model MIROC. In these hindcast experiments, initial ...conditions are obtained from an anomaly assimilation procedure using the observed oceanic temperature and salinity with prescribed natural and anthropogenic forcings on the basis of the historical data and future emission scenarios in the Intergovernmental Panel of Climate Change. Results of the multi-model ensemble in our hindcast experiments show that predictability of surface air temperature (SAT) anomalies on decadal timescales mostly originates from externally forced variability. Although the predictable component of internally generated variability has considerably smaller SAT variance than that of externally forced variability, ocean subsurface temperature variability has predictive skills over almost a decade, particularly in the North Pacific and the North Atlantic where dominant signals associated with Pacific decadal oscillation (PDO) and the Atlantic multidecadal oscillation (AMO) are observed. Initialization enhances the predictive skills of AMO and PDO indices and slightly improves those of global mean temperature anomalies. Improvement of these predictive skills in the multi-model ensemble is higher than that in a single-model ensemble.
► We evaluate the heat flux in an ice-ocean model with the ice thickness distribution. ► The conductive heat flux is underestimated when using grid-averaged thickness. ► The annual-mean bias is about ...5Wm2 in the Arctic and 2Wm2 in the Southern Ocean. ► This bias becomes smaller when using a harmonic mean thickness.
We evaluated the impact of subgrid-scale ice thickness distribution on the heat flux on and through sea ice in a numerical model. An ice-ocean coupled model with a subgrid-scale ice thickness distribution scheme, COCO4.5, is forced by an atmospheric climatology to simulate the present state of the sea ice and ocean. The modeled climatology reproduces the ice cover reasonably well with a realistic ice thickness distribution.
The heat flux on and through the sea ice is established using the grid-averaged sea-ice and snow-on-ice thickness from the results of the simulation. When the grid-averaged thickness is calculated as a weighted arithmetic mean, the conductive heat flux through the ice and snow is underestimated compared with that actually driving the model. This underestimation becomes smaller in magnitude when either a weighted harmonic mean or a weighted arithmetic mean with a modification based on the ratio of these two types of means is used. Rearrangement of the ice categories shows that the flux bias decreases with an increase in the number of categories. We also perform a sensitivity experiment in which the model is forced by the biased heat flux identified using the arithmetic mean of the ice thickness. A significant decrease in ice volume is found, notably in the Arctic Ocean. These results suggest that sea-ice models without an ice thickness distribution scheme underestimate the conductive heat flux through ice, and thereby the resultant sea-ice thickness, because the ice thickness from these models typically corresponds to the weighted arithmetic mean thickness.
To assess the skill of seasonal to inter-annual
predictions of the detrended sea ice extent in the Arctic Ocean (SIEAO)
and to clarify the underlying physical processes, we conducted ensemble
...hindcasts, started on 1 January, 1 April, 1 July and 1 October for
each year from 1980 to 2011, for lead times up to three years, using the
Model for Interdisciplinary Research on Climate (MIROC) version 5
initialised with the observed atmosphere and ocean anomalies and sea ice
concentration. Significant skill is found for the winter months: the
December SIEAO can be predicted up to 11 months ahead (anomaly
correlation coefficient is 0.42). This skill might be attributed to the
subsurface ocean heat content originating in the North Atlantic. A plausible
mechanism is as follows: the subsurface water flows into the Barents Sea
from spring to fall and emerges at the surface in winter by vertical mixing,
and eventually affects the sea ice variability there. Meanwhile, the
September SIEAO predictions are skillful for lead times of up to
two months, due to the persistence of sea ice in the Beaufort, Chukchi, and East
Siberian seas initialised in July, as suggested by previous studies.
In line with the experimental design for near-term climate prediction toward the 5th Assessment Report of the Intergovernmental Panel on Climate Change (IPCC-AR5) and the Coupled Model ...Intercomparison Project Phase 5 (CMIP5), we perform ensembles of initialized decadal hindcast experiments using two recent versions of the Model for Interdisciplinary Research On Climate (MIROC): MIROC4h (T213L56 AGCM and 1/6-1/4 deg. 48 level OGCM) and MIROC5 (T85L40 AGCM and 0.56-1.4 deg. 50 level OGCM). We analyze sets of 10-yearlong 9-ensemble hindcasts (3 members by MIROC4h and 6 members by MIROC5) with initialization every five years after 1961 and explore the predictability of decadal climate changes. The most predictable variation on decadal timescales is the global warming signal due to the favorable response of the models to external forcing. The results of these initialized hindcast experiments using MIROC5 validate our ability to enhance decadal predictability primarily through the initialization, particularly of the Pacific Decadal Oscillation (PDO) for a few years and of the Atlantic Multidecadal Oscillation (AMO) for almost a decade. The initialization has large impacts on the upper ocean temperature hindcasts over the mid-and high latitudes of the North Pacific and the high latitudes of the North Atlantic, where the PDO and AMO signals are observed to be strongest. In contribution to process and assessment studies in IPCC-AR5 and CMIP5, further analysis of our hindcast data (and near-term prediction data) using MIROC4h and MIROC5 is worthwhile. We note that the initialized hindcasts using MIROC4h have predictive skill inferior to the MIROC5 results and that at this stage, fully significant discussions may not be possible due to the small number of ensembles with limited computational resources.
Byline: Yoshiki Komuro (1,2), Hiroyasu Hasumi (1) Results from an ice-ocean coupled model are used to investigate the impact of long-term variability in sea ice transport at the Fram Strait on the ...intensity of the Atlantic deep circulation. An increase (or decrease) in sea ice transport through the Fram Strait leads to a stronger (or weaker) deep circulation in the Atlantic. Change in the sea ice transport is accompanied by a salinity anomaly in the surface layer of the Arctic Ocean. Such an anomaly could inversely affect the Atlantic circulation once it reaches deep water formation regions. If the Canadian Archipelago is closed, the anomaly is subsequently transported through the Fram Strait, and counters the initial changes in the Atlantic deep circulation. On the other hand, if the Canadian Archipelago is open, some of the anomaly is transported to the Canadian Archipelago, and the initial change in the Atlantic deep circulation persists. In the Arctic Ocean basin, the time scale and path of the salinity anomalys propagation depends on the large-scale flow at the surface of the Arctic Ocean. Our results suggest that the salinity anomaly transport and its propagation pathway out of the Arctic Ocean are important determinants of the role of sea ice transport variability through the Fram Strait in controlling the intensity of the Atlantic deep circulation. Author Affiliation: (1) Center for Climate System Research, University of Tokyo, General Research Building, 5-1-5 Kashiwanoha, Kashiwa-shi, Chiba, 277-8568, Japan (2) Frontier Research Center for Global Change, Japan Agency for Marine-Earth Science and Technology, 3173-25 Showamachi, Kanazawa-ku, Yokohama City, Kanagawa, 236-0001, Japan Article History: Registration Date: 13/03/2007 Received Date: 19/10/2006 Accepted Date: 06/03/2007 Online Date: 11/04/2007
Sea ice has a large impact on climatic system and its variability. A good reproducibility of the past state of the sea ice in global climate models will reduce uncertainty in future projection. Here, ...we present sea-ice simulations for new versions of atmosphere-ocean coupled general circulation models, the Model for Interdisciplinary Research on Climate version 4h (MIROC4h) and version 5 (MIROC5), and assess the reproducibility of the sea ice prior to the future projection. The horizontal resolution of MIROC4h is significantly high for a coupled climate model, although its sea-ice component is based on the previous version. MIROC5 employs some improved schemes including subgrid-scale ice thickness distribution. Hindcast simulations of twentieth-century climate by the new models are compared with observations and with the results of previous versions of MIROC. For the Northern Hemisphere, Arctic sea-ice simulations are improved in both MIROC4h and MIROC5 compared with previous models. MIROC5 generally agrees well with observational data, whereas in MIROC4h, the Arctic sea ice is smaller in summer extent and in thickness. Employment of the ice thickness distribution, which allows large heat exchange through subgrid-scale thin ice regardless of the grid-averaged thickness, and relatively high albedo parameters contribute to reproduce more realistic ice thickness in MIROC5 compared with that in MIROC4h. For the Southern Hemisphere, MIROC4h well reproduces the observed ice edge, especially in winter, while MIROC5 underestimates sea-ice extent. Both models indicate decreasing trends in Arctic sea ice in the late twentieth century. A heat budget analysis of the MIROC5 Arctic Ocean suggests that intensification of ice-albedo feedback accelerates the rate of Arctic ice decline.
In line with the experimental design for near-term climate prediction toward the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, we perform ensembles of initialized decadal ...hindcast and near-future projection using three versions of the coupled atmosphere-ocean model MIROC. In the present study, we explore interannual and multiyear predictability of tropical cyclone (TC) activity in the western North Pacific (WNP) using the initialized hindcasts and examine global warming impacts on TC activity in the near-future on the basis of near-future projections up to 2035. The hindcasts of year-to-year variation in TC number capture the observed values reasonably well. Moreover, interannual variability of TC genesis and occurrence frequency associated with the El Niño Southern Oscillation are found to be predictable, mainly through better prediction of sea surface temperature (SST) and large-scale vorticity anomalies in the lower troposphere. These results indicate that the models can reproduce the major basic mechanisms that link TC genesis with large-scale circulation. Skillful prediction of TC number is likely difficult on multiyear timescales, at least based on our hindcasts, but through initializations, the three-year-mean hindcasts from 1998 onward reasonably capture observed major characteristics of TC activity associated with Pacific climate shift during the late 1990s. Near-future projections (2016-2035) suggest significant reductions (approximately 14%) in TC number, particularly over the western part of the WNP, even under scenarios in which projected global warming is less prominent than that at the end of this century. This reduction is likely due to the suppression of large-scale lower tropospheric vorticity and relative humidity and the enhancement of vertical wind shear. The projected SST exhibits a more pronounced warming over the eastern tropical Pacific than over the western region and accompanies the weakening of Walker circulation via redistribution of tropical convection activity, which appears to be responsible for the change in the large-scale fields in the WNP.