Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “land surface temperature”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 973 records · Page 54

Classifying Microwave Radiometer Observations Over the Netherlands Into Dry, Shallow-, and Non-Shallow Precipitation Using A Random Forest Model

Spaceborne microwave radiometers represent an important component of the Global Precipitation Measurement (GPM) mission due to their frequent sampling of rain systems. Microwave radiometers measure microwave radiation (brightness temperatures, Tb), which can be converted into precipitation estimates with appropriate assumptions. However, detecting shallow precipitation systems using space-borne radiometers is challenging, especially over land, as their weak signals are hard to differentiate from those associated with dry conditions. This study uses a random forest model (RF) to classify microwave radiometer observations as dry, shallow, or nonshallow over the Netherlands - a region with varying surface conditions and frequent occurrence of shallow precipitation. The RF is trained on five years of data (2016-2020) and tested with two independent years (2015, 2021). The observations are classified using ground-based weather radar echo top heights. Various RF models are assessed, such as using only GPM’s Microwave Imager (GMI) Tb values as input features or including spatially aligned ERA-5 2-meter temperature and freezing level reanalysis and/or Dual Precipitation Radar (DPR) observations. Independent of the input features, the model performs best in summer and worst in winter. The model classifies observations from high-frequency channels (≥85 GHz) with lower Tb-values as non-shallow, higher values as dry, and those in between as shallow. Misclassified footprints exhibit radiometric characteristics corresponding to their assigned class. Case studies reveal dry observations misclassified as shallow are associated with lower Tb-values, likely resulting from the presence of ice particles in non-precipitating clouds. Shallow footprints misclassified as dry are likely related to the absence of ice particles.

Linda Bogerd↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

The Second Phase of the Global Land-Atmosphere Coupling Experiment (GLACE-2)

A major motivation for the study of the coupled land-atmosphere system is the idea that soil moisture anomalies may affect future meteorological variables through their effects on future surface energy and water budgets. If true, the accurate initialization of soil moisture in a subseasonal or seasonal forecast system may improve forecast skill, making the forecast products more valuable to society. The GLACE-2 project is examining, through a coordinated experiment using a wide variety of models, the degree to which subseasonal (out to two months) precipitation and air temperature forecasts improve through the realistic initialization of soil moisture. For the first time ever, a global consensus should emerge regarding the value of land initialization for forecasts, perhaps motivating national forecast centers to make full use of land moisture initialization in their operations

Koster, Randal D.↗

Understanding Volume Estimation Uncertainty of Lakes and Wetlands Using Satellites and Citizen Science

We studied variations in the volume of water stored in small lakes and wetlands using satellite remote sensing and lake water height data contributed by citizen scientists. A total of 94 water bodies across the globe were studied using satellite data in the optical and microwave wavelengths from Landsat 8, Sentinel-1, and Sentinel-2. The uncertainty in volume estimation as a function of geography and geophysical factors, such as cloud cover, precipitation, and water surface temperature, was studied. The key finding that emerged from this global study is that uncertainty is highest in regions with a distinct precipitation season, such as in the monsoon dominated South Asia or the Pacific Northwestern region of the USA. This uncertainty is further compounded when small lakes and wetlands are seasonal with alternating land use as a water body and agricultural land, such as the wetlands of Northeastern Bangladesh. On an average, 45% of studied lakes could be estimated of their volume change with a statistical significant uncertainty that is less than the expected volume in South Asia. In North America, this statistically significant uncertainty in volume estimation was found to be around 50% in lakes eastward of the 108th meridian with lowest uncertainty found in lakes along the East coast of the USA. The article provides a baseline for understanding the current state of the art in estimating volumetric change of lakes and wetlands using citizen science in anticipation of the recently launched Surface Water and Ocean Topography Mission.

Shahzaib Khan↗

Atmospheric Environments for Entry, Descent and Landing (EDL)

