Seawater carbonate chemistry, nutrient uptake and biological processes of coral Stylophora pistillata during experiments, 2011

The effects of ocean acidification and elevated seawater temperature on coral calcification and photosynthesis have been extensively investigated over the last two decades, whereas they are still unknown on nutrient uptake, despite their importance for coral energetics. We therefore studied the separate and combined impacts of increases in temperature and pCO2 on phosphate, ammonium, and nitrate uptake rates by the scleractinian coral S. pistillata. Three experiments were performed, during 10 days i) at three pHT conditions (8.1, 7.8, and 7.5) and normal temperature (26°C), ii) at three temperature conditions (26°, 29°C, and 33°C) and normal pHT(8.1), and iii) at three pHT conditions (8.1, 7.8, and 7.5) and elevated temperature (33°C). After 10 days of incubation, corals had not bleached, as protein, chlorophyll, and zooxanthellae contents were the same in all treatments. However, photosynthetic rates significantly decreased at 33°C, and were further reduced for the pHT 7.5. The photosynthetic efficiency of PSII was only decreased by elevated temperature. Nutrient uptake rates were not affected by a change in pH alone. Conversely, elevated temperature (33°C) alone induced an increase in phosphate uptake but a severe decrease in nitrate and ammonium uptake rates, even leading to a release of nitrogen into seawater. Combination of high temperature (33°C) and low pHT(7.5) resulted in a significant decrease in phosphate and nitrate uptake rates compared to control corals (26°C, pHT = 8.1). These results indicate that both inorganic nitrogen and phosphorus metabolism may be negatively affected by the cumulative effects of ocean warming and acidification.

(Table T1) Density, porosity, void ratio, diffusivity and permeability of ODP Site 185-1149 sediments

The subduction of oceanic plates regulates crustal growth, influences arc volcanism, and refertilizes the mantle. Continental growth occurs by subduction of crustal material (seawater components, marine sediments, and basaltic crust). The geochemical and physical evolution of the Earth's crust depends, in large part, on the fate of subducted material at convergent margins (Armstrong, 1968, doi:10.1029/RG006i002p00175; Karig and Kay, 1981, 10.1098/rsta.1981.0108). The crustal material on the downgoing plate is recycled to various levels in the subduction zone. The recycling process that takes place in the "Subduction Factory" is difficult to observe directly but is clearly illuminated using chemical tracers. Von Huene and Scholl (1991, doi:10.1029/91RG00969) and Plank and Langmuir (1993, doi:10.1038/362739a0) preliminarily calculated a large flux of subducted materials. By mass balancing the chemical tracers and measuring the fractionations that occur between them, the Subduction Factory work and the effect on the Earth's evolution can be estimated. In order to elucidate this mass balance, Ocean Drilling Program Leg 185 drilled two deepwater shales into the oceanic crust situated in the Mariana-Izu Trenches and recovered core samples of incoming oceanic crust. The calculations of mass circulation in the subduction zone, however, did not take into account the mass transfer properties within subducted oceanic crust, although the dewatering fluid and diffused ions may play an important role in various activities such as seismogeneity, serpentine diapiring, and arc volcanism. Thus, this paper focuses on the quantitative measurements of the physical and mass transfer properties of subducted oceanic crust.

Pollen flux off Cape Blanc

Fluxes of airborne freshwater diatoms (FD), phytoliths (PH), and pollen grains (PO) collected with sediment traps off Cape Blanc, northwest Africa, from 1988 till 1991 are presented. Both continental rainfall variations and wind mean strength and direction play a key role in the temporal fluctuations of the fluxes of eolian traces in the pelagic realm. Drier conditions in Northern Africa in 1987 could have preceded the high lithogenic input and moderate FD flux in 1988. The PH peak in summer 1988 was probably caused by increased wind velocity. Wetter rainy seasons of 1988/89 might have promoted a significant pollen production in summer 1989, and FD in late 1989 and early 1990, as well as contributed to the reduction of the lithogenic flux in 1989/90. Decreased fluxes of FD, PH and PO, and higher contribution of the 6-11 µm lithogenic fraction in 1991 would mainly reflect minor intensity and decreased amount of continental trade winds. Air-mass backward trajectories confirm that the Saharan Air Layer is predominantly involved in the spring/summer transport. Trade winds play a decisive role in the fall/winter months, but also contribute to the transport during late spring/summer. Origin of wind trajectories does not support a direct relationship between transporting wind-layers and material source areas in Northern Africa. High winter fluxes of eolian tracers and high amount of trade winds with continental origin in summer warn against a simplistic interpretation of the seasonal eolian signal preserved in the sediments off Cape Blanc, and the wind layer involved in its transport.

