Showing papers on "Pressure gradient published in 2007"

A Physical Relationship between Electron-Proton Temperature Equilibration and Mach Number in Fast Collisionless Shocks

Abstract: The analysis of Balmer-dominated optical spectra from nonradiative (adiabatic) SNRs has shown that the ratio of the electron to proton temperature at the blast wave is close to unity at v(s) approximate or smaller than 400 km s(-1) but declines sharply down to the minimum value of m(e) /m(p) dictated by the jump conditions at shock speeds exceeding 2000 km s(-1). We propose a physical model for the heating of electrons and ions in non-cosmic-ray-dominated, strong shocks (v(s) 400 km s(-1)) wherein the electrons are heated by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of electron heating over a wide range of shock speeds, producing a relationship (T(e)/T(p))sub 0 proportional to v(-2/s) (proportional to M(-2)) that is fully consistent with the observations.

Approach to an asymptotic state for zero pressure gradient turbulent boundary layers

Abstract: Flat plate turbulent boundary layers under zero pressure gradient at high Reynolds numbers are studied to reveal appropriate scale relations and asymptotic behaviour. Careful examination of the skin-friction coefficient results confirms the necessity for direct and independent measurement of wall shear stress. We find that many of the previously proposed empirical relations accurately describe the local Cf behaviour when modified and underpinned by the same experimental data. The variation of the integral parameter, H, shows consistent agreement between the experimental data and the relation from classical theory. In accordance with the classical theory, the ratio of D and d asymptotes to a constant. Then, the usefulness of the ratio of appropriately defined mean and turbulent time-scales to define and diagnose equilibrium flow is established. Next, the description of mean velocity profiles is revisited, and the validity of the logarithmic law is re-established using both the mean velocity profile and its diagnostic function. The wake parameter, P, is shown to reach an asymptotic value at the highest available experimental Reynolds numbers if correct values of logarithmic-law constants and an appropriate skin-friction estimate are used. The paper closes with a discussion of the Reynolds number trends of the outer velocity defect which are important to establish a consistent similarity theory and appropriate scaling.

Intrinsic rotation and electric field shear

Abstract: A novel mechanism for the generation and amplification of intrinsic rotation at the low-mode to high-mode transition is presented. The mechanism is one where the net parallel flow is accelerated by turbulence. A preferential direction of acceleration results from the breaking of k‖→−k‖ symmetry by sheared E×B flow. It is shown that the equilibrium pressure gradient contributes a piece of the parallel Reynolds stress, which is nonzero for vanishing parallel flow, and so can accelerate the plasma, driving net intrinsic rotation. Rotation drive, transport, and fluctuation dynamics are treated self-consistently.

Self-consistent high-Reynolds-number asymptotics for zero-pressure-gradient turbulent boundary layers

Abstract: The asymptotic behavior of mean velocity and integral parameters in flat plate turbulent boundary layers under zero pressure gradient are studied for Reynolds numbers approaching infinity. Using the classical two-layer approach of Millikan, Rotta, and Clauser with a logarithmic velocity profile in the overlap region between “inner” and “outer” layers, a fully self-consistent leading-order description of the mean velocity profile and all integral parameters is developed. It is shown that this description fits most high Reynolds number data, and in particular their Reynolds number dependence, exceedingly well; i.e., within experimental errors.

Predicting incipient motion, including the effect of turbulent pressure fluctuations in the bed

Abstract: [1] Incipient motion of particles on a river bed can, in principle, be understood and predicted from a balance of the forces acting on the particles. On a natural river bed the exposure of particles to the flow is variable. The effect of exposure on the initiation of motion, including the case of deep embedding, is studied in this paper. So far, the turbulence-induced lifting force has been derived exclusively from examinations of the surface flow. The understanding of destabilization processes of the riverbed is improved when turbulence-induced vertical pressure gradients in the upper sediment layer are included. This is particularly important for particles that lie in the lee of larger ones. An analytical solution for the critical Shields parameter of spherical particles is found, including the effects of pressure fluctuations in the bed, macroscale flow structures, average bed slope, and shallow flow. Existing laboratory data designed to test the effect of exposure and roughness conditions are in excellent agreement with the new model.

