Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data parameterized events”

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 163 records · Page 9

The Impact of the Evolving Satellite Data Record on Reanalysis Water and Energy Fluxes During the Past 30 Years

Retrospective analyses (reanalyses) use a fixed assimilation model to take diverse observations and synthesize consistent, time-dependent fields of state variables and fluxes (e.g. temperature, moisture, momentum, turbulent and radiative fluxes). Because they offer data sets of these quantities at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Of course, one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics. Another issue, perhaps less appreciated, is the fact that the discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (~ the last 30 years). In particular we highlight the effect on reanalyses of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We examine these data sets from two perspectives. The first is the ability to capture modes of variability that have coherent spatial structure (e.g. ENSO events and near-decadal coupling to SST changes) and how these modes are contained within trends in near global averages of key quantities. Secondly, we consider diagnostics that measure the consistency in energetic scaling in the hydrologic cycle, particularly the fractional changes in column-integrated water vapor versus precipitation as they are coupled to radiative flux constraints. These results will be discussed in the context of implications for science objectives and priorities of the NASA Energy and Water Cycle Study, NEWS.

Robertson, Franklin R.↗

Investigating Extreme Events in the NASA GEOS Multiscale Modeling System

Advances in computing capabilities and scientific development have come together to evolve general circulation models into multi-scale Earth system modeling tools. The Goddard Earth Observing System (GEOS) model is one such example of this evolution. The GEOS model is driven by the finite-volume cubed-sphere (FV3) non-hydrostatic dynamical core. Surrounding FV3 is a scale-aware physics package and data assimilation capability permitting multi-scale application of GEOS for sub-seasonal to seasonal climate prediction, medium range weather prediction, and global mesoscale modeling at convection allowing resolutions. GEOS also includes a comprehensive chemistry package representing a range of capabilities from basic chemistry and interactve aerosols and gaseous species, to carbon emissions and uptake, and complex ozone photochemistry. In this study, we apply the GEOS model to study the fidelity of these processes with increasing horizontal resolution and scale-aware processes in GEOS on extreme events. The GEOS model is run for 40-days beginning in August 2016 at three uniform global resolutions of 13-km (c768), 6-km (c1536) and 3-km (c3072) with 72 vertical levels up to 0.01mb. The model physics use the Grell-Freitas scale-aware convection scheme to dynamically reduce the role of parameterized deep convection as resolved scale processes in the model take over at higher resolutions. We will include high-resolution global emissions and fluxes of aerosols and carbon downscaled from recent satellite observations. We will compare these simulations with reanalyses and observations, focusing on rainfall, clouds and radiative forcing at hourly to monthly timescales. We will closely examine the probability distribution of precipitation intensities and radiative properties of clouds on a daily time scale. In addition, we will focus on extreme events, in particular the diurnal cycle of convection over the US and the frequency and physical nature of organized convection and heavy rain events across the globe.

Putman, William↗

Characterizing and Mitigating Intraday Variability: Reconstructing Source Structure in Accreting Black Holes with mm-VLBI

The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87* and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse (u, v) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87*. We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87*, finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.

79 ASTRONOMY AND ASTROPHYSICS↗

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗

Winter QPF Sensitivities to Snow Parameterizations and Comparisons to NASA CloudSat Observations

