MetUM-GOML: a near-globally coupled atmosphere–ocean-mixed-layer model

Well-resolved air–sea interactions are simulated in a new ocean mixed-layer, coupled configuration of the Met Office Unified Model (MetUM-GOML), comprising the MetUM coupled to the Multi-Column K Profile Parameterization ocean (MC-KPP). This is the first globally coupled system which provides a vertically resolved, high near-surface resolution ocean at comparable computational cost to running in atmosphere-only mode. As well as being computationally inexpensive, this modelling framework is adaptable– the independent MC-KPP columns can be applied selectively in space and time – and controllable – by using temperature and salinity corrections the model can be constrained to any ocean state. The framework provides a powerful research tool for process-based studies of the impact of air–sea interactions in the global climate system. MetUM simulations have been performed which separate the impact of introducing inter- annual variability in sea surface temperatures (SSTs) from the impact of having atmosphere–ocean feedbacks. The representation of key aspects of tropical and extratropical variability are used to assess the performance of these simulations. Coupling the MetUM to MC-KPP is shown, for example, to reduce tropical precipitation biases, improve the propagation of, and spectral power associated with, the Madden–Julian Oscillation and produce closer-to-observed patterns of springtime blocking activity over the Euro-Atlantic region.

JULES-crop: a parametrisation of crops in the Joint UK Land Environment Simulator

Studies of climate change impacts on the terrestrial biosphere have been completed without recognition of the integrated nature of the biosphere. Improved assessment of the impacts of climate change on food and water security requires the development and use of models not only representing each component but also their interactions. To meet this requirement the Joint UK Land Environment Simulator (JULES) land surface model has been modified to include a generic parametrisation of annual crops. The new model, JULES-crop, is described and evaluation at global and site levels for the four globally important crops; wheat, soybean, maize and rice. JULES-crop demonstrates skill in simulating the inter-annual variations of yield for maize and soybean at the global and country levels, and for wheat for major spring wheat producing countries. The impact of the new parametrisation, compared to the standard configuration, on the simulation of surface heat fluxes is largely an alteration of the partitioning between latent and sensible heat fluxes during the later part of the growing season. Further evaluation at the site level shows the model captures the seasonality of leaf area index, gross primary production and canopy height better than in the standard JULES. However, this does not lead to an improvement in the simulation of sensible and latent heat fluxes. The performance of JULES-crop from both an Earth system and crop yield model perspective is encouraging. However, more effort is needed to develop the parametrisation of the model for specific applications. Key future model developments identified include the introduction of processes such as irrigation and nitrogen limitation which will enable better representation of the spatial variability in yield.

Improved simulation of fire–vegetation interactions in the Land surface Processes and eXchanges dynamic global vegetation model (LPX-Mv1)

The Land surface Processes and eXchanges (LPX) model is a fire-enabled dynamic global vegetation model that performs well globally but has problems representing fire regimes and vegetative mix in savannas. Here we focus on improving the fire module. To improve the representation of ignitions, we introduced a reatment of lightning that allows the fraction of ground strikes to vary spatially and seasonally, realistically partitions strike distribution between wet and dry days, and varies the number of dry days with strikes. Fuel availability and moisture content were improved by implementing decomposition rates specific to individual plant functional types and litter classes, and litter drying rates driven by atmospheric water content. To improve water extraction by grasses, we use realistic plant-specific treatments of deep roots. To improve fire responses, we introduced adaptive bark thickness and post-fire resprouting for tropical and temperate broadleaf trees. All improvements are based on extensive analyses of relevant observational data sets. We test model performance for Australia, first evaluating parameterisations separately and then measuring overall behaviour against standard benchmarks. Changes to the lightning parameterisation produce a more realistic simulation of fires in southeastern and central Australia. Implementation of PFT-specific decomposition rates enhances performance in central Australia. Changes in fuel drying improve fire in northern Australia, while changes in rooting depth produce a more realistic simulation of fuel availability and structure in central and northern Australia. The introduction of adaptive bark thickness and resprouting produces more realistic fire regimes in Australian savannas. We also show that the model simulates biomass recovery rates consistent with observations from several different regions of the world characterised by resprouting vegetation. The new model (LPX-Mv1) produces an improved simulation of observed vegetation composition and mean annual burnt area, by 33 and 18% respectively compared to LPX.

