Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Ocean heat content variability and change in an ensemble of ocean reanalyses

Accurate knowledge of the location and magnitude of ocean heat content (OHC) variability and change is essential for understanding the processes that govern decadal variations in surface temperature, quantifying changes in the planetary energy budget, and developing constraints on the transient climate response to external forcings. We present an overview of the temporal and spatial characteristics of OHC variability and change as represented by an ensemble of dynamical and statistical ocean reanalyses (ORAs). Spatial maps of the 0–300 m layer show large regions of the Pacific and Indian Oceans where the interannual variability of the ensemble mean exceeds ensemble spread, indicating that OHC variations are well-constrained by the available observations over the period 1993–2009. At deeper levels, the ORAs are less well-constrained by observations with the largest differences across the ensemble mostly associated with areas of high eddy kinetic energy, such as the Southern Ocean and boundary current regions. Spatial patterns of OHC change for the period 1997–2009 show good agreement in the upper 300 m and are characterized by a strong dipole pattern in the Pacific Ocean. There is less agreement in the patterns of change at deeper levels, potentially linked to differences in the representation of ocean dynamics, such as water mass formation processes. However, the Atlantic and Southern Oceans are regions in which many ORAs show widespread warming below 700 m over the period 1997–2009. Annual time series of global and hemispheric OHC change for 0–700 m show the largest spread for the data sparse Southern Hemisphere and a number of ORAs seem to be subject to large initialization ‘shock’ over the first few years. In agreement with previous studies, a number of ORAs exhibit enhanced ocean heat uptake below 300 and 700 m during the mid-1990s or early 2000s. The ORA ensemble mean (±1 standard deviation) of rolling 5-year trends in full-depth OHC shows a relatively steady heat uptake of approximately 0.9 ± 0.8 W m−2 (expressed relative to Earth’s surface area) between 1995 and 2002, which reduces to about 0.2 ± 0.6 W m−2 between 2004 and 2006, in qualitative agreement with recent analysis of Earth’s energy imbalance. There is a marked reduction in the ensemble spread of OHC trends below 300 m as the Argo profiling float observations become available in the early 2000s. In general, we suggest that ORAs should be treated with caution when employed to understand past ocean warming trends—especially when considering the deeper ocean where there is little in the way of observational constraints. The current work emphasizes the need to better observe the deep ocean, both for providing observational constraints for future ocean state estimation efforts and also to develop improved models and data assimilation methods.

An ensemble of eddy-permitting global ocean reanalyses from the MyOcean project

A set of four eddy-permitting global ocean reanalyses produced in the framework of the MyOcean project have been compared over the altimetry period 1993–2011. The main differences among the reanalyses used here come from the data assimilation scheme implemented to control the ocean state by inserting reprocessed observations of sea surface temperature (SST), in situ temperature and salinity profiles, sea level anomaly and sea-ice concentration. A first objective of this work includes assessing the interannual variability and trends for a series of parameters, usually considered in the community as essential ocean variables: SST, sea surface salinity, temperature and salinity averaged over meaningful layers of the water column, sea level, transports across pre-defined sections, and sea ice parameters. The eddy-permitting nature of the global reanalyses allows also to estimate eddy kinetic energy. The results show that in general there is a good consistency between the different reanalyses. An intercomparison against experiments without data assimilation was done during the MyOcean project and we conclude that data assimilation is crucial for correctly simulating some quantities such as regional trends of sea level as well as the eddy kinetic energy. A second objective is to show that the ensemble mean of reanalyses can be evaluated as one single system regarding its reliability in reproducing the climate signals, where both variability and uncertainties are assessed through the ensemble spread and signal-to-noise ratio. The main advantage of having access to several reanalyses differing in the way data assimilation is performed is that it becomes possible to assess part of the total uncertainty. Given the fact that we use very similar ocean models and atmospheric forcing, we can conclude that the spread of the ensemble of reanalyses is mainly representative of our ability to gauge uncertainty in the assimilation methods. This uncertainty changes a lot from one ocean parameter to another, especially in global indices. However, despite several caveats in the design of the multi-system ensemble, the main conclusion from this study is that an eddy-permitting multi-system ensemble approach has become mature and our results provide a first step towards a systematic comparison of eddy-permitting global ocean reanalyses aimed at providing robust conclusions on the recent evolution of the oceanic state.

