Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsidence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

High-Resolution Simulation of Hurricane Bonnie (1998): Water Budget - Part 2

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, Scott A.↗

Interpretation of TOMS Observations of Tropical Tropospheric Ozone with a Global Model and In Situ Observations

We interpret the distribution of tropical tropospheric ozone columns (TTOCs) from the Total Ozone Mapping Spectrometer (TOMS) by using a global three-dimensional model of tropospheric chemistry (GEOS-CHEM) and additional information from in situ observations. The GEOS-CHEM TTOCs capture 44% of the variance of monthly mean TOMS TTOCs from the convective cloud differential method (CCD) with no global bias. Major discrepancies are found over northern Africa and south Asia where the TOMS TTOCs do not capture the seasonal enhancements from biomass burning found in the model and in aircraft observations. A characteristic feature of these northern topical enhancements, in contrast to southern tropical enhancements, is that they are driven by the lower troposphere where the sensitivity of TOMS is poor due to Rayleigh scattering. We develop an efficiency correction to the TOMS retrieval algorithm that accounts for the variability of ozone in the lower troposphere. This efficiency correction increases TTOC's over biomass burning regions by 3-5 Dobson units (DU) and decreases them by 2-5 DU over oceanic regions, improving the agreement between CCD TTOCs and in situ observations. Applying the correction to CCD TTOCs reduces by approximately DU the magnitude of the "tropical Atlantic paradox" [Thompson et al, 2000], i.e. the presence of a TTOC enhancement over the southern tropical Atlantic during the northern African biomass burning season in December-February. We reproduce the remainder of the paradox in the model and explain it by the combination of upper tropospheric ozone production from lightning NOx, peristent subsidence over the southern tropical Atlantic as part of the Walker circulation, and cross-equatorial transport of upper tropospheric ozone from northern midlatitudes in the African "westerly duct." These processes in the model can also account for the observed 13-17 DU persistent wave-1 pattern in TTOCs with a maximum above the tropical Atlantic and a minimum over the tropical Pacific during all seasons. The photochemical effects of mineral dust have only a minor role on the modeled distribution of TTOCs, including over northern Africa, due to multiple competing effects. The photochemical effects of mineral dust globally decease annual mean OH concentrations by 9%. A global lightning NOx source of 6 Tg N yr(sup -1) in the model produces a simulation that is most consistent with TOMS and in situ observations.

Martin, Randall V.↗

Outburst of Jupiter's synchrotron radiation after the impact of comet Shoemaker-Levy 9

Jupiter's nonthermal microwave emission, as measured by a global network of 11 radio telescopes, increased dramatically during the Shoemaker-Levy 9 impacts. The increase was wavelength-dependent, varying from approximately 10 percent at 70 to 90 centimeters to approximately 45 percent at 6 and 36 centimeters. The radio spectrum hardened (flattened toward shorter wavelengths) considerably during the week of impacts and continued to harden afterward. After the week of cometary impacts, the flux density began to subside at all wavelengths and was still declining 3 months later. Very Large Array and Australia Telescope images of the brightness distribution showed the enhancement to be localized in longitude and concentrated near the magnetic equator. The evidence therefore suggests that the increase in flux density was caused by a change in the resident particle population, for example, through an energization or spatial redistribution of the emitting particles.

Jupiter↗

Nonlinear simulations of Jupiter's 5-micron hot spots

Large-scale nonlinear simulations of Jupiter's 5-micron hot spots produce long-lived coherent structures that cause subsidence in local regions, explaining the low cloudiness and the dryness measured by the Galileo probe inside a hot spot. Like observed hot spots, the simulated coherent structures are equatorially confined, have periodic spacing, propagate west relative to the flow, are generally confined to one hemisphere, and have an anticyclonic gyre on their equatorward side. The southern edge of the simulated hot spots develops vertical shear of up to 70 meters per second in the eastward wind, which can explain the results of the Galileo probe Doppler wind experiment.

Jupiter↗

[Relationship Between Core Convective Structure and Intensity Change in Tropical Cyclones]

