This paper describes the development of a technically robust climate modelling system, HadGEM3, which couples the Met Office Unified Model atmosphere component, the NEMO ocean model and the Los Alamos sea ice model (CICE) using the OASIS coupler. Details of the coupling and technical solutions of the physical model (HadGEM3-AO) are documented, in addition to a description of the configurations of the individual submodels. The paper demonstrates that the implementation of the model has resulted in accurate conservation of heat and freshwater across the model components. The model performance in early versions of this climate model is briefly described to demonstrate that the results are scientifically credible. HadGEM3-AO is the basis for a number of modelling efforts outside of the Met Office, both within the UK and internationally. This documentation of the HadGEM3-AO system provides a detailed reference for developers of HadGEM3-based climate configurations

Abstract. We describe the HadGEM2 family of climate configurations of the Met Office Unified Model, MetUM. The concept of a model "family" comprises a range of specific model configurations incorporating different levels of complexity but with a common physical framework. The HadGEM2 family of configurations includes atmosphere and ocean components, with and without a vertical extension to include a well-resolved stratosphere, and an Earth-System (ES) component which includes dynamic vegetation, ocean biology and atmospheric chemistry. The HadGEM2 physical model includes improvements designed to address specific systematic errors encountered in the previous climate configuration, HadGEM1, namely Northern Hemisphere continental temperature biases and tropical sea surface temperature biases and poor variability. Targeting these biases was crucial in order that the ES configuration could represent important biogeochemical climate feedbacks. Detailed descriptions and evaluations of particular HadGEM2 family memCorrespondence to: G. M. Martin (gill.martin@metoffice.gov.uk) bers are included in a number of other publications, and the discussion here is limited to a summary of the overall performance using a set of model metrics which compare the way in which the various configurations simulate present-day climate and its variability.

[1] Measurements of the top-of-the-atmosphere outgoing longwave radiation (OLR) for July 2003 from Meteosat-7 are used to assess the performance of the numerical weather prediction version of the Met Office Unified Model. A significant difference is found over desert regions of northern Africa where the model emits too much OLR by up to 35 Wm À2 in the monthly mean. By cloud-screening the data we find an error of up to 50 Wm À2 associated with cloud-free areas, which suggests an error in the model surface temperature, surface emissivity, or atmospheric transmission. By building up a physical model of the radiative properties of mineral dust based on in situ, and surface-based and satellite remote sensing observations we show that the most plausible explanation for the discrepancy in OLR is due to the neglect of mineral dust in the model. The calculations suggest that mineral dust can exert a longwave radiative forcing by as much as 50 Wm À2 in the monthly mean for 1200 UTC in cloud-free regions, which accounts for the discrepancy between the model and the Meteosat-7 observations. This suggests that inclusion of the radiative effects of mineral dust will lead to a significant improvement in the radiation balance of numerical weather prediction models with subsequent improvements in performance.

This paper describes the development of a technically robust climate modelling system, HadGEM3, which couples the Met Office Unified Model atmosphere component, the NEMO ocean model and the Los Alamos sea ice model (CICE) using the OASIS coupler. Details of the coupling and technical solutions are documented in the paper in addition to a description of the configurations of the individual submodels. The paper demonstrates that the implementation of the model has resulted in accurate conservation of heat and freshwater across the model components. The model performance in early versions of this climate model is briefly described to demonstrate that the results are scientifically credible. HadGEM3 is the basis for a number of modelling efforts outside of the Met Office, both within the UK and internationally. This documentation of the HadGEM3 system provides a detailed reference for developers of HadGEM3-based climate configurations