The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations

We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0.

Including an ocean carbon cycle model into iLOVECLIM (v1.0)

The atmospheric carbon dioxide concentration plays a crucial role in the radiative balance and as such has a strong influence on the evolution of climate. Because of the numerous interactions between climate and the carbon cycle, it is necessary to include a model of the carbon cycle within a climate model to understand and simulate past and future changes of the carbon cycle. In particular, natural variations of atmospheric CO2 have happened in the past, while anthropogenic carbon emissions are likely to continue in the future. To study changes of the carbon cycle and climate on timescales of a few hundred to a few thousand years, we have included a simple carbon cycle model into the iLOVECLIM Earth System Model. In this study, we describe the ocean and terrestrial biosphere carbon cycle models and their performance relative to observational data. We focus on the main carbon cycle variables including the carbon isotope ratios δ13C and the Δ14C. We show that the model results are in good agreement with modern observations both at the surface and in the deep ocean for the main variables, in particular phosphates, dissolved inorganic carbon and the carbon isotopes.

The Met Office Global Coupled model 2.0 (GC2) configuration

The latest coupled configuration of the Met Office Unified Model (Global Coupled configuration 2, GC2) is presented. This paper documents the model components which make up the configuration (although the scientific description of these components is detailed elsewhere) and provides a description of the coupling between the components. The performance of GC2 in terms of its systematic errors is assessed using a variety of diagnostic techniques. The configuration is intended to be used by the Met Office and collaborating institutes across a range of timescales, with the seasonal forecast system (GloSea5) and climate projection system (HadGEM) being the initial users. In this paper GC2 is compared against the model currently used operationally in those two systems. Overall GC2 is shown to be an improvement on the configurations used currently, particularly in terms of modes of variability (e.g. mid-latitude and tropical cyclone intensities, the Madden–Julian Oscillation and El Niño Southern Oscillation). A number of outstanding errors are identified with the most significant being a considerable warm bias over the Southern Ocean and a dry precipitation bias in the Indian and West African summer monsoons. Research to address these is ongoing.

Ice-sheet configuration in the CMIP5/PMIP3 Last Glacial Maximum experiments

We describe the creation of a data set describing changes related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and elevation of ice-free land at the Last Glacial Maximum (LGM), which were used in LGM experiments conducted as part of the fifth phase of the Coupled Modelling Intercomparison Project (CMIP5) and the third phase of the Palaeoclimate Modelling Intercomparison Project (PMIP3). The CMIP5/PMIP3 data sets were created from reconstructions made by three different groups, which were all obtained using a model-inversion approach but differ in the assumptions used in the modelling and in the type of data used as constraints. The ice-sheet extent in the Northern Hemisphere (NH) does not vary substantially between the three individual data sources. The difference in the topography of the NH ice sheets is also moderate, and smaller than the differences between these reconstructions (and the resultant composite reconstruction) and ice-sheet reconstructions used in previous generations of PMIP. Only two of the individual reconstructions provide information for Antarctica. The discrepancy between these two reconstructions is larger than the difference for the NH ice sheets, although still less than the difference between the composite reconstruction and previous PMIP ice-sheet reconstructions. Although largely confined to the ice-covered regions, differences between the climate response to the individual LGM reconstructions extend over the North Atlantic Ocean and Northern Hemisphere continents, partly through atmospheric stationary waves. Differences between the climate response to the CMIP5/PMIP3 composite and any individual ice-sheet reconstruction are smaller than those between the CMIP5/PMIP3 composite and the ice sheet used in the last phase of PMIP (PMIP2).

A global empirical system for probabilistic seasonal climate prediction

