Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies.

Late Quaternary vegetation history of North Stradbroke Island, Queensland, eastern Australia

Currently there is a paucity of records of late Quaternary palaeoenvironmental variability available from the subtropics of Australia. The three continuous palaeoecological records presented here, from North Stradbroke Island, subtropical Queensland, assist in bridging this large spatial gap in the current state of knowledge. The dominance of arboreal taxa in the pollen records throughout the past >40,000 years is in contrast with the majority of records from temperate Australia, and indicates a positive moisture balance for North Stradbroke Island. The charcoal records show considerable inter-site variability indicating the importance of local-scale events on individual records, and highlighting the caution that needs to be applied when interpreting a single site as a regional record. The variability in the burning regimes is interpreted as being influenced by both climatic and human factors. Despite this inter-site variability, broad environmental trends are identifiable, with changes in the three records comparable with the OZ-INTIMATE climate synthesis for the last 35,000 years.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tides, extra-tropical storm surges and mean sea level

The occurrence of extreme water levels along low-lying, highly populated and/or developed coastlines can lead to considerable loss of life and billions of dollars of damage to coastal infrastructure. Therefore it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood management, engineering and future land-use planning. This ensures the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. This paper estimates for the first time present day extreme water level exceedence probabilities around the whole coastline of Australia. A high-resolution depth averaged hydrodynamic model has been configured for the Australian continental shelf region and has been forced with tidal levels from a global tidal model and meteorological fields from a global reanalysis to generate a 61-year hindcast of water levels. Output from this model has been successfully validated against measurements from 30 tide gauge sites. At each numeric coastal grid point, extreme value distributions have been fitted to the derived time series of annual maxima and the several largest water levels each year to estimate exceedence probabilities. This provides a reliable estimate of water level probabilities around southern Australia; a region mainly impacted by extra-tropical cyclones. However, as the meteorological forcing used only weakly includes the effects of tropical cyclones, extreme water level probabilities are underestimated around the western, northern and north-eastern Australian coastline. In a companion paper we build on the work presented here and more accurately include tropical cyclone-induced surges in the estimation of extreme water level. The multi-decadal hindcast generated here has been used primarily to estimate extreme water level exceedance probabilities but could be used more widely in the future for a variety of other research and practical applications.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tropical cyclone-induced storm surges

The incidence of major storm surges in the last decade have dramatically emphasized the immense destructive capabilities of extreme water level events, particularly when driven by severe tropical cyclones. Given this risk, it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood and erosion management, engineering and for future land-use planning and to ensure the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. Australia has a long history of coastal flooding from tropical cyclones. Using a novel integration of two modeling techniques, this paper provides the first estimates of present day extreme water level exceedance probabilities around the whole coastline of Australia, and the first estimates that combine the influence of astronomical tides, storm surges generated by both extra-tropical and tropical cyclones, and seasonal and inter-annual variations in mean sea level. Initially, an analysis of tide gauge records has been used to assess the characteristics of tropical cyclone-induced surges around Australia. However, given the dearth (temporal and spatial) of information around much of the coastline, and therefore the inability of these gauge records to adequately describe the regional climatology, an observationally based stochastic tropical cyclone model has been developed to synthetically extend the tropical cyclone record to 10,000 years. Wind and pressure fields derived for these synthetically generated events have then been used to drive a hydrodynamic model of the Australian continental shelf region with annual maximum water levels extracted to estimate exceedance probabilities around the coastline. To validate this methodology, selected historic storm surge events have been simulated and resultant storm surges compared with gauge records. Tropical cyclone induced exceedance probabilities have been combined with estimates derived from a 61-year water level hindcast described in a companion paper to give a single estimate of present day extreme water level probabilities around the whole coastline of Australia. Results of this work are freely available to coastal engineers, managers and researchers via a web-based tool (www.sealevelrise.info). The described methodology could be applied to other regions of the world, like the US east coast, that are subject to both extra-tropical and tropical cyclones.

Indian Ocean sea surface salinity variations in a coupled model

