Search NASA⌕ Search

SEARCH · Search NASA

Results for “errors total quality”

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.

102 records · Page 6

Remote Sensing of the Absorption Coefficients and Chlorophyll a Concentration in the U.S. Southern Middle Atlantic Bight from SeaWiFS and MODIS-Aqua

At present, satellite remote sensing of coastal water quality and constituent concentration is subject to large errors as compared to the capability of satellite sensors in oceanic waters. In this study, field measurements collected on a series of cruises within U.S. southern Middle Atlantic Bight (SMAB) were applied to improve retrievals of satellite ocean color products in order to examine the factors that regulate the bio-optical properties within the continental shelf waters of the SMAB. The first objective was to develop improvements in satellite retrievals of absorption coefficients of phytoplankton (a(sub ph)), colored dissolved organic matter (CDOM) (a(sub g)), non-pigmented particles (a(sub d)), and non-pigmented particles plus CDOM (a(sub dg)), and chlorophyll a concentration ([Chl_a]). Several algorithms were compared to derive constituent absorption coefficients from remote sensing reflectance (R(sub rs)) ratios. The validation match-ups showed that the mean absolute percent differences (MAPD) were typically less than 35%, although higher errors were found for a(sub d) retrievals. Seasonal and spatial variability of satellite-derived absorption coefficients and [Chl_a] was apparent and consistent with field data. CDOM is a major contributor to the bio-optical properties of the SMAB, accounting for 35-70% of total light absorption by particles plus CDOM at 443 nm, as compared to 30-45% for phytoplankton and 0-20% for non-pigmented particles. The overestimation of [Chl_a] from the operational satellite algorithms may be attributed to the strong CDOM absorption in this region. River discharge is important in controlling the bio-optical environment, but cannot explain all of the regional and seasonal variability of biogeochemical constituents in the SMAB.

Pan, Xiaoju↗

Revised estimates of NO 2 reductions during the COVID-19 lockdowns using updated TROPOMI NO 2 retrievals and model simulations

The TROPOspheric Monitoring Instrument (TROPOMI) observed unprecedented declines in NO 2 vertical column densities (VCD) over the world's most densely populated cities during the 2020 COVID-19 lockdowns. These favorable changes in NO 2 air quality were correlated with sharp reductions in traffic volume and economic activity during the lockdowns. In this comprehensive global study, we provide revised estimates of the declines in anthropogenic emissions for 36 megacities using a novel methodology for disentangling the anthropogenic emissions from the meteorological transport and natural variability. We further quantify the uncertainty associated with changes in the a priori profile shape information during the lockdowns due to reduced emissions. Satellite NO 2 retrieval techniques calculate an air mass factor that requires a priori NO 2 profile shape information representative of the local atmosphere. This information, which is typically obtained from a chemical transport model (CTM), was not available for the early studies. This study also accounts for the satellite sampling errors resulting from the selective sampling of non-cloudy scenes during the study period. For our analysis, we used CTM simulations that were generated with and without COVID-impacted emissions. We perform retrievals of tropospheric NO 2 columns with the NASA NO 2 algorithm, and then use observed and simulated data to disentangle the meteorological transport from the contribution due anthropogenic emissions. We found that the meteorological transport was most significant source of variability ranging between −35% and 22% of the change total tropospheric VCD. We also find that not accounting for changes in the a priori NO 2 profile shape information during the lockdowns resulted in systematic retrieval errors that were up to 12% of the estimated decline, and the elimination of cloud contaminated scenes resulted in sampling errors that in general ranged between varied ±15%.

NO2↗

Science Benefits of Onboard Spacecraft Navigation