Numerical modeling of self-channeling granular flows and of their levee-channel deposits

Abstract: When not laterally confined in valleys, pyroclastic flows create their own channel along the slope by selecting a given flowing width. Furthermore, the lobe-shaped deposits display a very specific morphology with high parallel lateral levees. A numerical model based on Saint Venant equations and the empirical variable friction coefficient proposed by Pouliquen and Forterre (2002) is used to simulate unconfined granular flow over an inclined plane with a constant supply. Numerical simulations successfully reproduce the self-channeling of the granular lobe and the levee-channel morphology in the deposits without having to take into account mixture concepts or polydispersity. Numerical simulations suggest that the quasi-static shoulders bordering the flow are created behind the front of the granular material by the rotation of the velocity field due to the balance between gravity, the two-dimensional pressure gradient, and friction. For a simplified hydrostatic model, competition between the decreasing friction coefficient and increasing surface gradient as the thickness decreases seems to play a key role in the dynamics of unconfined flows. The description of the other disregarded components of the stress tensor would be expected to change the balance of forces. The front's shape appears to be constant during propagation. The width of the flowing channel and the velocity of the material within it are almost steady and uniform. Numerical results suggest that measurement of the width and thickness of the central channel morphology in deposits in the field provides an estimate of the velocity and thickness during emplacement. Citation: Mangeney, A., F. Bouchut, N. Thomas, J. P. Vilotte, and M. O. Bristeau (2007), Numerical modeling of self-channeling granular flows and of their levee-channel deposits,

A wave propagation model of blood flow in large vessels using an approximate velocity profile function

Abstract: Lumped-parameter models (zero-dimensional) and wave-propagation models (one-dimensional) for pressure and flow in large vessels, as well as fully three-dimensional fluid–structure interaction models for pressure and velocity, can contribute valuably to answering physiological and patho-physiological questions that arise in the diagnostics and treatment of cardiovascular diseases. Lumped-parameter models are of importance mainly for the modelling of the complete cardiovascular system but provide little detail on local pressure and flow wave phenomena. Fully three-dimensional fluid–structure interaction models consume a large amount of computer time and must be provided with suitable boundary conditions that are often not known. One-dimensional wave-propagation models in the frequency and time domain are well suited to obtaining clinically relevant information on local pressure and flow waves travelling through the arterial system. They can also be used to provide boundary conditions for fully three-dimensional models, provided that they are defined in, or transferred to, the time domain. Most of the one-dimensional wave propagation models in the time domain described in the literature assume velocity profiles and therefore frictional forces to be in phase with the flow, whereas from exact solutions in the frequency domain a phase difference between the flow and the wall shear stress is known to exist. In this study an approximate velocity profile function more suitable for one-dimensional wave propagation is introduced and evaluated. It will be shown that this profile function provides first-order approximations for the wall shear stress and the nonlinear term in the momentum equation, as a function of local flow and pressure gradient in the time domain. The convective term as well as the approximate friction term are compared to their counterparts obtained from Womersley profiles and show good agreement in the complete range of the Womersley parameter α. In the limiting cases, for Womersley parameters α → 0 and α → ∞, they completely coincide. It is shown that in one-dimensional wave propagation, the friction term based on the newly introduced approximate profile function is important when considering pressure and flow wave propagation in intermediate-sized vessels.

Buoyancy-driven variable-density turbulence