The downshear reformation of Tropical Storm Gabrielle (2001) was investigated using radar reflectivity and lightning data that were nearly continuous in time, as well as frequent aircraft reconnaissance fights. Initially the storm was a marginal tropical storm in an environment with strong 850-200 hPa vertical wind shear of 12-13 m/s and an approaching upper tropospheric trough. Both the observed outflow and an adiabatic balance model calculation showed that the radial-vertical circulation increased with time as the trough approached. Convection was highly asymmetric, with almost all radar return located in one quadrant left of downshear in the storm. Reconnaissance data show that an intense mesovortex formed downshear of the original center. This vortex was located just south of: rather than within, a strong downshear left lightning outbreak, consistent with tilting of the horizontal vorticity associated with the vertical wind shear. The downshear mesovortex contained a 972 hPa minimum central pressure, 20 hPa lower than minimum pressure in the original vortex just three hours earlier. The mesovortex became the new center of the storm, but weakened somewhat prior to landfall. It is argued that dry air carried around the storm from the region of upshear subsidence, as well as the direct effects of the shear, prevented the reformed vortex from continuing to intense.

Source record↗

Mesoscale Aspects of the Downshear Reformation of a Tropical Cyclone

The downshear reformation of Tropical Storm Gabrielle (2001) was investigated using radar reflectivity and lightning data that were nearly continuous in time, as well as frequent aircraft reconnaissance fights. Initially the storm was a marginal tropical storm in an environment with strong 850-200 hPa vertical wind shear of 12-13 m/s and an approaching upper tropospheric trough. Both the observed outflow and an adiabatic balance model calculation showed that the radial-vertical circulation increased with time as the trough approached. Convection was highly asymmetric, with almost all radar return located in one quadrant left of downshear in the storm. Reconnaissance data show that an intense mesovortex formed downshear of the original center. This vortex was located just south of: rather than within, a strong downshear left lightning outbreak, consistent with tilting of the horizontal vorticity associated with the vertical wind shear. The downshear mesovortex contained a 972 hPa minimum central pressure, 20 hPa lower than minimum pressure in the original vortex just three hours earlier. The mesovortex became the new center of the storm, but weakened somewhat prior to landfall. It is argued that dry air carried around the storm from the region of upshear subsidence, as well as the direct effects of the shear, prevented the reformed vortex from continuing to intense.

Source record↗

Tropical Convective Responses to Microphysical and Radiative Processes: A Sensitivity Study With a 2D Cloud Resolving Model

Prognostic cloud schemes are increasingly used in weather and climate models in order to better treat cloud-radiation processes. Simplifications are often made in such schemes for computational efficiency, like the scheme being used in the National Centers for Environment Prediction models that excludes some microphysical processes and precipitation-radiation interaction. In this study, sensitivity tests with a 2D cloud resolving model are carried out to examine effects of the excluded microphysical processes and precipitation-radiation interaction on tropical thermodynamics and cloud properties. The model is integrated for 10 days with the imposed vertical velocity derived from the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment. The experiment excluding the depositional growth of snow from cloud ice shows anomalous growth of cloud ice and more than 20% increase of fractional cloud cover, indicating that the lack of the depositional snow growth causes unrealistically large mixing ratio of cloud ice. The experiment excluding the precipitation-radiation interaction displays a significant cooling and drying bias. The analysis of heat and moisture budgets shows that the simulation without the interaction produces more stable upper troposphere and more unstable mid and lower troposphere than does the simulation with the interaction. Thus, the suppressed growth of ice clouds in upper troposphere and stronger radiative cooling in mid and lower troposphere are responsible for the cooling bias, and less evaporation of rain associated with the large-scale subsidence induces the drying in mid and lower troposphere.

Li, Xiao-Fan↗

Effects of Cloud-Microphysics on Tropical Atmospheric Hydrologic Processes in the GEOS GCM