Preparing for episodes with risks of anomalous weather a month to a year ahead is an important challenge for governments, non-governmental organisations, and private companies and is dependent on the availability of reliable forecasts. The majority of operational seasonal forecasts are made using process-based dynamical models, which are complex, computationally challenging and prone to biases. Empirical forecast approaches built on statistical models to represent physical processes offer an alternative to dynamical systems and can provide either a benchmark for comparison or independent supplementary forecasts. Here, we present a simple empirical system based on multiple linear regression for producing probabilistic forecasts of seasonal surface air temperature and precipitation across the globe. The global CO2-equivalent concentration is taken as the primary predictor; subsequent predictors, including large-scale modes of variability in the climate system and local-scale information, are selected on the basis of their physical relationship with the predictand. The focus given to the climate change signal as a source of skill and the probabilistic nature of the forecasts produced constitute a novel approach to global empirical prediction. Hindcasts for the period 1961–2013 are validated against observations using deterministic (correlation of seasonal means) and probabilistic (continuous rank probability skill scores) metrics. Good skill is found in many regions, particularly for surface air temperature and most notably in much of Europe during the spring and summer seasons. For precipitation, skill is generally limited to regions with known El Niño–Southern Oscillation (ENSO) teleconnections. The system is used in a quasi-operational framework to generate empirical seasonal forecasts on a monthly basis.

SHIMMER (1.0): a novel mathematical model for microbial and biogeochemical dynamics in glacier forefield ecosystems

SHIMMER (Soil biogeocHemIcal Model for Microbial Ecosystem Response) is a new numerical modelling framework designed to simulate microbial dynamics and biogeochemical cycling during initial ecosystem development in glacier forefield soils. However, it is also transferable to other extreme ecosystem types (such as desert soils or the surface of glaciers). The rationale for model development arises from decades of empirical observations in glacier forefields, and enables a quantitative and process focussed approach. Here, we provide a detailed description of SHIMMER, test its performance in two case study forefields: the Damma Glacier (Switzerland) and the Athabasca Glacier (Canada) and analyse sensitivity to identify the most sensitive and unconstrained model parameters. Results show that the accumulation of microbial biomass is highly dependent on variation in microbial growth and death rate constants, Q10 values, the active fraction of microbial biomass and the reactivity of organic matter. The model correctly predicts the rapid accumulation of microbial biomass observed during the initial stages of succession in the forefields of both the case study systems. Primary production is responsible for the initial build-up of labile substrate that subsequently supports heterotrophic growth. However, allochthonous contributions of organic matter, and nitrogen fixation, are important in sustaining this productivity. The development and application of SHIMMER also highlights aspects of these systems that require further empirical research: quantifying nutrient budgets and biogeochemical rates, exploring seasonality and microbial growth and cell death. This will lead to increased understanding of how glacier forefields contribute to global biogeochemical cycling and climate under future ice retreat.

Evaluation of the Plant–Craig stochastic convection scheme (v2.0) in the ensemble forecasting system MOGREPS-R (24 km) based on the Unified Model (v7.3)

The Plant–Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant–Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant–Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

PMIP4-CMIP6: the contribution of the Paleoclimate Modelling Intercomparison Project to CMIP6

The goal of the Palaeoclimate Modelling Intercomparison Project (PMIP) is to understand the response of the climate system to changes in different climate forcings and to feedbacks. Through comparison with observations of the environmental impacts of these climate changes, or with climate reconstructions based on physical, chemical or biological records, PMIP also addresses the issue of how well state-of-the-art models simulate climate changes. Palaeoclimate states are radically different from those of the recent past documented by the instrumental record and thus provide an out-of-sample test of the models used for future climate projections and a way to assess whether they have the correct sensitivity to forcings and feedbacks. Five distinctly different periods have been selected as focus for the core palaeoclimate experiments that are designed to contribute to the objectives of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This manuscript describes the motivation for the choice of these periods and the design of the numerical experiments, with a focus upon their novel features compared to the experiments performed in previous phases of PMIP and CMIP as well as the benefits of common analyses of the models across multiple climate states. It also describes the information needed to document each experiment and the model outputs required for analysis and benchmarking.

Relações planta-animal num sistema de integração lavoura-pecuária

