Fingerprinting of the Atlantic meridional overturning circulation in climate models to aid in the design of proxy investigations

It is desirable to design proxy investigations that target regions where properties reconstructed from calibrated parameters potentially carry high-fidelity information concerning changes in large-scale climate systems. Numerical climate models can play an important role in this task, producing simulations that can be analyzed to produce spatial “fingerprints” of the expected response of various properties under a variety of different scenarios. We will introduce a new method of fingerprinting the Atlantic meridional overturning circulation (AMOC) that not only provides information concerning the sensitivity of the response at a given location to changes in the large-scale system, but also quantifies the linearity, monotonicity and symmetry of the response. In this way, locations that show high sensitivities to changes in the AMOC, but that exhibit, for example, strongly nonlinear behavior can be avoided during proxy investigations. To demonstrate the proposed approach we will use the example of the response of seawater temperatures to changes in the strength of the AMOC. We present results from an earth-system climate model which has been perturbed with an idealized freshwater forcing scenario in order to reduce the strength of the AMOC in a systematic manner. The seawater temperature anomalies that result from the freshwater forcing are quantified in terms of their sensitivity to the AMOC strength in addition to the linearity and monotonicity of their response. A first-order reversal curve (FORC) approach is employed to investigate and quantify the irreversibility of the temperature response to a slowing and recovering AMOC. Thus, FORCs allow the identification of areas that are unsuitable for proxy reconstructions because their temperature versus AMOC relationship lacks symmetry.     

Quantifying Arctic contributions to climate predictability in a regional coupled ocean-ice-atmosphere model

The relative importance of regional processes inside the Arctic climate system and the large scale atmospheric circulation for Arctic interannual climate variability has been estimated with the help of a regional Arctic coupled ocean-ice-atmosphere model. The study focuses on sea ice and surface climate during the 1980s and 1990s. Simulations agree reasonably well with observations. Correlations between the winter North Atlantic Oscillation index and the summer Arctic sea ice thickness and summer sea ice extent are found. Spread of sea ice extent within an ensemble of model runs can be associated with a surface pressure gradient between the Nordic Seas and the Kara Sea. Trends in the sea ice thickness field are widely significant and can formally be attributed to large scale forcing outside the Arctic model domain. Concerning predictability, results indicate that the variability generated by the external forcing is more important in most regions than the internally generated variability. However, both are in the same order of magnitude. Local areas such as the Northern Greenland coast together with Fram Straits and parts of the Greenland Sea show a strong importance of internally generated variability, which is associated with wind direction variability due to interaction with atmospheric dynamics on the Greenland ice sheet. High predictability of sea ice extent is supported by north-easterly winds from the Arctic Ocean to Scandinavia.     

An Observational Study of the Mesoscale Mistral Dynamics

We investigate the mesoscale dynamics of the mistral through the wind profiler observations of the MAP (autumn 1999) and ESCOMPTE (summer 2001) field campaigns. We show that the mistral wind field can dramatically change on a time scale less than 3 hours. Transitions from a deep to a shallow mistral are often observed at any season when the lower layers are stable. The variability, mainly attributed in summer to the mistral/land–sea breeze interactions on a 10-km scale, is highlighted by observations from the wind profiler network set up during ESCOMPTE. The interpretations of the dynamical mistral structure are performed through comparisons with existing basic theories. The linear theory of R. B. Smith [Advances in Geophysics, Vol. 31, 1989, Academic Press, 1–41] and the shallow water theory [Schär, C. and Smith, R. B.: 1993a, J. Atmos. Sci. 50, 1373–1400] give some complementary explanations for the deep-to-shallow transition especially for the MAP mistral event. The wave breaking process induces a low-level jet (LLJ) downstream of the Alps that degenerates into a mountain wake, which in turn provokes the cessation of the mistral downstream of the Alps. Both theories indicate that the flow splits around the Alps and results in a persistent LLJ at the exit of the Rhône valley. The LLJ is strengthened by the channelling effect of the Rhône valley that is more efficient for north-easterly than northerly upstream winds despite the north–south valley axis. Summer moderate and weak mistral episodes are influenced by land–sea breezes and convection over land that induce a very complex interaction that cannot be accurately described by the previous theories.     