Abstract: Buoyancy-generated motions in an unstably stratified medium composed of two incompressible miscible fluids with different densities, as occurs in the variable-density Rayleigh–Taylor instability, are examined using direct numerical simulations. The non-equilibrium homogeneous buoyantly driven problem is proposed as a unit problem for variable density turbulence to study: (i) the nature of variable density turbulence, (ii) the transition to turbulence and the generation of turbulence by the conversion of potential to kinetic energy; (iii) the role of non-Boussinesq effects; and (iv) a parameterization of the initial conditions by a static Reynolds number. Simulations are performed for Atwood numbers up to 0.5 with root mean square density up to 50% of the mean density and Schmidt numbers, 0.1 ≤ Sc ≤ 2. The benchmark problem has been designed to have the largest mass flux possible and is, in this configuration, the maximally unstable non-equilibrium flow possible. It is found that the mass flux, owing to its central role in the conversion of potential to kinetic energy, is probably the single most important dynamical quantity to predict in lower-dimensional models. Other primary findings include the evolution of the mean pressure gradient: during the non-Boussinesq portions of the flow, the evolution of the mean pressure gradient is non-hydrostatic (as opposed to a Boussinesq fluid) and is set by the evolution of the specific volume pressure gradient correlation. To obtain the numerical solution, a new pressure projection algorithm which treats the pressure step exactly, useful for simulations of non-solenoidal velocity flows, has been constructed.

Unsteady aerodynamic forces estimation on a square cylinder by TR-PIV

Abstract: The unsteady aerodynamic forces acting on a square cross-sectional cylinder are investigated by means of time-resolved particle image velocimetry (TR-PIV) at Reynolds number 4,900. The objective of the investigation is to prove the feasibility of non-intrusive force measurements around two-dimensional bodies. The PIV measurements performed at a rate of 1 kHz enable a time resolved (TR) description of the vortex shedding phenomenon occurring at 10 Hz and to follow the time evolution of vortex dominated wake. The instantaneous aerodynamic force coefficients are obtained from the integration of the force equations within a control volume enclosing the object. The required instantaneous pressure distribution is inferred making use of two physical models: Bernoulli relation is adopted in the potential slowly-evolving flow region; in the turbulent wake, the Navier–Stokes equations are invoked to determine the pressure gradient spatial distribution, which integrated in space yields the pressure distribution. The spatial acceleration field is directly obtained from the temporal difference of the time-filtered velocity field. For a choice of the control volume approximately one model height away from the surface the contributions to the aerodynamic forces coming from the different terms of the force equation are individually examined. The convective term dominates the unsteady lift forces whereas the pressure term prevails for the drag. The temporal evolution of C L returns a clear periodic pattern in phase with the vortex shedding at a frequency of 10.1 Hz (Strouhal number St = 0.128) with oscillation amplitude of 0.9, whereas C D barely shows periodicity. The measurement uncertainties associated to the evaluation of all the terms in the force equation and especially in relation to TR-PIV measurements are discussed.

Midlatitude Wind Stress–Sea Surface Temperature Coupling in the Vicinity of Oceanic Fronts

Abstract: The influences of strong gradients in sea surface temperature on near-surface cross-front winds are explored in a series of idealized numerical modeling experiments. The atmospheric model is the Naval Research Laboratory Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) model, which is fully coupled to the Regional Ocean Modeling System (ROMS) ocean model. A series of idealized, two-dimensional model calculations is carried out in which the wind blows from the warm-to-cold side or the cold-to-warm side of an initially prescribed ocean front. The evolution of the near-surface winds, boundary layer, and thermal structure is described, and the balances in the momentum equation are diagnosed. The changes in surface winds across the front are consistent with previous models and observations, showing a strong positive correlation with the sea surface temperature and boundary layer thickness. The coupling arises mainly as a result of changes in the flux Richardson number across the front, and the strength of the coupling coefficient grows quadratically with the strength of the cross-front geostrophic wind. The acceleration of the winds over warm water results primarily from the rapid change in turbulent mixing and the resulting unbalanced Coriolis force in the vicinity of the front. Much of the loss/gain of momentum perpendicular to the front in the upper and lower boundary layer results from acceleration/ deceleration of the flow parallel to the front via the Coriolis term. This mechanism is different from the previously suggested processes of downward mixing of momentum and adjustment to the horizontal pressure gradient, and is active for flows off the equator with sufficiently strong winds. Although the main focus of this work is on the midlatitude, strong wind regime, calculations at low latitudes and with weak winds show that the pressure gradient and turbulent mixing terms dominate the cross-front momentum budget, consistent with previous work.

