Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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 613 records · Page 34

Analysis of synthetic aperture radar imagery

The author has identified the following significant results. Average radar response for L-band like polarized system appeared to be related to the watershed runoff coefficients when the viewing angle was approximately 42 deg off nadir. Four requirements for radar systems used to verify applications of active microwave for water resources were identified: (1) first generation digital data will be required; (2) radar should be calibrated both internally and externally; (3) new systems should avoid radom use; and (4) images should be geometrically rectified prior to delivery to the user.

Blanchard, B. J.↗

Passive and active remote sensing of atmospheric precipitation

Both passive and active remote sensing of atmospheric precipitation are studied with the vector radiative transfer equations by making use of the Mie scattering phase functions and incorporating the raindrop-size distributions. For passive remote sensing, the Gaussian quadrature method is employed to solve for the brightness temperatures. For active remote sensing, an iterative approach carrying out to the second order in albedo is used to calculate for the bistatic coefficients, the backscattering cross sections/unit volume, and the interchannel cross talks. The calculated results are plotted as a function of rainfall rates and compared to various available experimental data. The theoretical model is easily applied to the remote sensing of aerosol particles, smoke, fog, and haze at infrared and visible frequencies.

Jin, Y. Q.↗

Mineralogical and geochemical anomalous data of the K-T boundary samples

Cretaceous-Tertiary boundary problem has been discussed previously from the geological research, mainly by fossil changes. Although geochemical bulk data of Ir anomaly suggest the extraterrestrial origin of the K-T boundary, the exact formation process discussed mainly by mineralogical and geochemical study has been started recently, together with noble gas contents. The K-T boundary sample at Kawaruppu River, Hokkaido was collected, in order to compare with the typical K-T boundary samples of Bubbio, Italy, Stevns Klint, Denmark, and El Kef, Tunisia. The experimental data of the silicas and calcites in these K-T boundary samples were obtained from the X-ray unit-cell dimension (i.e., density), ESR signal and total linear absorption coefficient, as well as He and Ne contents. The K-T boundary samples are usually complex mixture of the terrestrial activities after the K-T boundary event. The mineralogical and geochemical anomalous data indicate special terrestrial atmosphere at the K-T boundary formation probably induced by asteroid impact, followed the many various terrestrial activities (especially the strong role of sea-water mixture, compared with terrestrial highland impact and impact craters in the other earth-type planetary bodies).

Miura, Y.↗

Preliminary analysis of ERS-1 SAR for forest ecosystem studies

An image obtained by the C-band VV-polarized ERS-1 SAR is examined with respect to potential land applications. A scene obtained near noon on August 15, 1991, along the U.S.-Canadian border near Sault Ste. Marie is calibrated relative to an array of trihedral corner reflectors and active radar calibrators distributed across the swath. Extensive contemporaneous ground observations of forest stands are used to predict the radar backscattering coefficient sigma at the time of the SAR overpass using a first-order vector radiative transfer model. These predictions generally agree with the calibrated ERS-1 data to within 1 dB. It is demonstrated that the dynamic range of sigma is sufficient to perform limited discrimination of various forest and grassland communities even for a single-date observation. Furthermore, retrieval of near-surface soil moisture is feasible for grass-covered soils when plant biomass is less than 1 tonne/ha.

Dobson, M. G.↗

GPS-Based Reduced Dynamic Orbit Determination Using Accelerometer Data

