An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted

Influence of mass diffusion in sea ice dynamical models

The mixing of floes of different thickness caused by repeated deformation of the ice cover is modeled as diffusion, and the mass balance equation for sea ice accounting for mass diffusion is developed. The effect of deformational diffusion on the ice thickness balance is shown to reach 1% of the divergence effect, which describes ridging and lead formation. This means that with the same accuracy the mass balance equation can be written in terms of mean velocity rather than mean mass-weighted velocity, which one should correctly use for a multicomponent fluid such as sea ice with components identified by floe thickness. Mixing (diffusion) of sea ice also occurs because of turbulent variations in wind and ocean drags that are unresolved in models. Estimates of the importance of turbulent mass diffusion on the dynamic redistribution of ice thickness are determined using empirical data for the turbulent diffusivity. For long-time-scale prediction (≫5 days), where unresolved atmospheric motion may have a length scale on the order of the Arctic basin and the time scale is larger than the synoptic time scale of atmospheric events, turbulent mass diffusion can exceed 10% of the divergence effect. However, for short-time-scale prediction, for example, 5 days, the unresolved scales are on the order of 100 km, and turbulent diffusion is about 0.1% of the divergence effect. Because inertial effects are small in the dynamics of the sea ice pack, diffusive momentum transfer can be disregarded.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Evaluating adjusted forcing and model spread for historical and future scenarios in the CMIP5 generation of climate models

We utilize energy budget diagnostics from the Coupled Model Intercomparison Project phase 5 (CMIP5) to evaluate the models' climate forcing since preindustrial times employing an established regression technique. The climate forcing evaluated this way, termed the adjusted forcing (AF), includes a rapid adjustment term associated with cloud changes and other tropospheric and land-surface changes. We estimate a 2010 total anthropogenic and natural AF from CMIP5 models of 1.9 ± 0.9 W m−2 (5–95% range). The projected AF of the Representative Concentration Pathway simulations are lower than their expected radiative forcing (RF) in 2095 but agree well with efficacy weighted forcings from integrated assessment models. The smaller AF, compared to RF, is likely due to cloud adjustment. Multimodel time series of temperature change and AF from 1850 to 2100 have large intermodel spreads throughout the period. The intermodel spread of temperature change is principally driven by forcing differences in the present day and climate feedback differences in 2095, although forcing differences are still important for model spread at 2095. We find no significant relationship between the equilibrium climate sensitivity (ECS) of a model and its 2003 AF, in contrast to that found in older models where higher ECS models generally had less forcing. Given the large present-day model spread, there is no indication of any tendency by modelling groups to adjust their aerosol forcing in order to produce observed trends. Instead, some CMIP5 models have a relatively large positive forcing and overestimate the observed temperature change.

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

First-order turbulence closure for modelling complex canopy flows

Simple first-order closure remains an attractive way of formulating equations for complex canopy flows when the aim is to find analytic or simple numerical solutions to illustrate fundamental physical processes. Nevertheless, the limitations of such closures must be understood if the resulting models are to illuminate rather than mislead. We propose five conditions that first-order closures must satisfy then test two widely used closures against them. The first is the eddy diffusivity based on a mixing length. We discuss the origins of this approach, its use in simple canopy flows and extensions to more complex flows. We find that it satisfies most of the conditions and, because the reasons for its failures are well understood, it is a reliable methodology. The second is the velocity-squared closure that relates shear stress to the square of mean velocity. Again we discuss the origins of this closure and show that it is based on incorrect physical principles and fails to satisfy any of the five conditions in complex canopy flows; consequently its use can lead to actively misleading conclusions.

Upper-mantle seismic anisotropy from SKS splitting in the South American stable platform: A test of asthenospheric flow models beneath the lithosphere

Upper-mantle seismic anisotropy has been extensively used to infer both present and past deformation processes at lithospheric and asthenospheric depths. Analysis of shear-wave splitting (mainly from core-refracted SKS phases) provides information regarding upper-mantle anisotropy. We present average measurements of fast-polarization directions at 21 new sites in poorly sampled regions of intra-plate South America, such as northern and northeastern Brazil. Despite sparse data coverage for the South American stable platform, consistent orientations are observed over hundreds of kilometers. Over most of the continent, the fast-polarization direction tends to be close to the absolute plate motion direction given by the hotspot reference model HS3-NUVEL-1A. A previous global comparison of the SKS fast-polarization directions with flow models of the upper mantle showed relatively poor correlation on the continents, which was interpreted as evidence for a large contribution of ""frozen"" anisotropy in the lithosphere. For the South American plate, our data indicate that one of the reasons for the poor correlation may have been the relatively coarse model of lithospheric thicknesses. We suggest that improved models of upper-mantle flow that are based on more detailed lithospheric thicknesses in South America may help to explain most of the observed anisotropy patterns.

Deformation and geochronology of syntectonic granitoids emplaced in the Major Gercino Shear Zone, southeastern South America

