Progress in the Development and Application of Climate Ocean Models and Ocean-Atmosphere Coupled Models in China

A review is presented about the development and application of climate ocean models and oceanatmosphere coupled models developed in China as well as a review of climate variability and climate change studies performed with these models. While the history of model development is briefly reviewed, emphasis has been put on the achievements made in the last five years. Advances in model development are described along with a summary on scientific issues addressed by using these models. The focus of the review is the climate ocean models and the associated coupled models, including both global and regional models, developed at the Institute of Atmospheric Physics, Chinese Academy of Sciences. The progress of either coupled model development made by other institutions or climate modeling using internationally developed models also is reviewed.

A Regional Air-Sea Coupled Model and Its Application over East Asia in the Summer of 2000

A regional air-sea coupled model, comprising the Regional Integrated Environment Model System （RIEMS） and the Princeton Ocean Model （POM） was developed to simulate summer climate features over East Asia in 2000. The sensitivity of the model＇s behavior to the coupling time interval （CTI）, the causes of the sea surface temperature （SST） biases, and the role of air-sea interaction in the simulation of precipitation over China are investigated. Results show that the coupled model can basically produce the spatial pattern of SST, precipitation, and surface air temperature （SAT） with five different CTIs respectively. Also, using a CTI of 3, 6 or 12 hours tended to produce more successful simulations than if using 1 and 24 hours. Further analysis indicates that both a higher and lower coupling frequency result in larger model biases in air-sea heat flux exchanges, which might be responsible for the sensitivity of the coupled model＇s behavior to the CTI. Sensitivity experiments indicate that SST biases between the coupled and uncoupled POM occurring over the China coastal waters were due to the mismatch of the surface heat fluxes produced by the RIEMS with those required by the POM. In the coupled run, the air-sea feedbacks reduced the biases in surface heat fluxes, compared with the uncoupled RIEMS, consequently resulted in changes in thermal contrast over land and sea and led to a precipitation increase over South China and a decrease over North China. These results agree well observations in the summer of 2000.

A review of progress in coupled ocean-atmosphere model developments for ENSO studies in China

El Nino-Southern Oscillation(ENSO) is the strongest interannual signal that is producedby basinscale processes in the tropical Pacific,with significant effects on weather and climate worldwide.In the past,extensive and intensive international efforts have been devoted to coupled model developments for ENSO studies.A hierarchy of coupled ocean-atmo sphere models has been formulated;in terms of their complexity,they can be categorized into intermediate coupled models(ICMs),hybrid coupled models(HCMs),and fully coupled general circulation models(CGCMs).ENSO modeling has made significant progress over the past decades,reaching a stage where coupled models can now be used to successfully predict ENSO events 6 months to one year in advance.Meanwhile,ENSO exhibits great diversity and complexity as observed in nature,which still cannot be adequately captured by current state-of-the-art coupled models,presenting a challenge to ENSO modeling.We primarily reviewed the long-term efforts in ENSO modeling continually and steadily made at different institutions in China;some selected representative examples are presented here to review the current status of ENSO model developments and applications,which have been actively pursued with noticeable progress being made recently.As ENSO simulations are very sensitive to model formulations and process representations etc.,dedicated efforts have been devoted to ENSO model developments and improvements.Now,different ocean-atmosphere coupled models have been available in China,which exhibit good model performances and have already had a variety of applications to climate modeling,including the Coupled Model Intercomparison Project Phase 6(CMIP6).Nevertheless,large biases and uncertainties still exist in ENSO simulations and predictions,and there are clear rooms for their improvements,which are still an active area of researches and applications.Here,model performances of ENSO simulations are assessed in terms of advantages and disadvantages with these differently formulated coupled models,pinpointing to the areas where they need to be further improved for ENSO studies.These analyses provide valuable guidance for future improvements in ENSO simulations and predictions.

Primary Reasoning behind the Double ITCZ Phenomenon in a Coupled Ocean-Atmosphere General Circulation Model