Scientific measurements of atmospheric properties have been made by a wide variety of planetary flyby missions, orbiters, and landers. Although landers can make in-situ observations of near-surface atmospheric conditions (and can collect atmospheric data during their entry phase), the vast majority of data on planetary atmospheres has been collected by remote sensing techniques from flyby and orbiter spacecraft (and to some extent by Earth-based remote sensing). Many of these remote sensing observations (made over a variety of spectral ranges), consist of vertical profiles of atmospheric temperature as a function of atmospheric pressure level. While these measurements are of great interest to atmospheric scientists and modelers of planetary atmospheres, the primary interest for engineers designing entry descent and landing (EDL) systems is information about atmospheric density as a function of geometric altitude. Fortunately, as described in in this paper, it is possible to use a combination of the gas-law relation and the hydrostatic balance relation to convert temperature-versus-pressure, scientific observations into density-versus-altitude data for use in engineering applications. The following section provides a brief introduction to atmospheric thermodynamics, as well as constituents, and winds for EDL. It also gives methodology for using atmospheric information to do "back-of-the-envelope" calculations of various EDL aeroheating parameters, including peak deceleration rate ("g-load"), peak convective heat rate. and total heat load on EDL spacecraft thermal protection systems. Brief information is also provided about atmospheric variations and perturbations for EDL guidance and control issues, and atmospheric issues for EDL parachute systems. Subsequent sections give details of the atmospheric environments for five destinations for possible EDL missions: Venus. Earth. Mars, Saturn, and Titan. Specific atmospheric information is provided for these destinations, and example results are presented for the "back-of-the-envelope" calculations mentioned above.

Justus, Carl G.↗

Review of NASA antiskid braking research

NASA antiskid braking system research programs are reviewed. These programs include experimental studies of four antiskid systems on the Langley Landing Loads Track, flights tests with a DC-9 airplane, and computer simulation studies. Results from these research efforts include identification of factors contributing to degraded antiskid performance under adverse weather conditions, tire tread temperature measurements during antiskid braking on dry runway surfaces, and an assessment of the accuracy of various brake pressure-torque computer models. This information should lead to the development of better antiskid systems in the future.

Tanner, J. A.↗

Design of a global soil moisture initialization procedure for the simple biosphere model

Global soil moisture and land-surface evapotranspiration fields are computed using an analysis scheme based on the Simple Biosphere (SiB) soil-vegetation-atmosphere interaction model. The scheme is driven with observed precipitation, and potential evapotranspiration, where the potential evapotranspiration is computed following the surface air temperature-potential evapotranspiration regression of Thomthwaite (1948). The observed surface air temperature is corrected to reflect potential (zero soil moisture stress) conditions by letting the ratio of actual transpiration to potential transpiration be a function of normalized difference vegetation index (NDVI). Soil moisture, evapotranspiration, and runoff data are generated on a daily basis for a 10-year period, January 1979 through December 1988, using observed precipitation gridded at a 4 deg by 5 deg resolution.

Liston, G. E.↗

What is Going on Around Cloud Edges: A View from Above and Below

tend to produce the greatest precipitation enhancements. Further analysis of the response over the U.S. to the Pacific forcing highlights a number of noteworthy and to some extent unexpected results. These include a seasonal dependence of the precipitation response that is characterized by signal-to-noise ratios that peak in spring, and surface temperature signal-to-noise ratios that are both lower and show less agreement among the models than those found for the precipitation response. Another interesting result concerns what appears to be a substantially different character in the surface temperature response over the U.S. to the Pacific forcing by the only model examined here that was developed for use in numerical weather prediction. The response to the positive SST trend forcing pattern is an overall surface warming over the world's land areas with substantial regional variations that are in part reproduced in runs forced with a globally uniform SST trend forcing. The precipitation response to the trend forcing is weak in all the models.

Marshak, Alexander↗

Combined Dust Detection Algorithm by Using MODIS Infrared Channels over East Asia