The variability of the sea surface salinity (SSS) in the Indian Ocean is studied using a 100-year control simulation of the Community Climate System Model (CCSM 2.0). The monsoon-driven seasonal SSS pattern in the Indian Ocean, marked by low salinity in the east and high salinity in the west, is captured by the model. The model overestimates runoff int the Bay of Bengal due to higher rainfall over the Himalayan-Tibetan regions which drain into the Bay of Bengal through Ganga-Brahmaputra rivers. The outflow of low-salinity water from the Bay of Bengal is to strong in the model. Consequently, the model Indian Ocean SSS is about 1 less than that seen in the climatology. The seasonal Indian Ocean salt balance obtained from the model is consistent with the analysis from climatological data sets. During summer, the large freshwater input into the Bay of Bengal and its redistribution decide the spatial pattern of salinity tendency. During winter, horizontal advection is the dominant contributor to the tendency term. The interannual variability of the SSS in the Indian Ocean is about five times larger than that in coupled model simulations of the North Atlantic Ocean. Regions of large interannual standard deviations are located near river mouths in the Bay of Bengal and in the eastern equatorial Indian Ocean. Both freshwater input into the ocean and advection of this anomalous flux are responsible for the generation of these anomalies. The model simulates 20 significant Indian Ocean Dipole (IOD) events and during IOD years large salinity anomalies appear in the equatorial Indian Ocean. The anomalies exist as two zonal bands: negative salinity anomalies to the north of the equator and positive to the south. The SSS anomalies for the years in which IOD is not present and for ENSO years are much weaker than during IOD years. Significant interannual SSS anomalies appear in the Indian Ocean only during IOD years.

Tree-ring based hydroclimate reconstruction over northern Vietnam from Fokienia hodginsii: eighteenth century mega-drought and tropical Pacific influence

We present here the first statistically calibrated and verified tree-ring reconstruction of climate from continental Southeast Asia.The reconstructed variable is March-May (MAM) Palmer Drought Severity Index (PDSI) based on ring widths from 22 trees (42 radial cores) of rare and long-lived conifer, Fokienia hodginsii (Po Mu as locally called) from northern Vietnam. This is the first published tree ring chronology from Vietnam as well as the first for this species. Spanning 535 years, this is the longest cross-dated tree-ring series yet produced from continental Southeast Asia. Response analysis revealed that the annual growth of Fokienia at this site was mostly governed by soil moisture in the pre-monsoon season. The reconstruction passed the calibration-verification tests commonly used in dendroclimatology, and revealed two prominent periods of drought in the mid-eighteenth and late-nineteenth enturies. The former lasted nearly 30 years and was concurrent with a similar drought over northwestern Thailand inferred from teak rings, suggesting a ``mega-drought'' extending across Indochina in the eighteenth century. Both of our reconstructed droughts are consistent with the periods of warm sea surface temperature (SST)anomalies in the tropical Pacific. Spatial correlation analyses with global SST indicated that ENSO-like anomalies might play a role in modulating droughts over the region, with El Nio (warm) phases resulting in reduced rainfall. However, significant correlation was also seen with SST over the Indian Ocean and the north Pacific,suggesting that ENSO is not the only factor affecting the climate of the area. Spectral analyses revealed significant peaks in the range of 53.9-78.8 years as well as in the ENSO-variability range of 2.0 to 3.2 years.

The effect of horizontal resolution on simulation of very extreme US precipitation events in a global atmosphere model

We investigate the ability of a global atmospheric general circulation model (AGCM) to reproduce observed 20 year return values of the annual maximum daily precipitation totals over the continental United States as a function of horizontal resolution. We find that at the high resolutions enabled by contemporary supercomputers, the AGCM can produce values of comparable magnitude to high quality observations. However, at the resolutions typical of the coupled general circulation models used in the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, the precipitation return values are severely underestimated.

Quasi-periodic, global oscillations in sea level pressure on intraseasonal timescales