Este trabalho buscou estudar um sistema de integração lavoura-pecuária com diferentes alturas de pastos, no período de inverno, e seus reflexos sobre a dinâmica da pastagem e o desempenho animal. O experimento foi conduzido em uma pastagem de aveia + azevém manejada sob diferentes intensidades de pastejo. O delineamento experimental foi de blocos ao acaso com quatro tratamentos (10, 20, 30 e 40 cm de altura de manejo) e três repetições. Utilizaramse terneiros de corte de cruzamento industrial com idade e peso médio inicial de 10 meses e 210 kg, respectivamente. O método de pastejo foi contínuo com lotação variável. A adubação de base foi de 400 kg/ha de superfosfato simples e de 90 kg/ha de N em cobertura. As alturas do pasto afetaram a massa de forragem (MF), onde para cada cm de aumento na altura acima de 10 cm, houve incremento na matéria seca da pastagem em cerca de 86 kg/ha de MS. O aumento no ganho médio diário (GMD) foi condicionado pelo incremento na qualidade e/ou na quantidade de forragem disponível, e o modelo de resposta do GMD em relação às alturas do pasto, resultou em valores de 0,73 e 1,14 kg/animal/dia nos tratamentos de maior e menor GMD, respectivamente, que foram de 10 cm e 30 cm de altura. No rendimento de carcaça não houve diferença (P>0,05) entre os tratamentos, uma vez que todos os valores ficaram em torno de 51%. Quanto ao peso de carcaça quente e fria, peso de costilhar, escore de condição corporal e grau de acabamento, observou-se comportamento muito similar à evolução do ganho médio diário dos animais.

Desempenho de novilhas expostas à reprodução aos 14/15 meses de idade

Foram conduzidos dois experimentos para avaliar o comportamento reprodutivo de novilhas de corte pertencentes a quatro grupos genéticos (Hereford, ½ Nelore½ Hereford; ¼ Nelore ¾ Hereford; ½ Angus ½ Hereford) acasaladas aos 14/15 meses de idade. No experimento iniciado em 2001 foram utilizadas quatro alternativas de alimentação no outono/inverno: a) suplementação do campo nativo (CN) a 1,5% do peso vivo (PV) com ração comercial (RC) contendo 14% de proteína bruta (PB) e 68% de NDT (S68); b) suplementação do CN a 1,5% do PV com RC contendo 14% de PB e 75% de NDT (S75); c) suplementação em pastagem de azevém (Lolium multiflorum Lan), trevo branco (Trifolium repens) e cornichão (Lotus corniculatus) a 0,5% do PV de milho em grão (SPAS);d) confinamento a céu aberto com silagem de milho e 0,7% do PV de ração (CON). Após, em um só grupo, permaneceram em pastagem de azevém e aveia (Avena sativa) até o início do período reprodutivo. Durante o período de aplicação dos tratamentos alimentares o ganho médio diário (GMD) do tratamento SPAS e CON (0,755 vs 0,784 kg/dia) não apresentou diferença significativa (P>0,05), nem os S68 e S75 (0,511 vs 0,489 kg/dia). Durante o período conjunto em pastagem não foi determinada diferença significativa (P>0,05) no GMD entre S68, S75 e SPAS, nem entre S68, SPAS e CON. O peso ao início do período reprodutivo foi de 233,7 e 232,3 kg para S68 e S75 (P>0,05) e de 260,6 e 254,3 kg para SPAS e CON (P>0,05). As taxas de prenhez do SPAS e COM não apresentaram diferença entre si (62,16% e 53,86%), mas foram significativamente mais elevadas (P<0,05) em relação às S68 e S75 (20,51% e 25,64%) as quais também não diferiram entre si (P>0,05). As novilhas que conceberam foram mais pesadas, apresentaram maiores GMD de peso, condição corporal (CC) e escore trato reprodutivo (ETR) ao início do período reprodutivo em relação às que não conceberam. No experimento iniciado em 2002 foram utilizadas três alternativas de alimentação no outono/inverno: a) suplementação do CN com RC (SUR) a 1% do PV, contendo 14% de PB e 75% de NDT; b) suplementação do CN com farelo de arroz (SUFA) a 0,5% PV; c) pastejo contínuo em pastagem cultivada de azevém (PAST) e aveia. O tratamento PAST (0,478 kg/dia) apresentou GMD de peso mais elevado (P<0,05), seguido do tratamento SUR (0,326 kg/dia) superior (P<0,05) ao tratamento SUFA (0,100 kg/dia). Durante o período em conjunto em pastagem de azevém e aveia os tratamentos SUFA e PAST não apresentaram diferença significativa no GMD de peso (P>0,05) mas, foram superiores (P<0,05) ao tratamento SUR. O peso ao início do período reprodutivo foi de 265,4; 236,7 e 222,6 kg para os tratamentos PAST, SUR e SUFA, respectivamente (P<0,05). A taxa de prenhez do tratamento PAST (61,36%) foi superior (P<0,05) a dos tratamentos SUR e SUFA (20,0 e 22,73%; P>0,05). As novilhas que conceberam foram as mais pesadas e mais velhas, do início dos tratamentos alimentares até o final do período reprodutivo.