Primitive bodies (asteroids and comets), which have remained relatively unaltered since their formation, are important targets for scientific missions that seek to understand the evolution of the solar system. Often the first step is to fly by these bodies with robotic spacecraft. The key to maximizing data returns from these flybys is to determine the spacecraft trajectory relative to the target body-in short, navigate the spacecraft- with sufficient accuracy so that the target is guaranteed to be in the instruments' field of view. The most powerful navigation data in these scenarios are images taken by the spacecraft of the target against a known star field (onboard astrometry). Traditionally, the relative trajectory of the spacecraft must be estimated hours to days in advance using images collected by the spacecraft. This is because of (1)!the long round-trip light times between the spacecraft and the Earth and (2)!the time needed to downlink and process navigation data on the ground, make decisions based on the result, and build and uplink instrument pointing sequences from the results. The light time and processing time compromise navigation accuracy considerably, because there is not enough time to use more accurate data collected closer to the target-such data are more accurate because the angular capability of the onboard astrometry is essentially constant as the distance to the target decreases, resulting in better "plane-of- sky" knowledge of the target. Excellent examples of these timing limitations are high-speed comet encounters. Comets are difficult to observe up close; their orbits often limit scientists to brief, rapid flybys, and their coma further restricts viewers from seeing the nucleus in any detail, unless they can view the nucleus at close range. Comet nuclei details are typically discernable for much shorter durations than the roundtrip light time to Earth, so robotic spacecraft must be able to perform onboard navigation. This onboard navigation can be accomplished through a self- contained system that by eliminating light time restrictions dramatically improves the relative trajectory knowledge and control and subsequently increases the amount of quality data collected. Flybys are one-time events, so the system's underlying algorithms and software must be extremely robust. The autonomous software must also be able to cope with the unknown size, shape, and orientation of the previously unseen comet nucleus. Furthermore, algorithms must be reliable in the presence of imperfections and/or damage to onboard cameras accrued after many years of deep-space operations. The AutoNav operational flight software packages, developed by scientists at the Jet Propulsion Laboratory (JPL) under contract with NASA, meet all these requirements. They have been directly responsible for the successful encounters on all of NASA's close-up comet-imaging missions (see Figure !1). AutoNav is the only system to date that has autonomously tracked comet nuclei during encounters and performed autonomous interplanetary navigation. AutoNav has enabled five cometary flyby missions (Table!1) residing on four NASA spacecraft provided by three different spacecraft builders. Using this software, missions were able to process a combined total of nearly 1000 images previously unseen by humans. By eliminating the need to navigate spacecraft from Earth, the accuracy gained by AutoNav during flybys compared to ground-based navigation is about 1!order of magnitude in targeting and 2!orders of magnitude in time of flight. These benefits ensure that pointing errors do not compromise data gathered during flybys. In addition, these benefits can be applied to flybys of other solar system objects, flybys at much slower relative velocities, mosaic imaging campaigns, and other proximity activities (e.g., orbiting, hovering, and descent/ascent).

Autonomy↗

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗

Characterizing DebriSat Fragments: So Many Fragments, So Much Data, and So Little Time

To improve prediction accuracy, the DebriSat project was conceived by NASA and DoD to update existing standard break-up models. Updating standard break-up models require detailed fragment characteristics such as physical size, material properties, bulk density, and ballistic coefficient. For the DebriSat project, a representative modern LEO spacecraft was developed and subjected to a laboratory hypervelocity impact test and all generated fragments with at least one dimension greater than 2 mm are collected, characterized and archived. Since the beginning of the characterization phase of the DebriSat project, over 130,000 fragments have been collected and approximately 250,000 fragments are expected to be collected in total, a three-fold increase over the 85,000 fragments predicted by the current break-up model. The challenge throughout the project has been to ensure the integrity and accuracy of the characteristics of each fragment. To this end, the post hypervelocity-impact test activities, which include fragment collection, extraction, and characterization, have been designed to minimize handling of the fragments. The procedures for fragment collection, extraction, and characterization were painstakingly designed and implemented to maintain the post-impact state of the fragments, thus ensuring the integrity and accuracy of the characterization data. Each process is designed to expedite the accumulation of data, however, the need for speed is restrained by the need to protect the fragments. Methods to expedite the process such as parallel processing have been explored and implemented while continuing to maintain the highest integrity and value of the data. To minimize fragment handling, automated systems have been developed and implemented. Errors due to human inputs are also minimized by the use of these automated systems. This paper discusses the processes and challenges involved in the collection, extraction, and characterization of the fragments as well as the time required to complete the processes. The objective is to provide the orbital debris community an understanding of the scale of the effort required to generate and archive high quality data and metadata for each debris fragment 2 mm or larger generated by the DebriSat project.

