SEARCH · Engineering Papers

Results for “surface mass variations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A Scalable Compact Additively Manufactured Molten Salt to Supercritical Carbon Dioxide Heat Exchanger for Solar Thermal Application

Design of an additively manufactured molten salt (MS) to supercritical carbon dioxide (sCO 2 ) primary heat exchanger (PHE) for solar thermal power generation is presented. The PHE is designed to handle temperatures up to 720 °C on the MS side and an internal pressure of 200 bar on the sCO 2 side. In the core, MS flows through a three-dimensional periodic lattice network, while sCO 2 flows within pin arrays. The design includes integrated sCO 2 headers located within the MS flow, allowing for a counterflow design of the PHE. The sCO 2 headers are configured to enable uniform flow distribution into each sCO 2 plate while withstanding an internal pressure of 200 bar and minimizing obstruction to the flow of MS around it. The structural integrity of the design is verified on additively manufactured (AM) 316 stainless steel sub-scale specimens. An experimentally validated, correlation-based sectional PHE core thermofluidic model is developed to study the impact of flow and geometrical parameters on the PHE performance, with varied parameters including the mass flowrate, surface roughness, and PHE dimensions. A process-based cost model is used to determine the impact of parameter variation on build cost. The model results show that a heat exchanger with a power density of 18.6 MW/m 3 (including sCO 2 header volume) and effectiveness of 0.88 can be achieved at a heat capacity rate ratio of 0.8. As a result, the impact of design and AM machine parameters on the cost of the PHE are assessed.

14 SOLAR ENERGY↗

NASA NTRS1976

Initial results from the High Altitude Observatory white light coronagraph on Skylab - A progress report

Frequent periodic observations by the white-light coronagraph allow an examination of coronal variations over a broad range of temporal scales. Examples of the slowest and most rapid variations are presented. An example of extremely slow coronal variations is the gradual evolution - to a large equatorial streamer - in association with a marked decrease in solar activity, as the total magnetic flux in one hemisphere decreased. Another example is given of a long-lived quasi-stable coronal streamer, apparently associated with a stable filament channel; comparison of this streamer with coronal potential-magnetic-field computations show little correlation. Some results on coronal transients - the most rapid variations observed - are summarized. Characteristic masses and energies involved in mass-ejection transients, their temporal and spatial distributions, their associations with surface phenomena and possible interplanetary signatures, and their role in coronal evolution are briefly noted.

Macqueen, R. M.↗

NASA NTRS1980

Geoid anomalies in the vicinity of subduction zones

The regional geoid of the southwest Pacific is matched reasonably well by results from a model of the upper mantle density structure (including slabs) associated with subduction zones of the region. Estimates of the geoid are obtained from Geos-3 and Seasat radar altimeter data. These data are very well suited to the task of detecting intermediate wavelength (600-4000 km) geopotential variations. Actually, subducting slabs can be expected to produce primarily intermediate and longer wavelength variations. Gravimetric profiles across trench/island arc complexes resolve primarily short wavelengths. The model represents subducting slabs as thin surfaces of anomalous mass per unit area. These surfaces are positioned using published seismicity results which detail the configuration of the Benioff zones. Crustal effects are ignored. Effects due to the contrast between the young thermal lithosphere of the behind-arc regions (marginal basins) and the older lithosphere seaward of the trench are modelled. Results indicate that the New Hebrides slab possesses an average areal density anomaly of about 300,000 gm/sq cm. This is about three times that which is estimated for the Tonga-Kermadec slab. Additional modelling suggests that slabs worldwide may be an important source of large, long wavelength gravity highs; i.e., they may contribute substantially to geopotential power of harmonic degree as low as three or four up to twenty or more.

Mcadoo, D. C.↗

NASA NTRS2014

The Global Albedo of the Moon at 1064 nm from LOLA

