A statistical�Cdynamical scheme for reconstructing ocean forcing in the Atlantic. Part II: methodology, validation and application to high-resolution ocean models

A novel statistical?Cdynamical scheme has been developed to reconstruct the sea surface atmospheric variables necessary to force an ocean model. Multiple linear regressions are first built over a so-called learning period and over the entire Atlantic basin from the observed relationship between the surface wind conditions, or predictands, and the anomalous large scale atmospheric circulations, or predictors. The latter are estimated in the extratropics by 500?hPa geopotential height weather regimes and in the tropics by low-level wind classes. The transfer function further combined to an analog step is then used to reconstruct all the surface variables fields over 1958?C2002. We show that the proposed hybrid scheme is very skillful in reproducing the mean state, the seasonal cycle and the temporal evolution of all the surface ocean variables at interannual timescale. Deficiencies are found in the level of variance especially in the tropics. It is underestimated for 2-m temperature and humidity as well as for surface radiative fluxes in the interannual frequency band while it is slightly overestimated at higher frequency. Decomposition in empirical orthogonal function (EOF) shows that the spatial and temporal coherence of the forcing fields is however very well captured by the reconstruction method. For dynamical downscaling purposes, reconstructed fields are then interpolated and used to carry out a high-resolution oceanic simulation using the NATL4 (1/4°) model integrated over 1979?C2001. This simulation is compared to a reference experiment where the original observed forcing fields are prescribed instead. Mean states between the two experiments are virtually undistinguishable both in terms of surface fluxes and ocean dynamics estimated by the barotropic and the meridional overturning streamfunctions. The 3-dimensional variance of the simulated ocean is well preserved at interannual timescale both for temperature and salinity except in the tropics where it is underestimated. The main modes of interannual variability assessed through EOF are correctly reproduced for sea surface temperature, barotropic streamfunction and mixed layer depth both in terms of spatial structure and temporal evolution. Collectively, our results provide evidence that the statistical?Cdynamical scheme presented in this two-part study is an efficient and promising tool to infer oceanic changes (in particular those related to the wind-driven circulation) due to modifications in the large-scale atmospheric circulation. As a prerequisite, we have here validated the method for present-day climate; we encourage its use for climate change studies with some adaptations though.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere–ocean global coupled model. Part I: a mechanism governing the variability of ocean convection in a preindustrial experiment

A preindustrial climate experiment was conducted with the third version of the CNRM global atmosphere–ocean–sea ice coupled model (CNRM-CM3) for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). This experiment is used to investigate the main physical processes involved in the variability of the North Atlantic ocean convection and the induced variability of the Atlantic meridional overturning circulation (MOC). Three ocean convection sites are simulated, in the Labrador, Irminger and Greenland–Iceland–Norwegian (GIN) Seas in agreement with observations. A mechanism linking the variability of the Arctic sea ice cover and convection in the GIN Seas is highlighted. Contrary to previous suggested mechanisms, in CNRM-CM3 the latter is not modulated by the variability of freshwater export through Fram Strait. Instead, the variability of convection is mainly driven by the variability of the sea ice edge position in the Greenland Sea. In this area, the surface freshwater balance is dominated by the freshwater input due to the melting of sea ice. The ice edge position is modulated either by northwestward geostrophic current anomalies or by an intensification of northerly winds. In the model, stronger than average northerly winds force simultaneous intense convective events in the Irminger and GIN Seas. Convection interacts with the thermohaline circulation on timescales of 5–10 years, which translates into MOC anomalies propagating southward from the convection sites.     

Potential role of Atlantic Warm Pool-induced freshwater forcing in the Atlantic Meridional Overturning Circulation: ocean–sea ice model simulations

Recent studies have indicated that the multidecadal variations of the Atlantic Warm Pool (AWP) can induce a significant freshwater change in the tropical North Atlantic Ocean. In this paper, the potential effect of the AWP-induced freshwater flux on the Atlantic Meridional Overturning Circulation (AMOC) is studied by performing a series of ocean–sea ice model experiments. Our model experiments demonstrate that ocean response to the anomalous AWP-induced freshwater flux is primarily dominated by the basin-scale gyre circulation adjustments with a time scale of about two decades. The positive (negative) freshwater anomaly leads to an anticyclonic (cyclonic) circulation overlapping the subtropical gyre. This strengthens (weakens) the Gulf Stream and the recirculation in the interior ocean, thus increases warm (cold) water advection to the north and decreases cold (warm) water advection to the south, producing an upper ocean temperature dipole in the midlatitude. As the freshwater (salty water) is advected to the North Atlantic deep convection region, the AMOC and its associated northward heat transport gradually decreases (increases), which in turn lead to an inter-hemispheric SST seesaw. In the equilibrium state, a comma-shaped SST anomaly pattern develops in the extratropical region, with the largest amplitude over the subpolar region and an extension along the east side of the basin and into the subtropical North Atlantic. Based on our model experiments, we argue that the multidecadal AWP-induced freshwater flux can affect the AMOC, which plays a negative feedback role that acts to recover the AMOC after it is weakened or strengthened. The sensitivity of AMOC response to the AWP-induced freshwater forcing amplitude is also examined and discussed.     