Currently two gravity field satellite missions, CHAMP and GRACE, are equipped with high sensitivity electrostatic accelerometers, measuring the non-conservative forces acting on the spacecraft in three orthogonal directions. During the gravity field recovery these measurements help to separate gravitational and non-gravitational contributions in the observed orbit perturbations. For precise orbit determination purposes all these missions have a dual-frequency GPS receiver on board. The reduced dynamic technique combines the dense and accurate GPS observations with physical models of the forces acting on the spacecraft, complemented by empirical accelerations, which are stochastic parameters adjusted in the orbit determination process. When the spacecraft carries an accelerometer, these measured accelerations can be used to replace the models of the non-conservative forces, such as air drag and solar radiation pressure. This approach is implemented in a batch least-squares estimator of the GPS High Precision Orbit Determination Software Tools (GHOST), developed at DLR/GSOC and DEOS. It is extensively tested with data of the CHAMP and GRACE satellites. As accelerometer observations typically can be affected by an unknown scale factor and bias in each measurement direction, they require calibration during processing. Therefore the estimated state vector is augmented with six parameters: a scale and bias factor for the three axes. In order to converge efficiently to a good solution, reasonable a priori values for the bias factor are necessary. These are calculated by combining the mean value of the accelerometer observations with the mean value of the non-conservative force models and empirical accelerations, estimated when using these models. When replacing the non-conservative force models with accelerometer observations and still estimating empirical accelerations, a good orbit precision is achieved. 100 days of GRACE B data processing results in a mean orbit fit of a few centimeters with respect to high-quality JPL reference orbits. This shows a slightly better consistency compared to the case when using force models. A purely dynamic orbit, without estimating empirical accelerations thus only adjusting six state parameters and the bias and scale factors, gives an orbit fit for the GRACE B test case below the decimeter level. The in orbit calibrated accelerometer observations can be used to validate the modelled accelerations and estimated empirical accelerations computed with the GHOST tools. In along track direction they show the best resemblance, with a mean correlation coefficient of 93% for the same period. In radial and normal direction the correlation is smaller. During days of high solar activity the benefit of using accelerometer observations is clearly visible. The observations during these days show fluctuations which the modelled and empirical accelerations can not follow.

VanHelleputte, Tom↗

DSN Data Visualization Suite

The DSN Data Visualization Suite is a set of computer programs and reusable Application Programming Interfaces (APIs) that assist in the visualization and analysis of Deep Space Network (DSN) spacecraft-tracking data, which can include predicted and actual values of downlink frequencies, uplink frequencies, and antenna-pointing angles in various formats that can include tables of values and polynomial coefficients. The data can also include lists of antenna-pointing events, lists of antenna- limit events, and schedules of tracking activities. To date, analysis and correlation of these intricately related data before and after tracking have been difficult and time-consuming. The DSN Data Visualization Suite enables operators to quickly diagnose tracking-data problems before, during, and after tracking. The Suite provides interpolation on demand and plotting of DSN tracking data, correlation of all data on a given temporal point, and display of data with color coding configurable by users. The suite thereby enables rapid analysis of the data prior to transmission of the data to DSN control centers. At the control centers, the same suite enables operators to validate the data before committing the data to DSN subsystems. This software is also Web-enabled to afford its capabilities to international space agencies.

Bui, Bach X.↗

Creation of a Rapid High-Fidelity Aerodynamics Module for a Multidisciplinary Design Environment