Intrinsic pump-conduit behavior of lymphangions

Abstract: Lymphangions, segments of lymphatic vessels bounded by valves, have characteristics of both ventricles and arteries. They can act primarily like pumps when actively transporting lymph against a pressure gradient. They also can act as conduit vessels when passively transporting lymph down a pressure gradient. This duality has implications for clinical treatment of several types of edema, since the strategy to optimize lymph flow may depend on whether it is most beneficial for lymphangions to act as pumps or conduits. To address this duality, we employed a simple computational model of a contracting lymphangion, predicted the flows at both positive and negative axial pressure gradients, and validated the results with in vitro experiments on bovine mesenteric vessels. This model illustrates that contraction increases flow for normal axial pressure gradients. With edema, limb elevation, or external compression, however, the pressure gradient might reverse, and lymph may flow passively down a pressure gradient. In such cases, the valves may be forced open during the entire contraction cycle. The vessel thus acts as a conduit, and contraction has the effect of increasing resistance to passive flow, thus inhibiting flow rather than promoting it. This analysis may explain a possible physiological benefit of the observed flow-mediated inhibition of the lymphatic pump at high flow rates.

Heat transfer, pressure drop, and flow pattern recognition during condensation inside smooth, helical micro-fin, and herringbone tubes

Abstract: This paper presents a study of flow regimes, pressure drops, and heat transfer coefficients during refrigerant condensation inside a smooth, an 18° helical micro-fin, and a herringbone tubes. Experimental work was conducted for condensing refrigerants R-22, R-407C, and R-134a at an average saturation temperature of 40 °C with mass fluxes ranging from 400 to 800 kg m −2 s −1 , and with vapour qualities ranging from 0.85 to 0.95 at condenser inlet and from 0.05 to 0.15 at condenser outlet. These test conditions represent annular and intermittent (slug and plug) flow conditions. Results showed that transition from annular flow to intermittent flow, on average for the three refrigerants, occurred at a vapour quality of 0.49 for the smooth tube, 0.29 for the helical micro-fin tube, and 0.26 for the herringbone tube. These transition vapour qualities were also reflected in the pressure gradients, with the herringbone tube having the highest pressure gradient. The pressure gradients encountered in the herringbone tube were about 79% higher than that of the smooth tube and about 27% higher than that of the helical micro-fin tube. A widely used pressure drop correlation for condensation in helical micro-fin tubes was modified for the case of the herringbone tube. The modified correlation predicted the data within a 1% error with an absolute deviation of 7%. Heat transfer enhancement factors for the herringbone tube against the smooth tube were on average 70% higher while against the helical micro-fin tube it was 40% higher. A correlation for predicting heat transfer coefficients inside a helical micro-fin tube was modified for the herringbone tube. On average the correlation predicted the data to within 4% with an average standard deviation of 8%.

A One-dimensional Model of Blood Flow in Arteries with Friction and Convection Based on the Womersley Velocity Profile