Impact of resolving the diurnal cycle in an ocean–atmosphere GCM. Part 2: A diurnally coupled CGCM

Coupled ocean atmosphere general circulation models (GCM) are typically coupled once every 24 h, excluding the diurnal cycle from the upper ocean. Previous studies attempting to examine the role of the diurnal cycle of the upper ocean and particularly of diurnal SST variability have used models unable to resolve the processes of interest. In part 1 of this study a high vertical resolution ocean GCM configuration with modified physics was developed that could resolve the diurnal cycle in the upper ocean. In this study it is coupled every 3 h to atmospheric GCM to examine the sensitivity of the mean climate simulation and aspects of its variability to the inclusion of diurnal ocean-atmosphere coupling. The inclusion of the diurnal cycle leads to a tropics wide increase in mean sea surface temperature (SST), with the strongest signal being across the equatorial Pacific where the warming increases from 0.2°C in the central and western Pacific to over 0.3°C in the eastern equatorial Pacific. Much of this warming is shown to be a direct consequence of the rectification of daily mean SST by the diurnal variability of SST. The warming of the equatorial Pacific leads to a redistribution of precipitation from the Inter tropical convergence zone (ITCZ) toward the equator. In the western Pacific there is an increase in precipitation between Papa new guinea and 170°E of up to 1.2 mm/day, improving the simulation compared to climatology. Pacific sub tropical cells are increased in strength by about 10%, in line with results of part 1 of this study, due to the modification of the exchange of momentum between the equatorially divergent Ekman currents and the geostropic convergence at depth, effectively increasing the dynamical response of the tropical Pacific to zonal wind stresses. During the spring relaxation of the Pacific trade winds, a large diurnal cycle of SST increases the seasonal warming of the equatorial Pacific. When the trade winds then re-intensify, the increase in the dynamical response of the ocean leads to a stronger equatorial upwelling. These two processes both lead to stronger seasonal basin scale feedbacks in the coupled system, increasing the strength of the seasonal cycle of the tropical Pacific sector by around 10%. This means that the diurnal cycle in the upper ocean plays a part in the coupled feedbacks between ocean and atmosphere that maintain the basic state and the timing of the seasonal cycle of SST and trade winds in the tropical Pacific. The Madden–Julian Oscillation (MJO) is examined by use of a large scale MJO index, lag correlations and composites of events. The inclusion of the diurnal cycle leads to a reduction in overall MJO activity. Precipitation composites show that the MJO is stronger and more coherent when the diurnal cycle of coupling is resolved, with the propagation and different phases being far more distinct both locally and to larger lead times across the tropical Indo-Pacific. Part one of this study showed that that diurnal variability of SST is modulated by the MJO and therefore increases the intraseasonal SST response to the different phases of the MJO. Precipitation-based composites of SST variability confirm this increase in the coupled simulations. It is argued that including this has increased the thermodynamical coupling of the ocean and atmosphere on the timescale of the MJO (20–100 days), accounting for the improvement in the MJO strength and coherency seen in composites of precipitation and SST. These results show that the diurnal cycle of ocean–atmosphere interaction has profound impact on a range of up-scale variability in the tropical climate and as such, it is an important feature of the modelled climate system which is currently either neglected or poorly resolved in state of the art coupled models.     