In the traditional aerospace vehicle design process, each successive design phase is accompanied by an increment in the modeling fidelity of the disciplinary analyses being performed. This trend follows a corresponding shrinking of the design space as more and more design decisions are locked in. The correlated increase in knowledge about the design and decrease in design freedom occurs partly because increases in modeling fidelity are usually accompanied by significant increases in the computational expense of performing the analyses. When running high fidelity analyses, it is not usually feasible to explore a large number of variations, and so design space exploration is reserved for conceptual design, and higher fidelity analyses are run only once a specific point design has been selected to carry forward. The designs produced by this traditional process have been recognized as being limited by the uncertainty that is present early on due to the use of lower fidelity analyses. For example, uncertainty in aerodynamics predictions produces uncertainty in trajectory optimization, which can impact overall vehicle sizing. This effect can become more significant when trajectories are being shaped by active constraints. For example, if an optimal trajectory is running up against a normal load factor constraint, inaccuracies in the aerodynamic coefficient predictions can cause a feasible trajectory to be considered infeasible, or vice versa. For this reason, a trade must always be performed between the desired fidelity and the resources available. Apart from this trade between fidelity and computational expense, it is very desirable to use higher fidelity analyses earlier in the design process. A large body of work has been performed to this end, led by efforts in the area of surrogate modeling. In surrogate modeling, an up-front investment is made by running a high fidelity code over a Design of Experiments (DOE); once completed, the DOE data is used to create a surrogate model, which captures the relationships between input variables and responses into regression equations. Depending on the dimensionality of the problem and the fidelity of the code for which a surrogate model is being created, the initial DOE can itself be computationally prohibitive to run. Cokriging, a modeling approach from the field of geostatistics, provides a desirable compromise between computational expense and fidelity. To do this, cokriging leverages a large body of data generated by a low fidelity analysis, combines it with a smaller set of data from a higher fidelity analysis, and creates a kriging surrogate model with prediction fidelity approaching that of the higher fidelity analysis. When integrated into a multidisciplinary environment, a disciplinary analysis module employing cokriging can raise the analysis fidelity without drastically impacting the expense of design iterations. This is demonstrated through the creation of an aerodynamics analysis module in NASA’s OpenMDAO framework. Aerodynamic analyses including Missile DATCOM, APAS, and USM3D are leveraged to create high fidelity aerodynamics decks for parametric vehicle geometries, which are created in NASA’s Vehicle Sketch Pad (VSP). Several trade studies are performed to examine the achieved level of model fidelity, and the overall impact to vehicle design is quantified.

Aerodynamic Properties↗

A Fully Passive Daily Load-follow Operation in a Small PWR System

An innovative fully passive load-follow operation (PLFO) concept, based only on the negative reactivity feedbacks, is proposed by the authors. The new design concept requires no control rod adjustment or boron utilization for achieving load-follow operation in pressurized water reactors (PWRs). PLFO is motivated by the need for developing next-generation PWRs with enhanced power maneuvering capabilities, increased level of autonomism, and less dependence on active controllers. Enhanced power maneuvering capabilities are required to allow next-generation PWRs to actively contribute to electrical grid stability. Meanwhile, increased level of autonomism and less dependence on active controllers are required for safety and reliability of next-generation PWRs. However, developing an autonomous thermal reactor is challenging due to its strong thermal and xenon (Xe) reactivity feedback. Therefore, enhancing power maneuvering capabilities and simultaneously increasing the level of autonomy with less dependence on active controllers is very challenging in PWRs. However, PLFO is a promising solution. For the feasibility of PLFO, the coolant temperature coefficient (CTC) should be always clearly negative and this requires a soluble boron-free (SBF) coolant system. An SBF coolant system is also favorable in terms of reactor safety. Nevertheless, utilizing soluble boron in load-following operation is not favorable, as the change in boron concentration is a quite slow process. The power demand can be either regular (repeated pattern in a daily-basis) or irregular (varying demand). The feasibility of PLFO for regular demand was demonstrated by the authors in other works. This paper is dedicated for summarizing the feasibility of PLFO for irregular demand, and to point out the impact of Xe reactivity on PLFO for a varying power demand.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimizing Insulation Design for Transformers in Medium Voltage Power Conversion Systems