The sensitivity of tropical atmospheric hydrologic processes to cloud-microphysics is investigated using the NASA GEOS GCM. Results show that a faster autoconversion - rate produces more warm rain and less clouds at all levels. Fewer clouds enhances longwave cooling and reduces shortwave heating in the upper troposphere, while more warm rain produces increased condensation heating in the lower troposphere. This vertical heating differential destablizes the tropical atmosphere, producing a positive feedback resulting in more rain over the tropics. The feedback is maintained via a two-cell secondary circulation. The lower cell is capped by horizontal divergence and maximum cloud detrainment near the melting/freezing, with rising motion in the warm rain region connected to descending motion in the cold rain region. The upper cell is found above the freezing/melting level, with longwave-induced subsidence in the warm rain and dry regions, coupled to forced ascent in the deep convection region. The tropical large scale circulation is found to be very sensitive to the radiative-dynamic effects induced by changes in autoconversion rate. Reduced cloud-radiation processes feedback due to a faster autoconversion rate results in intermittent but more energetic eastward propagating Madden and Julian Oscillations (MJO). Conversely,-a slower autconversion rate, with increased cloud radiation produces MJO's with more realistic westward propagating transients, resembling a supercloud cluster structure. Results suggests that warm rain and associated low and mid level clouds, i.e., cumulus congestus, may play a critical role in regulating the time-intervals of deep convections and hence the fundamental time scales of the MJO.

Lau, K. M.↗

Impaired sympathetic vascular regulation in humans after acute dynamic exercise

1. The reduction in vascular resistance which accompanies acute dynamic exercise does not subside immediately during recovery, resulting in a post-exercise hypotension. This sustained vasodilatation suggests that sympathetic vascular regulation is altered after exercise. 2. Therefore, we assessed the baroreflex control of sympathetic outflow in response to arterial pressure changes, and transduction of sympathetic activity into vascular resistance during a sympatho-excitatory stimulus (isometric handgrip exercise) after either exercise (60 min cycling at 60% peak aerobic power (VO2,peak)) or sham treatment (60 min seated rest) in nine healthy subjects. 3. Both muscle sympathetic nerve activity and calf vascular resistance were reduced after exercise (-29.7 +/- 8.8 and -25.3 +/- 9.1%, both P < 0.05). The baroreflex relation between diastolic pressure and sympathetic outflow was shifted downward after exercise (post-exercise intercept, 218 +/- 38 total integrated activity (heartbeat)-1; post-sham intercept, 318 +/- 51 total integrated activity (heartbeat)-1, P < 0.05), indicating less sympathetic outflow across all diastolic pressures. Further, the relation between sympathetic activity and vascular resistance was attenuated after exercise (post-exercise slope, 0.0031 +/- 0.0007 units (total integrated activity)-1 min; post-sham slope, 0.0100 +/- 0.0033 units (total integrated activity)-1 min, P < 0.05), indicating less vasoconstriction with any increase in sympathetic activity. 4. Thus, both baroreflex control of sympathetic outflow and the transduction of sympathetic activity into vascular resistance are altered after dynamic exercise. We conclude that the vasodilation which underlies post-exercise hypotension results from both neural and vascular phenomena.

Non-NASA Center↗

Deflections in Lava Flow Directions Relative to Topography in the Tharsis Region of Mars: Indications of Post-Flow Tectonic Motion

High-resolution topographic data from the Mars Orbiter Laser Altimeter (MOLA), and imagery from the Mars Orbiter Camera (MOC) and the Thermal Emission Imaging System (THEMIS) allow for the first accurate assessment of lava flow directions relative to topographic slopes in the Tharsis region. Tharisis has long been recognized as the dominant tectonic and volcanic province on the planet, with a complex geologic history. In this study, lava flow directions on Daedalia Planum, Syria Planum, Tempe Terra, and near the Tharsis Montes are compared with MOLA topographic contours to look for deviations of flow directions from the local slope direction. The topographic deviations identified in this study are likely due to Tharsis tectonic deformation that has modified the regional topography subsequent to the emplacement of the flows, and can be used to model the mechanisms and magnitudes of relatively recent tectonism in the region. A similar approach was used to identify possible post-flow tectonic subsidence on the Snake River Plain in Idaho.

Chadwick, D. J.↗

Back pain during 6 degrees head-down tilt approximates that during actual microgravity