On the mechanisms in a tropical ocean–global atmosphere coupled general circulation model. Part I: mean state and the seasonal cycle

 The mechanisms responsible for the mean state and the seasonal and interannual variations of the coupled tropical Pacific-global atmosphere system are investigated by analyzing a thirty year simulation, where the LMD global atmospheric model and the LODYC tropical Pacific model are coupled using the delocalized physics method. No flux correction is needed over the tropical region. The coupled model reaches its regime state roughly after one year of integration in spite of the fact that the ocean is initialized from rest. Departures from the mean state are characterized by oscillations with dominant periodicites at annual, biennial and quadriennial time scales. In our model, equatorial sea surface temperature and wind stress fluctuations evolved in phase. In the Central Pacific during boreal autumn, the sea surface temperature is cold, the wind stress is strong, and the Inter Tropical Convergence Zone (ITCZ) is shifted northwards. The northward shift of the ITCZ enhances atmospheric and oceanic subsidence between the equator and the latitude of organized convention. In turn, the stronger oceanic subsidence reinforces equatorward convergence of water masses at the thermocline depth which, being not balanced by equatorial upwelling, deepens the equatorial thermocline. An equivalent view is that the deepening of the thermocline proceeds from the weakening of the meridional draining of near-surface equatorial waters. The inverse picture prevails during spring, when the equatorial sea surface temperatures are warm. Thus temperature anomalies tend to appear at the thermocline level, in phase opposition to the surface conditions. These subsurface temperature fluctuations propagate from the Central Pacific eastwards along the thermocline; when reaching the surface in the Eastern Pacific, they trigger the reversal of sea surface temperature anomalies. The whole oscillation is synchronized by the apparent meridional motion of the sun, through the seasonal oscillation of the ITCZ. This possible mechanism is partly supported by the observed seasonal reversal of vorticity between the equator and the ITCZ, and by observational evidence of eastward propagating subsurface temperature anomalies at the thermocline level. Received: 7 April 1997 / Accepted: 15 July 1998     

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 Nino–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 II of this study.     

Parameterization of Sheared Entrainment in a Well-developed CBL.Part Ⅱ:A Simple Model for Predicting the Growth Rate of the CBL

Following the parameterization of sheared entrainment obtained in the companion paper, Liu et al. (2016), the present study aims to further investigate the characteristics of entrainment, and develop a simple model for predicting the growth rate of a well-developed and sheared CBL. The relative stratification, defined as the ratio of the stratification in the free atmosphere to that in the entrainment zone, is found to be a function of entrainment flux ratio (A _e). This leads to a simple expression of the entrainment rate, in which A _e needs to be parameterized. According to the results in Liu et al. (2016), A _e can be simply expressed as the ratio of the convective velocity scale in the sheared CBL to that in the shear-free CBL. The parameterization of the convective velocity scale in the sheared CBL is obtained by analytically solving the bulk model with several assumptions and approximations. Results indicate that the entrainment process is influenced by the dynamic effect, the interaction between mean shear and environmental stratification, and one other term that includes the Coriolis effect. These three parameterizations constitute a simple model for predicting the growth rate of a well-developed and sheared CBL. This model is validated by outputs of LESs, and the results show that it performs satisfactorily. Compared with bulk models, this model does not need to solve a set of equations for the CBL. It is more convenient to apply in numerical models.     

Flow and transport modeling in the sea-breeze. Part I: A modifiedE — ∈ model with a non-equilibrium level 2.5 closure

In order to overcome shortcomings of the conventionalE — model and to develop a more general model applicable to the variety of atmospheric conditions observed, for example, during a land-sea breeze cycle, a modifiedE — model is proposed. The model is a simplified form of an algebraic stress model including wall proximity effects of Gibson and Launder (1978). The proposed model is similar to the popular Mellor and Yamada (1982) level 2.5 model but does not employ a local equilibrium assumption in the algebraic equations for Reynolds stresses. The resulting model shows a wider realization region under unstable conditions than the Mellor and Yamada model.The modified model is compared herein to observations, higher order closure simulations and large eddy simulations under neutral, stable and convective conditions. Various dissipation rate equations were employed and compared to understand their performance with the modified model. The modifiedE — model reproduced the observed behavior well under all conditions except near the base of an elevated inversion layer under convective conditions. The ability of the model to describe flow dynamics under a wide range of atmospheric stabilities suggests that the model can be used to describe the complicated diurnal behavior of the land-sea breeze circulation.     

Dynamical and Thermal Problems in Vortex Development and Movement. Part I: A PV–Q View

