Mechanisms underlying temperature extremes in Iberia: a Lagrangian perspective

The mechanisms underlying the occurrence of temperature extremes in Iberia are analysed considering a Lagrangian perspective of the atmospheric flow, using 6-hourly ERA-Interim reanalysis data for the years 1979–2012. Daily 2-m minimum temperatures below the 1st percentile and 2-m maximum temperatures above the 99th percentile at each grid point over Iberia are selected separately for winter and summer. Four categories of extremes are analysed using 10-d backward trajectories initialized at the extreme temperature grid points close to the surface: winter cold (WCE) and warm extremes (WWE), and summer cold (SCE) and warm extremes (SWE). Air masses leading to temperature extremes are first transported from the North Atlantic towards Europe for all categories. While there is a clear relation to large-scale circulation patterns in winter, the Iberian thermal low is important in summer. Along the trajectories, air mass characteristics are significantly modified through adiabatic warming (air parcel descent), upper-air radiative cooling and near-surface warming (surface heat fluxes and radiation). High residence times over continental areas, such as over northern-central Europe for WCE and, to a lesser extent, over Iberia for SWE, significantly enhance these air mass modifications. Near-surface diabatic warming is particularly striking for SWE. WCE and SWE are responsible for the most extreme conditions in a given year. For WWE and SCE, strong temperature advection associated with important meridional air mass transports are the main driving mechanisms, accompanied by comparatively minor changes in the air mass properties. These results permit a better understanding of mechanisms leading to temperature extremes in Iberia.

Influência do relevo e erosão na variabilidade espacial de um Latossolo em Jaboticabal (SP)

Com o objetivo de avaliar as relações entre formas de paisagem e erosão em um Latossolo Vermelho eutroférrico (LVef) por meio de técnicas geoestatísticas na Fazenda Santa Isabel, município de Jaboticabal (SP), identificaram-se modelos de paisagem côncavo e linear. Coletaram-se amostras de solo em ambas as formas de paisagem em uma malha regular espaçada de 50 x 50 m em sete transeções na profundidade de 0,00-0,20 m, totalizando 412 pontos em 93 ha. Determinou-se a erodibilidade dos solos pelo método indireto, granulometria, bem como o carbono orgânico para cada ponto da malha, sendo esses valores usados na estimativa do potencial natural de erosão (pne). Os dados obtidos das propriedades de solo, erodibilidade e potencial natural de erosão foram analisados por meio de estatística descritiva e geoestatística com a modelagem de semivariogramas, sendo este utilizado para a confecção de mapas de krigagem. Segundo os resultados, as propriedades do solo e do pne apresentaram maior variabilidade espacial na pedoforma côncava, apesar de esta forma de paisagem apresentar menores perdas de solo por erosão e menor variabilidade espacial da erodibilidade. Deste modo, conclui-se que a erosão não adicionou variabilidade espacial para as propriedades do solo na mesma magnitude do relevo.

Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas

As formas do relevo podem ser indicadores da variação dos atributos do solo, pois essa variabilidade é causada por pequenas alterações do declive que afetam os processos pedogenéticos bem como o transporte e o armazenamento de água no perfil do solo. O trabalho foi desenvolvido em Catanduva (SP), com o objetivo de caracterizar a variabilidade espacial de atributos do solo e fatores de erosão em diferentes pedoformas sob cultivo de cana-de-açúcar. de acordo com o modelo de Troeh classificou-se as formas do relevo em duas pedoformas, côncava e convexa. Com a utilização de um DGPS levantaram-se as cotas altimétricas, estabelecendo-se uma malha, com intervalos regulares de 50 m, com 270 pontos na pedoforma côncava e 353 pontos na pedoforma convexa, perfazendo um total de 623 pontos, coletados na profundidade de 0,0 - 0,2 m em uma área de 200 ha. em cada ponto da malha foram determinados os atributos químicos do solo, granulometria, espessura do solo e fatores de erosão locais, tais como erosividade (R), erodibilidade (K), fator topográfico (LS), uso e manejo (C), práticas conservacionistas (P), potencial natural de erosão (PNE), perda de solo (A) e risco de erosão (RE). Os dados foram avaliados primeiramente por uma análise estatística exploratória, calculando-se a média, mediana, variância, coeficiente de variação, coeficiente de assimetria, coeficiente de curtose e teste de normalidade. Posteriormente, a dependência espacial foi verificada por meio da técnica de geoestatística utilizando-se semivariogramas. As maiores perdas de solo, risco de erosão e potencial natural de erosão e menor espessura do solo ocorreram na pedoforma convexa, indicando forte dependência espacial com a forma do relevo. A pedoforma côncava proporcionou maior variabilidade espacial, demonstrando que a forma do relevo condiciona padrões diferenciados de variabilidade. A magnitude da variabilidade dos atributos do solo é mais influenciada pela forma do relevo que pela erosão. A espessura do horizonte A+E integrado com a forma do relevo é um indicador de processos erosivos para classe de Argissolos.