The Lunar Orbiter Laser Altimeter (LOLA) measures the backscattered energy of the returning altimetric laser pulse at its wavelength of 1064 nm, and these data are used to map the reflectivity of the Moon at zero-phase angle with a photometrically uniform data set. Global maps have been produced at 4 pixels per degree (about 8 kilometers at the equator) and 2 kilometers resolution within 20 deg latitude of each pole. The zero-phase geometry is insensitive to lunar topography, so these data enable characterization of subtle variations in lunar albedo, even at high latitudes where such measurements are not possible with the Sun as the illumination source. The geometric albedo of the Moon at 1064 nm was estimated from these data with absolute calibration derived from the Kaguya Multiband Imager and extrapolated to visual wavelengths. The LOLA estimates are within 2 sigma of historical measurements of geometric albedo. No consistent latitude-dependent variations in reflectance are observed, suggesting that solar wind does not dominate space weathering processes that modify lunar reflectance. The average normal albedo of the Moon is found to be much higher than that of Mercury consistent with prior measurements, but the normal albedo of the lunar maria is similar to that of Mercury suggesting a similar abundance of space weathering products. Regions within permanent shadow in the polar regions are found to be more reflective than polar surfaces that are sometimes illuminated. Limiting analysis to data with slopes less than 10 deg eliminates variations in reflectance due to mass wasting and shows a similar increased reflectivity within permanent polar shadow. Steep slopes within permanent shadow are also more reflective than similar slopes that experience at least some illumination. Water frost and a reduction in effectiveness of space weathering are offered as possible explanations for the increased reflectivity of permanent shadow; porosity is largely ruled out as the sole explanation. The south polar crater Shackleton is found to be among the most reflective craters in its size range globally but is not the most reflective, so mass wasting cannot be ruled out as a cause for the crater's anomalous reflectance. Models of the abundance of ice needed to account for the reflectance anomaly range from 3 to 14% by weight or area depending on assumptions regarding the effects of porosity on reflectance and whether ice is present as patches or is well mixed in the regolith. If differences in nanophase iron abundances are responsible for the anomaly, the permanently shadowed regions have between 50 and 80% the abundance of nanophase iron in mature lunar soil.

PSR↗

NASA NTRS2000

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.↗

DOE OSTI2022

Radial velocity variability and the evolution of hot subdwarf stars

Hot subdwarf stars represent a late and peculiar stage in the evolution of low-mass stars, since they are likely formed by close binary interactions. In this work, we perform a radial velocity (RV) variability study of a sample of 646 hot subdwarfs with multi-epoch radial velocities based on spectra from Sloan Digital Sky Survey (SDSS) and Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). The atmospheric parameters and RVs were taken from the literature. For stars with archival spectra but without literature values, we determined the parameters by fitting model atmospheres. In addition, we redetermined the atmospheric parameters and RVs for all the He-enriched sdO/Bs. This broad sample allowed us to study RV-variability as a function of the location in theT eff - log g- and T eff - log n (He)/n (H) diagrams in a statistically significant way. We used the fraction of RV-variable stars and the distribution of the maximum RV variations ΔRV max as diagnostics. Both indicators turned out to be quite inhomogeneous across the studied parameter ranges. A striking feature is the completely dissimilar behaviour of He-poor and He-rich hot subdwarfs. While the former have a high fraction of close binaries, almost no significant RV variations could be detected for the latter. This has led us to the conclusion that there is likely no evolutionary connection between these subtypes. On the other hand, intermediate He-rich- and extreme He-rich sdOB/Os are more likely to be related. Furthermore, we conclude that the vast majority of this population is formed via one or several binary merger channels. Hot subdwarfs with temperatures cooler than ~24 000 K tend to show fewer and smaller RV-variations. These objects might constitute a new subpopulation of binaries with longer periods and late-type or compact companions. The RV-variability properties of the extreme horizontal branch (EHB) and corresponding post-EHB populations of the He-poor hot subdwarfs match and confirm the predicted evolutionary connection between them. Stars found below the canonical EHB at somewhat higher surface gravities show large RV variations and a high RV variability fraction. These properties are consistent with most of them being low-mass EHB stars or progenitors of low-mass helium white dwarfs in close binaries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Tailoring Cu-Zr gradient nanoglass structures: Influence of nanoparticle size and cooling rates on glass-glass interfaces