Abstract: In this paper, we present a one-dimensional model for blood flow in arteries, without assuming an a priori shape for the velocity profile across an artery (Azer, Ph.D. thesis, Courant Institute, New York University, 2006). We combine the one-dimensional equations for conservation of mass and momentum with the Womersley model for the velocity profile in an iterative way. The pressure gradient of the one-dimensional model drives the Womersley equations, and the velocity profiles calculated then feed back into both the friction and nonlinear parts of the one-dimensional model. Besides enabling us to evaluate the friction correctly and also to use the velocity profile to correct the nonlinear terms, having the velocity profile available as output should be useful in a variety of applications. We present flow simulations using both structured trees and pure resistance models for the small arteries, and compare the resulting flow and pressure waves under various friction models. Moreover, we show how to couple the one-dimensional equations with the Taylor diffusion limit (Azer, Int J Heat Mass Transfer 2005;48:2735-40; Taylor, Proc R Soc Lond Ser A 1953;219:186-203) of the convection-diffusion equations to drive the concentration of a solute along an artery in time.

The forces exerted by aquatic suction feeders on their prey

Abstract: Successful prey capture by aquatic suction feeders depends on the ability of the predator to generate a flow of water external to the mouth that overcomes any movements and forces that the prey uses to resist the suction flow. Elucidating the nature and magnitude of these forces is a key to understanding what limits suction feeding performance. We identify three potential forces produced by the suction flow field: drag, acceleration reaction and the fluid pressure gradient. Using a mathematical model parametrized with empirical data from feeding bluegill, Lepomis macrochirus, we explore the relative magnitude of these forces under three encounter scenarios with a 5mm diameter, spherical prey: an immobile mid-water prey; a similar prey that executes an escape response; and a prey item that grips a substratum. Contrary to the almost exclusive emphasis on drag in the suction feeding literature, it made a minor contribution to the total forces in all three cases. In all three scenarios, the pressure gradient is the largest of the three forces. These results are important because previous researchers have emphasized drag and have not explicitly recognized a role for the pressure gradient force in suction feeding. The simulations suggest previously unrecognized mechanisms that suction feeders can use to enhance the forces that they exert, by increasing the steepness of the pressure gradient that the prey item is exposed to. This can be accomplished either by increasing the rate of increase in fluid velocity or by restricting the size of the mouth aperture, which creates a steeper spatial gradient in pressure.

Peristaltic transport of a third order fluid under the effect of a magnetic field

Abstract: The peristaltic transport of a third order fluid in a planar channel is considered. The fluid is electrically conducting by a transverse magnetic field. The perturbation solution is obtained using small Deborah number. Expressions of stream function, longitudinal velocity and pressure gradient valid for long wavelength are developed. Numerical integration is performed to analyze the effect of Hartman number on the pressure rise and frictional force. It is noted that both Hartman and Deborah numbers suppress the flow.

The Local Geometry of Gas Injection into Saturated Homogeneous Porous Media

Abstract: The injection of gases into liquid saturated porous media is of theoretical and practical interest (e.g., air sparging for the removal of volatile organic compounds from contaminated aquifer sediments). The influence of the rate of gas delivery and the vertical distance from the source are developed. The concept of a “near-injection region” is presented in which the pressure gradients exceed buoyant gradients and thus exhibits largely radial flow. The near-injection size is shown to have an area required to carry the injected gas flow under unit gradient. The parabolic movement of gas outside of this area which has often been observed is explained as reflecting the sum of many realizations of gas channels following random lateral movements as they precede upward independent of flux. These concepts are confirmed through comparison with published and experimental data of air injection into slabs consisting of saturated sands of a range of textures.

Buoyancy-driven crack propagation: the limit of large fracture toughness