The sea level pressure (SLP) variability in 30-60 day intraseasonal timescales is investigated using 25 years of reanalysis data addressing two issues. The first concerns the non-zero zonal mean component of SLP near the equator and its meridional connections, and the second concerns the fast eastward propagation (EP) speed of SLP compared to that of zonal wind. It is shown that the entire globe resonates with high amplitude wave activity during some periods which may last for few to several months, followed by lull periods of varying duration. SLP variations in the tropical belt are highly coherent from 25A degrees S to 25A degrees N, uncorrelated with variations in mid latitudes and again significantly correlated but with opposite phase around 60A degrees S and 65A degrees N. Near the equator (8A degrees S-8A degrees N), the zonal mean contributes significantly to the total variance in SLP, and after its removal, SLP shows a dominant zonal wavenumber one structure having a periodicity of 40 days and EP speeds comparable to that of zonal winds in the Indian Ocean. SLP from many of the atmospheric and coupled general circulation models show similar behaviour in the meridional direction although their propagation characteristics in the tropical belt differ widely.

Albedo enhancement of marine clouds to counteract global warming: impacts on the hydrological cycle

Recent studies have shown that changes in solar radiation affect the hydrological cycle more strongly than equivalent CO(2) changes for the same change in global mean surface temperature. Thus, solar radiation management ``geoengineering'' proposals to completely offset global mean temperature increases by reducing the amount of absorbed sunlight might be expected to slow the global water cycle and reduce runoff over land. However, proposed countering of global warming by increasing the albedo of marine clouds would reduce surface solar radiation only over the oceans. Here, for an idealized scenario, we analyze the response of temperature and the hydrological cycle to increased reflection by clouds over the ocean using an atmospheric general circulation model coupled to a mixed layer ocean model. When cloud droplets are reduced in size over all oceans uniformly to offset the temperature increase from a doubling of atmospheric CO(2), the global-mean precipitation and evaporation decreases by about 1.3% but runoff over land increases by 7.5% primarily due to increases over tropical land. In the model, more reflective marine clouds cool the atmospheric column over ocean. The result is a sinking motion over oceans and upward motion over land. We attribute the increased runoff over land to this increased upward motion over land when marine clouds are made more reflective. Our results suggest that, in contrast to other proposals to increase planetary albedo, offsetting mean global warming by reducing marine cloud droplet size does not necessarily lead to a drying, on average, of the continents. However, we note that the changes in precipitation, evaporation and P-E are dominated by small but significant areas, and given the highly idealized nature of this study, a more thorough and broader assessment would be required for proposals of altering marine cloud properties on a large scale.

Processes of 30-90 days sea surface temperature variability in the northern Indian Ocean during boreal summer

During summer, the northern Indian Ocean exhibits significant atmospheric intraseasonal variability associated with active and break phases of the monsoon in the 30-90 days band. In this paper, we investigate mechanisms of the Sea Surface Temperature (SST) signature of this atmospheric variability, using a combination of observational datasets and Ocean General Circulation Model sensitivity experiments. In addition to the previously-reported intraseasonal SST signature in the Bay of Bengal, observations show clear SST signals in the Arabian Sea related to the active/break cycle of the monsoon. As the atmospheric intraseasonal oscillation moves northward, SST variations appear first at the southern tip of India (day 0), then in the Somali upwelling region (day 10), northern Bay of Bengal (day 19) and finally in the Oman upwelling region (day 23). The Bay of Bengal and Oman signals are most clearly associated with the monsoon active/break index, whereas the relationship with signals near Somali upwelling and the southern tip of India is weaker. In agreement with previous studies, we find that heat flux variations drive most of the intraseasonal SST variability in the Bay of Bengal, both in our model (regression coefficient, 0.9, against similar to 0.25 for wind stress) and in observations (0.8 regression coefficient); similar to 60% of the heat flux variation is due do shortwave radiation and similar to 40% due to latent heat flux. On the other hand, both observations and model results indicate a prominent role of dynamical oceanic processes in the Arabian Sea. Wind-stress variations force about 70-100% of SST intraseasonal variations in the Arabian Sea, through modulation of oceanic processes (entrainment, mixing, Ekman pumping, lateral advection). Our similar to 100 km resolution model suggests that internal oceanic variability (i.e. eddies) contributes substantially to intraseasonal variability at small-scale in the Somali upwelling region, but does not contribute to large-scale intraseasonal SST variability due to its small spatial scale and random phase relation to the active-break monsoon cycle. The effect of oceanic eddies; however, remains to be explored at a higher spatial resolution.

Characteristics of Arabian Sea mini warm pool and Indian summer monsoon