This paper investigates the processes behind the double ITCZ phenomenon, a common problem in Coupled ocean-atmosphere General Circulation Models (CGCMs), using a CGCM—FGCM-0 (Flexible General Circulation Model, version 0). The double ITCZ mode develops rapidly during the ?rst two years of the integration and becomes a perennial phenomenon afterwards in the model. By way of Singular Value Decomposition (SVD) for SST, sea surface pressure, and sea surface wind, some air-sea interactions are analyzed. These interactions prompt the anomalous signals that appear at the beginning of the coupling to develop rapidly. There are two possible reasons, proved by sensitivity experiments: (1) the overestimated east-west gradient of SST in the equatorial Paci?c in the ocean spin-up process, and (2) the underestimated amount of low-level stratus over the Peruvian coast in CCM3 (the Community Climate Model, Version Three). The overestimated east-west gradient of SST brings the anomalous equatorial easterly. The anomalous easterly, a?ected by the Coriolis force in the Southern Hemisphere, turns into an anomalous westerly in a broad area south of the equator and is enhanced by atmospheric anomalous circulation due to the underestimated amount of low-level stratus over the Peruvian coast simulated by CCM3. The anomalous westerly leads to anomalous warm advection that makes the SST warm in the southeast Paci?c. The double ITCZ phenomenon in the CGCM is a result of a series of nonlocal and nonlinear adjustment processes in the coupled system, which can be traced to the uncoupled models, oceanic component, and atmospheric component. The zonal gradient of the equatorial SST is too large in the ocean component and the amount of low-level stratus over the Peruvian coast is too low in the atmosphere component.

Seasonal and Intraseasonal Variations of East Asian Summer Monsoon Precipitation Simulated by a Regional Air-Sea Coupled Model

The performance of a regional air-sea coupled model, comprising the Regional Integrated Environment Model System （RIEMS） and the Princeton Ocean Model （POM）, in simulating the seasonal and intraseasonal variations of East Asian summer monsoon （EASM） rainfall was investigated. Through comparisons of the model results among the coupled model, the uncoupled RIEMS, and observations, the impact of air-sea coupling on simulating the EASM was also evaluated. Results showed that the regional air sea coupled climate model performed better in simulating the spatial pattern of the precipitation climatology and produced more realistic variations of the EASM rainfall in terms of its amplitude and principal EOF modes. The coupled model also showed greater skill than the uncoupled RIEMS in reproducing the principal features of climatological intraseasonal oscillation （CISO） of EASM rainfall, including its dominant period, intensity, and northward propagation. Further analysis indicated that the improvements in the simulation of the EASM rainfall climatology and its seasonal variation in the coupled model were due to better simulation of the western North Pacific Subtropical High, while the improvements of CISO simulation were owing to the realistic phase relationship between the intraseasonal convection and the underlying SST resulting from the air-sea coupling.

Introduction to the Regional Coupled Model WRF4-LICOM: Performance and Model Intercomparison over the Western North Pacific

Regional coupled modeling is one of the frontiers of regional climate modeling,but intercomparison has not been well coordinated.In this study,a community regional climate model,WRF4,with a resolution of 15 km,was coupled with a high-resolution(0.1°)North Pacific Ocean model(LICOM_np).The performance of the regional coupled model,WRF4_LICOM,was compared to that of another regional coupled model,RegCM4_LICOM,which was a coupling of version 4 of the Regional Climate Model(RegCM4)with LICOM_np.The analysis focused on the 2005 western North Pacific summer monsoon rainfall.The results showed that the regional coupled models with either RegCM4 or WRF4 as their atmospheric model component simulated the broad features over the WNP reasonably well.Quantitative intercomparison of the regional coupled simulations exhibited different biases for different climate variables.RegCM4_LICOM exhibited smaller biases in its simulation of the averaged June–July–August SST and rainfall,while WRF4_LICOM better captured the tropical cyclone(TC)intensity,the percentage contributions of rainfall induced by TCs to the total rainfall,and the diurnal cycle of rainfall and stratiform percentages,especially over land areas.The different behaviors in rainfall simulated by the two models were partly ascribed to the behaviors in the simulated western North Pacific subtropical high(WNPSH).The stronger(weaker)WNPSH in WRF4_LICOM(RegCM4_LICOM)was driven by overestimated(underestimated)diabatic heating,which peaked at approximately 450 hPa over the region around the Philippines in association with different condensation–radiation processes.Coupling of WRF4 with LIOCM is a crucial step towards the development of the next generation of regional earth system models at the Chinese Academy of Sciences.

Coupling of a Regional Climate Model with a Crop Development Model and Evaluation of the Coupled Model across China

