Ensemble transform Kalman-Bucy filters

Two recent works have adapted the Kalman–Bucy filter into an ensemble setting. In the first formulation, the ensemble of perturbations is updated by the solution of an ordinary differential equation (ODE) in pseudo-time, while the mean is updated as in the standard Kalman filter. In the second formulation, the full ensemble is updated in the analysis step as the solution of single set of ODEs in pseudo-time. Neither requires matrix inversions except for the frequently diagonal observation error covariance. We analyse the behaviour of the ODEs involved in these formulations. We demonstrate that they stiffen for large magnitudes of the ratio of background error to observational error variance, and that using the integration scheme proposed in both formulations can lead to failure. A numerical integration scheme that is both stable and is not computationally expensive is proposed. We develop transform-based alternatives for these Bucy-type approaches so that the integrations are computed in ensemble space where the variables are weights (of dimension equal to the ensemble size) rather than model variables. Finally, the performance of our ensemble transform Kalman–Bucy implementations is evaluated using three models: the 3-variable Lorenz 1963 model, the 40-variable Lorenz 1996 model, and a medium complexity atmospheric general circulation model known as SPEEDY. The results from all three models are encouraging and warrant further exploration of these assimilation techniques.

Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum

We use a state-of-the-art ocean general circulation and biogeochemistry model to examine the impact of changes in ocean circulation and biogeochemistry in governing the change in ocean carbon-13 and atmospheric CO2 at the last glacial maximum (LGM). We examine 5 different realisations of the ocean's overturning circulation produced by a fully coupled atmosphere-ocean model under LGM forcing and suggested changes in the atmospheric deposition of iron and phytoplankton physiology at the LGM. Measured changes in carbon-13 and carbon-14, as well as a qualitative reconstruction of the change in ocean carbon export are used to evaluate the results. Overall, we find that while a reduction in ocean ventilation at the LGM is necessary to reproduce carbon-13 and carbon-14 observations, this circulation results in a low net sink for atmospheric CO2. In contrast, while biogeochemical processes contribute little to carbon isotopes, we propose that most of the change in atmospheric CO2 was due to such factors. However, the lesser role for circulation means that when all plausible factors are accounted for, most of the necessary CO2 change remains to be explained. This presents a serious challenge to our understanding of the mechanisms behind changes in the global carbon cycle during the geologic past.

Simulation of late-twenty-first-century changes in wintertime atmospheric circulation over Europe due to anthropogenic causes

Evidence is presented, based on an ensemble of climate change scenarios performed with a global general circulation model of the atmosphere with high horizontal resolution over Europe, to suggest that the end-of-century anthropogenic climate change over the North Atlantic--European region strongly projects onto the positive phase of the North Atlantic Oscillation during wintertime. It is reflected in a doubling of the residence frequency of the climate system in the associated circulation regime, in agreement with the nonlinear climate perspective. The strong increase in the amplitude of the response, compared to coarse-resolution coupled model studies, suggests that improved model representation of regional climate is needed to achieve more reliable projections of anthropogenic climate change on European climate.

Millennial climatic fluctuations are key to the structure of last glacial ecosystems

Whereas fossil evidence indicates extensive treeless vegetation and diverse grazing megafauna in Europe and northern Asia during the last glacial, experiments combining vegetation models and climate models have to-date simulated widespread persistence of trees. Resolving this conflict is key to understanding both last glacial ecosystems and extinction of most of the mega-herbivores. Using a dynamic vegetation model (DVM) we explored the implications of the differing climatic conditions generated by a general circulation model (GCM) in “normal” and “hosing” experiments. Whilst the former approximate interstadial conditions, the latter, designed to mimic Heinrich Events, approximate stadial conditions. The “hosing” experiments gave simulated European vegetation much closer in composition to that inferred from fossil evidence than did the “normal” experiments. Given the short duration of interstadials, and the rate at which forest cover expanded during the late-glacial and early Holocene, our results demonstrate the importance of millennial variability in determining the character of last glacial ecosystems.

Simple uncertainty frameworks for selecting weighting schemes and interpreting multimodel ensemble climate change experiments