Based on the Lagrangian change equation of vertical vorticity deduced from the equation of threedimensional Ertel potential vorticity(PV e),the development and movement of vortex are investigated from the view of potential vorticity and diabatic heating(PV-Q).It is demonstrated that the asymmetric distribution in the vortex of the non-uniform diabatic heating in both vertical and horizontal can lead to the vortex’s development and movement.The theoretical results are used to analyze the development and movement of a Tibetan Plateau(TP) vortex(TPV),which appeared over the TP,then slid down and moved eastward in late July 2008,resulting in heavy rainfall in Sichuan Province and along the middle and lower reaches of the Yangtze River.The relative contributions to the vertical vorticity development of the TPV are decomposed into three parts:the diabatic heating,the change in horizontal component of PV e(defined as PV 2),and the change in static stability θ z.The results show that in most cases,diabatic heating plays a leading role,followed by the change in PV 2,while the change of θ z usually has a negative impact in a stable atmosphere when the atmosphere becomes more stable,and has a positive contribution when the atmosphere approaches neutral stratification.The intensification of the TPV from 0600 to 1200 UTC 22 July 2008 is mainly due to the diabatic heating associated with the precipitation on the eastern side of the TPV when it uplifted on the up-slope of the northeastern edge of the Sichuan basin.The vertical gradient of diabatic heating makes positive(negative) PV e generation below(above) the maximum of diabatic heating;the positive PV e generation not only intensifies the low-level vortex but also enhances the vertical extent of the vortex as it uplifts.The change in PV e due to the horizontal gradient of diabatic heating depends on the vertical shear of horizontal wind that passes through the center of diabatic heating.The horizontal gradient of diabatic heating makes positive(negative) PV e generation on the right(left) side of the vertical shear of horizontal wind.The positive PV e generation on the right side of the vertical shear of horizontal wind not only intensifies the local vertical vorticity but also affects direction of movement of the TPV.These diagnostic results are in good agreement with the theoretic results developed from the PV-Q view.     

Proposed scheme for modeling of ocean equatorial currents in the phase of El Niño and La Niña: Implementation of the mesoscale turbulence theory

El Niño is a phenomenon of the catastrophic increase of surface temperature in the eastern part of the Pacific Ocean. It has a significant impact to weather of the American continent and western regions of the tropical Pacific, as well as on the weather and climate of entirely the Earth. Most important factors influencing El Niño are the wind, ocean currents and slope of the water surface (and temperature resulting from these factors) at the equator in the Pacific Ocean. The paper considers results of mathematical modeling of the equatorial Pacific Ocean currents in the El Niño and La Niña phases using the theory of mesoscale turbulence. This theory has been successfully tested in modeling of global circulation of atmosphere and ocean (Arsen’yev et al., 2010) and it has been able to calculate the ocean current changes at equator under changing external conditions. It is shown that the water currents at the equator have a four-tier vertical structure. The surface trade-wind current is located above the subsurface undercurrent, below which we observe the intermediate current, turning into the equatorial deep counter flow. When El Niño begins, the currents are rearranged, change signs and sometimes merge with each other. In the phase of maximum development of the phenomenon there is a two-tier structure: (1) surface current heading the American coast is underlain (below the depth of 440 m) by (2) deep equatorial current directed to the Indonesian coast. The theoretical calculations are compared with the physical observations of ocean currents in the El Niño and La Niña phases. The obtained results indicate that the proposed mathematical apparatus makes it possible to explain the set of physical observations in the Pacific Ocean.     

On the mechanisms in a tropical ocean–global atmosphere coupled general circulation model. Part II: interannual variability and its relation to the seasonal cycle

 The thirty year simulation of the coupled global atmosphere-tropical Pacific Ocean general circulation model of the Laboratoire de Métérologie Dynamique and the Laboratoire d’Océanographie Dynamique et de Climatologie presented in Part I is further investigated in order to understand the mechanisms of interannual variability. The model does simulate interannual events with ENSO characteristics; the dominant periodicity is quasi-biennial, though strong events are separated by four year intervals. The mechanism that is responsible for seasonal oscillations, identified in Part I, is also active in interannual variability with the difference that now the Western Pacific is dynamically involved. A warm interannual phase is associated with an equatorward shift of the ITCZ in the Western and Central Pacific. The coupling between the ITCZ and the ocean circulation is then responsible for the cooling of the equatorial subsurface by the draining mechanism. Cold subsurface temperature anomalies then propagate eastward along the mean equatorial thermocline. Upon reaching the Eastern Pacific where the mean thermocline is shallow, cold subsurface anomalies affect surface temperatures and reverse the phase of the oscillation. The preferred season for efficient eastward propagation of thermocline depth temperature anomalies is boreal autumn, when draining of equatorial waters towards higher latitudes is weaker than in spring by a factor of six. In that way, the annual cycle acts as a dam that synchronizes lower frequency oscillations. Received: 7 April 1997 / Accepted: 15 July 1998     

The Role of Land–sea Distribution and Orography in the Asian Monsoon. Part I: Land–sea Distribution

