Heat stress response for physiological traits in dairy and dual purpose cattle populations on phenotypic and genetic scales

The main objective of this thesis was to determine the potential impact of heat stress (HS) on physiological traits of lactating cows and semen quality of bulls kept in a temperate climate. The thesis is comprised of three studies. An innovative statistical modeling aspect common to all three studies was the application of random regression methodology (RRM) to study the phenotypic and genetic trajectory of traits in dependency of a continuous temperature humidity index (THI). In the first study, semen quality and quantity traits of 562 Holstein sires kept on an AI station in northwestern Germany were analyzed in the course of THI calculated from data obtained from the nearest weather station. Heat stress was identified based on a decline in semen quality and quantity parameters. The identified general HS threshold (THI = 60) and the thermoneutal zone (THI in the range from 50 to 60) for semen production were lower than detected in studies conducted in tropical and subtropical climates. Even though adult bulls were characterized by higher semen productivity compared to younger bulls, they responded with a stronger semen production loss during harsh environments. Heritabilities (low to moderate range) and additive genetic variances of semen characteristics varied with different levels of THI. Also, based on genetic correlations genotype, by environment interactions were detected. Taken together, these findings suggest the application of specific selection strategies for specific climate conditions. In the second study, the effect of the continuous environmental descriptor THI as measured inside the barns on rectal temperatures (RT), skin temperatures (ST), vaginal temperatures (VT), respiration rates (RR), and pulse rate (PR) of lactating Holstein Friesian (HF) and dual-purpose German black pied cattle (DSN) was analyzed. Increasing HS from THI 65 (threshold) to THI 86 (maximal THI) resulted in an increase of RT by 0.6 °C (DSN) and 1 °C (HF), ST by 3.5 °C (HF) and 8 °C (DSN), VT by 0.3 °C (DSN), and RR by 47 breaths / minute (DSN), and decreased PR by 7 beats / minute (DSN). The undesired effects of rising THI on physiological traits were most pronounced for cows with high levels of milk yield and milk constituents, cows in early days in milk and later parities, and during summer seasons in the year 2014. In the third study of this dissertation, the genetic components of the cow’s physiological responses to HS were investigated. Heat stress was deduced from indoor THI measurements, and physiological traits were recorded on native DSN cows and their genetically upgraded crosses with Holstein Friesian sires in two experimental herds from pasture-based production systems reflecting a harsh environment of the northern part of Germany. Although heritabilities were in a low range (from 0.018 to 0.072), alterations of heritabilities, repeatabilities, and genetic components in the course of THI justify the implementation of genetic evaluations including heat stress components. However, low repeatabilities indicate the necessity of using repeated records for measuring physiological traits in German cattle. Moderate EBV correlations between different trait combinations indicate the potential of selection for one trait to simultaneously improve the other physiological attributes. In conclusion, bulls of AI centers and lactating cows suffer from HS during more extreme weather conditions also in the temperate climate of Northern Germany. Monitoring physiological traits during warm and humid conditions could provide precious information for detection of appropriate times for implementation of cooling systems and changes in feeding and management strategies. Subsequently, the inclusion of these physiological traits with THI specific breeding values into overall breeding goals could contribute to improving cattle adaptability by selecting the optimal animal for extreme hot and humid conditions. Furthermore, the recording of meteorological data in close distance to the cow and visualizing the surface body temperature by infrared thermography techniques might be helpful for recognizing heat tolerance and adaptability in cattle.

An observing system simulation experiment for climate monitoring with GNSS radio occulation data: setup and testbed study