Medium-frequency transformers (MFTs) play a crucial role in medium-voltage (MV) solidstate transformer (SST) systems, particularly in extreme fast charging applications. Achieving partial discharge (PD)-free operation while maintaining high power density is a significant challenge due to the high electric field (E-field) stresses inherent in MV applications. This dissertation focuses on the insulation design and optimization of MFTs used in both the main power electronics circuits and auxiliary power supplies. The study begins with an overview of insulation testing methodologies, including high potential tests, basic insulation level tests, and PD tests, which are critical for evaluating MFT insulation reliability. Given the importance of PD-free operation for long-term reliability, particular emphasis is placed on understanding PD mechanisms, including void, corona, and surface discharge, and their mitigation strategies. A high voltage isolated auxiliary power supply is then introduced, utilizing a gapped transformer encapsulated in silicone gel. This design achieves PD-free insulation up to 18 kV RMS while maintaining low coupling capacitance to minimize common-mode current. The proposed solution ensures reliable operation in MV environments and offers a scalable approach for auxiliary power in cascaded SST architectures. To improve MFT insulation in main power conversion circuits, a novel structure is developed using polypropylene sheets and potting compounds to create a void-free air gap, effectively mitigating E-field intensity. A prototype transformer with this insulation structure is built and achieves PD-free operation up to 30 kV RMS. This design is experimentally validated in a resonant converter operating at 46 kW, demonstrating its feasibility for MV SST applications. Further optimization is implemented to enhance MFT performance for dual-active-bridge(DAB) converters by integrating a semiconductive shielding layer within the insulation structure. This shielding layer improves the magnetic coupling coefficient while effectively confining the E-field within high insulation materials, thereby reducing eddy current losses. The optimized MFT achieves PD-free operation at 12.6 kV RMS and is successfully tested in a DAB converter operating at 43 kW, which meets the insulation requirements for a 13.2 kV SST system. This dissertation advances MFT insulation design by introducing and experimentally validating novel approaches that improve high voltage insulation while optimizing magnetic coupling and manufacturability. The proposed insulation structures enable PD-free operation while minimizing insulation material usage and simplifying assembly, making them ideal for high power, high voltage applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Solar Activity Relations in Energetic Electron Events Measured By the MESSENGER Mission

Aims. We perform a statistical study of the relations between the properties of solar energetic electron (SEE) events measured by the MESSENGER mission from 2010 to 2015 and the parameters of the respective parent solar activity phenomena in order to identify the potential correlations between them. During the time of analysis, the MESSENGER heliocentric distance varied between 0.31 and 0.47 au. Methods. We used a published list of 61 SEE events measured by MESSENGER, which includes information on the near-relativistic electron peak intensities, the peak-intensity energy spectral indices, and the measured X-ray peak intensity of the flares related to the SEE events. Taking advantage of multi-viewpoint remote-sensing observations, we reconstructed, whenever possible, the associated coronal mass ejections (CMEs) and shock waves; and we determined the three-dimensional (3D) properties (location, speed, and width) of the CMEs and the maximum speed of the 3D CME-driven shocks in the corona. We used different methods (Spearman, Pearson, and a Bayesian approach, namely the Kelly method to linear regression) to estimate the correlation coefficients between the flare intensity, maximum speed at the apex of the CME-driven shock, CME speed at the apex, and CME width with the electron peak intensities and with the energy spectral indices. In this statistical study, we considered and addressed the limitations of the particle instrument on board MESSENGER (elevated background intensity level, anti-Sun pointing). Results. There is an asymmetry to the east in the range of connection angles (CAs) for which the SEE events present the highest peak intensities, where the CA is the longitudinal separation between the footpoint of the magnetic field connecting to the spacecraft and the flare location. Based on this asymmetry, we define a subsample of well-connected events as when −65° ≤ CA ≤ +33°. For the well-connected sample, we find moderate to strong correlations between the near-relativistic electron peak intensity and the 3D CME-driven shock maximum speed at the apex (Spearman: cc = 0.53 ± 0.05; Pearson: cc = 0.65 ± 0.04; Kelly: cc = 0.87 ± 0.20), the flare peak intensity (Spearman: cc = 0.63 ± 0.03; Pearson: cc = 0.59 ± 0.03; Kelly: cc = 0.74 ± 0.30), and the 3D CME speed at the apex (Spearman: cc = 0.50 ± 0.04; Pearson: cc = 0.46 ± 0.03; Kelly: cc = 0.60 ± 0.39). When including poorly connected events (full sample), the relations between the peak intensities and the solar-activity phenomena are blurred, showing lower correlation coefficients. Conclusions. Based on the comparison of the correlation coefficients presented in this study using near 0.4 au data, (1) both flare and shock-related processes may contribute to the acceleration of near relativistic electrons in large SEE events, in agreement with previous studies based on near 1 au data; and (2) the maximum speed of the CME-driven shock is a better parameter to investigate particle-acceleration-related mechanisms than the average CME speed, as suggested by the stronger correlation with the SEE peak intensities.