In this study, the CERES(Crop Estimation through Resource and Environment Synthesis) crop model was coupled with CLM3.5, the land module of the regional climate model RegCM4. The new coupled model was named RegCM4_CERES; and in this model, crop type was further divided into winter wheat, spring wheat, spring maize, summer maize, early rice, late rice,single rice, and other crop types based on each distribution fraction. The development of each crop sub-type was simulated by the corresponding crop model separately, with each planting and harvesting date. A simulation test using RegCM4_CERES was conducted across China from 1999 to 2008; a control test was also performed using the original RegCM4. Data on crop LAI(leaf area index), soil moisture at 10 cm depth, precipitation, and 2 m air temperature were collected to evaluate the performance of RegCM4_CERES. The evaluation provided comparison of single-station time series, regional distributions,seasonal variations, and statistical indices for RegCM4_CERES. The results revealed that the coupled model had an excellent ability to simulate the phonological changes and spatial variations in crops. The consideration of dynamic crop development in RegCM4_CERES corrected the wet bias of the original RegCM4 over North China and the cold bias over South China.However, the degree of improvement was minimal and the statistical indices for RegCM4_CERES were roughly the same as the original RegCM4.

SIMULATION OF SUMMER CLIMATE IN CHINA DURING 1997 AND 1998 USING A REGIONAL AIR-SEA COUPLED MODEL

Using the regional air-sea coupled climate model RegCM3-POM,a series of numerical experiments are performed to simulate the summer climate in 1997 and 1998 with different coupling time steps.The results show that the coupled model has good performance on the simulation of the summer sea surface temperature(SST) in 1997 and 1998,and the simulation results of CPL1(with the coupling time step at 1 hour) are similar to those of CPL6(with the coupling time step at 6 hours).The coupled model can well simulate SST differences between 1997 and 1998.As for the simulation of the drought in 1997 and the flood in 1998,the results of CPL6 are more accurate.The coupled model can well simulate the drought in 1997 over North China,and compared with the results of the atmosphere model RegCM3,the simulation ability of the coupled model is improved.The coupling model has better ability in the simulation of the circulation in the middle and low levels,and the water vapor transportation in the coupling model is reasonable in both 1997 and 1998.RegCM3(an uncoupled model) cannot correctly simulate the transportation path differences between 1997 and 1998,but the coupled model can simulate the differences well.

The Interannual Variability of Climate in a Coupled Ocean-Atmosphere Model

In this paper, the interannual variability simulated by the coupled ocean-atmosphere general circulation modelof the institute of Atmospheric Physics (IAP CGCM) in 40 year integrations is analyzed, and compared with that bythe corresponding IAP AGCM which uses the climatic sea surface temperature as the boundary condition in 25 yearintegrations.The mean climatic states of January and July simulated by IAP CGCM are in good agreement with that by IAPAGCM, i.e., no serious 'climate drift' occurs in the CGCM simulation. A comparison of the results from AGCM andCGCM indicates that the standard deviation of the monthly averaged sea level pressure simulated by IAP CGCM ismuch greater than that by IAP AGCM in tropical region. In addition, both Southern Oscillation (SO) and NorthAtlantic Oscillation (NAO) can be found in the CGCM simulation for January, but these two oscillations do not existin the AGCM simulation.The interannual variability of climate may be classified into two typest one is the variation of the annual mean,another is the variation of the annual amplitude. The ocean-atmosphere interaction mainly increases the first type ofvariability. By means of the rotated EOF, the most important patterns corresponding to the two types of interannualvariability are found to have different spatial and temporal characteristics.

THE IMPACT OF INITIAL FORCED WIND ON THE PREDICTABILITY OF THE ZEBIAK-CANE COUPLED OCEAN-ATMOSPHERE MODEL

With simultaneous observed sea surface temperature anomaly （SSTA）, the difference between NCEP/NCAR 925hPa reanalysis wind stress anomaly （NCEPWSA） and FSU wind stress anomaly （FSUWSA） is analyzed, and the prediction abilities of Zebiak-Cane coupled ocean-atmosphere model （ZC coupled model） with NCEPWSA and FSUWSA serving respectively as initialization wind are compared. The results are as follows. The distribution feature of NCEPWSA matches better with that of the observed SSTA than counterpart of FSUWSA both in 1980s and in 1990s; The ZC ocean model has a better skill under the forcing of NCEPWSA than that of FSUWSA, especially in 1990s. Meanwhile, the forecast abilities of the ZC coupled model in 1990s as well as in 1980s have been improved employing NCEPWSA as initialization wind instead of FSUWSA. Particularly, it succeeded in predicting 1997/1998 E1 Nino 6 to 8 months ahead, further analysis shows that on the antecedent and onset stages of the 1997/1998 E1 Nino event, the horizontal cold and warm distribution characteristics of the simulated SSTA from ZC ocean model, with NCEPWSA forcing compared to FSUWSA forcing, match better with counterparts of the corresponding observed SSTA, whereby providing better predication initialization conditions for ZC coupled model, which, in turn, is favorable to improve the forecast ability of the coupled model.