Laboratory study of the rates and products of the phototransformations of naphthalene adsorbed on samples of titanium dioxide,ferric oxide,muscovite, and fly ash

The photodegradation of naphthalene (NPH), chosen as a model of polynuclear aromatic pollutants, has been studied in the presence of a layer of four water-insoluble inorganic solids which can be found in the troposphere (TiO₂, Fe₂O₃, muscovite, and a fly ash sample). Direct photolysis of NPH is negligible at >340 nm. Dark adsorption of NPH on TiO₂ (mainly anatase, nonporous, 50 m² g^–1) at 293 K corresponds to a surface coverage ofca. 50% at equilibrium. Under these conditions (saturated surface), the stationary-state photocatalytic degradation reaches 0.4 molecule nm^–2 h^–1 (>340 nm, radiant fluxca. 22 mW cm^–2). Dioxygen is required and its partial pressure in air is such that the degradation is zero order in O₂. Water vapor markedly increases the rate. The other particulates have also an effect, less important than that of TiO₂, however quite noticeable with respect to surface area unit for the fly ash sample which contains 3.2% Fe₂O₃. Apart from 1,4-naphthoquinone, which is the main intermediate product in all cases, 2-naphthol, phthalide, phthaldialdehyde, phthalic acid, acetophenone, benzaldehyde, benzoic acid are also formed on dry TiO₂. Depending on their volatility, these compounds are transferred to the gas phase or remain principally adsorbed on the solid particles where they are further transformed. For instance, phthalic acid (or anhydride) and benzoic acid are generated from 1,4-naphthoquinone. Degradation mechanisms are briefly discussed.     

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

A copula-based multivariate analysis of Canadian RCM projected changes to flood characteristics for northeastern Canada

In the present work, climate change impacts on three spring (March–June) flood characteristics, i.e. peak, volume and duration, for 21 northeast Canadian basins are evaluated, based on Canadian regional climate model (CRCM) simulations. Conventional univariate frequency analysis for each flood characteristic and copula based bivariate frequency analysis for mutually correlated pairs of flood characteristics (i.e. peak–volume, peak–duration and volume–duration) are carried out. While univariate analysis is focused on return levels of selected return periods (5-, 20- and 50-year), the bivariate analysis is focused on the joint occurrence probabilities P1 and P2 of the three pairs of flood characteristics, where P1 is the probability of any one characteristic in a pair exceeding its threshold and P2 is the probability of both characteristics in a pair exceeding their respective thresholds at the same time. The performance of CRCM is assessed by comparing ERA40 (the European Centre for Medium-Range Weather Forecasts 40-year reanalysis) driven CRCM simulated flood statistics and univariate and bivariate frequency analysis results for the current 1970–1999 period with those observed at selected 16 gauging stations for the same time period. The Generalized Extreme Value distribution is selected as the marginal distribution for flood characteristics and the Clayton copula for developing bivariate distribution functions. The CRCM performs well in simulating mean, standard deviation, and 5-, 20- and 50-year return levels of flood characteristics. The joint occurrence probabilities are also simulated well by the CRCM. A five-member ensemble of the CRCM simulated streamflow for the current (1970–1999) and future (2041–2070) periods, driven by five different members of a Canadian Global Climate Model ensemble, are used in the assessment of projected changes, where future simulations correspond to A2 scenario. The results of projected changes, in general, indicate increases in the marginal values, i.e. return levels of flood characteristics, and the joint occurrence probabilities P1 and P2. It is found that the future marginal values of flood characteristics and P1 and P2 values corresponding to longer return periods will be affected more by anthropogenic climate change than those corresponding to shorter return periods but the former ones are subjected to higher uncertainties.     

Impacts of Indian Ocean SST biases on the Indian Monsoon: as simulated in a global coupled model

