Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coherent errors”

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.

333 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Planar Doppler Velocimetry for Large-Scale Wind Tunnel Testing

Recently, Planar Doppler Velocimetry (PDV) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike other optical means of measuring flow velocities, PDV is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDV have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDV in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDV capabilities for applications in large-scale wind tunnels that are intended primarily for the production testing of subsonic aircraft. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes and separations. To accomplish the objectives indicated, the PDV measurement process is first modeled and its performance evaluated computationally. The noise sources considered include those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with coherent illumination from pulsed lasers. The signal noise estimates are incorporated into the PDV signal analysis process and combined with computed scattering signals using a Mie scattering theory for polydisperse smoke particles. The relevant parameters incorporate a range of practical aerodynamic test conditions and facility sizes. The results define the optimum instrument configurations, show that the expected signal levels from a practical PDV system are sufficiently large to allow its useful application in large facilities, and show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing. Experimental studies using several experimental bench-top setups are then described that validate the physics of the PDV model and to calibrate its computed results. The validated model allows estimates of the uncertainties of PDV measurements and a complete definition of the PDV capabilities to be made with sufficient confidence to decide the viability of PDV for large-scale wind tunnel applications.

McKenzie, Robert L.↗
NASA NTRSDate not supplied

Navigation Doppler Lidar Signal Processing Architecture and Algorithms

A coherent Doppler lidar has been developed to address the need for a high-performance, compact, and cost-effective velocity and altitude sensor onboard landing vehicles. One of the core subsystems within NDL is its signal processing architecture. The primary function of this subsystem is to digitize the analog output of the homodyned signal received from the photo-diode associated with each telescope, process this signal in the frequency domain, and then turn this information into velocity and range measurements. A secondary function of this subsystem is to setup and control the internal components within NDL and to enable the instrument to interface with external systems by receiving commands, sending telemetry, and synchronizing its measurements with external clocks and pulses. This paper will give an overview of the signal processing architecture within NDL. This will include a discussion on its fundamental processing components and the core signal processing algorithms. Command, telemetry and internal component controls will also be discussed as will how the signal processing architecture has also evolved over the years of development. We will discuss provisions for robust operation including Triple-Mode Redundancy (TMR) within NDL’s Field Programmable Gate Array (FPGA) and storage elements, a fault-tolerant boot process, and Error Detection and Correction (EDAC) on the Static Random-Access Memories (SRAMs). Finally, we conclude this paper with how this architecture could continue to evolve.

Doppler Lidar↗
NASA NTRS2012

Measurement Via Optical Near-Nulling and Subaperture Stitching

A subaperture stitching interferometer system provides near-nulling of a subaperture wavefront reflected from an object of interest over a portion of a surface of the object. A variable optical element located in the radiation path adjustably provides near-nulling to facilitate stitching of subaperture interferograms, creating an interferogram representative of the entire surface of interest. This enables testing of aspheric surfaces without null optics customized for each surface prescription. The surface shapes of objects such as lenses and other precision components are often measured with interferometry. However, interferometers have a limited capture range, and thus the test wavefront cannot be too different from the reference or the interference cannot be analyzed. Furthermore, the performance of the interferometer is usually best when the test and reference wavefronts are nearly identical (referred to as a null condition). Thus, it is necessary when performing such measurements to correct for known variations in shape to ensure that unintended variations are within the capture range of the interferometer and accurately measured. This invention is a system for nearnulling within a subaperture stitching interferometer, although in principle, the concept can be employed by wavefront measuring gauges other than interferometers. The system employs a light source for providing coherent radiation of a subaperture extent. An object of interest is placed to modify the radiation (e.g., to reflect or pass the radiation), and a variable optical element is located to interact with, and nearly null, the affected radiation. A detector or imaging device is situated to obtain interference patterns in the modified radiation. Multiple subaperture interferograms are taken and are stitched, or joined, to provide an interferogram representative of the entire surface of the object of interest. The primary aspect of the invention is the use of adjustable corrective optics in the context of subaperture stitching near-nulling interferometry, wherein a complex surface is analyzed via multiple, separate, overlapping interferograms. For complex surfaces, the problem of managing the identification and placement of corrective optics becomes even more pronounced, to the extent that in most cases the null corrector optics are specific to the particular asphere prescription and no others (i.e. another asphere requires completely different null correction optics). In principle, the near-nulling technique does not require subaperture stitching at all. Building a near-null system that is practically useful relies on two key features: simplicity and universality. If the system is too complex, it will be difficult to calibrate and model its manufacturing errors, rendering it useless as a precision metrology tool and/or prohibitively expensive. If the system is not applicable to a wide range of test parts, then it does not provide significant value over conventional null-correction technology. Subaperture stitching enables simpler and more universal near-null systems to be effective, because a fraction of a surface is necessarily less complex than the whole surface (excepting the extreme case of a fractal surface description). The technique of near-nulling can significantly enhance aspheric subaperture stitching capability by allowing the interferometer to capture a wider range of aspheres. More over, subaperture stitching is essential to a truly effective near-nulling system, since looking at a fraction of the surface keeps the wavefront complexity within the capability of a relatively simple nearnull apparatus. Furthermore, by reducing the subaperture size, the complexity of the measured wavefront can be reduced until it is within the capability of the near-null design.