Future climate change projections are often derived from ensembles of simulations from multiple global circulation models using heuristic weighting schemes. This study provides a more rigorous justification for this by introducing a nested family of three simple analysis of variance frameworks. Statistical frameworks are essential in order to quantify the uncertainty associated with the estimate of the mean climate change response. The most general framework yields the “one model, one vote” weighting scheme often used in climate projection. However, a simpler additive framework is found to be preferable when the climate change response is not strongly model dependent. In such situations, the weighted multimodel mean may be interpreted as an estimate of the actual climate response, even in the presence of shared model biases. Statistical significance tests are derived to choose the most appropriate framework for specific multimodel ensemble data. The framework assumptions are explicit and can be checked using simple tests and graphical techniques. The frameworks can be used to test for evidence of nonzero climate response and to construct confidence intervals for the size of the response. The methodology is illustrated by application to North Atlantic storm track data from the Coupled Model Intercomparison Project phase 5 (CMIP5) multimodel ensemble. Despite large variations in the historical storm tracks, the cyclone frequency climate change response is not found to be model dependent over most of the region. This gives high confidence in the response estimates. Statistically significant decreases in cyclone frequency are found on the flanks of the North Atlantic storm track and in the Mediterranean basin.

Mid-Holocene land-surface conditions in northern Africa and the Arabian peninsula: a data set for the analysis of biogeophysical feedbacks in the climate system

Large changes in the extent of northern subtropical arid regions during the Holocene are attributed to orbitally forced variations in monsoon strength and have been implicated in the regulation of atmospheric trace gas concentrations on millenial timescales. Models that omit biogeophysical feedback, however, are unable to account for the full magnitude of African monsoon amplification and extension during the early to middle Holocene (˜9500–5000 years B.P.). A data set describing land-surface conditions 6000 years B.P. on a 1° × 1° grid across northern Africa and the Arabian Peninsula has been prepared from published maps and other sources of palaeoenvironmental data, with the primary aim of providing a realistic lower boundary condition for atmospheric general circulation model experiments similar to those performed in the Palaeoclimate Modelling Intercomparison Project. The data set includes information on the percentage of each grid cell occupied by specific vegetation types (steppe, savanna, xerophytic woods/scrub, tropical deciduous forest, and tropical montane evergreen forest), open water (lakes), and wetlands, plus information on the flow direction of major drainage channels for use in large-scale palaeohydrological modeling.

The composite-tendency RAW filter in semi-implicit integrations

The time discretization in weather and climate models introduces truncation errors that limit the accuracy of the simulations. Recent work has yielded a method for reducing the amplitude errors in leapfrog integrations from first-order to fifth-order. This improvement is achieved by replacing the Robert--Asselin filter with the RAW filter and using a linear combination of the unfiltered and filtered states to compute the tendency term. The purpose of the present paper is to apply the composite-tendency RAW-filtered leapfrog scheme to semi-implicit integrations. A theoretical analysis shows that the stability and accuracy are unaffected by the introduction of the implicitly treated mode. The scheme is tested in semi-implicit numerical integrations in both a simple nonlinear stiff system and a medium-complexity atmospheric general circulation model, and yields substantial improvements in both cases. We conclude that the composite-tendency RAW-filtered leapfrog scheme is suitable for use in semi-implicit integrations.

On the clustering of winter storm loss events over Germany

During the last decades, several windstorm series hit Europe leading to large aggregated losses. Such storm series are examples of serial clustering of extreme cyclones, presenting a considerable risk for the insurance industry. Clustering of events and return periods of storm series for Germany are quantified based on potential losses using empirical models. Two reanalysis data sets and observations from German weather stations are considered for 30 winters. Histograms of events exceeding selected return levels (1-, 2- and 5-year) are derived. Return periods of historical storm series are estimated based on the Poisson and the negative binomial distributions. Over 4000 years of general circulation model (GCM) simulations forced with current climate conditions are analysed to provide a better assessment of historical return periods. Estimations differ between distributions, for example 40 to 65 years for the 1990 series. For such less frequent series, estimates obtained with the Poisson distribution clearly deviate from empirical data. The negative binomial distribution provides better estimates, even though a sensitivity to return level and data set is identified. The consideration of GCM data permits a strong reduction of uncertainties. The present results support the importance of considering explicitly clustering of losses for an adequate risk assessment for economical applications.