The long-term stability, high accuracy, all-weather capability, high vertical resolution, and global coverage of Global Navigation Satellite System (GNSS) radio occultation (RO) suggests it as a promising tool for global monitoring of atmospheric temperature change. With the aim to investigate and quantify how well a GNSS RO observing system is able to detect climate trends, we are currently performing an (climate) observing system simulation experiment over the 25-year period 2001 to 2025, which involves quasi-realistic modeling of the neutral atmosphere and the ionosphere. We carried out two climate simulations with the general circulation model MAECHAM5 (Middle Atmosphere European Centre/Hamburg Model Version 5) of the MPI-M Hamburg, covering the period 2001–2025: One control run with natural variability only and one run also including anthropogenic forcings due to greenhouse gases, sulfate aerosols, and tropospheric ozone. On the basis of this, we perform quasi-realistic simulations of RO observables for a small GNSS receiver constellation (six satellites), state-of-the-art data processing for atmospheric profiles retrieval, and a statistical analysis of temperature trends in both the “observed” climatology and the “true” climatology. Here we describe the setup of the experiment and results from a test bed study conducted to obtain a basic set of realistic estimates of observational errors (instrument- and retrieval processing-related errors) and sampling errors (due to spatial-temporal undersampling). The test bed results, obtained for a typical summer season and compared to the climatic 2001–2025 trends from the MAECHAM5 simulation including anthropogenic forcing, were found encouraging for performing the full 25-year experiment. They indicated that observational and sampling errors (both contributing about 0.2 K) are consistent with recent estimates of these errors from real RO data and that they should be sufficiently small for monitoring expected temperature trends in the global atmosphere over the next 10 to 20 years in most regions of the upper troposphere and lower stratosphere (UTLS). Inspection of the MAECHAM5 trends in different RO-accessible atmospheric parameters (microwave refractivity and pressure/geopotential height in addition to temperature) indicates complementary climate change sensitivity in different regions of the UTLS so that optimized climate monitoring shall combine information from all climatic key variables retrievable from GNSS RO data.

Seasonal variations of the physical and optical characteristics of Saharan Dust: results from the Dust outflow and Deposition to the Ocean (DODO) experiment

North African dust is important for climate through its direct radiative effect on solar and terrestrial radiation and its role in the biogeochemical system. The Dust Outflow and Deposition to the Ocean project (DODO) aimed to characterize the physical and optical properties of airborne North African dust in two seasons and to use these observations to constrain model simulations, with the ultimate aim of being able to quantify the deposition of iron to the North Atlantic Ocean. The in situ properties of dust from airborne campaigns measured during February and August 2006, based at Dakar, Senegal, are presented here. Average values of the single scattering albedo (0.99, 0.98), mass specific extinction (0.85 m^2 g^-1 , 1.14 m^2 g^-1 ), asymmetry parameter (0.68, 0.68), and refractive index (1.53--0.0005i,1.53--0.0014i) for the accumulation mode were found to differ by varying degrees between the dry and wet season, respectively. It is hypothesized that these differences are due to different source regions and transport processes which also differ between the DODO campaigns. Elemental ratios of Ca/Al were found to differ between the dry and wet season (1.1 and 0.5, respectively). Differences in vertical profiles are found between seasons and between land and ocean locations and reflect the different dynamics of the seasons. Using measurements of the coarse mode size distribution and illustrative Mie calculations, the optical properties are found to be very sensitive to the presence and amount of coarse mode of mineral dust, and the importance of accurate measurements of the coarse mode of dust is highlighted.

Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional- and shear-wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse-echo method. The measurements were made both in vacuum-dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin-impregnated polished thin sections, X-ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional- and shear-wave velocities (V-p and V-s, respectively) decrease with increasing porosity and that V-p decreases approximately twice as fast as V-s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore-structure-dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot-Gassmann velocity values are greater than the measured velocity values due to the rock-fluid interaction. This is not accounted for in the Biot-Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time-average relationships overestimated the measured velocities even more than the Biot model.

Evaluating the century C model using long-term fertilizer trials in the Indo-Gangetic Plains, India

The GEFSOC Project developed a system for estimating soil carbon (C) stocks and changes at the national and sub-national scale. As part of the development of the system, the Century ecosystem model was evaluated for its ability to simulate soil organic C (SOC) changes in environmental conditions in the Indo-Gangetic Plains, India (IGP). Two long-term fertilizer trials (LTFT), with all necessary parameters needed to run Century, were used for this purpose: a jute (Corchorus capsularis L.), rice (Oryza sativa L.) and wheat (Triticum aestivum L.) trial at Barrackpore, West Bengal, and a rice-wheat trial at Ludhiana, Punjab. The trials represent two contrasting climates of the IGP, viz. semi-arid, dry with mean annual rainfall (MAR) of < 800 mm and humid with > 1600 turn. Both trials involved several different treatments with different organic and inorganic fertilizer inputs. In general, the model tended to overestimate treatment effects by approximately 15%. At the semi-arid site, modelled data simulated actual data reasonably well for all treatments, with the control and chemical N + farm yard manure showing the best agreement (RMSE = 7). At the humid site, Century performed less well. This could have been due to a range of factors including site history. During the study, Century was calibrated to simulate crop yields for the two sites considered using data from across the Indian IGP. However, further adjustments may improve model performance at these sites and others in the IGP. The availability of more longterm experimental data sets (especially those involving flooded lowland rice and triple cropping systems from the IGP) for testing and validation is critical to the application of the model's predictive capabilities for this area of the Indian sub-continent. (C) 2007 Elsevier B.V. All rights reserved.