In this study, the impact of the ocean–atmosphere coupling on the atmospheric mean state over the Indian Ocean and the Indian Summer Monsoon (ISM) is examined in the framework of the SINTEX-F2 coupled model through forced and coupled control simulations and several sensitivity coupled experiments. During boreal winter and spring, most of the Indian Ocean biases are common in forced and coupled simulations, suggesting that the errors originate from the atmospheric model, especially a dry islands bias in the Maritime Continent. During boreal summer, the air-sea coupling decreases the ISM rainfall over South India and the monsoon strength to realistic amplitude, but at the expense of important degradations of the rainfall and Sea Surface Temperature (SST) mean states in the Indian Ocean. Strong SST biases of opposite sign are observed over the western (WIO) and eastern (EIO) tropical Indian Ocean. Rainfall amounts over the ocean (land) are systematically higher (lower) in the northern hemisphere and the south equatorial Indian Ocean rainfall band is missing in the control coupled simulation. During boreal fall, positive dipole-like errors emerge in the mean state of the coupled model, with warm and wet (cold and dry) biases in the WIO (EIO), suggesting again a significant impact of the SST errors. The exact contributions and the distinct roles of these SST errors in the seasonal mean atmospheric state of the coupled model have been further assessed with two sensitivity coupled experiments, in which the SST biases are replaced by observed climatology either in the WIO (warm bias) or EIO (cold bias). The correction of the WIO warm bias leads to a global decrease of rainfall in the monsoon region, which confirms that the WIO is an important source of moisture for the ISM. On the other hand, the correction of the EIO cold bias leads to a global improvement of precipitation and circulation mean state during summer and fall. Nevertheless, all these improvements due to SST corrections seem drastically limited by the atmosphere intrinsic biases, including prominently the unimodal oceanic position of the ITCZ (Inter Tropical Convergence Zone) during summer and the enhanced westward wind stress along the equator during fall.     

Concentration profiles of particles settling in the neutral and stratified atmospheric boundary layer

An expression for the vertical equilibrium concentration profile of heavy particles, including the effects of canopy on the eddy diffusivity as well as corrections for atmospheric stability, is proposed. This expression is validated against measurements of vertical concentration profiles of corn pollen above a corn field. The fitted theoretical profiles show that particle settling is correctly accounted for. The sensitivity to variations in the turbulent Schmidt number, settling velocity and stability corrections are explicitly characterized. The importance of independent measurements of the surface flux of pollen in future experiments is noted.     

Projections of climate change impacts on central America tropical rainforest

Tropical rainforest plays an important role in the global carbon cycle, accounting for a large part of global net primary productivity and contributing to CO₂ sequestration. The objective of this work is to simulate potential changes in the rainforest biome in Central America subject to anthropogenic climate change under two emissions scenarios, RCP4.5 and RCP8.5. The use of a dynamic vegetation model and climate change scenarios is an approach to investigate, assess or anticipate how biomes respond to climate change. In this work, the Inland dynamic vegetation model was driven by the Eta regional climate model simulations. These simulations accept boundary conditions from HadGEM2-ES runs in the two emissions scenarios. The possible consequences of regional climate change on vegetation properties, such as biomass, net primary production and changes in forest extent and distribution, were investigated. The Inland model projections show reductions in tropical forest cover in both scenarios. The reduction of tropical forest cover is greater in RCP8.5. The Inland model projects biomass increases where tropical forest remains due to the CO₂ fertilization effect. The future distribution of predominant vegetation shows that some areas of tropical rainforest in Central America are replaced by savannah and grassland in RCP4.5. Inland projections under both RCP4.5 and RCP8.5 show a net primary productivity reduction trend due to significant tropical forest reduction, temperature increase, precipitation reduction and dry spell increments, despite the biomass increases in some areas of Costa Rica and Panama. This study may provide guidance to adaptation studies of climate change impacts on the tropical rainforests in Central America.     

A modified parameterization of flux-profile relationships in the surface layer using different roughness length values for heat and momentum

The surface-layer flux-profile formulae of Louis (1979), used in many atmospheric models, are modified in a simple way to allow for different values of the roughness lengths for heat and momentum. The modified set of formulae simplifies the calculation of surface-layer fluxes over most natural land surfaces, where the roughness length for momentum can be almost an order of magnitude greater than that for heat.