A STUDY ON THE IMPACTS OF LATENT HEAT PARAMETERIZATION SCHEME ON PREDICTION SKILL OF ENSO WITH A SIMPLE OCEAN-ATMOSPHERE COUPLED MODEL

This study revises Weare＇s latent heat parameterization scheme and conducts an associated theoretic analysis. The revised Weare＇s scheme is found to present potentially better results than Zebiak＇s scheme. The Zebiak-Cane coupled ocean-atmosphere model, initialized by the National Centers for Environmental Prediction and the National Center for Atmospheric Research （NCEP/NCAR） reanalysis of wind stress anomaly at 925 hPa, is referred to as the ZCW coupled model. The atmosphere models of the ZCW coupled model that use Zebiak＇s scheme and the revised Weare＇s scheme are referred to as the ZCW0 and ZCWN atmosphere models, respectively. The coupled ocean-atmosphere models that use Zebiak＇s scheme and the revised Weare＇s scheme are referred to as the ZCWoand ZCWN coupled models, respectively. The simulations between the ZCW0 and ZCWN atmosphere models and between the ZCW0 and ZCWN coupled models are analyzed. The results include： （1） The evolution of heat, meridional wind and divergence anomalies simulated by similar ZCW0 and ZCWN atmosphere models, although the magnitudes of the former are larger than those of the latter; （2） The prediction skill of the Nino3 index from 1982 to 1999 by the ZCWN coupled model shows improvement compared with those by the ZCW0 coupled model; （3） The analysis of E1 Nino events in 1982/1983, 1986/1987, and 1997/1998 and La Nifia events in 1984/1985, 1988/1989, and 1998/2000 suggests that the ZCWN coupled model is better than the ZCW0 coupled model in predicting warm event evolution and cold event generation. The results also show the disadvantage of the ZCWN coupled model for predicting E1 Nino.

A Hybrid Coupled Ocean-Atmosphere Model and ENSO Prediction Study

A hybrid coupled ocean-atmosphere model is designed, which consists of a global AGCM and a simple anomaly ocean model in the tropical Pacific. Retroactive experimental predictions initiated in each season from 1979 to 1994 are performed. Analyses indicate that: (1) The overall predictive capability of this model for SSTA over the central-eastern tropical Pacific can reach one year, and the error is not larger than 0.8 degrees C. (2) The prediction skill depends greatly on the season when forecasts start. However, the phenomenon of SPB (spring prediction barrier) is not found in the model. (3) The ensemble forecast method can effectively improve prediction results. A new initialization scheme is discussed.

Influence of the North American Dipole on ENSO onset as simulated by a coupled ocean-Atmosphere model

The North American Dipole(NAD)is a north-south seesaw pattern of sea level pressure anomalies over the western tropical North Atlantic and northeastern North America.Previous observational studies have demonstrated that the NAD can affect the outbreak of El Niño-Southern Oscillation(ENSO)events.The present study analyzed the NAD-ENSO relationship as simulated by a coupled ocean-atmosphere model-namely,the Flexible Global Ocean-Atmosphere-Land System model,gridpoint version 2(FGOALS-g2).Results indicated that the model can replicate a distinct dipole comprised of a low over northeastern North America and a high over the western tropical North Atlantic,which is the signature feature of the NAD.Further analysis verified that the winter NAD can initiate the central equatorial Pacific warming in the subsequent winter by effectively forcing an anticyclonic flow and sea surface temperature(SST)warming over the northeastern subtropical Pacific(NESP)during late winter or early spring.In addition,the probability of an El Niño event was increased by a factor of 1.8 in the assimilation experiment with the NAD.By comparison,the winter Northern Atlantic Oscillation had no significant impact on the occurrence of ENSO a year later owing to its failure to induce the SST and surface wind anomalies over the NESP.

A Simple Quasi-Geostrophic Coupled Ocean-Atmosphere Model