The Major Gercino Shear Zone is one of the NE-SW lineaments that separate the Neoproterozoic Dom Feliciano Belt, of Brazil and Uruguay, into two different domains: a northwestern supracrustal domain from a southeastern granitoid domain. The shear zone, striking NE, is composed of protomylonites to ultramylonites with mainly dextral kinematic indicators. In Santa Catarina State, southern Brazil, the shear zone is composed of two mylonite belts. The mylonites have mineral orientations produced under greenschist fades conditions at a high strain rate. Strong flattening and coaxial deformation indicate the transpressive character, while the role of pure shear is emphasized by the orientation of the mylonite belts in relation to the inferred stress field component. The quartz microstructures point out that different dynamic recrystallization regimes and crystal plasticity were the dominant mechanisms of deformation during the mylonitization process. Additionally, the fabrics suggest that the glide systems are activated for deformation conditions compatible with the metamorphism in the middle greenschist facies. Elongated granitoid intrusions belonging to two petrographically, geochemically and isotopically distinct rock associations occur between the two mylonite belts. The structures observed in the granites result from a deformation range from magmatic to solid-state conditions points to a continuum of magma straining during and just after its crystallization. Conventional U-Pb analysis of multi-crystal zircon fractions yielded essentially identical ages of 609 +/- 16 Ma and 614 +/- 2 Ma for the two granitic associations, and constrain the transpressive phase of the shear zone. K-Ar ages of biotites between 585 and 560 Ma record the slow cooling and uplift of the intrusions. Some K-Ar ages of micas in regional mylonites are similar, suggesting that thermo-tectonic activity was intense up to this time, probably related to the agglutination of the granite belt to the supracrustal belt NW of the MGSZ. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Influence diagnostics for linear models with first-order autoregressive elliptical errors

We introduce in this paper the class of linear models with first-order autoregressive elliptical errors. The score functions and the Fisher information matrices are derived for the parameters of interest and an iterative process is proposed for the parameter estimation. Some robustness aspects of the maximum likelihood estimates are discussed. The normal curvatures of local influence are also derived for some usual perturbation schemes whereas diagnostic graphics to assess the sensitivity of the maximum likelihood estimates are proposed. The methodology is applied to analyse the daily log excess return on the Microsoft whose empirical distributions appear to have AR(1) and heavy-tailed errors. (C) 2008 Elsevier B.V. All rights reserved.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

Short-term temporal changes of bare soil CO2 fluxes after tillage described by first-order decay models

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

Breeding value accuracy estimates for growth traits using random regression and multi-trait models in Nelore cattle

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

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows

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

Reavaliação de critérios estruturais na hidrogeologia de terrenos cristalinos, com ênfase na neotectônica e sensoriamento remoto

Hydrogeological prospecting in Northeast Brazil and in other crystalline terrains has been developed on the basis of structural and regional geology concepts that date back to the 50-60 decades and, as such, demand a natural re-evaluation and update. In this kind of terrain, the percolation and accumulation of ground water are controlled by fractures and other types of discontinuities, such as foliations and geological contacts that, through weathering, impart porosity and permeability to the rocks, allowing water flow and storage. Several factors should be considered in the process of locating water wells, as discussed in the literature. Among these, the kind of structures, fracture geometry (including aperture and connectivity) and their geological and chronological context. It is important to correlate fracture systems with the regional neotectonic framework. Fractures at low angle (sub parallel) with the principal stress axis (s1) are those which tend to open (actually they work as tension joints) and, in principle, would present major hydric potential; in the opposite side, fractures at high angle to s1 would behave as closed by a compressional component. Fractures diagonal to the compression and tension axes correspond to shear fractures and, due to their connectivity with second fractures, are also important in terms of hydric potential. Uplift followed by terrain denudation leads to decompression and a general tendency to open (aided by weathering processes) fractures and other rock discontinuities, at different orientations. Low angle fractures, formed in this context, are equally important to increase connectivity, collection of water and recharge of the aquifer systems. In a general way, an opening component (neotectonic or by terrain decompression) and several models to increase fracture connectivity correlate with a greater hydric potential of these structures. Together with parallel research, this thesis addresses models of ground water occurrence in crystalline terrains, either improving well established concepts like the (Riacho-Fenda model), but also stressing other possibilities, like the role of alluvium and paleo-regoliths (the Calha Elúvio-Aluvionar model) and of strongly altered, permo-porous zones placed at variable depths below the present surface, flanking several types of discontinuities, especially interconnected fracture arrays (the Bolsões de Intemperismo model). Different methodological approaches are also discussed in order to improve success rates in the location of water wells in crystalline terrains. In this methodological review, a number of case studies were selected in the eastern domain of the State of Rio Grande do Norte, involving the localities of Santa Cruz, Santo Antônio, Serrinha, Nova Cruz, Montanhas, Lagoa de Pedras and Lagoa Salgada. Besides the neotectonic analysis of brittle structures, this Thesis addresses the validation of remote sensing as a tool for ground water prospecting. Several techniques were tested in order to detect and select areas with higher potential for ground water accumulation, using Landsat 5-TM and RADARSAT images, besides conventional aerial photos. A number of filters were tested to emphasize lineaments in the images, improving their discrimination, to identify areas with higher overburden humidity, which could reflect subsurface water accumulation, as well as alluvium and other sedimentary covers that might act as recharge zones. The work started with a regional analysis with the orbital images, followed by analysis of aerial photos, up to a detailed structural study of rock exposures in the terrain. This last step involved the analysis of outcrops surrounding wells (in a ray of approximately 10 to 100 m) with distinct productivities, including dry examples. At the level required for detail, it was not possible to accomplish a statistical approach using the available well data catalogs, which lack the desired specific information. The methodology worked out in this Thesis must undergo a testing phase through location of new water wells. An increase in the success rates as desired will led to a further consolidation step with wider divulgation of the methodology to private companies and governmental agencies involved in ground water prospecting in crystalline terrains