Forbes, Greg↗
NASA NTRS1997

Measurement Capabilities of Planar Doppler Velocimetry in Large-Scale Wind Tunnels

Over the past few years, Planar Doppler Velocimetry (PDT) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike some other optical means of measuring flow velocities, PDT is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDT (also called Doppler Global Velocimetry by some authors) have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDT in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDT capabilities for applications in large-scale wind tunnels that are intended primarily for the testing of rotorcraft and subsonic aircraft at speeds typically less than 100 m/s. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region and with a field of view that is sufficient to encompass the entire region of interest. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes, separations, and rotor-tip vortices. To accomplish the objectives of these studies, the PDT measurement process was first modeled and its performance evaluated computationally. The noise sources considered included those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with the coherent illumination of aerosols from pulsed lasers. The signal noise estimates were incorporated into the PDT signal analysis process and combined with computed scattering signals using a Mae scattering theory for a distributed range of aerosol particle sizes. The results are used to define the necessary instrument configuration, to show that the expected signal levels from a practical PDT system are sufficiently large to allow its useful application in large facilities, and to show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing, Experimental studies using several bench-top setups are then described that validate the physics of the PDT model and demonstrate the ability to obtain accurate PDT measurements using procedures that are compatible with large-scale wind tunnel operations. The PDT measurement capabilities found in this study lead to the conclusion that PDT offers significant advantages compared to other means of measuring velocity fields in large-scale wind tunnels, including those operating at low speeds.

McKenzie, Robert L.↗
NASA NTRS2012

Absolute Position of Targets Measured Through a Chamber Window Using Lidar Metrology Systems