Methane in the atmosphere of the transiting hot Neptune GJ436b?

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5, and 8 μm obtained with the Infrared Array Camera on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, Hubble Space Telescope, and ground-based V, I, H, and Ks published observations, the range 0.5-10 μm can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the data set. Representative climate models were calculated by using a three-dimensional, pseudospectral general circulation model with idealized thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio-calculated, line list for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water, and other molecules. No clear evidence of carbon monoxide and carbon dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesized to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

A process-based analysis of ocean heat uptake in an AOGCM with an eddy-permitting ocean component

About 90% of the anthropogenic increase in heat stored in the climate system is found the oceans. Therefore it is relevant to understand the details of ocean heat uptake. Here we present a detailed, process-based analysis of ocean heat uptake (OHU) processes in HiGEM1.2, an atmosphere-ocean general circulation model (AOGCM) with an eddy-permitting ocean component of 1/3 degree resolution. Similarly to various other models, HiGEM1.2 shows that the global heat budget is dominated by a downward advection of heat compensated by upward isopycnal diffusion. Only in the upper tropical ocean do we find the classical balance between downward diapycnal diffusion and upward advection of heat. The upward isopycnal diffusion of heat is located mostly in the Southern Ocean, which thus dominates the global heat budget. We compare the responses to a 4xCO2 forcing and an enhancement of the windstress forcing in the Southern Ocean. This highlights the importance of regional processes for the global ocean heat uptake. These are mainly surface fluxes and convection in the high latitudes, and advection in the Southern Ocean mid-latitudes. Changes in diffusion are less important. In line with the CMIP5 models, HiGEM1.2 shows a band of strong OHU in the mid-latitude Southern Ocean in the 4xCO2 run, which is mostly advective. By contrast, in the high-latitude Southern Ocean regions it is the suppression of convection that leads to OHU. In the enhanced windstress run, convection is strengthened at high Southern latitudes, leading to heat loss, while the magnitude of the OHU in the Southern mid-latitudes is very similar to the 4xCO2 results. Remarkably, there is only very small global OHU in the enhanced windstress run. The wind stress forcing just leads to a redistribution of heat. We relate the ocean changes at high southern latitudes to the effect of climate change on the Antarctic Circumpolar Current (ACC). It weakens in the 4xCO2 run and strengthens in the wind stress run. The weakening is due to a narrowing of the ACC, caused by an expansion of the Weddell Gyre, and a flattening of the isopycnals, which are explained by a combination of the wind stress forcing and increased precipitation.

Climate change modelling and its roles to chinese crops yield

Climate has been changing in the last fifty years in China and will continue to change regardless any efforts for mitigation. Agriculture is a climate-dependent activity and highly sensitive to climate changes and climate variability. Understanding the interactions between climate change and agricultural production is essential for society stable development of China. The first mission is to fully understand how to predict future climate and link it with agriculture production system. In this paper, recent studies both domestic and international are reviewed in order to provide an overall image of the progress in climate change researches. The methods for climate change scenarios construction are introduced. The pivotal techniques linking crop model and climate models are systematically assessed and climate change impacts on Chinese crops yield among model results are summarized. The study found that simulated productions of grain crop inherit uncertainty from using different climate models, emission scenarios and the crops simulation models. Moreover, studies have different spatial resolutions, and methods for general circulation model (GCM) downscaling which increase the uncertainty for regional impacts assessment. However, the magnitude of change in crop production due to climate change (at 700 ppm CO2 eq correct) appears within ±10% for China in these assessments. In most literatures, the three cereal crop yields showed decline under climate change scenarios and only wheat in some region showed increase. Finally, the paper points out several gaps in current researches which need more studies to shorten the distance for objective recognizing the impacts of climate change on crops. The uncertainty for crop yield projection is associated with climate change scenarios, CO2 fertilization effects and adaptation options. Therefore, more studies on the fields such as free air CO2 enrichment experiment and practical adaptations implemented need to be carried out

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.