Arabian Sea Mini Warm Pool (ASMWP) is a part of the Indian Ocean Warm Pool and formed in the eastern Arabian Sea prior to the onset of the summer monsoon season. This warm pool attained its maximum intensity during the pre-monsoon season and dissipated with the commencement of summer monsoon. The main focus of the present work was on the triggering of the dissipation of this warm pool and its relation to the onset of summer monsoon over Kerala. This phenomenon was studied utilizing NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric and Research) re-analysis data, TRMM Micro wave Imager (TMI) and observational data. To define the ASMWP, sea surface temperature exceeding 30.25A degrees C was taken as the criteria. The warm pool attained its maximum dimension and intensity nearly 2 weeks prior to the onset of summer monsoon over Kerala. Interestingly, the warm pool started its dissipation immediately after attaining its maximum core temperature. This information can be included in the present numerical models to enhance the prediction capability. It was also found that the extent and intensity of the ASMWP varied depending on the type of monsoon i.e., excess, normal, and deficient monsoon. Maximum core temperature and wide coverage of the warm pool observed during the excess monsoon years compared to normal and deficient monsoon years. The study also revealed a strong relationship between the salinity in the eastern Arabian Sea and the nature of the monsoon.

Albedo enhancement over land to counteract global warming: impacts on hydrological cycle

A recent modelling study has shown that precipitation and runoff over land would increase when the reflectivity of marine clouds is increased to counter global warming. This implies that large scale albedo enhancement over land could lead to a decrease in runoff over land. In this study, we perform simulations using NCAR CAM3.1 that have implications for Solar Radiation Management geoengineering schemes that increase the albedo over land. We find that an increase in reflectivity over land that mitigates the global mean warming from a doubling of CO2 leads to a large residual warming in the southern hemisphere and cooling in the northern hemisphere since most of the land is located in northern hemisphere. Precipitation and runoff over land decrease by 13.4 and 22.3%, respectively, because of a large residual sinking motion over land triggered by albedo enhancement over land. Soil water content also declines when albedo over land is enhanced. The simulated magnitude of hydrological changes over land are much larger when compared to changes over oceans in the recent marine cloud albedo enhancement study since the radiative forcing over land needed (-8.2 W m(-2)) to counter global mean radiative forcing from a doubling of CO2 (3.3 W m(-2)) is approximately twice the forcing needed over the oceans (-4.2 W m(-2)). Our results imply that albedo enhancement over oceans produce climates closer to the unperturbed climate state than do albedo changes on land when the consequences on land hydrology are considered. Our study also has important implications for any intentional or unintentional large scale changes in land surface albedo such as deforestation/afforestation/reforestation, air pollution, and desert and urban albedo modification.

A study of vertical cloud structure of the Indian summer monsoon using CloudSat data

Precise specification of the vertical distribution of cloud optical properties is important to reduce the uncertainty in quantifying the radiative impacts of clouds. The new global observations of vertical profiles of clouds from the CloudSat mission provide opportunities to describe cloud structures and to improve parameterization of clouds in the weather and climate prediction models. In this study, four years (2007-2010) of observations of vertical structure of clouds from the CloudSat cloud profiling radar have been used to document the mean vertical structure of clouds associated with the Indian summer monsoon (ISM) and its intra-seasonal variability. Active and break monsoon spells associated with the intra-seasonal variability of ISM have been identified by an objective criterion. For the present analysis, we considered CloudSat derived column integrated cloud liquid and ice water, and vertically profiles of cloud liquid and ice water content. Over the South Asian monsoon region, deep convective clouds with large vertical extent (up to 14 km) and large values of cloud water and ice content are observed over the north Bay of Bengal. Deep clouds with large ice water content are also observed over north Arabian Sea and adjoining northwest India, along the west coast of India and the south equatorial Indian Ocean. The active monsoon spells are characterized by enhanced deep convection over the Bay of Bengal, west coast of India and northeast Arabian Sea and suppressed convection over the equatorial Indian Ocean. Over the Bay of Bengal, cloud liquid water content and ice water content is enhanced by similar to 90 and similar to 200 % respectively during the active spells. An interesting feature associated with the active spell is the vertical tilting structure of positive CLWC and CIWC anomalies over the Arabian Sea and the Bay of Bengal, which suggests a pre-conditioning process for the northward propagation of the boreal summer intra-seasonal variability. It is also observed that during the break spells, clouds are not completely suppressed over central India. Instead, clouds with smaller vertical extent (3-5 km) are observed due to the presence of a heat low type of circulation. The present results will be useful for validating the vertical structure of clouds in weather and climate prediction models.

