A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains

An analytical model of orographic gravity wave drag due to sheared flow past elliptical mountains is developed. The model extends the domain of applicability of the well-known Phillips model to wind profiles that vary relatively slowly in the vertical, so that they may be treated using a WKB approximation. The model illustrates how linear processes associated with wind profile shear and curvature affect the drag force exerted by the airflow on mountains, and how it is crucial to extend the WKB approximation to second order in the small perturbation parameter for these effects to be taken into account. For the simplest wind profiles, the normalized drag depends only on the Richardson number, Ri, of the flow at the surface and on the aspect ratio, γ, of the mountain. For a linear wind profile, the drag decreases as Ri decreases, and this variation is faster when the wind is across the mountain than when it is along the mountain. For a wind that rotates with height maintaining its magnitude, the drag generally increases as Ri decreases, by an amount depending on γ and on the incidence angle. The results from WKB theory are compared with exact linear results and also with results from a non-hydrostatic nonlinear numerical model, showing in general encouraging agreement, down to values of Ri of order one.

Resonant gravity-wave drag enhancement in linear stratified flow over mountains

High-drag states produced in stratified flow over a 2D ridge and an axisymmetric mountain are investigated using a linear, hydrostatic, analytical model. A wind profile is assumed where the background velocity is constant up to a height z1 and then decreases linearly, and the internal gravity-wave solutions are calculated exactly. In flow over a 2D ridge, the normalized surface drag is given by a closed-form analytical expression, while in flow over an axisymmetric mountain it is given by an expression involving a simple 1D integral. The drag is found to depend on two dimensionless parameters: a dimensionless height formed with z_1, and the Richardson number, Ri, in the shear layer. The drag oscillates as z_1 increases, with a period of half the hydrostatic vertical wavelength of the gravity waves. The amplitude of this modulation increases as Ri decreases. This behaviour is due to wave reflection at z_1. Drag maxima correspond to constructive interference of the upward- and downward-propagating waves in the region z < z_1, while drag minima correspond to destructive interference. The reflection coefficient at the interface z = z_1 increases as Ri decreases. The critical level, z_c, plays no role in the drag amplification. A preliminary numerical treatment of nonlinear effects is presented, where z_c appears to become more relevant, and flow over a 2D ridge qualitatively changes its character. But these effects, and their connection with linear theory, still need to be better understood.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part I: gravity wave forcing

Parameterization schemes for the drag due to atmospheric gravity waves are discussed and compared in the context of a simple one-dimensional model of the quasi-biennial oscillation (QBO). A number of fundamental issues are examined in detail, with the goal of providing a better understanding of the mechanism by which gravity wave drag can produce an equatorial zonal wind oscillation. The gravity wave–driven QBOs are compared with those obtained from a parameterization of equatorial planetary waves. In all gravity wave cases, it is seen that the inclusion of vertical diffusion is crucial for the descent of the shear zones and the development of the QBO. An important difference between the schemes for the two types of waves is that in the case of equatorial planetary waves, vertical diffusion is needed only at the lowest levels, while for the gravity wave drag schemes it must be included at all levels. The question of whether there is downward propagation of influence in the simulated QBOs is addressed. In the gravity wave drag schemes, the evolution of the wind at a given level depends on the wind above, as well as on the wind below. This is in contrast to the parameterization for the equatorial planetary waves in which there is downward propagation of phase only. The stability of a zero-wind initial state is examined, and it is determined that a small perturbation to such a state will amplify with time to the extent that a zonal wind oscillation is permitted.

Nonlinear saturation of baroclinic instability: part II: continuously stratified fluid

