Search NASA⌕ Search

SEARCH · Search NASA

Results for “forecast sensitivity to observations”

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 199 records · Page 11

Local Scale (3-M) Soil Moisture Mapping Using SMAP and Planet Superdove

A capability for mapping meter-level resolution soil moisture with frequent temporal sampling over large regions is essential for quantifying local-scale environmental heterogeneity and eco-hydrologic behavior. However, available surface soil moisture (SSM) products generally involve much coarser grain sizes ranging from 30 m to several 10s of kilometers. Hence a new method is proposed to estimate 3-m resolution SSM using a combination of multi-sensor fusion, machine- learning (ML) and Cumulative Distribution Function (CDF) matching approaches. This method established favorable SSM correspondence between 3-m pixels and overlying 9-km grid cells from overlapping Planet SuperDove (PSD) observations and NASA Soil Moisture Active-Passive (SMAP) mission products. The resulting 3-m SSM predictions showed improved accuracy by reducing ab- solute bias and RMSE by ~0.01 cm3/cm3 over the original SMAP data in relation to in-situ soil moisture measurements for the Australian Yanco region, while preserving the high sampling frequency (1-3 day global revisit) and sensitivity to surface wetness (R 0.865) from SMAP. Heterogeneous soil moisture distributions varying with vegetation biomass gradients and irrigation regimes were generally captured within a selected study area. Further algorithm refinement and implementation for regional applications will allow for improvement in water resources management, precision agriculture, and disaster forecasts and responses.

soil moisture↗

Development of a climatological data base to help forecast cloud cover conditions for shuttle landings at the Kennedy Space Center

The Space Shuttle is an extremely weather sensitive vehicle with very restrictive constraints for both launches and landings. The most important difference between Shuttle and normal aircraft landings is that the Shuttle has no go-around capability once it begins its decent into the earth's atmosphere. The de-orbit burn decision is generally made approximately 90 minutes before landing requiring a forecast with little room for error. Because of the Shuttle's rapid re-entry to earth, the pilot must be able to see all runway and visual navigation aids from high altitude to land the Shuttle. In addition, the heat resistant tiles which are used to protect the Shuttle during its re-entry into the earth's atmosphere are extremely sensitive to any type of precipitation. Extensive damage to these tiles could occur if the Shuttle passes through any cloud that contains precipitation size particles. To help guard against changing weather conditions or any type of weather problems that might occur prior to landing, flight rules have been developed as guidelines for all landings. Although the rules vary depending on the location of the landing (Kennedy Space Center or Edwards AFB), length of mission, and weight of vehicle, most of the rules can be condensed into 4 major groupings. These are: (1) Cloud ceilings should not be less than 3048 m (10,000 feet), (2) Visibility should not be less than 13 km (7 nm), (3) Cross-wind no greater than 5-8 m/s (10-15 knots); and (4) No showers or thunderstorms at or within 56 km (30 nm) of the Shuttle Landing Facility. This study consisted of developing a climatological database of the Shuttle Landing Facility (SLF) surface observations and performing an analysis of observed conditions one and two hours subsequent to given conditions at the SLF to help analyze the 0.2 cloud cover rule. Particular emphasis was placed on Shuttle landing weather violations and the amounts of cloud cover below 3048 m (10,000 ft.). This analysis has helped to determine the best and worst times to land the Shuttle at KSC. In addition, nomograms have been developed to help forecasters make cloud cover forecasts for End of Mission (EOM) and Return to Launch Site (RTLS) at KSC. Results of categorizing this data by month, season, time of day, and surface and upper-air wind direction are presented.

Atchison, M. Kevin↗

A Wrf-Chem Flash Rate Parameterization Scheme and LNO(x) Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