Abstract: We study steady vertical propagation of a crack filled with buoyant viscous fluid through an elastic solid with large effective fracture toughness. For a crack fed by a constant flux Q, a non-dimensional fracture toughness K=Kc /(3μQm 3/2)1/4 describes the relative magnitudes of resistance to fracture and resistance to viscous flow, where Kc is the dimensional fracture toughness, μ the fluid viscosity and m the elastic modulus. Even in the limit K ≫ 1, the rate of propagation is determined by viscous effects. In this limit the large fracture toughness requires the fluid behind the crack tip to form a large teardrop-shaped head of length O(K 2/3) and width O(K 4/3), which is fed by a much narrower tail. In the head, buoyancy is balanced by a hydrostatic pressure gradient with the viscous pressure gradient negligible except at the tip; in the tail, buoyancy is balanced by viscosity with elasticity also playing a role in a region within O(K 2/3) of the head. A narrow matching region of length O(K −2/5) and width O(K −4/15), termed the neck, connects the head and the tail. Scalings and asymptotic solutions for the three regions are derived and compared with full numerical solutions for K ≤ 3600 by analysing the integro-differential equation that couples lubrication flow in the crack to the elastic pressure gradient. Time-dependent numerical solutions for buoyancy-driven propagation of a constant-volume crack show a quasi-steady head and neck structure with a propagation rate that decreases like t −2/3 due to the dynamics of viscous flow in the draining tail.

Three-Dimensional Tidal Flow in an Elongated, Rotating Basin

Abstract: The three-dimensional tidal circulation in an elongated basin of arbitrary depth is described with a linear, constant-density model on the f plane. Rotation fundamentally alters the lateral flow, introducing a lateral recirculation comparable in magnitude to the axial flow, as long as friction is not too large. This circulation is due to the imbalance between the cross-channel sea level gradient, which is in near-geostrophic balance with the Coriolis acceleration associated with the vertically averaged axial flow, and the Coriolis acceleration associated with the vertically sheared axial flow. During flood condition, for example, the lateral Coriolis acceleration near the surface exceeds the pressure gradient, tending to accelerate the lateral flow, while the converse is true near the bottom. As a result, with rotation, fluid parcels tend to corkscrew into and out of the basin in a tidal period. The axial flow is only weakly modified by rotation. When friction is small, the axial velocity is uniform in each section, except in a narrow bottom boundary layer where it decreases to zero. The boundary layer thickness increases with friction, so that with moderate or large friction, axial velocities are sheared from bottom to surface. When friction is large, the local and Coriolis accelerations are both small and the dynamics are governed by a balance between friction and the pressure gradient.

Atmospheric Circulation Effects on Wind Speed Variability at Turbine Height

Abstract: Mean monthly wind speed at 70 m above ground level is investigated for 11 sites in Minnesota for the period 1995–2003. Wind speeds at these sites show significant spatial and temporal coherence, with prolonged periods of above- and below-normal values that can persist for as long as 12 months. Monthly variation in wind speed primarily is determined by the north–south pressure gradient, which captures between 22% and 47% of the variability (depending on the site). Regression on wind speed residuals (pressure gradient effects removed) shows that an additional 6%–15% of the variation can be related to the Arctic Oscillation (AO) and Nino-3.4 sea surface temperature (SST) anomalies. Wind speeds showed little correspondence with variation in the Pacific–North American (PNA) circulation index. The effect of the strong El Nino of 1997/98 on the wind speed time series was investigated by recomputing the regression equations with this period excluded. The north–south pressure gradient remains the primary ...

Diffusioosmotic flows in slit nanochannels.

Abstract: Diffusioosmotic flows of electrolyte solutions in slit nanochannels with homogeneous surface charges induced by electrolyte concentration gradients in the absence of externally applied pressure gradients and potential differences are investigated theoretically. A continuum mathematical model consisting of the strongly coupled Nernst-Planck equations for the ionic species' concentrations, the Poisson equation for the electric potential in the electrolyte solution, and the Navier-Stokes equations for the flow field is numerically solved simultaneously. The induced diffusioosmotic flow through the nanochannel is computed as functions of the externally imposed concentration gradient, the concentration of the electrolyte solution, and the surface charge density along the walls of the nanochannel. With the externally applied electrolyte concentration gradient, a strongly spatially dependent electric field and pressure gradient are induced within the nanochannel that, in turn, generate a spatially dependent diffusioosmotic flow. The diffusioosmotic flow is opposite to the applied concentration gradient for a relatively low bulk electrolyte concentration. However, the electrolyte solution flows from one end of the nanochannel with a higher electrolyte concentration to the other end with a lower electrolyte concentration when the bulk electrolyte concentration is relatively high. There is an optimal concentration gradient under which the flow rate attains the maximum. The induced flow is enhanced with the increase in the fixed surface charge along the wall of the nanochannel for a relatively low bulk electrolyte concentration.