Sun↗

In General, the Total Voltammetric Current from a Mixture of Redox-Active Substances will Not be the Sum of the Currents that Each Substance would Produce Independently at the Same Concentration as in the Mixture

At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.

Leventis, Nicholas↗

Translational diffusion in supercooled water at and near the glass transition temperature—136 K

The properties of amorphous solid water at and near the calorimetric glass transition temperature, T g , of 136 K have been debated for years. One hypothesis is that water turns into a “true” liquid at T g (i.e., it becomes ergodic) and exhibits all the characteristics of an ergodic liquid, including translational diffusion. A competing hypothesis is that only rotational motion becomes active at T g , while the “real” glass transition in water is at a considerably higher temperature. To address this dispute, we have investigated the diffusive mixing in nanoscale water films, with thicknesses up to ∼100 nm, using infrared (IR) spectroscopy. The experiments used films that were composed of at least 90% H 2 O with D 2 O making up the balance and were conducted under conditions where H/D exchange was essentially eliminated. Because the IR spectra of multilayer D 2 O films (e.g., thicknesses of ∼3–6 nm) embedded within thick H 2 O films are distinct from the spectrum of isolated D 2 O molecules within H 2 O, the diffusive mixing of (initially) isotopically layered water films could be followed as a function of annealing time and temperature. The results show that water films with total thicknesses ranging from ∼20 to 100 nm diffusively mixed prior to crystallization for temperatures between 120 and 144 K. The translational diffusion had an Arrhenius temperature dependence with an activation energy of 40.8 ± 3.5 kJ/mol, which indicates that water at and near T g is a strong liquid. The measured diffusion coefficient at 136 K is 6.25 ± 1.4 × 10 −21 m 2 /s.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Step-Down Test Procedure for Wavelet Shrinkage Using Bootstrapping

Wavelet thresholding (or shrinkage) attempts to remove the noises existing in the signals while preserving inherent pattern characteristics in the reconstruction of true signals. For data-denoising purpose, we present a new wavelet thresholding procedure which employs the step-down testing idea of identifying active contrasts in unreplicated fractional factorial experiments. The proposed method employs bootstrapping methods to a step-down test for thresholding wavelet coefficients. By introducing the concept of a false discovery error rate in testing wavelet coefficients, we shrink the wavelet coefficients with p -values higher than the error rate. The error rate controls the expected proportion of wrongly accepted coefficients among chosen wavelet coefficients. Bootstrap samples are used to approximate the p -value for computational efficiency. We also present some guidelines for selecting the values of hyper-parameters which affect the performance in the step-down thresholding procedure. Based on some common testing signals and an air-conditioner sounds example, the comparison of our proposed procedure with other thresholding methods in the literature is performed. The analytical results show that the proposed procedure has a potential in data-denoising and data-reduction in a variety of signal reconstruction applications.

42 ENGINEERING↗

Documentation of GEMASS entry to touchdown simulation

The entry-to-touchdown space shuttle orbiter simulation capability incorporated into the GEMASS subprogram 33 (3-DOF) is documented. A digital autopilot interfaces between GEMASS and the guidance. Vehicle attitude is determined by use of the ability of GEMASS to integrate differential equations in addition to the equations of motion. Vehicle aerodynamic characteristics are obtained from an aerodynamic data tape and control surface deflections required to trim the vehicle, and trimmed aerodynamic coefficients are determined internally. Several indicators to allow evaluation of subsystem performance are included; as is the ability of the user to activate any of several error dispersion sources. The performance of the simulation compares well with more sophisticated simulations.