The study of gradient nanoglasses (GNGs) has gained attention due to their unique mechanical properties and potential applications in advanced materials. This study employs molecular dynamics simulations to synthesize a GNG using Cu-Zr metallic glass nanoparticles (NPs) sized from 3 to 15 nm. The NPs were produced by melting and quenching metallic clusters at a relatively slow quench rate of 10 9 K/s. The synthesis of GNG is elucidated along with the characterization of its heterogeneous metallic glass nanostructure. A seamless GNG structure is formed through cold compression of Cu 64 Zr 36 amorphous NPs of varying sizes. The influence of NP size on the GNG structure is investigated, utilizing deeply relaxed NPs, which exhibit a characteristic Cu segregation pattern on their surfaces. The results highlight an increase in structural heterogeneity due to heterogeneous mass transport and the development of local composition and density variations caused by Cu segregation at glass-glass interfaces (GGIs). A reduction in NP size is correlated with decreased Cu atomic displacements and local density at GGIs, suggesting that larger NPs may produce stronger GGIs. This research presents a novel methodology for synthesizing heterogeneous metallic glasses, demonstrating the capacity to control and customize nanostructure heterogeneity through the manipulation of NP sizes and cooling rates. Furthermore, these findings enhance our understanding of structural evolution during nanoglass synthesis and lay the foundation for further exploration in nanomaterial synthesis and characterization.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Integration of the fundamental knowledge on solvent-packing interactions into the multiscale framework for column scale design and optimization

The interfacial area, also known as the effective mass transfer area, is a key factor for determining the mass transfer for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Most of the widely used empirical and semi-empirical models for interfacial area were derived indirectly through absorption mass transfer with simplifications based on fast chemical kinetics. This report presents the comprehensive unique multiscale approach to develop a surrogate model for effective mass transfer area in structured packed columns that accounts local hydrodynamics as well as variation in physical properties, and changes in solid surface characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1991

Frequency effects on the stability of a journal bearing for periodic loading

The stability of a journal bearing is numerically predicted when a unidirectional periodic external load is applied. The analysis is performed using a cavitation algorithm, which mimics the Jakobsson-Floberg and Olsson (JFO) theory by accounting for the mass balance through the complete bearing. Hence, the history of the film is taken into consideration. The loading pattern is taken to be sinusoidal and the frequency of the load cycle is varied. The results are compared with the predictions using Reynolds boundary conditions for both film rupture and reformation. With such comparisons, the need for accurately predicting the cavitation regions for complex loading patterns is clearly demonstrated. For a particular frequency of loading, the effects of mass, amplitude of load variation and frequency of journal speed are also investigated. The journal trajectories, transient variations in fluid film forces, net surface velocity and minimum film thickness, and pressure profiles are also presented.

Vijayaraghavan, D.↗

NASA NTRS1992

Microcraters on Mars: Evidence of past climatic variations

On Mars, we do not expect to see the full spectrum of crater sizes that we see on the Moon, because the Martian atmosphere decelerates smaller incoming meteorites. Mars climate models predict thirty-fold variations in Mars' atmospheric pressure due to large amplitude quasi-periodic variations in Mars' obliquity. More recent studies have shown that uncertainties in Mars' moment of inertia may have resulted in an underestimation of this range. In this study, we find that potential variations in the mass of the Martian atmosphere should have dramatic effects on the production rates of centimeter sized craters on the surface of Mars, and suggest that observations of the densities of Martian microcraters would provide an important validation of the astronomical theory for Martian climate change.

Vasavada, A. R.↗

NASA NTRS1986

Development of a variational SEASAT data analysis technique

Oceans are data-sparse areas in terms of conventional weather observations. The surface pressure field obtained solely by analyzing the conventional weather data is not expected to possess high accuracy. On the other hand, in entering asynoptic data such as satellite-derived temperature soundings into an atmospheric prediction system, an improved surface analysis is crucial for obtaining more accurate weather predictions because the mass distribution of the entire atmosphere will be better represented in the system as a result of the more accurate surface pressure field. In order to obtain improved surface pressure analyses over the oceans, a variational adjustment technique was developed to help blend the densely distributed surface wind data derived from the SEASAT-A radar observations into the sparsely distributed conventional pressure data. A simple marine boundary layer scheme employed in the adjustment technique was discussed. In addition, a few aspects of the current technique were determined by numerical experiments.

Sasaki, Y. K.↗

NASA NTRS2007

Evidence for Interannual to Decadal Variations in Hadley and Walker Circulations and Links to Water and Energy Fluxes