Astronauts often experience back pain during spaceflight. Retrospectively, Wing et al. (14) found that during spaceflight, 14 of 19 Shuttle crewmembers experienced back pain, which they described as dull (62%), localized to the lower back (50%), and with an intensity of 2 on a 5-point scale. Further, the spine lengthens 4-7 cm in microgravity. Our objective was to compare back pain and spinal lengthening (body height increase) during simulated microgravity (6 degrees head-down tilt, HDT) with the same parameters during actual microgravity. Eight male subjects completed a modified McGill pain questionnaire with intensity graded from zero (no pain) to five (intense and incapacitating pain) each day at 7:00 pm during 2 d pre-HDT control, 16 d HDT, and 1 d post-HDT recovery periods. Also, the subjects' heights were measured each day while supine (control and recovery) and during HDT. Back pain increased from zero (pre-tilt control period) to 2.3 +/- 0.4 at days 1 to 3 of HDT, and was categorized as dull and/or burning pain in subjects' lower backs. Only 2 subjects reported any pain after day 9 of HDT and during recovery. Heights increased 2.1 +/- 0.5 cm by day 3 of HDT and remained at that level until the end of the HDT period. Although spinal lengthening in space is greater than that during HDT, the HDT model approximates the level, type, distribution, and time course of back pain associated with actual microgravity. In the HDT model, pain subsides in intensity when spinal lengthening stops.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

VA/Q distribution during heavy exercise and recovery in humans: implications for pulmonary edema

Ventilation-perfusion (VA/Q) inequality has been shown to increase with exercise. Potential mechanisms for this increase include nonuniform pulmonary vasoconstriction, ventilatory time constant inequality, reduced large airway gas mixing, and development of interstitial pulmonary edema. We hypothesized that persistence of VA/Q mismatch after ventilation and cardiac output subside during recovery would be consistent with edema; however, rapid resolution would suggest mechanisms related to changes in ventilation and blood flow per se. Thirteen healthy males performed near-maximal cycle ergometry at an inspiratory PO2 of 91 Torr (because hypoxia accentuates VA/Q mismatch on exercise). Cardiorespiratory variables and inert gas elimination patterns were measured at rest, during exercise, and between 2 and 30 min of recovery. Two profiles of VA/Q distribution behavior emerged during heavy exercise: in group 1 an increase in VA/Q mismatch (log SDQ of 0.35 +/- 0.02 at rest and 0.44 +/- 0.02 at exercise; P less than 0.05, n = 7) and in group 2 no change in VA/Q mismatch (n = 6). There were no differences in anthropometric data, work rate, O2 uptake, or ventilation during heavy exercise between groups. Group 1 demonstrated significantly greater VA/Q inequality, lower vital capacity, and higher forced expiratory flow at 25-75% of forced vital capacity for the first 20 min during recovery than group 2. Cardiac index was higher in group 1 both during heavy exercise and 4 and 6 min postexercise. However, both ventilation and cardiac output returned toward baseline values more rapidly than did VA/Q relationships. Arterial pH was lower in group 1 during exercise and recovery. We conclude that greater VA/Q inequality in group 1 and its persistence during recovery are consistent with the hypothesis that edema occurs and contributes to the increase in VA/Q inequality during exercise. This is supported by observation of greater blood flows and acidosis and, presumably therefore, higher pulmonary vascular pressures in such subjects.

NASA Discipline Cardiopulmonary↗

M.I.T./Canadian vestibular experiments on the Spacelab-1 mission: 4. Space motion sickness: symptoms, stimuli, and predictability