Quasi-biennial oscillation and tracer distributions in a coupled chemistry-climate model

[ 1] We have used a fully coupled chemistry-climate model (CCM), which generates its own wind and temperature quasi-biennial oscillation (QBO), to study the effect of coupling on the QBO and to examine the QBO signals in stratospheric trace gases, particularly ozone. Radiative coupling of the interactive chemistry to the underlying general circulation model tends to prolong the QBO period and to increase the QBO amplitude in the equatorial zonal wind in the lower and middle stratosphere. The model ozone QBO agrees well with Stratospheric Aerosol and Gas Experiment II and Total Ozone Mapping Spectrometer satellite observations in terms of vertical and latitudinal structure. The model captures the ozone QBO phase change near 28 km over the equator and the column phase change near +/- 15 degrees latitude. Diagnosis of the model chemical terms shows that variations in NOx are the main chemical driver of the O-3 QBO around 35 km, i.e., above the O-3 phase change.

Overview of observations from the RADAGAST experiment in Niamey, Niger: 2. Radiative fluxes and divergences

Broadband shortwave and longwave radiative fluxes observed both at the surface and from space during the Radiative Atmospheric Divergence using ARM Mobile Facility, GERB data and AMMA Stations (RADAGAST) experiment in Niamey, Niger, in 2006 are presented. The surface fluxes were measured by the Atmospheric Radiation Measurement (ARM) Program Mobile Facility (AMF) at Niamey airport, while the fluxes at the top of the atmosphere (TOA) are from the Geostationary Earth Radiation Budget (GERB) instrument on the Meteosat-8 satellite. The data are analyzed as daily averages, in order to minimize sampling differences between the surface and top of atmosphere instruments, while retaining the synoptic and seasonal changes that are the main focus of this study. A cloud mask is used to identify days with cloud versus those with predominantly clear skies. The influence of temperature, water vapor, aerosols, and clouds is investigated. Aerosols are ubiquitous throughout the year and have a significant impact on both the shortwave and longwave fluxes. The large and systematic seasonal changes in temperature and column integrated water vapor (CWV) through the dry and wet seasons are found to exert strong influences on the longwave fluxes. These influences are often in opposition to each other, because the highest temperatures occur at the end of the dry season when the CWV is lowest, while in the wet season the lowest temperatures are associated with the highest values of CWV. Apart from aerosols, the shortwave fluxes are also affected by clouds and by the seasonal changes in CWV. The fluxes are combined to provide estimates of the divergence of radiation across the atmosphere throughout 2006. The longwave divergence shows a relatively small variation through the year, because of a partial compensation between the seasonal variations in the outgoing longwave radiation (OLR) and surface net longwave radiation. A simple model of the greenhouse effect is used to interpret this result in terms of the dependence of the normalized greenhouse effect at the TOA and of the effective emissivity of the atmosphere at the surface on the CWV. It is shown that, as the CWV increases, the atmosphere loses longwave energy to the surface with about the same increasing efficiency with which it traps the OLR. When combined with the changes in temperature, this maintains the atmospheric longwave divergence within the narrow range that is observed. The shortwave divergence is mainly determined by the CWV and aerosol loadings and the effect of clouds is much smaller than on the component fluxes.

The potential to narrow uncertainty in projections of stratospheric ozone over the 21st century

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an "ensemble of opportunity" of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21st century, up-to and after the time when ozone concentrations return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

Chemical composition and in vitro fermentation of tannin-rich tree fruits

Dry and mature tree fruits are a potential source of protein for goats in the semi-arid areas of southern Africa, but their chemical composition and feeding value is largely unknown. This study presents the chemical composition and in vitro fermentation of indehiscent whole fruits and separated seed and hull fractions from Acacia nilotica, Acacia erubescens, Acacia sieberiana, Acacia erioloba, Piliostigma thonningii and Dichrostachys cinerea trees. Results indicate that the N contents of whole fruits ranged between 13.5 g/kg DM (A. nilotica) and 27.1 g/kg DM (A. erubescens). Seeds had a higher N content than hulls for all tree species. A. nilotica, D. cinerea and P thonningii fruits had high levels of extractable phenolics (758, 458 and 299 g/kg DM, respectively). Soluble phenolics (SPh) and ytterbium precipitable phenolics (YbPh) levels were negatively correlated to in vitro gas production but positively correlated to in vitro organic matter degradability (iOMD). Partition factors for whole fruits at 48 h ranged between 3.6 mg/ml for A. erioloba and 7.8 mg/ml for A. nilotica. Seeds of A. erioloba, A. erubescens and P thonningii were consistently fermented more efficiently throughout the incubation period compared to their whole fruits or hulls. Estimating in vitro degradability of phenolic-rich substrates through filtration procedures can give erroneous results due to the loss of soluble phenolics, which are not necessarily degradable. The feeding value of fruits from D. cinerea and A. nilotica tree species may be reduced due to the presence of high levels of phenolics. (C) 2007 Elsevier B.V. All rights reserved.