Lidar is a useful tool for taking metrology measurements without the need for physical contact with the parts under test. Lidar instruments are aimed at a target using azimuth and elevation stages, then focus a beam of coherent, frequency modulated laser energy onto the target, such as the surface of a mechanical structure. Energy from the reflected beam is mixed with an optical reference signal that travels in a fiber path internal to the instrument, and the range to the target is calculated based on the difference in the frequency of the returned and reference signals. In cases when the parts are in extreme environments, additional steps need to be taken to separate the operator and lidar from that environment. A model has been developed that accurately reduces the lidar data to an absolute position and accounts for the three media in the testbed air, fused silica, and vacuum but the approach can be adapted for any environment or material. The accuracy of laser metrology measurements depends upon knowing the parameters of the media through which the measurement beam travels. Under normal conditions, this means knowledge of the temperature, pressure, and humidity of the air in the measurement volume. In the past, chamber windows have been used to separate the measuring device from the extreme environment within the chamber and still permit optical measurement, but, so far, only relative changes have been diagnosed. The ability to make accurate measurements through a window presents a challenge as there are a number of factors to consider. In the case of the lidar, the window will increase the time-of-flight of the laser beam causing a ranging error, and refract the direction of the beam causing angular positioning errors. In addition, differences in pressure, temperature, and humidity on each side of the window will cause slight atmospheric index changes and induce deformation and a refractive index gradient within the window. Also, since the window is a dispersive media, the effect of both phase and group indices have to be considered. Taking all these factors into account, a method was developed to measure targets through multiple regions of different materials and produce results that are absolute measurements of target position in three-dimensional space, rather than simply relative position. The environment in which the lidar measurements are taken must be broken down into separate regions of interest and each region solved for separately. In this case, there were three regions of interest: air, fused silica, and vacuum. The angular position of the target inside the chamber is solved using only phase index and phase velocity, while the ranging effects due to travel from air to glass to vacuum/air are solved with group index and group velocity. When all parameters are solved simultaneously, an absolute knowledge of the position of each target within an environmental chamber can be derived. Novel features of this innovation include measuring absolute position of targets through multiple dispersive and non-dispersive media, deconstruction of lidar raw data from a commercial off-the-shelf unit into reworkable parameters, and use of group velocities to reduce range data. Measurement of structures within a vacuum chamber or other harsh environment, such as a furnace, may now be measured as easily as if they were in an ambient laboratory. This analysis permits transformation of the raw data into absolute spatial units (e.g., mm). This technique has also been extended to laser tracker, theodolite, and cathetometer measurements through refractive media.

Kubalak, David↗
NASA NTRS2016

The New Sun-Sky-Lunar Cimel CE318-T Multiband Photometer - A Comprehensive Performance Evaluation

This paper presents the new photometer CE318-T, able to perform daytime and night-time photometric measurements using the sun and the moon as light source. Therefore,this new device permits a complete cycle of diurnal aerosol and water vapour measurements valuable to enhance atmospheric monitoring to be extracted. In this study wehave found significantly higher precision of triplets when comparing the CE318-T master instrument and the Cimel AErosol RObotic NET work (AERONET) master (CE318-AERONET) triplets as a result of the new CE318-T tracking system. Regarding the instrument calibration, two new methodologies to transfer the calibration from a reference instrument using only daytime measurements (Sun Ratio and Sun-Moon gain factor techniques) are presented and discussed. These methods allow the reduction of the previous complexities inherent to nocturnal calibration. A quantitative estimation of CE318-T AOD uncertainty by means of error propagation theory during daytime revealed AOD uncertainties (u(sup D)(sub AOD)) for Langley-calibrated instruments similar to the expected values for other reference instruments (0.002-0.009). We have also found u(sup D)(sub AOD) values similar to the values reported in sun photometry for field instruments (approximately 0.015). In the case of the night-time period, the CE318-T-estimated standard combined uncertainty (u(sup N)(sub AOD)) is dependent not only on the calibration technique but also on illumination conditions and the instrumental noise. These values range from 0.011-0.018 for Lunar Langley-calibrated instruments to 0.012-0.021 for instruments calibrated using the Sun Ratio technique. In the case of moon-calibrated instruments using the Sun-Moon gain factor method and sun calibrated using the Langley technique, we found u(sup N)(sub AOD) ranging from 0.016 to 0.017 (up to 0.019 in 440 nm channel), not dependent on any lunar irradiance model. A subsequent performance evaluation including CE318-T and collocated measurements from independent reference instruments has served to assess the CE318-T performance as well as to confirm its estimated uncertainty. Daytime AOD evaluation, performed at Izana station from March to June 2014, encompassed measurements from a reference CE318-T, a CE318-AERONET master instrument, a Precision Filter Radiometer (PFR) and a Precision Spectroradiometer (PSR) prototype, reporting low AOD discrepancies between the four instruments (up to 0.006). The nocturnal AOD evaluation was performed using CE318-T- and starphotometer-collocated measurements and also by means of a day/night coherence transition test using the CE318-T master instrument and the CE318 daytime data from the CE318-AERONET master instrument. Results showed low discrepancies with the star photometer at 870 and 500 nm channels(less tna or equal to 0.013) and differences with AERONET daytime data (1 h after and before sunset and sunrise) in agreement with the estimated u(sup N)(sub AOD) values at all illumination conditions in the case of channels within the visible spectral range, and only for high moon's illumination conditions in the case of near infrared channels. Perceptible water vapour (PWV) validation showed a good agreement between CE318-T and Global Navigation Satellite System (GNSS) PWV values for all illumination conditions, within the expected precision for sun photometry. Finally, two case studies have been included to highlight the ability of the new CE318-T to capture the diurnal cycle of aerosols and water vapour as well as short-term atmospheric variations, critical for climate studies.