An estimate of equilibrium sensitivity of global terrestrial carbon cycle using NCAR CCSM4

Increasing concentrations of atmospheric CO2 influence climate, terrestrial biosphere productivity and ecosystem carbon storage through its radiative, physiological and fertilization effects. In this paper, we quantify these effects for a doubling of CO2 using a low resolution configuration of the coupled model NCAR CCSM4. In contrast to previous coupled climate-carbon modeling studies, we focus on the near-equilibrium response of the terrestrial carbon cycle. For a doubling of CO2, the radiative effect on the physical climate system causes global mean surface air temperature to increase by 2.14 K, whereas the physiological and fertilization on the land biosphere effects cause a warming of 0.22 K, suggesting that these later effects increase global warming by about 10 % as found in many recent studies. The CO2-fertilization leads to total ecosystem carbon gain of 371 Gt-C (28 %) while the radiative effect causes a loss of 131 Gt-C (10 %) indicating that climate warming damps the fertilization-induced carbon uptake over land. Our model-based estimate for the maximum potential terrestrial carbon uptake resulting from a doubling of atmospheric CO2 concentration (285-570 ppm) is only 242 Gt-C. This highlights the limited storage capacity of the terrestrial carbon reservoir. We also find that the terrestrial carbon storage sensitivity to changes in CO2 and temperature have been estimated to be lower in previous transient simulations because of lags in the climate-carbon system. Our model simulations indicate that the time scale of terrestrial carbon cycle response is greater than 500 years for CO2-fertilization and about 200 years for temperature perturbations. We also find that dynamic changes in vegetation amplify the terrestrial carbon storage sensitivity relative to a static vegetation case: because of changes in tree cover, changes in total ecosystem carbon for CO2-direct and climate effects are amplified by 88 and 72 %, respectively, in simulations with dynamic vegetation when compared to static vegetation simulations.

Variability of cloud liquid water and ice over South Asia from TMI estimates

In this study, the Tropical Rainfall Measurement Mission based Microwave Imager estimates (2A12) have been used to compare and contrast the characteristics of cloud liquid water and ice over the Indian land region and the ocean surrounding it, during the premonsoon (May) and monsoon (June-September) seasons. Based on the spatial homogeneity of rainfall, we have selected five regions for our study (three over ocean, two over land). Comparison across three ocean regions suggests that the cloud liquid water (CLW) over the orographically influenced Arabian Sea (close to the Indian west coast) behaves differently from the CLW over a trapped ocean (Bay of Bengal) or an open ocean (equatorial Indian Ocean). Specifically, the Arabian Sea region shows higher liquid water for a lower range of rainfall, whereas the Bay of Bengal and the equatorial Indian Ocean show higher liquid water for a higher range of rainfall. Apart from geographic differences, we also documented seasonal differences by comparing CLW profiles between monsoon and premonsoon periods, as well as between early and peak phases of the monsoon. We find that the CLW during the lean periods of rainfall (May or June) is higher than during the peak and late monsoon season (July-September) for raining clouds. As active and break phases are important signatures of the monsoon progression, we also analysed the differences in CLW during various phases of the monsoon, namely, active, break, active-to-break and break-to-active transition phases. We find that the cloud liquid water content during the break-to-active transition phase is significantly higher than during the active-to-break transition phase over central India. We speculate that this could be attributed to higher amount of aerosol loading over this region during the break phase. We lend credence to this aerosol-CLW/rain association by comparing the central Indian CLW with that over southeast Asia (where the aerosol loading is significantly smaller) and find that in the latter region, there are no significant differences in CLW during the different phases of the monsoon. While our hypothesis needs to be further investigated with numerical models, the results presented in this study can potentially serve as a good benchmark in evaluating the performance of cloud resolving models over the Indian region.