Space sickness symptoms were observed by 4 specially trained observers on Spacelab-1. Three reported persistent symptoms, and vomited repeatedly during the first and/or second day of flight. Head movements on all axes were provocative, particularly in pitch and roll. Head acceleration data recorded from 2 symptomatic crewmen showed that after several hours of physical activity in orbit, symptoms appeared, and thereafter both crewmen were compelled to limit head movements. Firm body contact with motionless surfaces helped alleviate symptoms. When crewmembers floated into unfamiliar body orientations in the cabin, inherent ambiguities in static visual orientation cues sometimes produced spatial reorientation episodes which were also provocative. Symptoms largely resembled those of other forms of prolonged motion sickness, superimposed upon other symptoms attributable to fluid shift. All 4 eventually used anti-motion sickness drugs. When they did, vomiting frequency was reduced. By the 4th day, symptoms subsided, and head accelerations again increased in magnitude and variability. Sickness intensity in orbit was not predicted by statistically concordant results of 6 acute preflight susceptibility tests. However, results from a longer duration preflight prism goggles test showed an apparent correlation. All subjects were asymptomatic making head movements in parabolic flight 4 days after the mission, but not 1 year later. Overall, results support the view that space sickness is a motion sickness.

Non-NASA Center↗

Computational Analysis of Gravitational Effects in Low-Density Gas Jets

This study deals with the computational analysis of buoyancy-induced instability in the nearfield of an isothermal helium jet injected into quiescent ambient air environment. Laminar, axisymmetric, unsteady flow conditions were considered for the analysis. The transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum were solved using a staggered grid finite volume method. The jet Richardson numbers of 1.5 and 0.018 were considered to encompass both buoyant and inertial jet flow regimes. Buoyancy effects were isolated by initiating computations in Earth gravity and subsequently, reducing gravity to simulate the microgravity conditions. Computed results concur with experimental observations that the periodic flow oscillations observed in Earth gravity subside in microgravity.

Satti, Rajani P.↗

Schlieren Measurements of Buoyancy Effects on Flow Transition in Low-Density Gas Jets

The transition from laminar to turbulent flow in helium jets discharged into air was studied using Rainbow Schlieren Deflectometry technique. In particular, the effects of buoyancy on jet oscillations and flow transition length were considered. Experiments to simulate microgravity were conducted in the 2.2s drop tower at NASA Glenn Research Center. The jet Reynolds numbers varied from 800 to1200 and the jet Richardson numbers ranged between 0.01 and 0.004. Schlieren images revealed substantial variations in the flow structure during the drop. Fast Fourier Transform (FFT) analysis of the data obtained in Earth gravity experiments revealed the existence of a discrete oscillating frequency in the transition region, which matched the frequency in the upstream laminar regime. In microgravity, the transition occurred farther downstream indicating laminarization of the jet in the absence of buoyancy. The amplitude of jet oscillations was reduced by up to an order of magnitude in microgravity. Results suggest that jet oscillations were buoyancy induced and that the brief microgravity period may not be sufficient for the oscillations to completely subside.

Pasumarthi, Kasyap S.↗

Support of NASA ADR/ Cross-Enterprise NRA Advanced Adiabatic Demagnetization Refrigerators for Continuous Cooling from 10K to 50mK, Development of a Heat Switch

Mechanical heat switches are used in conjunction with sorption refrigerators, adiabatic demagnetization refrigerators and for other cryogenic tasks including the pre-cooling cryogenic systems. They use a mechanical actuator which closes Au plated Cu jaws on an Au plated Cu bar. The thermal conductance in the closed position is essentially independent of the area of the jaws and proportional to the force applied. It varies linearly with T. It is approximately 10mW/K for 200 N at 1.5K. In some applications, the heat switch can be driven from outside the cryostat by a rotating rod and a screw. Such heat switches are available commercially from several sources. In other applications, including systems for space, it is desirable to drive the switch using a cold linear motor, or solenoid. Superconducting windings are used at temperatures s 4.2K to minimize power dissipation, but are not appropriate for pre-cooling a system at higher temperatures. This project was intended to improve the design of solenoid activated mechanical heat switches and to provide such switches as required to support the development of Advanced Adiabatic Demagnetization Refrigerators for Continuous Cooling from 10 K to 50 mK at GSFC. By the time funding began in 5/1/01, the immediate need for mechanical heat switches at GSFC had subsided but, at the same time, the opportunity had arisen to improve the design of mechanical heat switching by incorporating a "latching solenoid". In this device, the solenoid current is required only for changing the state of the switch and not during the whole time that the switch is closed.

Richards, Paul L.↗