Efeito da suplementação protéica sobre o desempenho de novilhos de corte em pastagem de forrageira tropical

Foi realizado um experimento para avaliar o efeito da suplementação com sal proteinado e sal mineralizado sobre o desempenho de novilhos de corte em pastejo de capim estrela (Cynodon nlemfuensis Vanderyst) (CE) diferido. Os suplementos avaliados foram produtos comerciais disponíveis no mercado: sal mineralizado (SM); sal proteinado (SP) e sem suplementação (SS). O experimento teve uma duração de 84 dias e foram utilizados 36 novilhos Braford com um peso médio de 267 kg, distribuídos aleatoriamente em seis piquetes com área de 4 ha, num delineamento completamente casualizado. A pastagem de CE apresentou valores médios de proteína bruta (PB), fibra insolúvel em detergente neutro (FDN) e digestibilidade in vitro da matéria orgânica (DIVMO) de 7,6, 78,7 e 45,3 % respectivamente e não sendo observadas diferenças significativas entre tratamentos (p>0,05). Entretanto, observou-se que as condições climáticas observadas durante o experimento podem ter afetado os parâmetros de avaliação da qualidade e disponibilidade da pastagem. O ganho médio diário (GMD) apresentado pelos animais não mostrou diferença significativa (p>0,05) entre os diferentes tratamentos, sendo de 0,009 kg/dia para SM, 0,104 kg/dia para SP e 0,198 kg/dia para SS. Nas condições em que foi conduzido o presente experimento, a utilização de sal mineralizado e de sal proteinado não melhorou o ganho de peso por animal e por área em relação aos novilhos que não receberam suplementação.

Alternativas de manejo para os campos de cima da serra

Os experimentos, foram conduzidos em uma área de campo natural, localizada na Fazenda Boa Vista, no Município de São José dos Ausentes, região dos Campos de Cima da Serra, Rio Grande do Sul. Os objetivos foram descrever o ambiente de pastagem natural da região, determinar a capacidade de suporte destas pastagens e avaliar alternativas de manejo para o inverno. As diferentes ofertas de forragem proporcionaram aumentos lineares no massa de forragem da pastagem, na taxa de acúmulo de matéria seca e no percentual de material morto. Nos parâmetros de resposta animal, as ofertas de forragem geraram modelos lineares para carga animal (kg/ha), GMD (kg/animal/dia), para a avaliação de produção por área foi obtida resposta quadrática, sendo o maior valor 69,7 kg/ha/ano. Para o levantamento florístico constatou-se que o campo é composto basicamente por gramíneas cespitosas estivais, sendo Andropogon lateralis a espécie mais abundante. Após a aplicação de quatro níveis de oferta de forragem (5, 9, 13 e 17% do PV), observou-se um aumento na participação de cyperaceas, principalmente Bulbostylis sphaerocephala. Para a avaliação dos suplementos, observou-se diferença entre os tratamentos, sendo o pastejo temporário superior em GMD à ração e às misturas múltiplas. Verificou-se que a aplicação de calcário em superfície é recomendável e soluciona os problemas decorrentes dos altos teores de Al trocável. A introdução de espécies sobre o campo natural da região é indicada.