The Deep Convective Clouds and Chemistry (DC3) field campaign in 2012 provided a plethora of aircraft and ground-based observations (e.g., trace gases, lightning and radar) to study deep convective storms, their convective transport of trace gases, and associated lightning occurrence and production of nitrogen oxides (NOx). Based on the measurements taken of the 29-30 May 2012 Oklahoma thunderstorm, an analysis against a Weather Research and Forecasting Chemistry (WRF-Chem) model simulation of the same event at 3-km horizontal resolution was performed. One of the main objectives was to include various flash rate parameterization schemes (FRPSs) in the model and identify which scheme(s) best captured the flash rates observed by the National Lightning Detection Network (NLDN) and Oklahoma Lightning Mapping Array (LMA). The comparison indicates how well the schemes predicted the timing, location, and number of lightning flashes. The FRPSs implemented in the model were based on the simulated thunderstorms physical features, such as maximum vertical velocity, cloud top height, and updraft volume. Adjustment factors were added to each FRPS to best capture the observed flash trend and a sensitivity study was performed to compare the range in model-simulated lightning-generated nitrogen oxides (LNOx) generated by each FRPS over the storms lifetime. Based on the best FRPS, model-simulated LNOx was compared against aircraft measured NOx. The trace gas analysis, along with the increased detail in the model specification of the vertical distribution of lightning flashes as suggested by the LMA data, provide guidance in determining the scenario of NO production per intracloud and cloud-to-ground flash that best matches the NOx mixing ratios observed by the aircraft.

Cloud physics↗

A WRF-Chem Flash Rate Parameterization Scheme and LNOx Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

The Deep Convective Clouds and Chemistry (DC3) field campaign in 2012 provided a plethora of aircraft and ground-based observations (e.g., trace gases, lightning and radar) to study deep convective storms, their convective transport of trace gases, and associated lightning occurrence and production of nitrogen oxides (NOx). Based on the measurements taken of the 29-30 May 2012 Oklahoma thunderstorm, an analysis against a Weather Research and Forecasting Chemistry (WRF-Chem) model simulation of the same event at 3-km horizontal resolution was performed. One of the main objectives was to include various flash rate parameterization schemes (FRPSs) in the model and identify which scheme(s) best captured the flash rates observed by the National Lightning Detection Network (NLDN) and Oklahoma Lightning Mapping Array (LMA). The comparison indicates how well the schemes predicted the timing, location, and number of lightning flashes. The FRPSs implemented in the model were based on the simulated thunderstorms physical features, such as maximum vertical velocity, cloud top height, and updraft volume. Adjustment factors were applied to each FRPS to best capture the observed flash trend and a sensitivity study was performed to compare the range in model-simulated lightning-generated nitrogen oxides (LNOx) generated by each FRPS over the storms lifetime. Based on the best FRPS, model-simulated LNOx was compared against aircraft measured NOx. The trace gas analysis, along with the increased detail in the model specification of the vertical distribution of lightning flashes as suggested by the LMA data, provide guidance in determining the scenario of NO production per intracloud and cloud-to-ground flash that best matches the NOx mixing ratios observed by the aircraft.

Atmospheric chemistry↗

A Mission to Observe Ice in Clouds from Space

To date there have been multiple satellite missions to observe and retrieve cloud top properties and the liquid in, and precipitation from, clouds. There are currently a few missions that attempt to measure cloud ice properties as a byproduct of other observations. However, we do not yet quantitatively understand the processes that control the water budget of the upper troposphere where ice is the predominant phase, and how these processes are linked to precipitation processes and the radiative energy budget. The ice in clouds either melts into rain or is detrained, and persists, as cirrus clouds affecting the hydrological and energy cycle, respectively. Fully modeling the Earth's climate and improving weather and climate forecasts requires accurate satellite measurements of various cloud properties at the temporal and spatial scales of cloud processes. The uncertainty in knowledge of these ice characteristics is reflected in the large discrepancies in model simulations of the upper tropospheric water budget. Model simulations are sensitive to the partition of ice between precipitation and outflow processes, i.e., to the parameterization of ice clouds and ice processes. This presentation will describe the Submillimeter-wave InfraRed Ice Cloud Experiment (SIRICE) concept, a satellite mission designed to acquire global Earth radiance measurements in the infrared and submillimeter-wave region (183-874 GHz). If successful, this mission will bridge the measurement gap between microwave sounders and shorter-wavelength infrared and visible sensors. The brightness temperatures at submillimeter-wave frequencies are especially sensitive to cirrus ice particle sizes (because they are comparable to the wavelength). This allows for more accurate ice water path estimates when multiple channels are used to probe into the cloud layers. Further, submillimeter wavelengths offer simplicity in the retrieval algorithms because they do not probe into the liquid and near surface portions of clouds, thus requiring only one term of the radiative transfer equation (ice scattering) to relate brightness temperatures to ice. Scientific justification and the SIRICE approach to measuring ice water path and particle size that span a range encompassing both the hydrologically active and radiatively active components of cloud systems will be presented.