Shiotani, B.↗

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗

Improving Global Analysis and Short-Range Forecast Using Rainfall and Moisture Observations Derived from TRMM and SSM/I Passive Microwave Instruments

The Global Precipitation Mission, a satellite project under consideration as a follow-on to the Tropical Rainfall Measuring Mission (TRMM) by the National Aeronautics and Space Agency (NASA) in the United States, the National Space Development Agency (NASDA) in Japan, and other international partners, comprises an improved TRMM-like satellite and a constellation of 8 satellites carrying passive microwave radiometers to provide global rainfall measurements at 3-hour intervals. The success of this concept relies on the merits of rainfall estimates derived from passive microwave radiometers. This article offers a proof-of-concept demonstration of the benefits of using, rainfall and total precipitable water (TPW) information derived from such instruments in global data assimilation with observations from the TRMM Microwave Imager (TMI) and 2 Special Sensor Microwave/Imager (SSM/I) instruments. Global analyses that optimally combine observations from diverse sources with physical models of atmospheric and land processes can provide a comprehensive description of the climate systems. Currently, such data analyses contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. We show that assimilating the 6-h averaged TMI and SSM/I surface rainrate and TPW retrievals improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the upper tropospheric moisture in the analysis produced by the Goddard Earth Observing System (GEOS) Data Assimilation System, as verified against radiation measurements by the Clouds and the Earth's Radiant Energy System (CERES) instrument and brightness temperature observations by the TIROS Operational Vertical Sounder (TOVS) instruments. Typically, rainfall assimilation improves clouds and radiation in areas of active convection, as well as the latent heating and large-scale motions in the tropics, while TPW assimilation leads to reduced moisture biases and improved radiative fluxes in clear-sky regions. Ensemble forecasts initialized with analyses that incorporate TMI and SSM/I rainfall and TPW data also yield better short-range predictions of geopotential heights, winds, and precipitation in the tropics. This study offers a compelling illustration of the potential of using rainfall and TPW information derived from passive microwave instruments to significantly improve the quality of 4-dimensional global datasets for climate analysis and weather forecasting applications.

Hou, Arthur Y.↗

Convective Weather Forecast Accuracy Analysis at Center and Sector Levels

This paper presents a detailed convective forecast accuracy analysis at center and sector levels. The study is aimed to provide more meaningful forecast verification measures to aviation community, as well as to obtain useful information leading to the improvements in the weather translation capacity models. In general, the vast majority of forecast verification efforts over past decades have been on the calculation of traditional standard verification measure scores over forecast and observation data analyses onto grids. These verification measures based on the binary classification have been applied in quality assurance of weather forecast products at the national level for many years. Our research focuses on the forecast at the center and sector levels. We calculate the standard forecast verification measure scores for en-route air traffic centers and sectors first, followed by conducting the forecast validation analysis and related verification measures for weather intensities and locations at centers and sectors levels. An approach to improve the prediction of sector weather coverage by multiple sector forecasts is then developed. The weather severe intensity assessment was carried out by using the correlations between forecast and actual weather observation airspace coverage. The weather forecast accuracy on horizontal location was assessed by examining the forecast errors. The improvement in prediction of weather coverage was determined by the correlation between actual sector weather coverage and prediction. observed and forecasted Convective Weather Avoidance Model (CWAM) data collected from June to September in 2007. CWAM zero-minute forecast data with aircraft avoidance probability of 60% and 80% are used as the actual weather observation. All forecast measurements are based on 30-minute, 60- minute, 90-minute, and 120-minute forecasts with the same avoidance probabilities. The forecast accuracy analysis for times under one-hour showed that the errors in intensity and location for center forecast are relatively low. For example, 1-hour forecast intensity and horizontal location errors for ZDC center were about 0.12 and 0.13. However, the correlation between sector 1-hour forecast and actual weather coverage was weak, for sector ZDC32, about 32% of the total variation of observation weather intensity was unexplained by forecast; the sector horizontal location error was about 0.10. The paper also introduces an approach to estimate the sector three-dimensional actual weather coverage by using multiple sector forecasts, which turned out to produce better predictions. Using Multiple Linear Regression (MLR) model for this approach, the correlations between actual observation and the multiple sector forecast model prediction improved by several percents at 95% confidence level in comparison with single sector forecast.

