A regional sea ice-ocean coupled model for the Arctic Ocean was developed, based on the MITgcm ocean circulation model and classical Hibler79 type two categorythermodynamics-dynamics sea ice model. The sea ice dynamics and thermodynamicswere considered based on Viscous-Plastic (VP) and Winton three-layer models, respectively. A detailed configuration of coupled model has been introduced. Special attention has been paid to the model grid setup, subgrid paramerization, ice-ocean coupling and open boundary treatment. The coupled model was then applied and two test run examples were presented. The first model run was a climatology simulation with 10 years (1992?002) averaged NCAR/NCEP reanalysis data as atmospheric forcing. The second model run was a seasonal simulation for the period of 1992?007. The atmospheric forcing was daily NCAR/NCEP reanalysis. The climatology simulation captured the general pattern of the sea ice thickness distribution of the Arctic, i.e., the thickest sea ice is situated around the CanadaArchipelago and the north coast of the Greenland. For the second model run, themodeled September Sea ice extent anomaly from 1992?007 was highly correlated with the observations, with a linear correlation coefficient of 0.88. Theminimum of the Arctic sea ice area in the September of 2007 was unprecedented. The modeled sea ice area and extent for this minimum was overestimated relative to the observations. However, it captured the general pattern of the sea ice retreat.