Ackerman, S.↗

Diagnosing the Sensitivity of Local Land-Atmosphere Coupling via the Soil Moisture-Boundary Layer Interaction

The inherent coupled nature of earth s energy and water cycles places significant importance on the proper representation and diagnosis of land atmosphere (LA) interactions in hydrometeorological prediction models. However, the precise nature of the soil moisture precipitation relationship at the local scale is largely determined by a series of nonlinear processes and feedbacks that are difficult to quantify. To quantify the strength of the local LA coupling (LoCo), this process chain must be considered both in full and as individual components through their relationships and sensitivities. To address this, recent modeling and diagnostic studies have been extended to 1) quantify the processes governing LoCo utilizing the thermodynamic properties of mixing diagrams, and 2) diagnose the sensitivity of coupled systems, including clouds and moist processes, to perturbations in soil moisture. This work employs NASA s Land Information System (LIS) coupled to the Weather Research and Forecasting (WRF) mesoscale model and simulations performed over the U.S. Southern Great Plains. The behavior of different planetary boundary layers (PBL) and land surface scheme couplings in LIS WRF are examined in the context of the evolution of thermodynamic quantities that link the surface soil moisture condition to the PBL regime, clouds, and precipitation. Specifically, the tendency toward saturation in the PBL is quantified by the lifting condensation level (LCL) deficit and addressed as a function of time and space. The sensitivity of the LCL deficit to the soil moisture condition is indicative of the strength of LoCo, where both positive and negative feedbacks can be identified. Overall, this methodology can be applied to any model or observations and is a crucial step toward improved evaluation and quantification of LoCo within models, particularly given the advent of next-generation satellite measurements of PBL and land surface properties along with advances in data assimilation schemes.

Santanello, Joseph A., Jr.↗

Analysis of Multiple Precipitation Products and Preliminary Assessment of Their Impact on Global Land Data Assimilation System (GLDAS) Land Surface States

Land surface models (LSMs) are computer programs, similar to weather and climate prediction models, which simulate the stocks and fluxes of water (including soil moisture, snow, evaporation, and runoff) and energy (including the temperature of and sensible heat released from the soil) after they arrive on the land surface as precipitation and sunlight. It is not currently possible to measure all of the variables of interest everywhere on Earth with sufficient accuracy and space-time resolution. Hence LSMs have been developed to integrate the available observations with our understanding of the physical processes involved, using powerful computers, in order to map these stocks and fluxes as they change in time. The maps are used to improve weather forecasts, support water resources and agricultural applications, and study the Earth's water cycle and climate variability. NASA's Global Land Data Assimilation System (GLDAS) project facilitates testing of several different LSMs with a variety of input datasets (e.g., precipitation, plant type). Precipitation is arguably the most important input to LSMs. Many precipitation datasets have been produced using satellite and rain gauge observations and weather forecast models. In this study, seven different global precipitation datasets were evaluated over the United States, where dense rain gauge networks contribute to reliable precipitation maps. We then used the seven datasets as inputs to GLDAS simulations, so that we could diagnose their impacts on output stocks and fluxes of water. In terms of totals, the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) had the closest agreement with the US rain gauge dataset for all seasons except winter. The CMAP precipitation was also the most closely correlated in time with the rain gauge data during spring, fall, and winter, while the satellitebased estimates performed best in summer. The GLDAS simulations revealed that modeled soil moisture is highly sensitive to precipitation, with differences in spring and summer as large as 45% depending on the choice of precipitation input.