Rigorous upper bounds are derived that limit the finite-amplitude growth of arbitrary nonzonal disturbances to an unstable baroclinic zonal flow in a continuously stratified, quasi-geostrophic, semi-infinite fluid. Bounds are obtained bath on the depth-integrated eddy potential enstrophy and on the eddy available potential energy (APE) at the ground. The method used to derive the bounds is essentially analogous to that used in Part I of this study for the two-layer model: it relies on the existence of a nonlinear Liapunov (normed) stability theorem, which is a finite-amplitude generalization of the Charney-Stern theorem. As in Part I, the bounds are valid both for conservative (unforced, inviscid) flow, as well as for forced-dissipative flow when the dissipation is proportional to the potential vorticity in the interior, and to the potential temperature at the ground. The character of the results depends on the dimensionless external parameter γ = f02ξ/β0N2H, where ξ is the maximum vertical shear of the zonal wind, H is the density scale height, and the other symbols have their usual meaning. When γ ≫ 1, corresponding to “deep” unstable modes (vertical scale ≈H), the bound on the eddy potential enstrophy is just the total potential enstrophy in the system; but when γ≪1, corresponding to ‘shallow’ unstable modes (vertical scale ≈γH), the eddy potential enstrophy can be bounded well below the total amount available in the system. In neither case can the bound on the eddy APE prevent a complete neutralization of the surface temperature gradient which is in accord with numerical experience. For the special case of the Charney model of baroclinic instability, and in the limit of infinitesimal initial eddy disturbance amplitude, the bound states that the dimensionless eddy potential enstrophy cannot exceed (γ + 1)2/24&gamma2h when γ ≥ 1, or 1/6;&gammah when γ ≤ 1; here h = HN/f0L is the dimensionless scale height and L is the width of the channel. These bounds are very similar to (though of course generally larger than) ad hoc estimates based on baroclinic-adjustment arguments. The possibility of using these kinds of bounds for eddy-amplitude closure in a transient-eddy parameterization scheme is also discussed.

Mean motions induced by baroclinic instability in a jet

A study is made of the zonal-mean motions induced by a growing baroclinic wave in several contexts, under the framework of three different analysis schemes: the conventional Eulerian mean (EM), the transformed Eulerian mean (TEM), and the generalized Lagrangian mean (GLM). The effect of meridional shear in the initial jet on these induced mean motions is considered by treating the instability problem in the context of the two-layer model. The conceptual simplicity of the TEM formulation is shown to be useful in diagnosing the dynamics of instability, much as it has been found helpful in many problems of wave, mean-flow interaction. In addition, it is found that the TEM vertical velocity is a very good indicator of the GLM vertical velocity. However, the GLM meridional velocity is always convergent towards the centre of instability activity, and is not at all well represented by the nondivergent TEM meridional velocity. In comparing the results with Uryu's (1979) calculation of the GLM circulation induced by a growing Eady wave, it is found that the inclusion of meridional jet shear in the present work leads to some strikingly different effects in the GLM zonal wind acceleration. In the case of pure baroclinic instability treated by Uryu, the Eulerian and Stokes accelerations nearly cancel each other in the centre of the channel, leaving a weak Lagrangian acceleration opposed to the Eulerian one. In the more general case of mixed baroclinic-barotropic instability, however, the Eulerian and Stokes accelerations can reinforce one another, leading to a very strong Lagrangian zonal wind

Lower-stratospheric radiative damping and polar-night jet oscillation events

The effect of stratospheric radiative damping time scales on stratospheric variability and on stratosphere–troposphere coupling is investigated in a simplified global circulation model by modifying the vertical profile of radiative damping in the stratosphere while holding it fixed in the troposphere. Perpetual-January conditions are imposed, with sinusoidal topography of zonal wavenumber 1 or 2. The depth and duration of the simulated sudden stratospheric warmings closely track the lower-stratospheric radiative time scales. Simulations with the most realistic profiles of radiative damping exhibit extended time-scale recoveries analogous to polar-night jet oscillation (PJO) events, which are observed to follow sufficiently deep stratospheric warmings. These events are characterized by weak lower-stratospheric winds and enhanced stability near the tropopause, which persist for up to 3 months following the initial warming. They are obtained with both wave-1 and wave-2 topography. Planetary-scale Eliassen–Palm (EP) fluxes entering the vortex are also suppressed, which is in agreement with observed PJO events. Consistent with previous studies, the tropospheric jets shift equatorward in response to the warmings. The duration of the shift is closely correlated with the period of enhanced stability. The magnitude of the shift in these runs, however, is sensitive only to the zonal wavenumber of the topography. Although the shift is sustained primarily by synoptic-scale eddies, the net effect of the topographic form drag and the planetary-scale fluxes is not negligible; they damp the surface wind response but enhance the vertical shear. The tropospheric response may also reduce the generation of planetary waves, further extending the stratospheric dynamical time scales.

Modelling the reorientation of sea-ice faults as the wind changes direction