Noninvasive measurement of intra-aneurysmal pressure and flow pattern using phase contrast with vastly undersampled isotropic projection imaging.

Abstract: BACKGROUND AND PURPOSE: Currently, more reliable parameters to predict the risk of aneurysmal rupture are needed. Intra-aneurysmal pressure gradients and flow maps could provide additional information regarding the risk of rupture. Our hypothesis was that phase contrast with vastly undersampled isotropic projection reconstruction (PC-VIPR), a novel 3D MR imaging sequence, could accurately assess intra-aneurysmal pressure gradients in a canine aneurysmal model when compared with invasive measurements. MATERIALS AND METHODS: A total of 13 surgically created aneurysms in 8 canines were included in this study. Pressure measurements were performed in the parent vessel, aneurysm neck, and 5 regions within the aneurysmal sac with a microcatheter. PC-VIPR sequence was used to obtain cardiac-gated velocity measurements in a region covering the entire aneurysm. The velocity and pressure gradient maps derived from the PC-VIPR data were then coregistered with the anatomic DSA images and compared with catheter measurements. RESULTS: In 7 of the bifurcation aneurysms, the velocity flow maps demonstrated a recirculation flow pattern with a small neck-to-dome pressure gradient (mean, +0.5 mm Hg). In 1 bifurcation aneurysm, a flow jet extending from the neck to the dome with significantly greater pressure gradient (+50.2 mm Hg) was observed. All sidewall aneurysms had low flow in the sac with intermediate pressure gradients (mean, +8.3 mm Hg). High statistical correlation existed between PC-VIPR aneurysmal pressures and microcatheter pressure measurements ( R = 0.82, P CONCLUSION: PC-VIPR can provide anatomic as well as noninvasive quantitative and qualitative hemodynamic information in the canine aneurysm model. The PC-VIPR intra-aneurysmal pressure measurements correlated well with catheter measurements.

Remote sensing of surface pressure on Mars with the Mars Express/OMEGA spectrometer: 1. Retrieval method

Abstract: [1] Observing and analyzing the variations of pressure on the surface of a planet is essential to understand the dynamics of its atmosphere. On Mars the absorption by atmospheric CO2 of the solar light reflected on the surface allows us to measure the surface pressure by remote sensing. We use the imaging spectrometer OMEGA aboard Mars Express, which provides an excellent signal to noise ratio and the ability to produce maps of surface pressure with a resolution ranging from 400 m to a few kilometers. Surface pressure is measured by fitting spectra of the CO2 absorption band centered at 2 μm. To process the hundreds of thousands of pixels present in each OMEGA image, we have developed a fast and accurate algorithm based on a line-by-line radiative transfer model which includes scattering and absorption by dust aerosols. In each pixel the temperature profile, the dust opacity, and the surface spectrum are carefully determined from the OMEGA data set or from other sources to maximize the accuracy of the retrieval. We estimate the 1-σ relative error to be around 7 Pa in bright regions and about 10 Pa in darker regions, with a possible systematic bias on the absolute pressure lower than 30 Pa (4%). The method is first tested by comparing an OMEGA pressure retrieval obtained over the Viking Lander 1 (VL1) landing site with in situ measurements recorded 30 years ago by the VL1 barometer. The retrievals are further validated using a surface pressure predictor which combines the VL1 pressure records with the MOLA topography and meteorological pressure gradients simulated with a General Circulation Model. A good agreement is obtained. In particular, OMEGA is able to monitor the seasonal variations of the surface pressure in Isidis Planitia. Such a tool can be applied to detect meteorological phenomena, as described by Spiga et al. (2007).