Gottschalck, Jon↗

Section on Observed Impacts on El Nino

Agricultural applications of El Nino forecasts are already underway in some countries and need to be evaluated or re-evaluated. For example, in Peru, El Nino forecasts have been incorporated into national planning for the agricultural sector, and areas planted with rice and cotton (cotton being the more drought-tolerant crop) are adjusted accordingly. How well are this and other such programs working? Such evaluations will contribute to the governmental and intergovernmental institutions, including the Inter-American Institute for Global Change Research and the US National Ocean and Atmospheric Agency that are fostering programs to aid the effective use of forecasts. As El Nino climate forecasting grows out of the research mode into operational mode, the research focus shifts to include the design of appropriate modes of utilization. Awareness of and sensitivity to the costs of prediction errors also grow. For example, one major forecasting model failed to predict the very large El Nino event of 1997, when Pacific sea-surface temperatures were the highest on record. Although simple correlations between El Nino events and crop yields may be suggestive, more sophisticated work is needed to understand the subtleties of the interplay among the global climate system, regional climate patterns, and local agricultural systems. Honesty about the limitations of an forecast is essential, especially when human livelihoods are at stake. An end-to-end analysis links tools and expertise from the full sequence of ENSO cause-and-effect processes. Representatives from many disciplines are needed to achieve insights, e.g, oceanographers and atmospheric scientists who predict El Nino events, climatologists who drive global climate models with sea-surface temperature predictions, agronomists who translate regional climate connections in to crop yield forecasts, and economists who analyze market adjustments to the vagaries of climate and determine the value of climate forecasts. Methods include historical studies to understand past patterns and to test hindcasts of the prediction tools, crop modeling, spatial analysis and remote sensing. This research involves expanding, deepening, and applying the understanding of physical climate to the fields of agronomy and social science; and the reciprocal understanding of crop growth and farm economics to climatology. Delivery of a regional climate forecast with no information about how the climate forecast was derived limits its effectiveness. Explanation of a region's major climate driving forces helps to place a seasonal forecast in context. Then, a useful approach is to show historical responses to previous El Nino events, and projections, with uncertainty intervals, of crop response from dynamic process crop growth models. Regional ID forecasts should be updated with real-time weather conditions. Since every El Nino event is different, it is important to track, report and advise on each new event as it unfolds. The stability of human enterprises depends on understanding both the potentialities and the limits of predictability. Farmers rely on past experience to anticipate and respond to fluctuations in the biophysical systems on which their livelihoods depend. Now scientists are improving their ability to predict some major elements of climate variability. The improvements in the reliability of El Nino forecasts are encouraging, but seasonal forecasts for agriculture are not, and will probably never be completely infallible, due to the chaotic nature of the climate system. Uncertainties proliferate as we extend beyond Pacific sea-surface temperatures to climate teleconnections and agricultural outcomes. The goal of this research is to shed as a clear light as possible on these inherent uncertainties and thus to contribute to the development of appropriate responses to El Nino and other seasonal forecasts for a range of stakeholders, which, ultimately, includes food consumers everywhere.

Rosenzweig, Cynthia↗

Forecasting disease risk for increased epidemic preparedness in public health

Emerging infectious diseases pose a growing threat to human populations. Many of the world's epidemic diseases (particularly those transmitted by intermediate hosts) are known to be highly sensitive to long-term changes in climate and short-term fluctuations in the weather. The application of environmental data to the study of disease offers the capability to demonstrate vector-environment relationships and potentially forecast the risk of disease outbreaks or epidemics. Accurate disease forecasting models would markedly improve epidemic prevention and control capabilities. This chapter examines the potential for epidemic forecasting and discusses the issues associated with the development of global networks for surveillance and prediction. Existing global systems for epidemic preparedness focus on disease surveillance using either expert knowledge or statistical modelling of disease activity and thresholds to identify times and areas of risk. Predictive health information systems would use monitored environmental variables, linked to a disease system, to be observed and provide prior information of outbreaks. The components and varieties of forecasting systems are discussed with selected examples, along with issues relating to further development.