Steady increases in computing power have allowed for numerical weather prediction models to be initialized and run at high spatial resolution, permitting a transition from larger scale parameterizations of the effects of clouds and precipitation to the simulation of specific microphysical processes and hydrometeor size distributions. Although still relatively coarse in comparison to true cloud resolving models, these high resolution forecasts (on the order of 4 km or less) have demonstrated value in the prediction of severe storm mode and evolution and are being explored for use in winter weather events . Several single-moment bulk water microphysics schemes are available within the latest release of the Weather Research and Forecast (WRF) model suite, including the NASA Goddard Cumulus Ensemble, which incorporate some assumptions in the size distribution of a small number of hydrometeor classes in order to predict their evolution, advection and precipitation within the forecast domain. Although many of these schemes produce similar forecasts of events on the synoptic scale, there are often significant details regarding precipitation and cloud cover, as well as the distribution of water mass among the constituent hydrometeor classes. Unfortunately, validating data for cloud resolving model simulations are sparse. Field campaigns require in-cloud measurements of hydrometeors from aircraft in coordination with extensive and coincident ground based measurements. Radar remote sensing is utilized to detect the spatial coverage and structure of precipitation. Here, two radar systems characterize the structure of winter precipitation for comparison to equivalent features within a forecast model: a 3 GHz, Weather Surveillance Radar-1988 Doppler (WSR-88D) based in Omaha, Nebraska, and the 94 GHz NASA CloudSat Cloud Profiling Radar, a spaceborne instrument and member of the afternoon or "A-Train" of polar orbiting satellites tasked with cataloguing global cloud characteristics. Each system provides a unique perspective. The WSR-88D operates in a surveillance mode, sampling cloud volumes of Rayleigh scatterers where reflectivity is proportional to the sixth moment of the size distribution of equivalent spheres. The CloudSat radar provides enhanced sensitivity to smaller cloud ice crystals aloft, as well as consistent vertical profiles along each orbit. However, CloudSat reflectivity signatures are complicated somewhat by resonant Mie scattering effects and significant attenuation in the presence of cloud or rain water. Here, both radar systems are applied to a case of light to moderate snowfall within the warm frontal zone of a cold season, synoptic scale storm. Radars allow for an evaluation of the accuracy of a single-moment scheme in replicating precipitation structures, based on the bulk statistical properties of precipitation as suggested by reflectivity signatures.

Molthan, Andrew↗

Rationalizing Burned Carbon with Carbon Monoxide Exported from South America