Waibel, R. H.↗

Diffusion of triglycine sulfate in aqueous solution

The diffusion coefficient of triglycine sulfate (TGS) in water was measured for several concentrations over a temperature range of 25 to 55 C. The activation energy for diffusion obtained from these measurements was 4180 cal/mol. No concentration dependence was seen. The maximum difference in D for the various ionic species present was determined by Raman spectroscopy to be about 5 percent.

Kroes, R. L.↗

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

Performance of MODIS Thermal Emissive Bands On-orbit Calibration Algorithms

Two nearly identical copies of the Moderate Resolution Imaging Spectroradiometer (MODIS) are currently operated on-board the Terra and Aqua spacecrafts, launched in December 1999 and May 2002, respectively. Together, they have produced an unprecedented amount of science data products, which are widely used for the studies of changes in the Earth's system of land, oceans, and atmosphere. MODIS is a cross-track scanning radiometer, which uses a two-sided scan mirror and collects data continuously over a wide scan angle range (+/-55 degree relative to the instrument nadir) each scan of 1.47 seconds. It has 36 spectral bands with wavelengths ranging from visible (VIS) to long-wave infrared (LWIR). MODIS bands 1-19 and 26 are the reflective solar bands (RSB) and bands 20-25 and 27-36 are the thermal emissive bands (TEB). MODIS was developed and designed with improvements made over its heritage sensors (such as AVHRR and Landsat) and, in particular, with more stringent calibration requirements. Because of this, MODIS was built with a set of state-of-art on-board calibrators (OBC), which include a solar diffuser (SD), a solar diffuser stability monitor (SDSM), a blackbody (BB), a spectroradiometric calibration assembly (SRCA), and a space view (SV) port. With the exception of view angle differences, MODIS OBC measurements and the Earth View (EV) observations are made via the same optical path. MODIS TEB have a total of 160 individual TEB detectors (10 per band), which are located on two cold focal plane assemblies (CFPA). For nominal on-orbit operation, the CFPA temperature is controlled at 83K via a passive radiative cooler. For the TEB, the calibration requirements at specified typical scene radiances are less than or equal to 1% with an exception for the fire detection (low gain) band. MODIS TEB on-orbit calibration is performed on a scan-by-scan basis using a quadratic calibration algorithm, and data collected from sensor responses to the onboard BB and SV. The BB serves as the thermal calibration source and the SV provides measurements for the sensor's background and offsets. MODIS on-board BB is a v-grooved plate with its temperature measured using 12 platinum resistive thermistors (PRT) uniformly embedded in the BB substrate. All the BB thermistors were characterized pre-launch with reference to the NIST temperature standards. Unlike typical BB operations in many heritage sensors, which have no temperature control capability, the MODIS on-board BB can be operated at any temperatures between instrument ambient (about 270K) and 315K and can also be varied continuously within this range. This feature has significantly enhanced the MODIS' capability of tracking and updating the TEB nonlinear calibration coefficients over its entire mission. Following a brief description of MODIS TEB on-orbit calibration methodologies and its onboard BB operational activities, this paper provides a comprehensive performance assessment of MODIS TEB quadratic calibration algorithm. It examines the scan-by-scan, orbit-by-orbit, daily, and seasonal variations of detector responses and associated impact due changes in the CFPA and instrument temperatures. Specifically, this paper will analyze the contribution by each individual thermal emissive source term (BB, scan cavity, and scan mirror), the impact on the Level 1 B data product quality due to pre-launch and on-orbit calibration uncertainties. A comparison of Terra and Aqua TEB on-orbit performance, lessons learned, and suggestions for future improvements will also be made.

Xiong, Xiaoxiong↗