Abstract.A modification to the standard bending-angle correction used in GPS radio occultation (GPS-RO) is proposed. The modified approach should reduce systematic residual ionospheric errors in GPS radio occultation climatologies. A new second-order term is introduced in order to account for a known source of systematic error, which is generally neglected. The new term has the form κ(a) × (α L1 (a) − α L2 (a)) 2 , where a is the impact parameter and (α L1 , α L2 ) are the L1 and L2 bending angles, respectively. The variable κ is a weak function of the impact parameter, a, but it does depend on a priori ionospheric information. The theoretical basis of the new term is examined. The sensitivity of κ to the assumed ionospheric parameters is investigated in one-dimensional simulations, and it is shown that κ 10-20 rad −1 . We note that the current implicit assumption is κ = 0, and this is probably adequate for numerical weather prediction applications. However, the uncertainty in κ should be included in the uncertainty estimates for the geophysical climatologies produced from GPS-RO measurements. The limitations of the new ionospheric correction when applied to CHAMP (Challenging Minisatellite Payload) measurements are noted. These arise because of the assumption that the refractive index is unity at the satellite, made when deriving bending angles from the Doppler shift values.

Abstract. In this study, a new model was explored which corrects for higher order ionospheric residuals in Global Positioning System (GPS) radio occultation (RO) data. Recently, the theoretical basis of this new "residual ionospheric error model" has been outlined (Healy and Culverwell, 2015). The method was tested in simulations with a one-dimensional model ionosphere.The proposed new model for computing the residual ionospheric error is the product of two factors, one of which expresses its variation from profile to profile and from time to time in terms of measurable quantities (the L1 and L2 bending angles), while the other describes the weak variation with altitude. A simple integral expression for the residual error (Vorob'ev and Krasil'nikova, 1994) has been shown to be in excellent numerical agreement with the exact value, for a simple Chapman layer ionosphere. In this case, the "altitudinal" element of the residual error varies (decreases) by no more than about 25 % between ∼ 10 and ∼ 100 km for physically reasonable Chapman layer parameters. For other simple model ionospheres the integral can be evaluated exactly, and results are in reasonable agreement with those of an equivalent Chapman layer.In this follow-up study the overall objective was to explore the validity of the new residual ionospheric error model for more detailed simulations, based on modeling through a complex three-dimensional ionosphere.The simulation study was set up, simulating day and night GPS RO profiles for the period of a solar cycle with and without an ionosphere. The residual ionospheric error was studied, the new error model was tested, and temporal and spatial variations of the model were investigated. The model performed well in the simulation study, capturing the temporal variability of the ionospheric residual. Although it was not possible, due to high noise of the simulated bending-angle profiles at mid-to high latitudes, to perform a thorough latitudinal investigation of the performance of the model, first positive and encouraging results were found at low latitudes. Furthermore, first application tests of the model on the data showed a reduction in temperature level of the ionospheric residual at 40 km from about −2.2 to −0.2 K.

We describe a new methodology for comparing satellite radiation budget data with a numerical weather prediction (NWP) model. This is applied to data from the Geostationary Earth Radiation Budget (GERB) instrument on Meteosat‐8. The methodology brings together, in near‐real time, GERB broadband shortwave and longwave fluxes with simulations based on analyses produced by the Met Office global NWP model. Results for the period May 2003 to February 2005 illustrate the progressive improvements in the data products as various initial problems were resolved. In most areas the comparisons reveal systematic errors in the model's representation of surface properties and clouds, which are discussed elsewhere. However, for clear‐sky regions over the oceans the model simulations are believed to be sufficiently accurate to allow the quality of the GERB fluxes themselves to be assessed and any changes in time of the performance of the instrument to be identified. Using model and radiosonde profiles of temperature and humidity as input to a single‐column version of the model's radiation code, we conduct sensitivity experiments which provide estimates of the expected model errors over the ocean of about ±5–10 W m−2 in clear‐sky outgoing longwave radiation (OLR) and ±0.01 in clear‐sky albedo. For the more recent data the differences between the observed and modeled OLR and albedo are well within these error estimates. The close agreement between the observed and modeled values, particularly for the most recent period, illustrates the value of the methodology. It also contributes to the validation of the GERB products and increases confidence in the quality of the data, prior to their release.

Sum.m~ry. A the?ry of compositional viscous fingering with no adjustable parameters that reduces to the Todd and Longstaff model for mIscIble floo?s IS prese~ted. The t~~ory gi~es excellent predictive agreement with simulation results for a wide range of recovery ?rocesses. The hlgh-r~s~lutlOn compOSItIonal SImulations used to validate the theory are the first to resolve viscous fingering adequately III flow other than mIscIble flow.