We present several estimates cross-checking the fluxes of carbon to the atmosphere from burning, comparing models that are based on simple land-surface parameterizations and atmospheric transport dynamics. Both estimates made by NASA Ames and USP modeling techniques are quite high compared to some detailed satellite/land-use studies of emissions. The flux of carbon liberated to the atmosphere via biomass burning is important for several reasons. This flux is a fundamental statistic for the parameterization of the large-scale flux of gases controlling the reactive greenhouse gases methane and ozone. Similarly, it is central to the estimation of the translocation of nitrogen and pyrodenitrification in the tropics. Thirdly, CO2 emitted from rainforest clearing contributes directly to carbon lost from the rainforest system as it contributes to greenhouse gas forcing. While CO2 from pasturage, agriculture, etc, is considered to be reabsorbed seasonally, and so "off budget" for the carbon cycle, it must also be accounted. CO2 anomalies related to daily weather and interannual climatic variation are strong enough to perturb our scientific perception of long-term carbon storage trends. We compare fluxes deduced from land-use statistics (originally, W.M. Hao) and from satellite hot pixels (A. Setzer) with atmospheric fluxes determined by the mesoscale/continental scale models RAMS and MM5, and point to some new work with highly resolved global models (the NASA Data Assimilation Office's GEOS4). Our simulations are tied to events, so that measured tracers like CO tie the models directly to the burning and meteorology of a specific period. We point out a particular sensitivity in estimates based on CO, and indicate how analysis of CO2 along with other biomass-burning tracers may lead to an improved multi-species estimator of carbon burned.

Chatfield, R.↗

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year−1 (1Tg = 1012 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH4 year−1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of −34 to −46 Tg CH4 year−1, we estimated that the global net land methane emissions are 149–176 Tg CH4 year−1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

wetland methane emission↗

Climate Induced Spillover and Implications for U.S. Security.

Developing nations incur a greater risk to climate change than the developed world due to poorly managed human/natural resources, unreliable infrastructure and brittle governing/economic institutions. These vulnerabilities often give rise to a climate induced “domino effect” of reduced natural resource production-leading to economic hardship, social unrest, and humanitarian crises. Integral to this cascading set of events is increased human migration, leading to the “spillover” of impacts to adjoining areas with even broader impact on global markets and security. Given the complexity of factors influencing human migration and the resultant spill-over effect, quantitative tools are needed to aid policy analysis. Toward this need, a series of migration models were developed along with a system dynamics model of the spillover effect. The migration decision models were structured according to two interacting paths, one that captured long-term “chronic” impacts related to protracted deteriorating quality of life and a second focused on short-term “acute” impacts of disaster and/or conflict. Chronic migration dynamics were modeled for two different cases; one that looked only at emigration but at a national level for the entire world; and a second that looked at both emigration and immigration but focused on a single nation. Model parameterization for each of the migration models was accomplished through regression analysis using decadal data spanning the period 1960-2010. A similar approach was taken with acute migration dynamics except regression analysis utilized annual data sets limited to a shorter time horizon (2001-2013). The system dynamics spillover model was organized around two broad modules, one simulating the decision dynamics of migration and a second module that treats the changing environmental conditions that influence the migration decision. The environmental module informs the migration decision, endogenously simulating interactions/changes in the economy, labor, population, conflict, water, and food. A regional model focused on Mali in western Africa was used as a test case to demonstrate the efficacy of the model.

Tidwell, Vincent C.↗

The persistent shadow of the supermassive black hole of M87. II. Model comparisons and theoretical interpretations

The Event Horizon Telescope (EHT) observation of M87∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Spinning Spacecraft Attitude Estimation Using Markley Variables: Filter Implementation And Results

Attitude estimation is often more difficult for spinning spacecraft than for three-axis stabilized platforms due to the need to follow rapidly-varying state vector elements and the lack of three-axis rate measurements from gyros. The estimation problem simplifies when torques are negligible and nutation has damped out, but the general case requires a sequential filter with dynamics propagation. This paper describes the implementation and test results for an extended Kalman filter for spinning spacecraft attitude and rate estimation based on a novel set of variables suggested in a paper by Markley [AAS93-3301 (referred to hereafter as Markley variables). Markley has demonstrated that the new set of variables provides a superior parameterization for numerical integration of the attitude dynamics for spinning or momentum-biased spacecraft. The advantage is that the Markley variables have fewer rapidly-varying elements than other representations such as the attitude quaternion and rate vector. A filter based on these variables was expected to show improved performance due to the more accurate numerical state propagation. However, for a variety of test cases, it has been found that the new filter, as currently implemented, does not perform significantly better than a quaternion-based filter that was developed and tested in parallel. This paper reviews the mathematical background for a filter based on Markley variables. It also describes some features of the implementation and presents test results. The test cases are based on a mission using magnetometer and Sun sensor data and gyro measurements on two axes normal to the spin axis. The orbit and attitude scenarios and spacecraft parameters are modeled after one of the THEMIS (Time History of Events and Macroscale Interactions during Substorms) probes. Several tests are presented that demonstrate the filter accuracy and convergence properties. The tests include torque-free motion with various nutation angles, large constant-torque attitude slews, sensor misalignments, large initial attitude and rate errors, and cases with low data frequency. It is found that the convergence is rapid, the radius of convergence is large, and the results are reasonably accurate even in the presence of unmodeled perturbations.

Sedlak, Joseph E.↗

Ozone Production by Corona Discharges During a Convective Event in DISCOVER-AQ Houston

An ozonesonde launched near electrically active convection in Houston, TX on 5 September 2013 during the NASA DISCOVER-AQ project measured a large enhancement of ozone throughout the troposphere. A separate ozonesonde was launched from Smith Point, TX (approx. 58 km southeast of the Houston site) at approximately the same time as the launch from Houston and did not measure that enhancement. Furthermore, ozone profiles for the descent of both sondes agreed well with the ascending Smith Point profile, suggesting a highly localized event in both space and time in which an anomalously large enhancement of 70 - 100 ppbv appeared in the ascending Houston ozonesonde data. Compared to literature values, such an enhancement appears to be the largest observed to date. Potential sources of the localized ozone enhancement such as entrainment of urban or biomass burning emissions, downward transport from the stratosphere, photochemical production from lightning NO(sub x), and direct ozone production from corona discharges were investigated using model simulations. We conclude that the most likely explanation for the large ozone enhancement is direct ozone production by corona discharges. Integrating the enhancement seen in the Houston ozone profile and using the number of electrical discharges detected by the NLDN (or HLMA), we estimate a production of 2.48 x 10(exp. 28) molecules of ozone per flash which falls within the range of previously recorded values (9.89 x 10(exp. 26) - 9.82 x 10)exp. 28) molecules of ozone per flash). Since there is currently no parameterization for the direct production of ozone from corona discharges we propose the implementation of an equation into a chemical transport model. Ultimately, additional work is needed to further understand the occurrence and impact of corona discharges on tropospheric chemistry on short and long timescales.

ozone↗

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year -1 (1Tg = 10 12 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH 4 year -1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of -34 to -46 Tg CH 4 year -1 , we estimated that the global net land methane emissions are 149–176 Tg CH 4 year -1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

58 GEOSCIENCES↗

Birds Shed RNA-Viruses According to the Pareto Principle

A major challenge in disease ecology is to understand the role of individual variation of infection load on disease transmission dynamics and how this influences the evolution of resistance or tolerance mechanisms. Such information will improve our capacity to understand, predict, and mitigate pathogen-associated disease in all organisms. In many hostpathogen systems, particularly macroparasites and sexually transmitted diseases, it has been found that approximately 20% of the population is responsible for approximately 80% of the transmission events. Although host contact rates can account for some of this pattern, pathogen transmission dynamics also depend upon host infectiousness, an area that has received relatively little attention. Therefore, we conducted a meta-analysis of pathogen shedding rates of 24 host (avian) – pathogen (RNA-virus) studies, including 17 bird species and five important zoonotic viruses. We determined that viral count data followed the Weibull distribution, the mean Gini coefficient (an index of inequality) was 0.687 (0.036 SEM), and that 22.0% (0.90 SEM) of the birds shed 80% of the virus across all studies, suggesting an adherence of viral shedding counts to the Pareto Principle. The relative position of a bird in a distribution of viral counts was affected by factors extrinsic to the host, such as exposure to corticosterone and to a lesser extent reduced food availability, but not to intrinsic host factors including age, sex, and migratory status. These data provide a quantitative view of heterogeneous virus shedding in birds that may be used to better parameterize epidemiological models and understand transmission dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

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↗

Quantified Event Automata: Towards Expressive and Efficient Runtime Monitors

Runtime verification is the process of checking a property on a trace of events produced by the execution of a computational system. Runtime verification techniques have recently focused on parametric specifications where events take data values as parameters. These techniques exist on a spectrum inhabited by both efficient and expressive techniques. These characteristics are usually shown to be conflicting - in state-of-the-art solutions, efficiency is obtained at the cost of loss of expressiveness and vice-versa. To seek a solution to this conflict we explore a new point on the spectrum by defining an alternative runtime verification approach.We introduce a new formalism for concisely capturing expressive specifications with parameters. Our technique is more expressive than the currently most efficient techniques while at the same time allowing for optimizations.

state machines↗

Simulation of 1986 South China Sea Monsoon with a Regional Climate Model

A Regional Land-Atmosphere Climate Simulation System (RELACS) project is being developed at NASA Goddard Space Flight Center. One of the major goals of RELACS is to use a regional scale model with improved physical processes and in particular land-related processes, to understand the role of the land surface and its interaction with convection and radiation as well as the water/energy cycles in the IndoChina/South China Sea (SCS) region. The Penn State/NCAR MM5 atmospheric modeling system, a state of the art atmospheric numerical model designed to simulate regional weather and climate, has been successfully coupled to the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model. The original MM5 model (without PLACE) includes the option for either a simple slab soil model or a five-layer soil model (MRF) in which the soil moisture availability evolves over time. However, the MM5 soil models do not include the effects of vegetation, and thus important physical processes such as evapotranspiration and interception are precluded. The PLACE model incorporates vegetation type and has been shown in international comparisons to accurately predict evapotranspiration and runoff over a wide variety of land surfaces. The coupling of MM5 and PLACE creates a numerical modeling system with the potential to more realistically simulate atmosphere and land surface processes including land-sea interaction, regional circulations such as monsoons, and flash flood events. In addition, the Penn State/NCAR MM5 atmospheric modeling system has been: (1) coupled to the Goddard Ice Microphysical scheme; (2) coupled to a turbulent kinetic energy (TKE) scheme; (3) modified to ensure cloud budget balance; and (4) incorporated initialization with the Goddard EOS data sets at NASA/Goddard Laboratory for Atmospheres. The improved MM5 with two nested domains (60 and 20 km horizontal resolution) was used to simulate convective activity over IndoChina and the South China Sea, during the monsoon season, from May 6 to May 20, 1986. The model results captured several dominant observed features, such as twin cyclones, a depression system over the Bay of Bengal, strong south-westerly winds over IndoChina before and during the on-set of convection over the SCS, and a vortex over the SCS. Two additional MM5 runs with different land process models, Blackadar and MRF, were performed, and their results are compared to the run with PLACE. The preliminary results indicate that the MM5 results using PLACE and Blackadar are in very good agreement, but the results using MRF do not contain the south-westerly wind over IndoChina prior to the on-set of convection over the SCS.

Tao, W. -K.↗