Barreto, Africa↗
NASA NTRSDate not supplied

Ultra-Low Noise Laser and Optical Frequency Comb-Based Timing System for the Event Horizon Explorer (EHE) Mission

Very Long Baseline Interferometry (VLBI) requires a highly stable time reference and synchronization system to maintain coherence between observations recorded independently at multiple stations – for example, the array of ground-based radio telescopes that comprise the Event Horizon Telescope (EHT). Instability in the time reference or inaccuracy in the synchronization impacts phase conference and reduces the performance of the VLBI system. For space-based telescopes like the Event Horizon Explorer (EHE) mission, the precision timing system must meet the size, weight, and power (SWaP) constraints of the spacecraft platform and be capable of operating in the space environment. The EHE mission is a proposed extension to the EHT that can improve the angular and time resolution of observations by enabling longer interferometric baselines than are possible on the Earth and sampling a wide range of Fourier spatial frequencies throughout the orbital motion. In this effort, we demonstrate the performance of one option for the EHE mission: the use of a space-qualified, ultra-low noise laser developed as part of the Laser Interferometer Space Antenna (LISA) mission as the timing reference, and an optical frequency comb to transfer the stability of this laser to the microwave regime for instrumentation use. We describe the implementation of the microwave down-conversion, in which the LISA cavity-stabilized laser is locked to a self-referenced optical frequency comb to divide the optical frequency down to 100 MHz. The phase noise of the 100 MHz signal is measured and validated using a phase noise analyzer that is referenced to a separate laboratory ultra-stable laser system. We present the results of this experiment, which demonstrates that the performance of this system meets the EHE requirement: a relative frequency error of 1e-14 from 1 to 30 seconds.

Hannah Tomio↗
NASA NTRS1999

Search for Pulsations from a Nearby Millisecond Pulsar and Wasilewski 49: Mirror for a Hidden Seyfert 1 Nucleus