A new dust detection algorithm is developed by combining the results of multiple dust detectionmethods using IR channels onboard the MODerate resolution Imaging Spectroradiometer (MODIS). Brightness Temperature Difference (BTD) between two wavelength channels has been used widely in previous dust detection methods. However, BTDmethods have limitations in identifying the offset values of the BTDto discriminate clear-sky areas. The current algorithm overcomes the disadvantages of previous dust detection methods by considering the Brightness Temperature Ratio (BTR) values of the dual wavelength channels with 30-day composite, the optical properties of the dust particles, the variability of surface properties, and the cloud contamination. Therefore, the current algorithm shows improvements in detecting the dust loaded region over land during daytime. Finally, the confidence index of the current dust algorithm is shown in 10 × 10 pixels of the MODIS observations. From January to June, 2006, the results of the current algorithm are within 64 to 81% of those found using the fine mode fraction (FMF) and aerosol index (AI) from the MODIS and Ozone Monitoring Instrument (OMI). The agreement between the results of the current algorithm and the OMI AI over the non-polluted land also ranges from 60 to 67% to avoid errors due to the anthropogenic aerosol. In addition, the developed algorithm shows statistically significant results at four AErosol RObotic NETwork (AERONET) sites in East Asia.

East Asia↗

Restoration of the Apollo Heat Flow Experiments Metadata

Geothermal heat flow probes were deployed on the Apollo 15 and 17 missions as part of the Apollo Lunar Surface Experiments Package (ALSEP). At each landing site, the astronauts drilled 2 holes, 10-m apart, and installed a probe in each. The holes were 1- and 1.5-m deep at the Apollo 15 site and 2.5-m deep at the Apollo 17 sites. The probes monitored surface temperature and subsurface temperatures at different depths. At the Apollo 15 site, the monitoring continued from July 1971 to January 1977. At the Apollo 17 site, it did from December 1972 to September 1977. Based on the observations made through December 1974, Marcus Langseth, the principal investigator of the heat flow experiments (HFE), determined the thermal conductivity of the lunar regolith by mathematically modeling how the seasonal temperature fluctuation propagated down through the regolith. He also determined the temperature unaffected by diurnal and seasonal thermal waves of the regolith at different depths, which yielded the geothermal gradient. By multiplying the thermal gradient and the thermal conductivity, Langseth obtained the endogenic heat flow of the Moon as 21 mW/m(exp 2) at Site 15 and 16 mW/m(exp 2) at Site 17.

HFE metadata↗

Microwave Polarized Signatures Generated within Cloud Systems: SSM/I Observations Interpreted with Radiative Transfer Simulations

Special Sensor Microwave /Imager (SSM/I) observations in cloud systems are studied over the tropics. Over optically thick cloud systems, presence of polarized signatures at 37 and 85 GHz is evidenced and analyzed with the help of cloud top temperature and optical thickness extracted from visible and IR satellite observations. Scattering signatures at 85 GHz (TbV(85) less than or = 250 K) are associated with polarization differences greater than or = 6 K, approx. 50%, of the time over ocean and approx. 40% over land. In addition. over thick clouds the polarization difference at 37 GHz is rarely negligible. The polarization differences at 37 and 85 GHz do not stem from the surface but are generated in regions of relatively homogeneous clouds having high liquid water content. To interpret the observations, a radiative transfer model that includes the scattering by non-spherical particles is developed. based on the T-matrix approach and using the doubling and adding method. In addition to handling randomly and perfectly oriented particles, this model can also simulate the effect of partial orientation of the hydrometeors. Microwave brightness temperatures are simulated at SSM/I frequencies and are compared with the observations. Polarization differences of approx. 2 K can be simulated at 37 GHz over a rain layer, even using spherical drops. The polarization difference is larger for oriented non-spherical particles. The 85 GHz simulations are very sensitive to the ice phase of the cloud. Simulations with spherical particles or with randomly oriented non-spherical ice particles cannot replicate the observed polarization differences. However, with partially oriented non-spherical particles, the observed polarized signatures at 85 GHz are explained, and the sensitivity of the scattering characteristics to the particle size, asphericity, and orientation is analyzed. Implications on rain and ice retrievals are discussed.

Prigent, Catherine↗

Groundwater Recharge Estimated by Land Surface Models: An Evaluation in the Conterminous U.S.