The quasi-geostrophic atmospheric and oceanic equations of momentum and thermodynamics with dissipation factors are used to create a simple coupled ocean-atmosphere model describing the large-scale shallow-water motion. We discuss the ocean-atmosphere coupling effect in mid-high and low latitudes separately and analyze characteristics of which the oscillatory periods of coupled low-frequency modes (ocean mode) vary with the coupling frequency and latitudinal number. This can interpret the correlation between low-frequency oscillation and ocean-atmosphere interaction. Then from the dispersion curves of atmosphere and ocean, we reveal effect of the coupling strength on the propagation of Rossby waves. The convection mechanism between the two modes is also discussed in view of the slowly varying wave train.The results show that Newtonian cooling and Rayleigh friction play a stable rule in oceanic Rossby waves, the period of coupled low-frequency mode grows with the increment of the coupling frequency. The larger the latitudinal number is, the more rapidly it grows. When the coupling frequency tends to critical value, the oceanic Rossby waves become static. When the ocean-atmosphere coupling strength grows to some degree, the propagation of oceanic Rossby waves will become opposite to its original direction. One part of the oceanic Rossby waves is converted into atmospheric Rossby waves, the energy conversion coefficient is also solved out.

Projected Regional 1.50℃and 2.00℃Warming Threshold-crossing Time Worldwide Using the CMIP6 Models

The Paris Agreement aims to limit global warming to well below 2.00℃and pursue efforts to limit the temperature increase to 1.50℃.However,the response of climate change to unbalanced global warming is affected by spatial and temporal sensitivities.To better understand the regional warming response to global warming at 1.50℃and 2.00℃,we detected the 1.50℃and 2.00℃warming threshold-crossing time(WTT)above pre-industrial levels globally using the Coupled Model Intercomparison Project phase 6(CMIP6)models.Our findings indicate that the 1.50℃or 2.00℃WTT differs substantially worldwide.The warming rate of land would be approximately 1.35–1.46 times that of the ocean between 60°N–60°S in 2015–2100.Consequently,the land would experience a 1.50℃(2.00℃)warming at least 10–20 yr earlier than the time when the global mean near-surface air temperature reaches 1.50℃(2.00℃)WTT.Meanwhile,the Southern Ocean between 0°and 60°S considerably slows down the global 1.50℃and 2.00℃WTT.In 2040–2060,over 98.70%(77.50%),99.70%(89.30%),99.80%(93.40%),and 100.00%(98.00%)of the land will have warmed by over 1.50℃(2.00℃)under SSP(Shared Socioeconomic Pathway)1–2.6,SSP2-4.5,SSP3-7.0,and SSP5-8.5,respectively.We conclude that regional 1.50℃(2.00℃)WTT should be fully considered,especially in vulnerable high-latitude and high-altitude regions.

Does regional air–sea coupling improve the simulation of the summer monsoon over the western North Pacific in the WRF4 model?

A new regional coupled ocean–atmosphere model,WRF4-LICOM,was used to investigate the impacts of regional air–sea coupling on the simulation of the western North Pacific summer monsoon(WNPSM),with a focus on the normal WNPSM year 2005.Compared to WRF4,WRF4-LICOM improved the simulation of the summer mean monsoon rainfall,circulations,sea surface net heat fluxes,and propagations of the daily rainband over the WNP.The major differences between the models were found over the northern South China Sea and east of the Philippines.The warmer SST reduced the gross moist stability of the atmosphere and increased the upward latent heat flux,and then drove local ascending anomalies,which led to the increase of rainfall in WRF4-LICOM.The resultant enhanced atmospheric heating drove a low-level anomalous cyclone to its northwest,which reduced the simulated circulation biases in the stand-alone WRF4 model.The local observed daily SST over the WNP was a response to the overlying summer monsoon.In the WRF4 model,the modeled atmosphere exhibited passive response to the underlying daily SST anomalies.With the inclusion of regional air–sea coupling,the simulated daily SST–rainfall relationship was significantly improved.WRF4-LICOM is recommended for future dynamical downscaling of simulations and projections over this region.

Atmospheric and oceanic responses to Super Typhoon Mangkhut in the South China Sea:a coupled CROCO-WRF simulation