A number of AGCM simulations were performed by including various land–sea distributions (LSDs), such as meridional LSDs, zonal LSDs, tropical large-scale LSDs, and subcontinental-scale LSDs, to identify their effects on the Asian monsoon. In seven meridional LSD experiments with the continent/ocean located to the north/south of a certain latitude, the LSDs remain identical except the southern coastline is varied from 40 ° to 4 ° N in intervals of 5.6° . In the experiments with the coastline located to the n...     

A Study on the Mesoscale Convective Systems during the Summer Monsoon Onset over the South China Sea in 1998Part I: Analysis of Large-Scale Fields for Occurrence and Development of the Mesoscale Convective Systems

1. Introduction China is located in the East Asian monsoon re- gion. Every year's weather and climate in this region is deeply affected by the monsoon activities. Es- pecially, during flooding season (May to September), the summer monsoon controls large-scale precipitation patterns, the movement of seasonal rain belt and oc- currence of drought/flood disasters. The Asian mon- soon can be divided into two systems (Tao and Chen, 1987). As a major component and its unique location, the South …     

Comparison Between Large-Eddy Simulation and Reynolds-Averaged Navier–Stokes Computations for the MUST Field Experiment. Part I: Study of the Flow for an Incident Wind Directed Perpendicularly to the Front Array of Containers

The large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) methodologies are used to simulate the air flow inside the container’s array geometry of the Mock Urban Setting Test (MUST) field experiment. Both tools are assessed and compared in a configuration for which the incident wind direction is perpendicular to the front array. The assessment is carried out against available wind-tunnel data. Effects of including small geometrical irregularities present in the experiments are analysed by considering LES and RANS calculations on two geometries: an idealized one with a perfect alignment and an identical shape of the containers, and a second one including the small irregularities considered in the experiment. These effects are assessed in terms of the local time-mean average and as well in terms of spatial average properties (relevant in atmospheric modelling) given for the velocity and turbulent fields. The structural flow properties obtained using LES and RANS are also compared. The inclusion of geometrical irregularities is found significant on the local time-mean flow properties, in particular the repeated flow patterns encountered in a perfect regular geometry is broken. LES and RANS provide close results for the local mean streamwise velocity profiles and shear-stress profiles, however the LES predictions are closer to the experimental values for the local vertical mean velocity. When considering the spatial average flow properties, the effects of geometrical irregularities are found insignificant and LES and RANS provide similar results.     

On the Use of the Jesuit Order Private Correspondence Records in Climate Reconstructions: A Case Study from Castille (Spain) for 1634–1648 A.D.

This work analyses the climatic information content of more than 1000 letters belonging to the private correspondence of the Jesuit order during the period 1634–1648. The information obtained mainly originates in Castille, and it was codified to obtain quantitative indices representative of the evolution of the thermal and rainfall anomalies. The results show that intense rainfall and cold air waves predominated during the study period, in accordance with other results on the climatic evolution of the Iberian Peninsula during the Little Ice Age.     

A Study on the Mesoscale Convective Systems during the Summer Monsoon Onset over the South China Sea in 1998 Part Ⅰ: Analysis of Large-Scale Fields for Occurrence and Development of the Mesoscale Convective Systems

The summer monsoon onset over the northern South China Sea (SCS) in May 16-20, 1998 was characterized by the abrupt onset of mesoscale convective activities and rapid increase of precipitation. The possible mechanism for formation of the mesoscale convective systems (MCSs) and related rain belts were revealed through discussing their forming physical conditions under the large-scale background: (1) The high pseudo-equivalent potential temperature and the convective instability in the lower troposphere, the low-level southwesterly confluence and the high-level divergence over South China and the northern SCS provided the favorable large-scale thermodynamic and dynamic conditions for development of MCSs. The southwesterly flow from the Bay of Bengal (BOB) interacted with that to the western flank of the subtropical high, which constituted the major moisture channels, thus bringing about deep wet layers and strong moisture convergence;(2) triggered by several cold troughs from high and mid latitudes, the convectively unstable energy was released and the convective activities over the northern SCS broke out abruptly;(3)analysis of retrieved precipitation based on the dual-Doppler radar during South China Sea Monsoon Experiment (SCSMEX) indicated that active convection influenced by the monsoon trough and corresponding wind shear line organized and formed continually some mesoscale convective rainbelts. During May 15-19,about 12 precipitation processes with 6-12-hour life span or more were observed;and (4) under the favorable synoptic conditions, establishment of the monsoon trough and shear line in the low levels, as well as production and development of mesoscale low vortex were all necessary conditions for the formation and maintenance of MCSs.