A discrete-element model of sea ice is used to study how a 90° change in wind direction alters the pattern of faults generated through mechanical failure of the ice. The sea-ice domain is 400km in size and consists of polygonal floes obtained through a Voronoi tessellation. Initially the floes are frozen together through viscous–elastic joints that can break under sufficient compressive, tensile and shear deformation. A constant wind-stress gradient is applied until the initially frozen ice pack is broken into roughly diamond-shaped aggregates, with crack angles determined by wing-crack formation. Then partial refreezing of the cracks delineating the aggregates is modelled through reduction of their length by a particular fraction, the ice pack deformation is neglected and the wind stress is rotated by 90°. New cracks form, delineating aggregates with a different orientation. Our results show the new crack orientation depends on the refrozen fraction of the initial faults: as this fraction increases, the new cracks gradually rotate to the new wind direction, reaching 90° for fully refrozen faults. Such reorientation is determined by a competition between new cracks forming at a preferential angle determined by the wing-crack theory and at old cracks oriented at a less favourable angle but having higher stresses due to shorter contacts across the joints

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.

Vertical ion flow in the polar ionosphere

Outflowing ions from the polar ionosphere fall into two categories: the classical polar wind and the suprathermal ion flows. The flows in both these categories vary a great deal with altitude. The classical polar wind is supersonic at high altitude: at ∼3 RE geocentric, the observed polar wind is H+ dominated and has a Mach number of 2.5–5.1. At 400–600 km, thermal and suprathermal upward O+ ion fluxes frequently occur at the poleward edge of the nightside auroral oval during magnetically active times. Above 500 km, ions are accelerated transverse to the local geomagnetic field. At 1400 km, transversely accelerated ions are frequently observed in winter nights but rarely appear in the summer. In the dayside cleft above ∼2000 km, ions of all species are transversely heated and upwell with significant number and heat fluxes, forming a cleft ion fountain as they convect across the polar cap. Upwelling ions are observed most (least) frequently in the summer (winter). At yet higher altitudes, energetic (>10 eV to several kiloelectron volts) upflowing H+ and O+ ions are frequently observed, their active time occurrence frequency being as high as 0.7 at auroral latitudes and 0.3 in the polar cap. Their composition, intensity, and angular characteristics vary quantitatively with solar activity, being O+ dominant and more intense near solar maximum. Their resulting ion outflow is dominated by ions below 1 keV and reaches 3.5×10^26 O+ and 7×10^25 H+ ions s^{−1} at magnetically active times (Kp≥5) near solar maximum. In comparison, the estimated polar wind ion outflow at times of moderate solar activity is 7×10^25H+ and 4×10^24 He+ ions s^{−1}. The estimated <10-eV cleft ion fountain flow is 3.8×10^25 O+ and 8.6×10^23 H+ ions s^{−1} near solar maximum.

Correlation between vertical misfits and stresses transmitted to implants from metal frameworks

An inappropriate prosthetic fit could cause stress over the interface implant/bone. The objective of this study was to compare stresses transmitted to implants from frameworks cast using different materials and to investigate a possible correlation between vertical misfits and these stresses. Fifteen one-piece cast frameworks simulating bars for fixed prosthesis in a model with five implants were fabricated and arranged into three different groups according to the material used for casting: CP Ti (commercially pure titanium), Co-Cr (cobalt-chromium) or Ni-Cr-Ti (nickel-chromium-titanium) alloys. Each framework was installed over the metal model with all screws tightened to a 10 N cm torque and then, vertical misfits were measured using an optical microscope. The stresses transmitted to implants were measured using quantitative photoelastic analysis in values of maximum shear stress (T), when each framework was tightened to the photoelastic model to a 10 N cm standardized torque. Stress data were statistically analyzed using one-way ANOVA and Tukey`s test and correlation tests were performed using Pearson`s rank correlation (alpha = 0.05). Mean and standard deviation values of vertical misfit are presented for CP Ti (22.40 +/- 9.05 mu m), Co-Cr (66.41 +/- 35.47 mu m) and Ni-Cr-Ti (32.20 +/- 24.47 mu m). Stresses generated by Co-Cr alloy (tau = 7.70 +/- 2.16 kPa) were significantly higher than those generated by CP Ti (tau = 5.86 +/- 1.55 kPa, p = 0.018) and Ni-Cr-Ti alloy (tau =5.74 +/- 3.05 kPa, p = 0.011), which were similar (p = 0.982). Correlations between vertical misfits and stresses around the implants were not significant as for any evaluated materials. (C) 2011 Elsevier Ltd. All rights reserved.