The South China Sea(SCS)is the largest marginal sea in the Northwest Pacific Ocean,and it encounters frequent typhoons.The atmosphere and ocean will create significant thermal and dynamic responses during the intense disturbance caused by typhoons.However,these responses have not been thoroughly investigated owing to the complicated marine environment.According to the satellite data,the SCS Basin was observed to have a strong sea surface temperature(SST)response to Typhoon Mangkhut,resulting in widespread SST cooling.A coupled model was used to investigate the atmospheric and oceanic responses to Typhoon Mangkhut.Best-track data,satellite SST,and ARGO measurements show that the coupled WRF-CROCO simulation displays better track,intensity,SST,temperature,and salinity profiles than those of the WRF-only simulation.Results show that the typhoon induced rightward intensifications in wind speed,ocean current,and SST.The following are some remarkable atmosphere and ocean responses:(1)the SST below the inner-core region is cooled by 1℃,resulting in a 37%-44%decrease in wet enthalpy,and the central pressure is increased by~9 hPa.Therefore,the changes in SST below the innercore region of the SCS Basin have a significant impact on air-sea fluxes under high-wind conditions;(2)the ocean boundary layer analysis shows that near-inertial oscillations on the right side of the typhoon track and a strong inertial current up to~2.28 m/s in the upper ocean were observed,which resonated with the local wind and flow field on the right side and induced strong SST cooling;(3)a decrease in SST decreased the moist static energy of the typhoon boundary layer,thereby weakening the typhoon’s intensity.The difference in equivalent potential temperature and sea surface pressure have a good correlation,indicating that the influence of moist static energy on typhoon intensity cannot be overlooked.

A Hybrid Coupled Model for the Pacific Ocean–Atmosphere System.Part I: Description and Basic Performance

A hybrid coupled model （HCM） is constructed for El Nifio-Southern Oscillation （ENSO）-related modeling studies over almost the entire Pacific basin.An ocean general circulation model is coupled to a statistical atmospheric model for interannual wind stress anomalies to represent their dominant coupling with sea surface temperatures.In addition,various relevant forcing and feedback processes exist in the region and can affect ENSO in a significant way; their effects are simply represented using historical data and are incorporated into the HCM,including stochastic forcing of atmospheric winds,and feedbacks associated with freshwater flux,ocean biology-induced heating （OBH）,and tropical instability waves （TIWs）.In addition to its computational efficiency,the advantages of making use of such an HCM enable these related forcing and feedback processes to be represented individually or collectively,allowing their modulating effects on ENSO to be examined in a clean and clear way.In this paper,examples are given to illustrate the ability of the HCM to depict the mean ocean state,the circulation pathways connecting the subtropics and tropics in the western Pacific,and interannual variability associated with ENSO.As satellite data are taken to parameterize processes that are not explicitly represented in the HCM,this work also demonstrates an innovative method of using remotely sensed data for climate modeling.Further model applications related with ENSO modulations by extratropical influences and by various forcings and feedbacks will be presented in Part Ⅱ of this study.

Does Industry-University-Research Cooperation Matter? An Analysis of Its Coupling Effect on Regional Innovation and Economic Development

The dislocation between regional innovation and economic development directly influences the economic effect of regional innovation. However, no in-depth researches have been made on how to solve this problem. Using data from Henan Province, China,employing geographical detector technology, this paper focuses on testing whether the industry-university-research cooperation can contribute to coordinating the relation between regional innovation and economic development. It is shown that: 1) the industry-universityresearch cooperation in Henan Province is increasing gradually, and the network presents a core-edge structure, and the coupling degree between regional innovation and economic development is spatially unbalanced, which is similar to the spatial distribution of the intensity of industry-university-research cooperation;2) as an important approach to effectively connect scientific researches with market demands, the industry-university-research cooperation can help form an interactive, interconnected, coupled and coordinated virtuous relation between regional innovation and economic development. Compared with the cooperation between organizations of the same type and the separate innovation of organizations, the improvement of the industry-university-research cooperation level can better coordinate the relation between regional innovation and economic development;3) the cooperative innovation model between enterprises and universities can better promote the coupling between regional innovation and economic development, compared with many industryuniversity-research cooperation models. For underdeveloped areas lacking local knowledge base, industry-university-research cooperation should be considered as a long-term development strategy, especially using the knowledge sources of external universities and scientific research institutions to enhance innovation capability and achieve economic growth.

Regional earth system modeling：review and future directions

The authors review recent advances in the development of coupled Regional Earth System Models （RESMs）,a field that is still in its early stages.To date,coupled regional atmosphere-ocean-sea ice,atmosphere-aerosol and atmosphere-biosphere models have been developed,but they have been applied onlyto limited regional settings.Much more work is thus needed to assess their transferability to a wide range of settings.Future challenges in regional climate modeling are identified,including the development of fully coupled RESMs encompassing not only atmosphere,ocean,cryosphere,biosphere,chemosphere,but also the human component in a fully interactive way.