Estimating diffuse recharge of precipitation is fundamental to assessing groundwater sustainability. Diffuse recharge is also the process through which climate and climate change directly affect groundwater. In this study, we evaluated diffuse recharge over the conterminous U.S. simulated by a suite of land surface models (LSMs) that were forced using a common set of meteorological input data. Simulated annual recharge exhibited spatial patterns that were similar among the LSMs, with the highest values in the eastern U.S. and Pacific Northwest. However, the magnitudes of annual recharge varied significantly among the models and were associated with differences in simulated ET, runoff and snow. Evaluation against two independent datasets did not answer the question of whether the ensemble mean performs the best, due to inconsistency between those datasets. The amplitude and timing of seasonal maximum recharge differed among the models, influenced strongly by model physics governing deep soil moisture drainage rates and, in cold regions, snowmelt. Evaluation using in situ soil moisture observations suggested that true recharge peaks 1-3 months later than simulated recharge, indicating systematic biases in simulating deep soil moisture. However, recharge from lateral flows and through preferential flows cannot be inferred from soil moisture data, and the seasonal cycle of simulated groundwater storage actually compared well with in situ groundwater observations. Long-term trends in recharge were not consistently correlated with either precipitation trends or temperature trends. This study highlights the need to employ dynamic flow models in LSMs, among other improvements, to enable more accurate simulation of recharge.

snowmelt↗

Lunar Dust Contamination Effects on Lunar Base Thermal Control Systems

Many studies have been conducted to develop a thermal control system that can operate under the extreme thermal environments found on the lunar surface. While these proposed heat rejection systems use different methods to reject heat, each system contains a similar component, a thermal radiator system. These studies have always considered pristine thermal control system components and have overlooked the possible deleterious effects of lunar dust contamination. Since lunar dust has a high emissivity and absorptivity (greater than 0.9) and is opaque, dust accumulation on a surface should radically alter its optical properties and therefore alter its thermal response compared to ideal conditions. In addition, the non-specular nature of the dust particles will alter the performance of systems that employ specular surfaces to enhance heat rejection. To date, few studies have examined the effect of dust deposition on the normal control system components. These studies only focused on a single heat rejection or photovoltaic system. These studies did show that lunar dust accumulations alter the optical properties of any lunar base hardware, which in turn affects component temperatures, and heat rejection. Therefore, a new study was conducted to determine the effect of lunar dust contamination on heat rejection systems. For this study, a previously developed dust deposition model was incorporated into the Thermal Synthesizer System (TSS) model. This modeling scheme incorporates the original method of predicting dust accumulation due to vehicle landings by assuming that the thin dust layer can be treated as a semitransparent surface slightly above and in thermal contact with the pristine surface. The results of this study showed that even small amounts of dust deposits can radically alter the performance of the heat rejection systems. Furthermore. this study indicates that heat rejection systems be either located far from any landing sites or be protected from dust producing mechanisms.

Keller, John R.↗

Seismoball: A Small Europa Orbiter Drop-Off Probe for Early Exploration of the Europan Surface

Recent magnetometry data received from Galileo indicate that the most likely explanation for the magnetic signature there is indeed a global conducting layer below the surface. This conducting layer is well matched by a salty, mineral rich strata beneath the Europan ice crust or a salt water ocean. Galileo imaging results show a variety of terrain types thought to contain young material; for example, lineaments, chaotic terrain, and eruption features. Additionally, Galileo images have shown indications of areas of up-welling where subsurface material periodically gets pushed to the surface due to the forces of fracturing, butting, and refreezing of the ice sheet. While Europa Orbiter will provide close-flyby high resolution images, as well as magnetometry, spectroscopy and other remote sensing data of the surface, it will not be able to provide essential engineering data like surface hardness and surface ice structure needed to support eventual landed missions. Additionally, ice chemical composition at microscopic scales can only be studied in detail through in situ instrumentation. Seismoball is a small probe designed to be injected into a surface intersect orbit around Europa. Using small reverse thrusters, the probe will be capable of nulling the high horizontal injection velocity as it approaches the 2 km surface injection altitude, thus allowing it to fall to the surface at an impact velocity of < 100m/sec (much less than the DS-2 impact velocities). The external breakaway thruster structure and crushable exterior shell absorb the impact energy while allowing the science instrument suite to remain intact. JPL has already started analyzing the entry dynamics and designing/building a small, low mass probe which will withstand the impact g-forces and fit as a 'carry-on' on board the Europa Orbiter. The probe will carry a suite of 5-6 micro-instruments for imaging the surface (both microscopic and far-field), surface and shallow subsurface ice temperatures, surface hardness, crustal dynamics and periodicity, and compositional chemistry. If selected, this flight development activity will provide a unique science opportunity and adjunct to the primary Orbiter science mission. The final flight system will be designed to accommodate orbiter mass, volume, and power interface constraints, as well as entry dynamics, g-load mitigation, and arbitrary landing orientation.