Variabilidade espacial de atributos químicos do solo sob diferentes usos e manejos

O uso e manejo do solo e da cultura são importantes condicionadores da variabilidade de atributos do solo. O trabalho foi desenvolvido em Selvíria (MS), com o objetivo de avaliar a variabilidade espacial do pH, potássio (K), cálcio (Ca), magnésio (Mg) e saturação por bases (V) em Latossolo Vermelho sob diferentes usos e manejos. Os solos foram amostrados em malha, com intervalos regulares de 2 m, perfazendo o total de 64 pontos, nas profundidades de 0,0-0,1 e 0,1-0,2 m, nas seguintes áreas: vegetação natural (Cerrado), plantio direto, plantio convencional e pastagem. As maiores variabilidades, medidas por meio do coeficiente de variação, foram observadas para K, Mg e Ca; o pH apresentou o menor coeficiente de variação nos diferentes usos e manejo do solo, e o atributo V, coeficiente de variação médio. Os sistemas preparo convencional e pastagem apresentaram os menores alcances quando comparado aos sistemas Cerrado e plantio direto.

Simulation of within field variability of corn yield with Ceres-Maize model

The CERES-Maize model was used to estimate the spatial variability in corn (Zea mays L.) yield for 1995 and 1996 using data measured on soil profiles located on a 30.5 m grid within a 3.9 ha field in Michigan. The model was calibrated for one grid profile for the 1995 and then used to simulate corn yield for all grid points for the 2 yrs. For the calibration for 1995, the model predicted corn yield within 2%. For 1995, the model predicted yield variability very well (r(2) = 0.85), producing similar yield maps with differences generally within +/- 300 kg ha(-1). For 1996, the model predicted low grain yields (1167 kg ha(-1)) compared with measured (8928 kg ha(-1)) because the model does not account for horizontal water movement within the landscape or water contributions from a water table. Under nonlimiting water conditions, the model performed well (average of 8717 vs. 8948 kg ha(-1)) but under-estimated the measured yield variability.

Fractal dimension and anisotropy of soil CO2 emission in an agricultural field during fallow

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelagem de distorções entre realizações de referenciais geodésicos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Comparação de métodos de mapeamento da distribuição espacial da infestação de plantas daninhas

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Agricultura de precisão em argissolo com variação nas formas de relevo sob o cultivo de cana-de-açúcar

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica do processo de secagem de madeiras

Pós-graduação em Engenharia Mecânica - FEG

High-resolution electrophoretic simulations: performance characteristics of one-dimensional simulators

Three comprehensive one-dimensional simulators were used on the same PC to simulate the dynamics of different electrophoretic configurations, including two migrating hybrid boundaries, an isotachophoretic boundary and the zone electrophoretic separation of ten monovalent anions. Two simulators, SIMUL5 and GENTRANS, use a uniform grid, while SPRESSO uses a dynamic adaptive grid. The simulators differ in the way components are handled. SIMUL5 and SPRESSO feature one equation for all components, whereas GENTRANS is based on the use of separate modules for the different types of monovalent components, a module for multivalent components and a module for proteins. The code for multivalent components is executed more slowly compared to those for monovalent components. Furthermore, with SIMUL5, the computational time interval becomes smaller when it is operated with a reduced calculation space that features moving borders, whereas GENTRANS offers the possibility of using data smoothing (removal of negative concentrations), which can avoid numerical oscillations and speed up a simulation. SPRESSO with its adaptive grid could be employed to simulate the same configurations with smaller numbers of grid points and thus is faster in certain but not all cases. The data reveal that simulations featuring a large number of monovalent components distributed such that a high mesh is required throughout a large proportion of the column are fastest executed with GENTRANS.