Interannual Variations in Tropical Upper-Tropospheric Humidity: Understanding Tropical Convective and Dynamical Processes

Uncertainty remains as to what extent variability in mid to upper tropospheric moisture, especially over the tropics, behaves as constant relative humidity during interannual climate variations associated with ENSO. Systematic variations in HIRS 6.7 micron and MLS 205 GHz suggest that dry subtropical regions evolving during warm SST events depress relative humidity, but the interpretation of these events is still uncertain. Additional specific concerns have to do with regional signatures of convective processes: How does the origin of dry air in the eastern subtropical N. Pacific differ in ENSO warm versus cold years? The dynamics of Rossby wave forcing by convective heating, subtropical jet stream dynamics, and dynamics driven subsidence all come into play here. How variations in precipitating ice hydrometeors from tropical anvils relate to variations in UTH is also a subject of debate? Do variations in precipitating ice, cloud cover and water vapor behavior show any support for the Iris-hypothesis mechanism? Here we examine historical records of SSM/T-2 data to gain a better physical understanding of the effects of deep convective moisture sources and dynamically-induced vertical circulations on UTH. These high frequency microwave measurements (183.3 GHz) take advantage of far less sensitivity to cloud hydrometeors than the 6.7 micron data to yield a record of upper tropospheric relative humidity. Furthermore, signatures of precipitating ice from these channels facilitate comparisons to TRMM hydrometeors detected by radar. In analyzing these observations, we isolate water vapor and temperature change components that affect brightness temperatures and the inferred relative humidity. Trajectory modeling is also used to understand interannual humidity anomalies in terms of outflow fbm convective regions and history of diabatically-driven sinking which modifies relative humidity.

Robertson, Franklin R.↗

Airport Pricing Strategies

Conventional economic wisdom suggests that congestion pricing would be an appropriate response to cope with the growing congestion levels currently experienced at many airports. Several characteristics of aviation markets, however, may make naive congestion prices equal to the value of marginal travel delays a non-optimal response. This paper has developed a model of airport pricing that captures a number of these features. The model in particular reflects that airlines typically have market power and are engaged in oligopolistic competition at different sub-markets; that part of external travel delays that aircraft impose are internal to an operator and hence should not be accounted for in congestion tolls. We presented an analytical treatment for a simple bi-nodal symmetric network, which through the use of 'hyper-networks' would be readily applicable to dynamic problems (in discrete time) such as peak - off-peak differences, and some numerical exercises for the same symmetric network, which was only designed to illustrate the possible comparative static impacts of tolling, in addition to marginal equilibrium conditions as could be derived for the general model specification. Some main conclusions are that second-best optimal tolls are typically lower than what would be suggested by congestion costs alone and may even be negative, and that the toll as derived by Brueckner (2002) may not lead to an increase in total welfare. While Brueckner (2002) has made clear that congestion tolls on airports may be smaller than expected when congestion costs among aircraft are internal for a firm, our analysis adds to this that a further downward adjustment may be in order due to market power. The presence of market power (which causes prices to exceed marginal costs) may cause the pure congestion toll to be suboptimal, because the resulting decrease in demand is too high (the pure congestion tall does not take into account the decrease in consumer surplus). The various downward adjustments in welfare maximizing tolls may well cause the optimal values of these to be negative. Insofar as subsidization is considered unacceptable for whichever reason, our results warn that the most efficient among the non-negative tolls may actually be a zero toll; the pure congestion toll may actually decrease welfare compared to the base case. The model in this paper contains a few simplifying assumptions that may be relaxed in future work. Load factors and aircraft capacity are fixed in this model for simplicity. In a more advanced version of this model, load factors and aircraft capacity can be endogenized. This makes the derivation of the optimality conditions far more complicated, but it should be feasible in a numerical experiment. One can also add a fourth layer to the model, describing the airport's optimization problem. For example, the airport can maximize profits under a cost recovery constraint. The model then deals with interactions between four types of agents. No distinction is made between peak and off-peak traffic in this paper. Finally, the results of the numerical exercise in this paper need to be checked against an asymmetric equilibrium.

Pels, Eric↗