Disease Outbreaks↗

Regional Model Nesting Within GFS Daily Forecasts Over West Africa

The study uses the RM3, the regional climate model at the Center for Climate Systems Research of Columbia University and the NASA/Goddard Institute for Space Studies (CCSR/GISS). The paper evaluates 30 48-hour RM3 weather forecasts over West Africa during September 2006 made on a 0.5 grid nested within 1 Global Forecast System (GFS) global forecasts. September 2006 was the Special Observing Period #3 of the African Monsoon Multidisciplinary Analysis (AMMA). Archived GFS initial conditions and lateral boundary conditions for the simulations from the US National Weather Service, National Oceanographic and Atmospheric Administration were interpolated four times daily. Results for precipitation forecasts are validated against Tropical Rainfall Measurement Mission (TRMM) satellite estimates and data from the Famine Early Warning System (FEWS), which includes rain gauge measurements, and forecasts of circulation are compared to reanalysis 2. Performance statistics for the precipitation forecasts include bias, root-mean-square errors and spatial correlation coefficients. The nested regional model forecasts are compared to GFS forecasts to gauge whether nesting provides additional realistic information. They are also compared to RM3 simulations driven by reanalysis 2, representing high potential skill forecasts, to gauge the sensitivity of results to lateral boundary conditions. Nested RM3/GFS forecasts generate excessive moisture advection toward West Africa, which in turn causes prodigious amounts of model precipitation. This problem is corrected by empirical adjustments in the preparation of lateral boundary conditions and initial conditions. The resulting modified simulations improve on the GFS precipitation forecasts, achieving time-space correlations with TRMM of 0.77 on the first day and 0.63 on the second day. One realtime RM3/GFS precipitation forecast made at and posted by the African Centre of Meteorological Application for Development (ACMAD) in Niamey, Niger is shown.

Druyan, Leonard M.↗

Impact of Optimized land Surface Parameters on the Land-Atmosphere Coupling in WRF Simulations of Dry and Wet Extremes

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty module in NASA's Land Information System (LIS-OPT), whereby parameter sets are calibrated in the Noah land surface model and classified according to the land cover and soil type mapping of the observations and the full domain. The impact of the calibrated parameters on the a) spin up of land surface states used as initial conditions, and b) heat and moisture fluxes of the coupled (LIS-WRF) simulations are then assessed in terms of ambient weather, PBL budgets, and precipitation along with L-A coupling diagnostics. In addition, the sensitivity of this approach to the period of calibration (dry, wet, normal) is investigated. Finally, tradeoffs of computational tractability and scientific validity (e.g.,. relating to the representation of the spatial dependence of parameters) and the feasibility of calibrating to multiple observational datasets are also discussed.

Kumar, Sujay↗

Assessing the Usefulness of AIRS Radiance Observations in a 4D-Var Assimilation Scheme Using the Penn State/NCAR Mesoscale Model Version 5 (MM5) and a Stand Alone Radiative Transfer Algorithm (SARTA)

Short-term precipitation forecast skill remains poor despite further advances in numerical modeling and data analysis. Utilizing indirect observations, such as the radiance observations obtained from the Atmospheric Infrared Sounder (AIRS), in a 4D-Var assimilation scheme may lead to a rapid improvement in the skill of 12-h QPF. In order to investigate the usefulness of AIRS observations, the adjoint code of a radiative transfer model (SARTA) is developed and linked to MM5 4D-Var system. A statistical comparison of the AIRS radiances with those calculated using the linked MMS-SARTA model is carried out to assess errors in the radiative transfer model, MM5 and AIRS data, and to develop a quality control scheme. Outliers and any bias that may exist in radiative transfer model, MM5 and AIRS data are identified. An adjoint sensitivity study, utilizing the SARTA adjoint model, is then done to identify the relative sensitivities of each spectral channel to such parameters as atmospheric temperature, pressure, and mixing ratio. Another sensitivity study is done using the MM5 alone to determine the sensitivity of the 12-h QPF to model parameters such as atmospheric temperature, pressure, and mixing ratio under various weather regimes (i.e. summertime convection, winter snowfall). Results from these two types of adjoint sensitivity studies are used to determine which AIRS channels contribute most to 12-h QPF under differing weather regimes. This information is useful when conducting 4D-Var experiments using a subset of the 2378 AIRS spectral channels. Finally, 4D-Var experiments are done to analyze the impact of the AIRS observations on the quality of the 12-h QPF results. This project falls under the data assimilation general session topic as specified by the conference organizers.