Explaining patterns of avian diversity and endemicity: climate and biomes of southern Africa over the last 140,000 years

Aim Test hypotheses that present biodiversity and endemic species richness are related to climatic stability and/or biome persistence.Location Africa south of 15° S. Methods Seventy eight HadCM3 general circulation model palaeoclimate experiments spanning the last 140,000 years, plus a pre-industrial experiment,were used to calculate measures of climatic variability for 0.5° grid cells. Models were fitted relating distributions of the nine biomes of South Africa,Lesotho and Swaziland to present climate. These models were used to simulate potential past biome distribution and extent for the 78 palaeoclimate experiments, and three measures of biome persistence. Climatic response surfaces were fitted for 690 bird species regularly breeding in the region and used to simulate present species richness for cells of the 0.5° grid. Species richness was evaluated for residents, mobile species (nomadic or partially/altitudinally migrant within the region), and intra-African migrants, and also separately for endemic/near-endemic (hereafter ‘endemic’) species as a whole and those associated with each biome. Our hypotheses were tested by analysing correlations between species richness and climatic variability or biome persistence. Results The magnitude of climatic variability showed clear spatial patterns. Marked changes in biome distributions and extents were projected, although limited areas of persistence were projected for some biomes. Overall species richness was not correlated with climatic variability, although richness of mobile species showed a weak negative correlation. Endemic species richness was significantly negatively correlated with climatic variability. Strongest correlations, however, were positive correlations between biome persistence and richness of endemics associated with individual biomes. Main conclusions Low climatic variability, and especially a degree of stability enabling biome persistence, is strongly correlated with species richness of birds endemic to southern Africa. This probably principally reflects reduced extinction risk for these species where the biome to which they are adapted persisted

IPCC global coupled model simulations of the South America monsoon system

This study examines the variability of the South America monsoon system (SAMS) over tropical South America (SA). The onset, end, and total rainfall during the summer monsoon are investigated using precipitation pentad estimates from the global precipitation climatology project (GPCP) 1979-2006. Likewise, the variability of SAMS characteristics is examined in ten Intergovernmental Panel on Climate Change (IPCC) global coupled climate models in the twentieth century (1981-2000) and in a future scenario of global change (A1B) (2081-2100). It is shown that most IPCC models misrepresent the intertropical convergence zone and therefore do not capture the actual annual cycle of precipitation over the Amazon and northwest SA. Most models can correctly represent the spatiotemporal variability of the annual cycle of precipitation in central and eastern Brazil such as the correct phase of dry and wet seasons, onset dates, duration of rainy season and total accumulated precipitation during the summer monsoon for the twentieth century runs. Nevertheless, poor representation of the total monsoonal precipitation over the Amazon and northeast Brazil is observed in a large majority of the models. Overall, MI-ROC3.2-hires, MIROC3.2-medres and MRI-CGCM3.2.3 show the most realistic representation of SAMS`s characteristics such as onset, duration, total monsoonal precipitation, and its interannual variability. On the other hand, ECHAM5, GFDL-CM2.0 and GFDL-CM2.1 have the least realistic representation of the same characteristics. For the A1B scenario the most coherent feature observed in the IPCC models is a reduction in precipitation over central-eastern Brazil during the summer monsoon, comparatively with the present climate. The IPCC models do not indicate statistically significant changes in SAMS onset and demise dates for the same scenario.

Transformed symmetric models

For the first time, we introduce a class of transformed symmetric models to extend the Box and Cox models to more general symmetric models. The new class of models includes all symmetric continuous distributions with a possible non-linear structure for the mean and enables the fitting of a wide range of models to several data types. The proposed methods offer more flexible alternatives to Box-Cox or other existing procedures. We derive a very simple iterative process for fitting these models by maximum likelihood, whereas a direct unconditional maximization would be more difficult. We give simple formulae to estimate the parameter that indexes the transformation of the response variable and the moments of the original dependent variable which generalize previous published results. We discuss inference on the model parameters. The usefulness of the new class of models is illustrated in one application to a real dataset.