In vitro biological activity of tannins from Acacia and other tree fruits: Correlations with colorimetric and gravimetric phenolic assays

This Study was designed to investigate impact of tannins on in vitro ruminal fermentation parameters as well as relationships between concentration and in vitro biological activity of tannins present in tree fruits. Dry and mature fruits of known phenolic content harvested from Acacia nilotica, A. erubescens, A. erioloba, A. sieberiana, Piliostigima thonningii and Dichrostachys cinerea tree species were fermented with rumen fluid in vitro with or without polyethylene glycol (PEG). Correlation between in vitro biological activity and phenolic concentration was determined. Polyethylene glycol inclusion increased Cumulative gas production from all fruit substrates. The largest Increase (225%) after 48 h incubation was observed in D. cinerea fruits while the least (12.7%) increase was observed in A. erubescens fruits. Organic matter degradability (48 h) was increased by PEG inclusion for all tree species except A. erubescens and P. thonningii. For D. cinerea fruits, colorimetric assays were poorly correlated to Increases In gas production due to PEG treatment. Ytterbium precipitable phenolics (YbPh) were also poorly correlated with response to PEG for A. erioloba and P. thonningii fruits. However, YbPh were strongly and positively correlated to the increase In Cumulative gas production due to PEG for A. erubescens and A. nilotica. Folin-Ciocalteau assayed phenolics (SPh) were not correlated to response to PEG in P. thonningii and A. sieberiana. It was Concluded that the PEG effect oil in vitro fermentation was closely related to some measures of phenolic concentration but the relationships varied with tree species.

Design and analysis of climate model experiments for the efficient estimation of anthropogenic signals

Presented herein is an experimental design that allows the effects of several radiative forcing factors on climate to be estimated as precisely as possible from a limited suite of atmosphere-only general circulation model (GCM) integrations. The forcings include the combined effect of observed changes in sea surface temperatures, sea ice extent, stratospheric (volcanic) aerosols, and solar output, plus the individual effects of several anthropogenic forcings. A single linear statistical model is used to estimate the forcing effects, each of which is represented by its global mean radiative forcing. The strong colinearity in time between the various anthropogenic forcings provides a technical problem that is overcome through the design of the experiment. This design uses every combination of anthropogenic forcing rather than having a few highly replicated ensembles, which is more commonly used in climate studies. Not only is this design highly efficient for a given number of integrations, but it also allows the estimation of (nonadditive) interactions between pairs of anthropogenic forcings. The simulated land surface air temperature changes since 1871 have been analyzed. The changes in natural and oceanic forcing, which itself contains some forcing from anthropogenic and natural influences, have the most influence. For the global mean, increasing greenhouse gases and the indirect aerosol effect had the largest anthropogenic effects. It was also found that an interaction between these two anthropogenic effects in the atmosphere-only GCM exists. This interaction is similar in magnitude to the individual effects of changing tropospheric and stratospheric ozone concentrations or to the direct (sulfate) aerosol effect. Various diagnostics are used to evaluate the fit of the statistical model. For the global mean, this shows that the land temperature response is proportional to the global mean radiative forcing, reinforcing the use of radiative forcing as a measure of climate change. The diagnostic tests also show that the linear model was suitable for analyses of land surface air temperature at each GCM grid point. Therefore, the linear model provides precise estimates of the space time signals for all forcing factors under consideration. For simulated 50-hPa temperatures, results show that tropospheric ozone increases have contributed to stratospheric cooling over the twentieth century almost as much as changes in well-mixed greenhouse gases.

Mozzarella-type curd made from buffalo, cows’ and ultrafiltered cows’ milk. 2. Physicochemical properties, curd yield and quality, casein fractions & micelles size