Carrier, M.↗

Timing the Flames: Geostationary Satellite Detection of Diurnally Shifting Stubble Burning in Northwestern India

Post-monsoon open-field stubble burning in northwestern (NW) India—a key agricultural region known as the “breadbasket”—is a longstanding practice used to clear fields. Satellite observations spanning over two decades have revealed significant upward trends in crop production, vegetative greenness, and the frequency of post-harvest fires, with this last contributing to hazardous air quality during the peak burning season (mid-October to mid-November). Since 2022, thermal anomaly data from Aqua-MODIS and SNPP-VIIRS sensors have shown a sharp decline in reported fire events—an observation that contrasts starkly with the concurrent rise in regional aerosol loading detected from space. This apparent discrepancy became particularly pronounced in 2024–2025, prompting a closer examination using high-temporal-resolution imagery from the Advanced Meteorological Imager (AMI) on the geostationary satellite GEO-KOMPSAT-2A. These observations revealed a clear spike in fire-related signals occurring around and after 4:00 p.m. local time, i.e., outside the typical noon to 2:00 p.m. detection window of the MODIS and VIIRS. A fire detection algorithm exploiting the fire-sensitive shortwave-infrared 3.8 μm signal and its contrast to 11.2 μm infrared observations is designed to adopt AMI observations and applied to its multi-year observations (2019–2025). The resulting fire dataset unambiguously shows a gradual shift in stubble burning activity toward the late afternoon hours beginning in 2022 which is underreported by polar-orbiting satellites. The orbital drift of NASA’s MODIS sensor on the Aqua platform allows detection of some of the gradually shifting fires during afternoon hours, but the MODIS still misses a large number of fires occurring around and after 4 p.m. The AMI’s relatively coarse spatial resolution (~4 km), a consequence of its slant viewing geometry over NW India, imposes inherent limitations on quantifying the full extent of fire occurrences. The operational air quality forecasting models currently assimilate satellite fire detections predominantly captured during early afternoon overpasses of the MODIS and VIIRS. The temporal shift in fire activity complicates such forecast, leading to a substantial underestimation of emissions. Intense stubble burning and the resulting air pollution highlight the need for effective crop residue management practices for mitigating the frequency of open biomass burning and thereby reducing episodic degradation of air quality and its associated public health and economic impacts.

post-monsoon stubble burning; northwestern India;↗

TPSAS-NF1676L-20351-DND

Characterization of atmospheric variability is crucial to the design of the Entry, Descent, and Landing (EDL) activities at Mars. The distribution of dust in the Martian atmosphere has a significant impact on the atmospheric structure and variability. Information obtained through reconstruction of EDL data from past missions can be used in studying the performance of a model used in preflight atmospheric characterization for future missions. We use the Mars Weather Research and Forecasting (MarsWRF) General Circulation Model to simulate the atmospheric temperature during the landing time periods of Mars Exploration Rovers MER-A and MER-B in 2004. The multi-scale capability of MarsWRF is useful for conducting high- resolution nested simulation centered at the landing site. This version of MarsWRF includes a high spectral accuracy k-distribution radiative transfer model. We use the IR dust opacity profile data from the Thermal Emission spectrometer (TES) limb scan data as well as the TES Nadir data to prescribe the dust opacity distribution. Data from MY 26, which corresponds to the year of MER landing, and data from a climatology based on the median of 4 years (MY24 to MY27) of TES data are used. Model temperature profiles are compared with the reconstructions from EDL data and TES observations. Sensitivity of the model simulation to variations in the input opacity data are shown.

Murali Natarajan↗

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↗