Seawater carbonate chemistry, sample density and Strongylocentrotus purpuratus size, filtering and respiration rate during experiments, 2011

Anthropogenic CO2 emissions are acidifying the world's oceans. A growing body of evidence is showing that ocean acidification impacts growth and developmental rates of marine invertebrates. Here we test the impact of elevated seawater pCO2 (129 Pa, 1271 µatm) on early development, larval metabolic and feeding rates in a marine model organism, the sea urchin Strongylocentrotus purpuratus. Growth and development was assessed by measuring total body length, body rod length, postoral rod length and posterolateral rod length. Comparing these parameters between treatments suggests that larvae suffer from a developmental delay (by ca. 8%) rather than from the previously postulated reductions in size at comparable developmental stages. Further, we found maximum increases in respiration rates of + 100 % under elevated pCO2, while body length corrected feeding rates did not differ between larvae from both treatments. Calculating scope for growth illustrates that larvae raised under high pCO2 spent an average of 39 to 45% of the available energy for somatic growth, while control larvae could allocate between 78 and 80% of the available energy into growth processes. Our results highlight the importance of defining a standard frame of reference when comparing a given parameter between treatments, as observed differences can be easily due to comparison of different larval ages with their specific set of biological characters.

222Radon flux map for Europe in netCDF format

A high-resolution 222Radon (222Rn) flux map for Europe was developed, based on a parameterization of 222Rn production and transport in the soil. The 222Rn exhalation rate is parameterized based on soil properties, uranium content, and modelled soil moisture from two different land-surface reanalysis data sets. Spatial variations in exhalation rates are primarily determined by the uranium content of the soil, but also influenced by soil texture and local water table depth. Temporal variations are related to soil moisture variations as the molecular diffusion in the unsaturated soil zone depends on available air-filled pore space. Monthly 222Rn exhalation rates from European soils were calculated with a nominal spatial resolution of 0.083° x 0.083°. The two realizations of the 222Rn flux map, based on the different soil moisture data sets, both realistically reproduce the observed seasonality in the fluxes but yield considerable differences for absolute flux values. The mean 222Rn flux from soils in Europe is estimated to be 10 mBq/m**2/s (ERA-Interim/Land soil moisture) or 15 mBq/m**2/s (GLDAS-Noah soil moisture) for the period 2006-2010. The 222Rn flux maps for Europe are available for the application in atmospheric transport studies, e.g to evaluate the performance of atmospheric transport models.

Oxygen Utilization and Downward Carbon Flux in an Oxygen-Depleted Eddy in the Eastern Tropical North Atlantic

The occurrence of mesoscale eddies that develop suboxic environments at shallow depth (about 40-100 m) has recently been reported for the eastern tropical North Atlantic (ETNA). Their hydrographic structure suggests that the water mass inside the eddy is well isolated from ambient waters supporting the development of severe near-surface oxygen deficits. So far, hydrographic and biogeochemical characterization of these eddies was limited to a few autonomous surveys, with the use of moorings, under water gliders and profiling floats. In this study we present results from the first dedicated biogeochemical survey of one of these eddies conducted in March 2014 near the Cape Verde Ocean Observatory (CVOO). During the survey the eddy core showed oxygen concentrations as low as 5 µmol kg-1 with a pH of around 7.6 at approximately 100 m depth. Correspondingly, the aragonite saturation level dropped to 1 at the same depth, thereby creating unfavorable conditions for calcifying organisms. To our knowledge, such enhanced acidity within near-surface waters has never been reported before for the open Atlantic Ocean. Vertical distributions of particulate organic matter and dissolved organic matter (POM and DOM), generally showed elevated concentrations in the surface mixed layer (0-70 m), with DOM also accumulating beneath the oxygen minimum. With the use of reference data from the upwelling region where these eddies are formed, the oxygen utilization rate was calculated by determining oxygen consumption through the remineralization of organic matter. Inside the core, we found these rates were almost 1 order of magnitude higher (apparent oxygen utilization rate (aOUR); 0.26 µmol kg-1 day-1) than typical values for the open North Atlantic. Computed downward fluxes for particulate organic carbon (POC), were around 0.19 to 0.23 g C m-2 day-1 at 100 m depth, clearly exceeding fluxes typical for an oligotrophic open-ocean setting. The observations support the view that the oxygen-depleted eddies can be viewed as isolated, westwards propagating upwelling systems of their own, thereby represent re-occurring alien biogeochemical environments in the ETNA.

Air-sea fluxes of CO2 and CH4 from the Penlee Point Atmospheric Observatory on the south-west coast of the UK