Tamppari, L.↗

The Mauna Loa sulfur flow as an analog to secondary sulfur flows on Io

Independent evidence suggests that both sulfur and silicate materials exist on the surface of Io. Spectral data indicate the presence of sulfur compounds, some of which are suggested to be of fumarolic origin. Morphological evidence and inferences of the physical properties of some land-forms suggest that silicate volcanism has occurred, which would involve temperatures 650 C or greater. Because the liquidus of sulfur is only about 115 C, it is likely that sulfur in close proximity to 'hot spots' or to active silicate volcanic areas on Io would be melted and mobilized as flows. The Mauna Loa sulfur flow may serve as an analog for such flows, as it consists of fumarolic sulfur that was melted as a consequence of a basaltic eruption and produced a small flow superimposed on silicate lavas.

Greeley, R.↗

Fusion Crust and the Measurement of Surface Ages of Antarctic Ordinary Chondrites

Natural thermoluminescence (TL) reflects radiation exposure and storage temperature. Meteorites generally exhibit thermoluminescence acquired during their long exposure to galactic cosmic rays in space. During atmospheric passage, temperatures are high enough to completely drain the TL, in the first mm of material under the fusion crust. We therefore refer to this surface layer as "fusion crust" although it does include some unmelted material just below the crust. When the meteorite lands on earth this drained layer will begin to build up natural TL once again due to radiation from cosmic rays and internal radionuclides. Cosmic ray annual dose is estimated to be between 0.04 and 0.06 rad/yr on the earth's surface in Antarctica while the internal radionuclides contribute only about 0.01 rad/yr. Therefore the total annual dose received by the meteorite while it is on the surface is between 0.05 and 0.07 rad/yr. If the meteorite is buried deeply in the ice it is effectively shielded from most cosmic rays and thus only internal radioactivity contributes to the annual dose.

Akridge, Jannette M. C.↗

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↗

Does Aerosol Weaken or Strengthen the South Asian Monsoon?

Aerosols are known to have the ability to block off solar radiation reaching the earth surface, causing it to cool - the so-called solar dimming (SDM) effect. In the Asian monsoon region, the SDM effect by aerosol can produce differential cooling at the surface reducing the meridional thermal contrast between land and ocean, leading to a weakening of the monsoon (Ramanathan et al. 2005). On the other hand, absorbing aerosols such as black carbon and dust, when forced up against the steep slopes of the southern Tibetan Plateau can produce upper tropospheric heating, and induce convection-dynamic feedback leading to an advance of the rainy season over northern India and an enhancement of the South Asian monsoon through the "Elevated Heat Pump" (EHP) effect (Lau et al. 2006). In this paper, we present modeling results showing that in a coupled ocean-atmosphere-land system in which concentrations of greenhouse gases are kept constant, the response of the South Asian monsoon to dust and black carbon forcing is the net result of the two opposing effects of SDM and EHP. For the South Asian monsoon, if the increasing upper tropospheric thermal contrast between the Tibetan Plateau and region to the south spurred by the EHP overwhelms the reduction in surface temperature contrast due to SDM, the monsoon strengthens. Otherwise, the monsoon weakens. Preliminary observations are consistent with the above findings. We find that the two effects are strongly scale dependent. On interannual and shorter time scales, the EHP effect appears to dominate in the early summer season (May-June). On decadal or longer time scales, the SDM dominates for the mature monsoon (July-August). Better understanding the physical mechanisms underlying the SDM and the EHP effects, the local emission and transport of aerosols from surrounding deserts and arid-regions, and their interaction with monsoon water cycle dynamics are important in providing better prediction and assessment of climate change impacts on precipitation of the Asian monsoon land regions.

Lau, William K. M.↗