The Sensitivity of WRF Daily Summertime Simulations over West Africa to Alternative Parameterizations. Part 1: African Wave Circulation

The performance of the NCAR Weather Research and Forecasting Model (WRF) as a West African regional-atmospheric model is evaluated. The study tests the sensitivity of WRF-simulated vorticity maxima associated with African easterly waves to 64 combinations of alternative parameterizations in a series of simulations in September. In all, 104 simulations of 12-day duration during 11 consecutive years are examined. The 64 combinations combine WRF parameterizations of cumulus convection, radiation transfer, surface hydrology, and PBL physics. Simulated daily and mean circulation results are validated against NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA) and NCEP/Department of Energy Global Reanalysis 2. Precipitation is considered in a second part of this two-part paper. A wide range of 700-hPa vorticity validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve correlations against reanalysis of 0.40-0.60 and realistic amplitudes of spatiotemporal variability for the 2006 focus year while a parallel-benchmark simulation by the NASA Regional Model-3 (RM3) achieves higher correlations, but less realistic spatiotemporal variability. The largest favorable impact on WRF-vorticity validation is achieved by selecting the Grell-Devenyi cumulus convection scheme, resulting in higher correlations against reanalysis than simulations using the Kain-Fritch convection. Other parameterizations have less-obvious impact, although WRF configurations incorporating one surface model and PBL scheme consistently performed poorly. A comparison of reanalysis circulation against two NASA radiosonde stations confirms that both reanalyses represent observations well enough to validate the WRF results. Validation statistics for optimized WRF configurations simulating the parallel period during 10 additional years are less favorable than for 2006.

variability↗

The GOES-R Lightning Mapper Sensor

The Lightning Mapper Sensor on GOES-R builds on previous measurements of lightning from low earth orbit by the OTD (Optical Transient Detector) and LIS (Lightning Imaging Sensor) sensors. Unlike observations from low earth orbit, the GOES-R platform will allow continuous monitoring of lightning activity over the Continental United States and southern Canada, Central and South America, and portions of the Atlantic and Pacific Oceans. The LMS will detect total (cloud-to-ground and intracloud) lightning at storm scale resolution (approx. 8 km) using a highly sensitive Charge Coupled Device (CCD) detector array. Discrimination between lightning optical transients and a bright sunlit background scene is accomplished by employing spectral, spatial, and temporal filtering along with a background subtraction technique. The result is 24 hour detection capability of total lightning. These total lightning observations can be made available to users within about 20 seconds. Research indicates a number of ways that total lightning observations from LMS could benefit operational activities, including 1) potential increases in lead times and reduced false alarms for severe thunderstorm and tornado Warnings, 2) improved routing of &rail around thunderstorms, 3) support for spacecraft launches and landings, 4) improved ability to monitor tropical cyclone intensity, 5) ability to monitor thunderstorm intensification/weakening during radar outages or where radar coverage is poor, 6) better identification of deep convection for the initialization of numerical prediction models, 7) improved forest fire forecasts, 8) identification of convective initiation, 9) identification of heavy convective snowfall, and 10) enhanced temporal resolution of storm evolution (1 minute) than is available from radar observations. Total lightning data has been used in an operational environment since July 2003 at the Huntsville, Alabama National Weather Service office. Total lightning measurements are obtained by the North Alabama Lightning Mapping Array (LMA) and have successfully been used in warning decisions. Every 2 minutes, total lightning counts in 2 km by 2 km horizontal, 1 km vertical grids are available to forecasters on an AWIPS (Advanced Weather Interactive Processing System) workstation. Experience with the LMA total lightning data is used to illustrate the potential use of LMS data that would be available to forecasters across the US. This abstract is for submission as a presentation to the National Weather Association Annual Meeting to be held 16-21 October 2004 in Portland, OR. This abstract will be published in the conference proceedings.

Buechler, Dennis↗

Evaluation of Model Microphysics Within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is ~0.25 meters per second too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were ~0.25 meters per second too slow, while the SBU-YLIN was 0.25 to 0.5 meters per second too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