Excitation and charge transfer in H-H+ collisions at 5-80 keV and application to astrophysical shocks

In astrophysical regimes where the collisional excitation of hydrogen atoms is relevant, the cross-sections for the interactions of hydrogen atoms with electrons and protons are necessary for calculating line profiles and intensities. In particular, at relative velocities exceeding ∼1000 km s−1, collisional excitation by protons dominates over that by electrons. Surprisingly, the H–H+ cross-sections at these velocities do not exist for atomic levels of n≥ 4, forcing researchers to utilize extrapolation via inaccurate scaling laws. In this study, we present a faster and improved algorithm for computing cross-sections for the H–H+ collisional system, including excitation and charge transfer to the n≥ 2 levels of the hydrogen atom. We develop a code named BDSCX which directly solves the Schrödinger equation with variable (but non-adaptive) resolution and utilizes a hybrid spatial-Fourier grid. Our novel hybrid grid reduces the number of grid points needed from ∼4000n6 (for a ‘brute force’, Cartesian grid) to ∼2000n4 and speeds up the computation by a factor of ∼50 for calculations going up to n= 4. We present (l, m)-resolved results for charge transfer and excitation final states for n= 2–4 and for projectile energies of 5–80 keV, as well as fitting functions for the cross-sections. The ability to accurately compute H–H+ cross-sections to n= 4 allows us to calculate the Balmer decrement, the ratio of Hα to Hβ line intensities. We find that the Balmer decrement starts to increase beyond its largely constant value of 2–3 below 10 keV, reaching values of 4–5 at 5 keV, thus complicating its use as a diagnostic of dust extinction when fast (∼1000 km s−1) shocks are impinging upon the ambient interstellar medium.

Antarctic temperature changes during the last millennium: evaluation of simulations and reconstructions

Temperature changes in Antarctica over the last millennium are investigated using proxy records, a set of simulations driven by natural and anthropogenic forcings and one simulation with data assimilation. Over Antarctica, a long term cooling trend in annual mean is simulated during the period 1000–1850. The main contributor to this cooling trend is the volcanic forcing, astronomical forcing playing a dominant role at seasonal timescale. Since 1850, all the models produce an Antarctic warming in response to the increase in greenhouse gas concentrations. We present a composite of Antarctic temperature, calculated by averaging seven temperature records derived from isotope measurements in ice cores. This simple approach is supported by the coherency displayed between model results at these data grid points and Antarctic mean temperature. The composite shows a weak multi-centennial cooling trend during the pre-industrial period and a warming after 1850 that is broadly consistent with model results. In both data and simulations, large regional variations are superimposed on this common signal, at decadal to centennial timescales. The model results appear spatially more consistent than ice core records. We conclude that more records are needed to resolve the complex spatial distribution of Antarctic temperature variations during the last millennium.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Optimization of multidimensional cross-section tables for few-group core calculations

Multigroup diffusion codes for three dimensional LWR core analysis use as input data pre-generated homogenized few group cross sections and discontinuity factors for certain combinations of state variables, such as temperatures or densities. The simplest way of compiling those data are tabulated libraries, where a grid covering the domain of state variables is defined and the homogenized cross sections are computed at the grid points. Then, during the core calculation, an interpolation algorithm is used to compute the cross sections from the table values. Since interpolation errors depend on the distance between the grid points, a determined refinement of the mesh is required to reach a target accuracy, which could lead to large data storage volume and a large number of lattice transport calculations. In this paper, a simple and effective procedure to optimize the distribution of grid points for tabulated libraries is presented. Optimality is considered in the sense of building a non-uniform point distribution with the minimum number of grid points for each state variable satisfying a given target accuracy in k-effective. The procedure consists of determining the sensitivity coefficients of k-effective to cross sections using perturbation theory; and estimating the interpolation errors committed with different mesh steps for each state variable. These results allow evaluating the influence of interpolation errors of each cross section on k-effective for any combination of state variables, and estimating the optimal distance between grid points.