Five studies are reported in this final report. The recently discovered 5.3 ms pulsar J1012+5307 at a distance of 520 pc is in an area of the sky which is particularly deficient in absorbing gas. The column density along the line of sight is less than 7.5 x 10(exp 19) CM(exp -2) which facilitates soft X-ray observations. Halpern reported a possible ROSAT Position Sensitive Proportional Counter (PSPC) detection of the pulsar in a serendipitous, off-axis observation. We have now confirmed the X-ray emission of PSR J1012+,5307 in a 23 ksec observation with the ROSAT High Resolution Imager (HRI). A point source is detected within 3" of the radio position. Its count rate of 1.6 +/- 0.3 x 10(exp -3) s(exp -1) corresponds to an unabsorbed 0. 1-2.4 keV flux of 6.4 x 10(exp -14) ergs cm(exp -2) s(exp -1), similar to that reported previously. This counts-to-flux conversion is valid for N(sub H) = 5 x 10(exp 19) cm(exp -2), and either a power-law spectrum of photon index 2.5 or a blackbody of kT = 0.1 keV. The implied X-ray luminosity of 2.0 x 10(exp 30) ergs s(exp -1) is 5 X 10(exp -4) of the pulsar's spin-down power dot-E, and similar to that of the nearest millisecond pulsar J0437-4715, which is nearly a twin of J1012+5307 in P dot-E. We subjected the 37 photons (and 13 background counts) within the source region to a pulsar search, but no evidence for pulsation was found. The pulsar apparently emits over a large fraction of its rotation cycle, and the absence of sharp modulation can be taken as evidence for surface thermal emission, as is favored for PSR J0437-4715, rather than magnetospheric X-ray emission which is apparent in the sharp pulses of the much more energetic millisecond pulsar B1821-24. A further test of this interpretation will be made with a longer ROSAT observation, which will increase the number of photons collected by a factor of 5, and permit a more sensitive examination of the light curve for modulation due to emission from heated polar caps. If found, such modulation will be further evidence that surface reheating by the impact of particles accelerated along open fiel;d lines operates in these approximately 10(exp -9) year old pulsars. In a second study, a new AM Her star serendipitously in a 25 day observation was detected with the EUVE satellite. A coherent period of 85.82 min is present in the EUVE Deep Survey imager light curve of this source. A spectroscopic optical identification is made with a 19th magnitude blue star that has H and He emission lines, and broad cyclotron humps typical of a magnetic cataclysmic variable. A lower limit to the polar magnetic field of 50 MG is estimated from the spacing of the cyclotron harmonics. EUVE J0425.6-5714 is also detected in archival ROSAT HRI observations spanning two months, and its stable and highly structured light curve permits us to fit a coherent ephemeris linking the ROSAT and EUVE data over a 1.3 yr gap. The derived period is 85.82107 +/- 0.00020 min, and the ephemeris should be accurate to 0.1 cycles until the year 2005. A narrow but partial X-ray eclipse suggests that this object belongs to the group of AM Her stars whose viewing geometry is such that the accretion stream periodically occults the soft X-ray emitting accretion spot on the surface of the white dwarf. A non-detection of hard X-rays from ASCA observations that are contemporaneous with the ROSAT HRI shows that the soft X-rays must dominate by at least an order of magnitude, which is consistent with a known trend among AM Her stars with large magnetic field. This object should not be confused with the Seyfert galaxy IH 0419-577 (= LB 1727), another X-ray/EUV source which lies only 3'95 away, and was the principal target of these monitoring observations. In a third report; the identity of the persistent high-energy (> 100 MeV) gamma-ray sources in the Galaxy, still largely a mystery is investigated. The second installment of the EGRET (2EG) lists a total of 128 sources, of which 51 are likely or possibly identified with AGNs, five with rotation-powered pulsars, and one is the LMC. There are 71 unidentified sources, of which 33, or almost half, lie in the narrow band of absolute value of b < 10 degrees along the Galactic plane. This excess of low-latitude sources must, therefore, constitute a Galactic population that is either similar to the already identified gamma-ray pulsars, or an entirely new class of gamma-ray emitters associated with the disk population. We are continuing our program, begun in A06, that is aimed at intermediate-latitude sources, arguing that X-ray detection of them is the most plausible method of identifying the Galactic population. The sources at high latitude must statistically be mostly AGNS, and are more straightforwardly identified through radio and optical means. For the fourth report, identification of the EGRET Source the identity of the persistent, high-energy gamma-ray sources in the Galactic plane is a mystery. The most likely scenario is a population of middle-aged pulsars, many of which could be radio quiet like Geminga. We have an ongoing program of ROSAT, VLA, and optical observations of selected EGRET error circles at intermediate Galactic latitude. For one of these fields, at (l, b) = (106 degrees, +3), our complete census of X-ray and radio sources, reveals a remarkable association between a radio shell with unique properties, and a complete census of X-ray and radio sources reveals a remarkable association between a radio shell with unique properties, and a compact X-ray source. Further observations are needed to determine whether or not this source has a hard X-ray spectrum like that of other gamma ray pulsars and, ideally, to find its pulsations. A fifth report reviews, Wasilewski 49, an interacting pair of Seyfert galaxies at z = 0.063, one of which contains a hidden Seyfert I nucleus as evidenced by broad wings on its Balmer lines.

Halpern, Jules P.↗
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