Southern Ocean isopycnal mixing and ventilation changes driven by winds

Observed and predicted changes in the strength of the westerly winds blowing over the Southern Ocean have motivated a number of studies of the response of the Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation (MOC) to wind perturbations and led to the discovery of the``eddy-compensation" regime, wherein the MOC becomes insensitive to wind changes. In addition to the MOC, tracer transport also depends on mixing processes. Here we show, in a high-resolution process model, that isopycnal mixing by mesoscale eddies is strongly dependent on the wind strength. This dependence can be explained by mixing-length theory and is driven by increases in eddy kinetic energy; the mixing length does not change strongly in our simulation. Simulation of a passive ventilation tracer (analogous to CFCs or anthropogenic CO$_2$) demonstrates that variations in tracer uptake across experiments are dominated by changes in isopycnal mixing, rather than changes in the MOC. We argue that, to properly understand tracer uptake under different wind-forcing scenarios, the sensitivity of isopycnal mixing to winds must be accounted for.

STUDY OF THE EQUATORIAL ATLANTIC OCEAN MIXING LAYER USING A ONE-DIMENSIONAL TURBULENCE MODEL

The General Ocean Turbulence Model (GOTM) is applied to the diagnostic turbulence field of the mixing layer (ML) over the equatorial region of the Atlantic Ocean. Two situations were investigated: rainy and dry seasons, defined, respectively, by the presence of the intertropical convergence zone and by its northward displacement. Simulations were carried out using data from a PIRATA buoy located on the equator at 23 degrees W to compute surface turbulent fluxes and from the NASA/GEWEX Surface Radiation Budget Project to close the surface radiation balance. A data assimilation scheme was used as a surrogate for the physical effects not present in the one-dimensional model. In the rainy season, results show that the ML is shallower due to the weaker surface stress and stronger stable stratification; the maximum ML depth reached during this season is around 15 m, with an averaged diurnal variation of 7 m depth. In the dry season, the stronger surface stress and the enhanced surface heat balance components enable higher mechanical production of turbulent kinetic energy and, at night, the buoyancy acts also enhancing turbulence in the first meters of depth, characterizing a deeper ML, reaching around 60 m and presenting an average diurnal variation of 30 m.

An analysis of the environmental energetics associated with the transition of the first South Atlantic hurricane

This study presents the first analysis of the energetics associated with a hybrid cyclone`s transition in the Southern Hemisphere, Hurricane Catarina ( March 2004). Catarina has earned a place in history as the first documented South Atlantic hurricane, but its unusual tropical transition is still poorly understood. Here we show that Catarina`s transition was preceded by marked environmental changes in the Lorenz energy cycle, with an abrupt shift from a baroclinic to a predominantly barotropic state. Such changes help to explain the unusual vortex`s growth until its transition was completed. Although the vortex`s energy flux is not explicitly calculated, a likely mechanism linking the environmental energetics with Catarina is the extraction of eddy kinetic energy from horizontal momentum and heat transfers within the through component of the blocking. The results advance the understanding of this rare event and suggest that the technique has a great potential to study transitioning systems in general.

Weak nonmesonic decay spectra of hypernuclei

We compute one- and two-nucleon kinetic-energy spectra and opening-angle distributions for the nonmesonic weak decay of several hypernuclei, and compare our results with some recent data. The decay dynamics is described by transition potentials of the one-meson-exchange type, and the nuclear structure aspects by two versions of the independent-particle shell model (IPSM). In version IPSM-a, the bole states are treated as stationary, while in version IPSM-b the deep-hole ones are considered to be quasi-stationary and are described by Breit-Wigner distributions.

On the dynamic behaviour of the "click" mechanism in dipteran flight