Mass and energy transports associated with the Hadley and Walker circulations are important components of the earth s climate system and are strongly linked to hydrologic processes. Interannual to decadal variation in these flows likely signify a combination of natural climate noise as well as a response to anthropgenic forcing. There remains considerable uncertainty in quantifying variations in these flows. Evidence in the surface pressure record supports a weakening of the Walker circulation over the Pacific in recent decades. Conversely the NCEP / NCAR and ERA 40 reanalyses indicate that the Hadley circulation has increased in strength over the last two decades, though these analyses depict significantly different mass circulation changes. Interestingly, the NCEP - II / DOE reanalysis contains essentially no Hadley circulation changes. Most climate model integrations anticipate a weakening of both tropical circulations associated with stronger static stability. Clearly there is much uncertainty not only with the mass transports, but also how they are linked to water and energy balance of the planet through variations in turbulent heat and radiative fluxes and horizontal exports / imports of energy. Here we examine heat and water budget variations from a number of reanalysis products and focus on the linear and nonlinear response of ENSO warm and cold events as opportunities to study budget variations over the past 15-20 years. Our analysis addresses such questions as To what extent do Hadley and Walker Cell variations compensate each other on mass and energy transport? Do static stability adjustments appear to constrain fractional precipitation response vs. fractional water vapor response? We appeal to constraints offered by GPCP precipitation, SSWI ocean evaporation estimates, and ISCCP-FD radiative fluxes, and other satellite data sets to interpret and confirm reanalysis-based diagnostics. Using our findings we also attempt to place in context the recent findings that tropical ocean evaporation increased by order 5% or more during the 1990s, reconciling this with GPCP precipitation variations.

Robertson, Franklin↗

NASA NTRS1995

Development of a Small Area Sniffer

The aim of this project is to develop and implement a sniffer that is capable of measuring the mass flow rate of air through a small area of pinholes whose diameters are on the magnitude of thousandths of an inch. The sniffer is used to scan a strip of a leading edge panel, which is being used in a hybrid laminar flow control experiment, in order to survey the variations in the amount of air that passes through the porous surface at different locations. Spanwise scans are taken at different chord locations by increasing the pressure in a control volume that is connected to the sniffer head, and recording the drop in pressure as the air is allowed to flow through the tiny holes. This information is used to obtain the mass flow through the structure. More importantly, the deviations from the mean flow rate are found and used to determine whether there are any significant variations in the flow rate from one area to the next. The preliminary results show little deviation in the spanwise direction. These results are important when dealing with the location and amount of suction that will be applied to the leading edge in the active laminar flow control experiment.

Meade, Laurie A.↗

NASA NTRS1986

Systematics of the CHON and other light-element particle populations in Comet Halley

Based on chemical signatures measured by the PIA experiment during the Giotto flyby of comet Halley, particle classifications were designated. In addition to silicate-like grains and particles of mixed (cosmic) composition, there appear to be several light-element rich populations, including the CHON, (H,C), (H,C,O), and (H,C,N) particle types. These compositional classes are further distinguished by differences in mass distributions, a density indicator, and variations in relative abundance within the coma. These particle populations are evidence for chemical heterogeneity in the surface of the cometary nucleus. Particles found mainly in the inner coma may be volatile icy grains. Most of the N of the comet may be found in up to three different populations of grains; one or more of these may be responsible for the observation of cyanojets.

Clark, Benton↗

NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗

NASA NTRS1994

Photosynthetic fractionation of 13C and concentrations of dissolved CO2 in the central equatorial Pacific during the last 255,000 years

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C [is equivalent to epsilon(p)] in a central equatorial Pacific sediment core that spans the last approximately 255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic composition of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon(p), derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 [is equivalent to c(e)] based on a new, site-specific calibration of the relationship between epsilon(p) and c(e). The calibration was based on reassessment of existing epsilon(p) versus c(e) data, which support a physiologically based model in which epsilon(p) is inversely related to c(e). Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index U(K/37). Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon(p) and 1/c(e). These are discussed in detail and it is concluded that the observed record of epsilon(p) most probably reflects significant variations in delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from approximately 110 microatmospheres during glacial intervals (ocean > atmosphere) to approximately 60 microatmospheres during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p<0.01) inverse correlations of delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

NASA Discipline Exobiology↗

NASA NTRS1972

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗

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↗