Cambrian mafic to felsic magmatism and its connections with transcurrent shear zones of the Borborema Province (NE Brazil): Implications for the late assembly of the West Gondwana

New U-Pb (SHRIMP) and (40)Ar/(39)Ar isotopic data of igneous rocks and mylonites of the Borborema Province (NE Brazil) show that a wide range of tectonothermal events affected the province during the transition from the Precambrian to the Cambrian. Concordant zircon U-Pb data constrained the crystallization age of mafic stocks, mafic to felsic dikes and granite batholiths between 548 and 533 Ma. These bodies were emplaced in a regional strain field combining extension and dextral shearing. The ductile shear deformation overprinted an older basement fabric to develop a low- to medium metamorphic grade vertical mylonite belt that cut the province in the E-W direction. Magnetic fabrics of the Cambrian batholiths determined by anisotropy of magnetic susceptibility are consistent with syntectonic emplacement. The magmatic pulses and shear deformation would have supplied enough heat to reset the synkinematic micas of mylonites to yield (40)Ar/(39)Ar plateau cooling ages between ca. 550 and 510 Ma. These results provide evidence that emplacement of Early Cambrian mafic and felsic magmas were accompanied by regional-scale shear deformations, probably in the consequence of late collisions along the West Gondwana margin. (C) 2010 Published by Elsevier B.V.

Vertical displacement during late-collisional escape tectonics (Brasiliano Orogeny) in the Ribeira Belt, São Paulo State, Brazil

During the Brasiliano-Pan-African Orogeny, West Gondwana formed by collisional processes around the Sao Francisco-Congo Craton. The Ribeira belt, in southeastern Brazil, resulted from northwestward collision (650-600 Ma), followed by large-scale northeast-southwest dextral strike-slip shear movements related to late-collisional escape tectonics (ca 600 Ma).In São Paulo State, three groups, also interpreted as terranes, are recognised in the Ribeira Belt, the Embu, Itapira and Sao Rogue Groups. The Embu and Itapira Groups are formed of sillimanite-gneisses, schists and migmatites intruded by Neoproterozoic calc-alkaline granitoids, all thrusted northwestward. The Sao Rogue Group is composed of metasediments and metavolcanics in greenschist-facies. Its deformation indicates a transpressional regime associated with tectonic escape. Sub-alkaline granites were emplaced in shallow levels during this regime. Microstructural studies along the Itu, Moreiras and Taxaquara Shear Zones demonstrate the coexistence of horizontal and Vertical displacement components during the transpressional regime. The vertical component is regarded as responsible for the lateral juxtaposition of different crustal levels. (C) 1999 Elsevier B.V. Limited. All rights reserved.

Avaliação do desempenho de um aeromotor de eixo vertical

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

The influence of large-scale circulation, transient and local processes on the seasonal circulation of the Eastern Brazilian Shelf, 13 degrees S

The circulation at the Eastern Brazilian Shelf (EBS), near 13 degrees S, is discussed in terms of the currents and hydrography, associating large-scale circulation, transient and local processes to establish a regional picture of the EBS circulation. The results show that the circulation within the continental shelf and slope region is strongly affected by the seasonal changes in the wind field and mesa/large-scale circulation. Transient processes associated to the passage of Cold Front systems or meso-scale activity and the presence of a local canyon add more complexity to the system. During the austral spring and summer seasons, the prevailing upwelling favorable winds blowing from E-NE were responsible for driving southwestward shelf currents. The interaction with the Western Boundary Current (the Brazil Current), especially during summer, was significant and a considerable vertical shear in the velocity field was observed at the outer shelf. The passage of a Cold Front system during the springtime caused a complete reversal of the mean flow and contributed to the deepening of the Mixed Layer Depth (MLD). In addition, the presence of Salvador Canyon, subject to an upwelling favorable boundary current, enhanced the upwelling system, when compared to the upwelling observed at the adjacent shelf. During the austral autumn and winter seasons the prevailing downwelling favorable winds blowing from the SE acted to total reverse the shelf circulation, resulting in a northeastward flow. The passage of a strong Cold Front, during the autumn season, contributed not only to the strengthening of the flow but also to the deepening of the MLD. The presence of the Salvador Canyon, when subject to a downwelling favorable boundary current, caused an intensification of the downwelling process. Interestingly, the alongshore velocity at the shelf region adjacent to the head of the canyon was less affected when compared to the upwelling situation.