Wang, Yao↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

Development of Jet Noise Power Spectral Laws

High-quality jet noise spectral data measured at the Aero-Acoustic Propulsion Laboratory (AAPL) at NASA Glenn is used to develop jet noise scaling laws. A FORTRAN algorithm was written that provides detailed spectral prediction of component jet noise at user-specified conditions. The model generates quick estimates of the jet mixing noise and the broadband shock-associated noise (BBSN) in single-stream, axis-symmetric jets within a wide range of nozzle operating conditions. Shock noise is emitted when supersonic jets exit a nozzle at imperfectly expanded conditions. A successful scaling of the BBSN allows for this noise component to be predicted in both convergent and convergent-divergent nozzles. Configurations considered in this study consisted of convergent and convergent- divergent nozzles. Velocity exponents for the jet mixing noise were evaluated as a function of observer angle and jet temperature. Similar intensity laws were developed for the broadband shock-associated noise in supersonic jets. A computer program called sJet was developed that provides a quick estimate of component noise in single-stream jets at a wide range of operating conditions. A number of features have been incorporated into the data bank and subsequent scaling in order to improve jet noise predictions. Measurements have been converted to a lossless format. Set points have been carefully selected to minimize the instability-related noise at small aft angles. Regression parameters have been scrutinized for error bounds at each angle. Screech-related amplification noise has been kept to a minimum to ensure that the velocity exponents for the jet mixing noise remain free of amplifications. A shock-noise-intensity scaling has been developed independent of the nozzle design point. The computer program provides detailed narrow-band spectral predictions for component noise (mixing noise and shock associated noise), as well as the total noise. Although the methodology is confined to single streams, efforts are underway to generate a data bank and algorithm applicable to dual-stream jets. Shock-associated noise in high-powered jets such as military aircraft can benefit from these predictions.

Khavaran, Abbas↗

GRB 091127: The Cooling Break Race on Magnetic Fuel

Using high-quality, broad-band afterglow data for GRB 091127, we investigate the validity of the synchrotron fireball model for gamma-ray bursts, and infer physical parameters of the ultra-relativistic outflow. Methods. We used multi-wavelength (NIR to X-ray) follow-up observations obtained with GROND simultaneously in the g' r' t' i' z' JH filters and the XRT onboard the Swift satellite in the 0.3 to 10 keY energy range. The resulting afterglow light curve is of excellent accuracy with relative photometric errors as low as 1 %, and the spectral energy distribution (SED) is well-sampled over 5 decades in energy. These data present one of the most comprehensive observing campaigns for a single GRB afterglow and allow us to test several proposed emission models and outflow characteristics in unprecedented detail. Results. Both the multi-color light curve and the broad-band SED of the afterglow of GRB 091127 show evidence of a cooling break moving from high to lower energies. The early light curve is well described by a broken power-law, where the initial decay in the optical/NlR wavelength range is considerably flatter than at X-rays. Detailed fitting of the time-resolved SED shows that the break is very smooth with a sharpness index of 2.2 +/- 0.2, and evolves towards lower frequencies as a power-law with index -1.23 +/- 0.06. These are the first accurate and contemporaneous measurements of both the sharpness of the spectral break and its time evolution. Conclusions. The measured evolution of the cooling break (V(sub c) varies as t(sup -1.2) is not consistent with the predictions of the standard model, wherein V(sub c) varies as t(sup -05) is expected. A possible explanation for the observed behavior is a time dependence of the microphysical parameters, in particular the fraction of the total energy in the magnetic field epsilon(sub Beta). This conclusion provides further evidence that the standard fireball model is too simplistic, and time-dependent micro-physical parameters may be required to model the growing number of well-sampled afterglow light curves.

Filgas, R.↗

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↗