In this paper, the dynamic behaviour of the "click" mechanism is analysed. A more accurate model is used than in the past, in which the limits of movement due to the geometry of the flight mechanism are imposed. Moreover, the effects of different damping models are investigated. In previous work, the damping model was assumed to be of the linear viscous type for simplicity, but it is likely that the damping due to drag forces is nonlinear. Accordingly, a model of damping in which the damping force is proportional to the square of the velocity is used, and the results are compared with the simpler model of linear viscous damping. Because of the complexity of the model an analytical approach is not possible so the problem has been cast in terms of non-dimensional variables and solved numerically. The peak kinetic energy of the wing root per energy input in one cycle is chosen to study the effectiveness of the "click" mechanism compared with a linear resonant mechanism. It is shown that, the "click" mechanism has distinct advantages when it is driven below its resonant frequency. When the damping is quadratic, there are some further advantages compared to when the damping is linear and viscous, provided that the amplitude of the excitation force is large enough to avoid the erratic behaviour of the mechanism that occurs for small forces. (C) 2011 Elsevier Ltd. All rights reserved.

NUMERICAL SIMULATION of CONVECTION-DIFFUSION PROBLEMS BY THE CONTROL-VOLUME-BASED FINITE-ELEMENT METHOD

This work presents a numerical study of the tri-dimensional convection-diffusion equation by the control-volume-based on finite-element method using quadratic hexahedral elements. Considering that the equation governing this problem in its main variable may represent several properties, including temperature, turbulent kinetic energy, viscous dissipation rate of the turbulent kinetic energy, specific dissipation rate of the turbulent kinetic energy, or even the concentration of a contaminant in a given medium, among others, the wide applicability of this problem is thus evidenced. Three cases of temperature distributions will be studied specifically in this work, in addition to one case of pollutant dispersion upon analysis of the concentration of a contaminant in a fixed flow point. Some comparisons will be carried out against works found in the open literature, while others will be done according to each phenomenon characteristics.

LEAST ACTION PRINCIPLE and THE INCOMPRESSIBLE EULER EQUATIONS WITH VARIABLE DENSITY

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

The Fermi pseudo-potential in one dimension

Wu and Yu recently examined point interactions in one dimension in the form of the Fermi pseudo-potential. on the other hand there are point interactions in the form of self-adjoint extensions (SAEs) of the kinetic energy operator. We examine the relationship between the point interactions in these two forms in the one-channel and two-channel cases. In the one-channel case the pseudo-potential leads to the standard three-parameter family of SAEs. In the two-channel case the pseudo-potential furnishes a ten-parameter family of SAEs.

PT-invariant point interactions in one dimension

By using Wu and Yu's pseudo-potential, we construct point interactions in one dimension that are complex but conform to space-time reflection (PT) invariance. The resulting point interactions are equivalent to those obtained by Albeverio, Fei and Kurasov as self-adjoint extensions of the kinetic energy operator.

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

Weak nonmesonic decay spectra of hypernuclei

We compute one- and two-nucleon kinetic-energy spectra and opening-angle distributions for the nonmesonic weak decay of several hypernuclei, and compare our results with some recent data. The decay dynamics is described by transition potentials of the one-meson-exchange type, and the nuclear structure aspects by two versions of the independent-particle shell model (IPSM). In version IPSM-a, the bole states are treated as stationary, while in version IPSM-b the deep-hole ones are considered to be quasi-stationary and are described by Breit-Wigner distributions.

Dynamical properties of a particle in a wave packet: Scaling invariance and boundary crisis

Some dynamical properties present in a problem concerning the acceleration of particles in a wave packet are studied. The dynamics of the model is described in terms of a two-dimensional area preserving map. We show that the phase space is mixed in the sense that there are regular and chaotic regions coexisting. We use a connection with the standard map in order to find the position of the first invariant spanning curve which borders the chaotic sea. We find that the position of the first invariant spanning curve increases as a power of the control parameter with the exponent 2/3. The standard deviation of the kinetic energy of an ensemble of initial conditions obeys a power law as a function of time, and saturates after some crossover. Scaling formalism is used in order to characterise the chaotic region close to the transition from integrability to nonintegrability and a relationship between the power law exponents is derived. The formalism can be applied in many different systems with mixed phase space. Then, dissipation is introduced into the model and therefore the property of area preservation is broken, and consequently attractors are observed. We show that after a small change of the dissipation, the chaotic attractor as well as its basin of attraction are destroyed, thus leading the system to experience a boundary crisis. The transient after the crisis follows a power law with exponent -2. (C) 2011 Elsevier Ltd. All rights reserved.