Buffalo milk contains (40–60 %) more protein, fat and calcium than cows’ milk. These constituents were enhanced by ultrafiltration (UF) of cows’ milk to give a product with similar levels to those found in the buffalo milk. Mozzarella-type curd was made from buffalo, cows’ and UF cows’ milk to compare the overall curd yield and quality. The curd yield on both dry and wet weight basis, curd moisture content and overall curd fat retention were found to be higher in the UF cows’ milk than for either the buffalo or the cows’ milk preparations. The minimum whey fat losses occurred in the UF cows’ curd when compared to the cows’ and the buffalo curd. The whey protein losses were found to be higher in the UF cows’ curd than those for the buffalo and the cows’ curds. The total mineral content of the curd was also higher in the UF cows’ milk than that found in either the buffalo or the cows’ milk. SEM micrographs showed that casein micelles sizes were different in the two different types of milk. Casein micelles were also observed to be deformed in the UF cows’ milk samples. UF cows’ milk contained higher amounts of both the αs1- and αs2-casein moieties than either the buffalo or the cows’ milk. Buffalo milk was found to contain a higher concentration of β-casein than either the UF cows’ or untreated cows’ milk samples. Gel strength was found to be higher in the resultant buffalo curd than for curds made from either native cows’ milk or those made from UF cows’ milk. The mineral distribution was also different in the three different types of bovine milk, measured by energy-dispersive X-ray (EDX) analysis. Differences in the curd quality observed between the buffalo and the cows’ milk appear to result from the differences in casein composition and overall micelle structure, rather than casein concentration alone.

Evaluation of water and energy budgets in regional climate models applied over Europe

This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.

Global and regional temperature-change potentials for near-term climate forcers

We examine the climate effects of the emissions of near-term climate forcers (NTCFs) from 4 continental regions (East Asia, Europe, North America and South Asia) using radiative forcing from the task force on hemispheric transport of air pollution source-receptor global chemical transport model simulations. These simulations model the transport of 3 aerosol species (sulphate, particulate organic matter and black carbon) and 4 ozone precursors (methane, nitric oxides (NOx), volatile organic compounds and carbon monoxide). From the equilibrium radiative forcing results we calculate global climate metrics, global warming potentials (GWPs) and global temperature change potentials (GTPs) and show how these depend on emission region, and can vary as functions of time. For the aerosol species, the GWP(100) values are −37±12, −46±20, and 350±200 for SO2, POM and BC respectively for the direct effects only. The corresponding GTP(100) values are −5.2±2.4, −6.5±3.5, and 50±33. This analysis is further extended by examining the temperature-change impacts in 4 latitude bands. This shows that the latitudinal pattern of the temperature response to emissions of the NTCFs does not directly follow the pattern of the diagnosed radiative forcing. For instance temperatures in the Arctic latitudes are particularly sensitive to NTCF emissions in the northern mid-latitudes. At the 100-yr time horizon the ARTPs show NOx emissions can have a warming effect in the northern mid and high latitudes, but cooling in the tropics and Southern Hemisphere. The northern mid-latitude temperature response to northern mid-latitude emissions of most NTCFs is approximately twice as large as would be implied by the global average.

Final warming of the Southern Hemisphere polar vortex in high- and low-top CMIP5 models

The final warming date of the polar vortex is a key component of Southern Hemisphere stratospheric and tropospheric variability in spring and summer. We examine the effect of external forcings on Southern Hemisphere final warming date, and the sensitivity of any projected changes to model representation of the stratosphere. Final warming date is calculated using a temperature-based diagnostic for ensembles of high- and low-top CMIP5 models, under the CMIP5 historical, RCP4.5, and RCP8.5 forcing scenarios. The final warming date in the models is generally too late in comparison with those from reanalyses: around two weeks too late in the low-top ensemble, and around one week too late in the high-top ensemble. Ensemble Empirical Mode Decomposition (EEMD) is used to analyse past and future change in final warming date. Both the low- and high-top ensemble show characteristic behaviour expected in response to changes in greenhouse gas and stratospheric ozone concentrations. In both ensembles, under both scenarios, an increase in final warming date is seen between 1850 and 2100, with the latest dates occurring in the early twenty-first century, associated with the minimum in stratospheric ozone concentrations in this period. However, this response is more pronounced in the high-top ensemble. The high-top models show a delay in final warming date in RCP8.5 that is not produced by the low-top models, which are shown to be less responsive to greenhouse gas forcing. This suggests that it may be necessary to use stratosphere resolving models to accurately predict Southern Hemisphere surface climate change.