We present air–sea fluxes of carbon dioxide (CO2), methane (CH4), momentum, and sensible heat measured by the eddy covariance method from the recently established Penlee Point Atmospheric Observatory (PPAO) on the south-west coast of the United Kingdom. Measurements from the south-westerly direction (open water sector) were made at three different sampling heights (approximately 15, 18, and 27m above mean sea level, a.m.s.l.), each from a different period during 2014–2015. At sampling heights ≥18ma.m.s.l., measured fluxes of momentum and sensible heat demonstrate reasonable (≤ ±20% in the mean) agreement with transfer rates over the open ocean. This confirms the suitability of PPAO for air–sea exchange measurements in shelf regions. Covariance air–sea CO2 fluxes demonstrate high temporal variability. Air-to-sea transport of CO2 declined from spring to summer in both years, coinciding with the breakdown of the spring phytoplankton bloom. We report, to the best of our knowledge, the first successful eddy covariance measurements of CH4 emissions from a marine environment. Higher sea-to-air CH4 fluxes were observed during rising tides (20±3; 38±3; 29±6 μmolem-2 d-1 at 15, 18, 27ma.m.s.l.) than during falling tides (14±2; 22±2; 21±5 μmolem-2 d-1), consistent with an elevated CH4 source from an estuarine outflow driven by local tidal circulation. These fluxes are a few times higher than the predicted CH4 emissions over the open ocean and are significantly lower than estimates from other aquatic CH4 hotspots (e.g. polar regions, freshwater). Finally, we found the detection limit of the air–sea CH4 flux by eddy covariance to be 20 μmolem-2 d-1 over hourly timescales (4 μmolem-2 d-1 over 24 h).

Air-sea fluxes of CO2 and CH4 from the Penlee Point Atmospheric Observatory on the south-west coast of the UK

We present air–sea fluxes of carbon dioxide (CO2), methane (CH4), momentum, and sensible heat measured by the eddy covariance method from the recently established Penlee Point Atmospheric Observatory (PPAO) on the south-west coast of the United Kingdom. Measurements from the south-westerly direction (open water sector) were made at three different sampling heights (approximately 15, 18, and 27m above mean sea level, a.m.s.l.), each from a different period during 2014–2015. At sampling heights ≥18ma.m.s.l., measured fluxes of momentum and sensible heat demonstrate reasonable (≤ ±20% in the mean) agreement with transfer rates over the open ocean. This confirms the suitability of PPAO for air–sea exchange measurements in shelf regions. Covariance air–sea CO2 fluxes demonstrate high temporal variability. Air-to-sea transport of CO2 declined from spring to summer in both years, coinciding with the breakdown of the spring phytoplankton bloom. We report, to the best of our knowledge, the first successful eddy covariance measurements of CH4 emissions from a marine environment. Higher sea-to-air CH4 fluxes were observed during rising tides (20±3; 38±3; 29±6 μmolem-2 d-1 at 15, 18, 27ma.m.s.l.) than during falling tides (14±2; 22±2; 21±5 μmolem-2 d-1), consistent with an elevated CH4 source from an estuarine outflow driven by local tidal circulation. These fluxes are a few times higher than the predicted CH4 emissions over the open ocean and are significantly lower than estimates from other aquatic CH4 hotspots (e.g. polar regions, freshwater). Finally, we found the detection limit of the air–sea CH4 flux by eddy covariance to be 20 μmolem-2 d-1 over hourly timescales (4 μmolem-2 d-1 over 24 h).

A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity

Models of the air-sea transfer velocity of gases may be either empirical or mechanistic. Extrapolations of empirical models to an unmeasured gas or to another water temperature can be erroneous if the basis of that extrapolation is flawed. This issue is readily demonstrated for the most well-known empirical gas transfer velocity models where the influence of bubble-mediated transfer, which can vary between gases, is not explicitly accounted for. Mechanistic models are hindered by an incomplete knowledge of the mechanisms of air-sea gas transfer. We describe a hybrid model that incorporates a simple mechanistic view—strictly enforcing a distinction between direct and bubble-mediated transfer—but also uses parameterizations based on data from eddy flux measurements of dimethyl sulphide (DMS) to calibrate the model together with dual tracer results to evaluate the model. This model underpins simple algorithms that can be easily applied within schemes to calculate local, regional, or global air-sea fluxes of gases.

A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity

Models of the air-sea transfer velocity of gases may be either empirical or mechanistic. Extrapolations of empirical models to an unmeasured gas or to another water temperature can be erroneous if the basis of that extrapolation is flawed. This issue is readily demonstrated for the most well-known empirical gas transfer velocity models where the influence of bubble-mediated transfer, which can vary between gases, is not explicitly accounted for. Mechanistic models are hindered by an incomplete knowledge of the mechanisms of air-sea gas transfer. We describe a hybrid model that incorporates a simple mechanistic view—strictly enforcing a distinction between direct and bubble-mediated transfer—but also uses parameterizations based on data from eddy flux measurements of dimethyl sulphide (DMS) to calibrate the model together with dual tracer results to evaluate the model. This model underpins simple algorithms that can be easily applied within schemes to calculate local, regional, or global air-sea fluxes of gases.

Design and Flux-Weakening Control of an Interior Permanent Magnet Synchronous Motor for Electric Vehicles

Permanent magnet synchronous motors (PMSMs) provide a competitive technology for EV traction drives owing to their high power density and high efficiency. In this paper, three types of interior PMSMs with different PM arrangements are modeled by the finite element method (FEM). For a given amount of permanent magnet materials, the V-shape interior PMSM is found better than the U-shape and the conventional rotor topologies for EV traction drives. Then the V-shape interior PMSM is further analyzed with the effects of stator slot opening and the permanent magnet pole chamfering on cogging torque and output torque performance. A vector-controlled flux-weakening method is developed and simulated in Matlab to expand the motor speed range for EV drive system. The results show good dynamic and steady-state performance with a capability of expanding speed up to four times of the rated. A prototype of the V-shape interior PMSM is also manufactured and tested to validate the numerical models built by the FEM.

Multiphase modeling of melting / solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Flame-propagation rates in ammonia-air combustion at high pressure

The reluctance of ammonia to be ignited easily and the necessity to advance the spark for optimum performance in combustion engines led to this investigation. Ionization gap techniques showed that long induction times and slow flame speed were both contributors to the observed phenomena. The most important factor in determining combustion characteristics of ammonia was found to be the extent of predissociation prior to attempted ignition. There was evidence from combustion gas analysis that NO + was the probable most abundant ion contributing to succesful application of the ionization gap techniques, and that the NO + was probably produced in the pyrolysis of ammonia rather than by equilibrium reactions. © 1967 Combustion Institute.

Post-shock thermochemistry in hypervelocity CO2 and air flow

This work represents ongoing efforts to study high-enthalpy carbon dioxide flows in anticipation of the upcoming Mars Science Laboratory (MSL) and future missions to the red planet. The work is motivated by observed anomalies between experimental and numerical studies in hypervelocity impulse facilities for high enthalpy carbon dioxide flows. In this work, experiments are conducted in the Hypervelocity Expansion Tube (HET) which, by virtue of its flow acceleration process, exhibits minimal freestream dissociation in comparison to reflected shock tunnels. This simplifies the comparison with computational result as freestream dissociation and considerable thermochemical excitation can be neglected. Shock shapes of the MSL aeroshell and spherical geometries are compared with numerical simulations incorporating detailed CO2 thermochemical modeling. The shock stand-off distance has been identified in the past as sensitive to the thermochemical state and as such, is used here as an experimental measurable for comparison with CFD and two different theoretical models. It is seen that models based upon binary scaling assumptions are not applicable for the low-density, small-scale conditions of the current work. Mars Science Laboratory shock shapes at zero angle of attack are also in good agreement with available data from the LENS X expansion tunnel facility, confi rming results are facility-independent for the same type of flow acceleration, and indicating that the flow velocity is a suitable first-order matching parameter for comparative testing. In an e ffort to address surface chemistry issues arising from high-enthalpy carbon dioxide ground-test based experiments, spherical stagnation point and aeroshell heat transfer distributions are also compared with simulation. Very good agreement between experiment and CFD is seen for all shock shapes and heat transfer distributions fall within the non-catalytic and super-catalytic solutions. We also examine spatial temperature profiles in the non-equilibrium relaxation region behind a stationary shock wave in a hypervelocity air Mach 7.42 freestream. The normal shock wave is established through a Mach reflection from an opposing wedge arrangement. Schlieren images confirm that the shock con guration is steady and the location is repeatable. Emission spectroscopy is used to identify dissociated species and to make vibrational temperature measurements using both the nitric oxide and the hydroxyl radical A-X band sequences. Temperature measurements are presented at selected locations behind the normal shock. LIFBASE is used as the simulation spectrum software for OH temperature-fitting, however the need to access higher vibrational and rotational levels for NO leads to the use of an in-house developed algorithm. For NO, results demonstrate the contribution of higher vibrational and rotational levels to the spectra at the conditions of this study. Very good agreement is achieved between the experimentally measured NO vibrational temperatures and calculations performed using an existing state-resolved, three-dimensional forced harmonic oscillator thermochemical model. The measured NO A-X vibrational temperatures are significantly higher than the OH A-X temperatures.