Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global sensitivity analysis”

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 541 records · Page 30

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗

A double peak in the seasonality of California's photosynthesis as observed from space

Solar-induced chlorophyll fluorescence (SIF) has been shown to be a powerful proxy for photosynthesis and gross primary productivity (GPP). The recently launched TROPOspheric Monitoring Instrument (TROPOMI) features the required spectral resolution and signal-to-noise ratio to retrieve SIF from space. Here, we present a downscaling method to obtain 500 m spatial resolution SIF over California. We report daily values based on a 14 d window. TROPOMI SIF data show a strong correspondence with daily GPP estimates at AmeriFlux sites across multiple ecosystems in California. We find a linear relationship between SIF and GPP that is largely invariant across ecosystems with an intercept that is not significantly different from zero. Measurements of SIF from TROPOMI agree with MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation indices – the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation index (NIR$_v$) – at annual timescales but indicate different temporal dynamics at monthly and daily timescales. TROPOMI SIF data show a double peak in the seasonality of photosynthesis, a feature that is not present in the MODIS vegetation indices. The different seasonality in the vegetation indices may be due to a clear-sky bias in the vegetation indices, whereas previous work has shown SIF to have a low sensitivity to clouds and to detect the downregulation of photosynthesis even when plants appear green. We further decompose the spatiotemporal patterns in the SIF data based on land cover. The double peak in the seasonality of California's photosynthesis is due to two processes that are out of phase: grasses, chaparral, and oak savanna ecosystems show an April maximum, while evergreen forests peak in June. An empirical orthogonal function (EOF) analysis corroborates the phase offset and spatial patterns driving the double peak. The EOF analysis further indicates that two spatiotemporal patterns explain 84 % of the variability in the SIF data. Results shown here are promising for obtaining global GPP at sub-kilometer spatial scales and identifying the processes driving carbon uptake.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of Spaceborne and Ground Radar Comparison Results to Data Analysis Methods and Constraints

With the availability of active weather radar observations from space from the Precipitation Radar (PR) on board the Tropical Rainfall Measuring Mission (TR.MM) satellite, numerous studies have been performed comparing PR reflectivity and derived rain rates to similar observations from ground-based weather radars (GR). These studies have used a variety of algorithms to compute matching PR and GR volumes for comparison. Most studies have used a fixed 3-dimensional Cartesian grid centered on the ground radar, onto which the PR and GR data are interpolated using a proprietary approach and/or commonly available GR analysis software (e.g., SPRINT, REORDER). Other studies have focused on the intersection of the PR and GR viewing geometries either explicitly or using a hybrid of the fixed grid and PR/GR common fields of view. For the Dual-Frequency Precipitation Radar (DPR) of the upcoming Global Precipitation Measurement (GPM) mission, a prototype DPR/GR comparison algorithm based on similar TRMM PR data has been developed that defines the common volumes in terms of the geometric intersection of PR and GR rays, where smoothing of the PR and GR data are minimized and no interpolation is performed. The PR and GR volume-averaged reflectivity values of each sample volume are accompanied by descriptive metadata, for attributes including the variability and maximum of the reflectivity within the sample volume, and the fraction of range gates in the sample average having reflectivity values above an adjustable detection threshold (typically taken to be 18 dBZ for the PR). Sample volumes are further characterized by rain type (Stratiform or Convective), proximity to the melting layer, underlying surface (land/water/mixed), and the time difference between the PR and GR observations. The mean reflectivity differences between the PR and GR can differ between data sets produced by the different analysis methods; and for the GPM prototype, by the type of constraints and categorization applied to the data. In this paper, we will show results comparing the 3-D gridded analysis "black box" approach to the GPM prototype volume-matching approach, using matching TRMM PR and WSR-88D ground radar data. The affects of applying data constraints and data categorizations on the volume-matched data to the results will be shown, and explanations of the differences in terms of data and analysis algorithm characteristics will be presented. Implications of the differences to the determination of PR/DPR calibration differences and use of ground radar data to evaluate the PR and DPR attenuation correction algorithms will be discussed.

Morris, Kenneth R.↗

Intraseasonal Variations in Tropical Energy Balance: Relevance to Climate Sensitivity?

Intraseasonal variability of deep convection represents a fundamental mode of organization for tropical convection. While most studies of intraseasonal oscillations (ISOs) have focused on the spatial propagation and dynamics of convectively coupled circulations, here we examine the projection of ISOs on the tropically-averaged heat and moisture budget. One unresolved question concerns the degree to which observable variations in the "fast" processes (e.g. convection, radiative / turbulent fluxes) can inform our understanding of feedback mechanisms operable in the context of climate change. Our analysis use daily data from satellite observations, the Modern Era analysis for Research and Applications (MERRA), and other model integrations to address these questions: (i) How are tropospheric temperature variations related to that tropical deep convection and the associated ice cloud fractional amount (ICF), ice water path (IWP), and properties of warmer liquid clouds? (ii) What role does moisture transport play vis-a-vis ocean latent heat flux in enabling the evolution of deep convection to mediate PBL - free atmospheric temperature equilibration? (iii) What affect do convectively generated upper-tropospheric clouds have on the TOA radiation budget? Our methodology is similar to that of Spencer et al., (2007 GRL ) whereby a composite time series of various quantities over 60+ ISO events is built using tropical mean tropospheric temperature signal as a reference to which the variables are related at various lag times (from -30 to +30 days). The area of interest encompasses the global oceans between 20oN/S. The increase of convective precipitation cannot be sustained by evaporation within the domain, implying strong moisture transports into the tropical ocean area. The decrease in net TOA radiation that develops after the peak in deep convective rainfall, is part of the response that constitutes a "discharge" / "recharge" mechanism that facilitates tropical heat balance maintenance on these time scales. However, water vapor and hydrologic scaling relationships for this mode of variability cast doubt on the utility of ISO variations as proxies for climate sensitivity response to external radiatively forced (e.g. greenhouse gas-induced) climate change.

Robertson, Franklin R.↗

Whole-soil warming leads to substantial soil carbon emission in an alpine grassland

The sensitivity of soil organic carbon (SOC) decomposition in seasonally frozen soils, such as alpine ecosystems, to climate warming is a major uncertainty in global carbon cycling. Here we measure soil CO 2 emission during four years (2018–2021) from the whole-soil warming experiment (4 °C for the top 1 m) in an alpine grassland ecosystem. We find that whole-soil warming stimulates total and SOC-derived CO 2 efflux by 26% and 37%, respectively, but has a minor effect on root-derived CO 2 efflux. Moreover, experimental warming only promotes total soil CO 2 efflux by 7-8% on average in the meta-analysis across all grasslands or alpine grasslands globally (none of these experiments were whole-soil warming). We show that whole-soil warming has a much stronger effect on soil carbon emission in the alpine grassland ecosystem than what was reported in previous warming experiments, most of which only heat surface soils.

54 ENVIRONMENTAL SCIENCES↗

The LUVOIR Extreme Coronagraph for Living Planetary Systems (ECLIPS) II. Performance Evaluation, Aberration Sensitivity Analysis and Exoplanet Detection Simulations

Future space missions such as the Large UV-Optical-Infrared Surveyor (LUVOIR) and the Habitable Exoplanet Observatory (HabEx) require coronagraphs with active wavefront control to suppress starlight to discover and characterize habitable exoplanets. The Extreme Coronagraph for Living Planetary Systems (ECLIPS) is the coronagraph instrument on the LUVOIR Surveyor mission concept, an 8–15m segmented telescope. ECLIPS is split into three channels: UV (200 to 400 nm), optical (400 nm to 850 nm), and NIR (850 nm to 2.0 microns), with each channel equipped with two deformable mirrors for wavefront control, a suite of coronagraph masks, a low-order/out-of-band wavefront sensor, and separate science imagers and spectrographs. The Apodized Pupil Lyot Coronagraph (APLC) and the Vector Vortex Coronagraph (VVC) are the baselined mask technologies for ECLIPS to enable the required 10 -10 contrast for observations in the habitable zones of nearby stars. Their performance depends on active wavefront sensing and control, as well as metrology subsystems to compensate for aberrations induced by segment errors (piston and tip/tilt, among others), secondary mirror misalignment, and global low-order wavefront errors. Here we present the latest results of the simulation of these effects for the two technologies and discuss the achieved contrast for exoplanet detection and characterization after closed-loop wavefront estimation and control algorithms have been applied. Finally, we show simulated observations using high-fidelity spatial and spectral input models of complete planetary systems generated with the Haystacks code framework.

Roser Juanola-Parramon↗

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗

Atomic oxygen in the Martian thermosphere

The Mariner 9 Ultraviolet Spectrometer (UVS) made extensive observations of air-glow emissions from the thermosphere of Mars throughout the nominal mission (November 1971 - February 1972), during late summer in the southern hemisphere. Limb and disc measurements of the 130 nm triplet emission from thermospheric atomic oxygen were modelled by Strickland et al. Recently, the thermospheric general circulation models (TGCMs) developed for the Earth and Venus have been applied to Mars; we refer to it as the MTGCM. Our analysis shows that the oxygen mixing ratio is the fundamental unknown controlling the 130 nm brightness. Our radiative transport calculation shows that the emergent intensity at 130 nm is not very sensitive to variations in thermospheric temperature. The pattern of diurnal variation derived from our analysis is roughly the same as Strickland et al. although with somewhat lower values for the O mixing ratio. The main reasons for this difference are the more important role played by the photoelectron source in our model, and the somewhat larger 130 nm solar flux; thus, we require less oxygen to match the observed brightnesses. Strickland et al. also found that the OI 130 nm emission on Mars is correlated with solar activity. We find that the correlation is virtually non-existent during the early orbits when the planet was covered with a thick global dust storm, but later orbits, during the clearing of the storm, show a persistent correlation.

Stewart, A. I. F.↗

Quantifying Fractional Ground Cover on the Climate Sensitive High Plains Using AVIRIS and Landsat TM Data

The High Plains is an economically important and climatologically sensitive region of the United States and Canada. The High Plains contain 100,000 sq km of Holocene sand dunes and sand sheets that are currently stabilized by natural vegetation. Droughts and the larger threat of global warming are climate phenomena that could cause depletion of natural vegetation and make this region susceptible to sand dune reactivation. This thesis is part of a larger study that is assessing the effect of climate variability on the natural vegetation that covers the High Plains using Landsat 5 and Landsat 7 data. The question this thesis addresses is how can fractional vegetation cover be mapped with the Landsat instruments using linear spectral mixture analysis and to what accuracy. The method discussed in this thesis made use of a high spatial and spectral resolution sensor called AVIRIS (Airborne Visible and Infrared Imaging Spectrometer) and field measurements to test vegetation mapping in three Landsat 7 sub-scenes. Near-simultaneous AVIRIS images near Ft. Morgan, Colorado and near Logan, New Mexico were acquired on July 10, 1999 and September 30, 1999, respectively. The AVIRIS flights preceded Landsat 7 overpasses by approximately one hour. These data provided the opportunity to test spectral mixture algorithms with AVIRIS and to use these data to constrain the multispectral mixed pixels of Landsat 7. The comparisons of mixture analysis between the two instruments showed that AVIRIS endmembers can be used to unmix Landsat 7 data with good estimates of soil cover, and reasonable estimates of non-photosynthetic vegetation and green vegetation. Landsat 7 derived image endmembers correlate with AVIRIS fractions, but the error is relatively large and does not give a precise estimate of cover.

Warner, Amanda Susan↗

Using Ensemble Short-Term Initialized Coupled NASA GEOS5 Climate Model Integrations to Study Convective Bias Growth

The representation of convective processes, particularly deep convection in the tropics, remains a persistent problem in climate models. In fact structural biases in the distribution of tropical rainfall in the CMIP5 models is hardly different than that of the CMIP3 versions. Given that regional climate change at higher latitudes is sensitive to the configuration of tropical forcing, this persistent bias is a major issue for the credibility of climate change projections. In this study we use model output from integrations of the NASA Global Earth Observing System Five (GEOS5) climate modeling system to study the evolution of biases in the location and intensity of convective processes. We take advantage of a series of hindcast experiments done in support of the US North American Multi-Model Ensemble (NMME) initiative. For these experiments a nine-month forecast using a coupled model configuration is made approximately every five days over the past 30 years. Each forecast is started with an updated analysis of the ocean, atmosphere and land states. For a given calendar month we have approximately 180 forecasts with daily means of various quantities. These forecasts can be averaged to essentially remove "weather scales" and highlight systematic errors as they evolve. Our primary question is to ask how the spatial structure of daily mean precipitation over the tropics evolves from the initial state and what physical processes are involved. Errors in parameterized convection, various water and energy fluxes and the divergent circulation are found to set up on fast time scales (order five days) compared to errors in the ocean, although SST changes can be non-negligible over that time. For the month of June the difference between forecast day five versus day zero precipitation looks quite similar to the difference between the June precipitation climatology and that from the Global Precipitation Climatology Project (GPCP). We focus much of our analysis on the influence of SST gradients, associated PBL baroclinicity enabled by turbulent mixing, the ensuing PBL moisture convergence, and how changes in these processes relate to convective precipitation bias growth over this short period.

Cohen, Charlie↗

First simultaneous global QCD analysis of dihadron fragmentation functions and transversity parton distribution functions

We perform a comprehensive study within quantum chromodynamics (QCD) of dihadron observables in electron-positron annihilation, semi-inclusive deep-inelastic scattering, and proton-proton collisions, including recent cross section data from Belle and azimuthal asymmetries from STAR. We extract simultaneously for the first time π + π − dihadron fragmentation functions (DiFFs) and the nucleon transversity distributions for up and down quarks as well as antiquarks. For the transversity distributions we impose their small- x asymptotic behavior and the Soffer bound. In addition, we utilize a new definition of DiFFs that has a number density interpretation to then calculate expectation values for the dihadron invariant mass and momentum fraction. Furthermore, we investigate the compatibility of our transversity results with those from single-hadron fragmentation (from a transverse momentum dependent/collinear twist-3 framework) and the nucleon tensor charges computed in lattice QCD. We find a universal nature to all of this available information. Future measurements of dihadron production can significantly further this research, especially, as we show, those that are sensitive to the region of large parton momentum fractions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 1: Aft Skirt analysis

Using the ANSYS finite element program, a global model of the aft skirt and a detailed nonlinear model of the failure region was made. The analysis confirmed the area of failure in both STA-2B and STA-3 tests as the forging heat affected zone (HAZ) at the aft ring centerline. The highest hoop strain in the HAZ occurs in this area. However, the analysis does not predict failure as defined by ultimate elongation of the material equal to 3.5 percent total strain. The analysis correlates well with the strain gage data from both the Wyle influence test of the original design aft sjirt and the STA-3 test of the redesigned aft skirt. it is suggested that the sensitivity of the failure area material strength and stress/strain state to material properties and therefore to small manufacturing or processing variables is the most likely cause of failure below the expected material ultimate properties.

Berry, David M.↗

MinMap: An imaging spectrometer for high resolution compositional mapping of the Moon

MinMap has been selected by the Lunar Scout program to characterize and map the mineral composition of the Moon. The instrument will be built as a collaborative effort between Brown University, SETS Technology Inc., and Ball Aerospace Corp. MinMap is a visible to near-infrared imaging spectrometer that contains 192 spectral channels from 0.35-2.4 microns with signal to noise greater than 200 and 256 cross-track spatial elements. The spectrometer design has a 6 deg field of view (FOV) and utilizes grating dispersive elements and two dimensional detectors (no moving parts). An 'image cube' of data is produced that contains two dimensions of spatial information and one dimension of spectral information. All spectral channels and cross-track spatial elements are recorded simultaneously with spacecraft motion scanning the second spatial dimension. The high spectral resolution and continuous spectral range of MinMap are designed to measure the diagnostic absorption features of principal lunar minerals and their lithologic mixtures. Since the optical properties of lunar materials change in a regular manner upon exposure to the space environment, this spectral range is also quite sensitive to variations in exposure history (soil maturity). Nominal measurement stragegy is to obtain full global data of the Moon at 180 m/pixel from a 450 km polar orbit during the first month or two of operation. A 100 km orbit is anticipated for the remaining part of a 1 year mission allowing higher resolution data (approx. 80 m/pixel) to be obtained for targeted regions. MinMap exceeds LE x SWG's measurement recommendations and will provide the highest spatial resolution compositional map of lunar rocks and soils currently planned for orbital missions. Since all spectral channels are co-registered and obtained simultaneously, 'image cube' data swaths will be available for analysis almost immediately.

Pieters, C. M.↗

Climate Model Evaluation using New Datasets from the Clouds and the Earth's Radiant Energy System (CERES)

There are some in the science community who believe that the response of the climate system to anthropogenic radiative forcing is unpredictable and we should therefore call off the quest . The key limitation in climate predictability is associated with cloud feedback. Narrowing the uncertainty in cloud feedback (and therefore climate sensitivity) requires optimal use of the best available observations to evaluate and improve climate model processes and constrain climate model simulations over longer time scales. The Clouds and the Earth s Radiant Energy System (CERES) is a satellite-based program that provides global cloud, aerosol and radiative flux observations for improving our understanding of cloud-aerosol-radiation feedbacks in the Earth s climate system. CERES is the successor to the Earth Radiation Budget Experiment (ERBE), which has widely been used to evaluate climate models both at short time scales (e.g., process studies) and at decadal time scales. A CERES instrument flew on the TRMM satellite and captured the dramatic 1998 El Nino, and four other CERES instruments are currently flying aboard the Terra and Aqua platforms. Plans are underway to fly the remaining copy of CERES on the upcoming NPP spacecraft (mid-2010 launch date). Every aspect of CERES represents a significant improvement over ERBE. While both CERES and ERBE measure broadband radiation, CERES calibration is a factor of 2 better than ERBE. In order to improve the characterization of clouds and aerosols within a CERES footprint, we use coincident higher-resolution imager observations (VIRS, MODIS or VIIRS) to provide a consistent cloud-aerosol-radiation dataset at climate accuracy. Improved radiative fluxes are obtained by using new CERES-derived Angular Distribution Models (ADMs) for converting measured radiances to fluxes. CERES radiative fluxes are a factor of 2 more accurate than ERBE overall, but the improvement by cloud type and at high latitudes can be as high as a factor of 5. Diurnal cycles are explicitly resolved by merging geostationary satellite observations with CERES and MODIS. Atmospheric state data are provided from a frozen version of the Global Modeling and Assimilation Office- Data Assimilation System at the NASA Goddard Space Flight Center. In addition to improving the accuracy of top-of-atmosphere (TOA) radiative fluxes, CERES also produces radiative fluxes at the surface and at several levels in the atmosphere using radiative transfer modeling, constrained at the TOA by CERES (ERBE was limited to the TOA). In all, CERES uses 11 instruments on 7 spacecraft all integrated to obtain climate accuracy in TOA to surface fluxes. This presentation will provide an overview of several new CERES datasets of interest to the climate community (including a new adjusted TOA flux dataset constrained by estimates of heat storage in the Earth system), show direct comparisons between CERES ad ERBE, and provide a detailed error analysis of CERES fluxes at various time and space scales. We discuss how observations can be used to reduce uncertainties in cloud feedback and climate sensitivity and strongly argue why we should NOT "call off the quest".

Loeb, Norman G.↗

Extreme Lake-Effect Snow from a GPM Microwave Imager Perspective: Observational Analysis and Precipitation Retrieval Evaluation

This study focuses on the ability of the Global Precipitation Measurement (GPM) passive microwave sensors to detect and provide quantitative precipitation estimates (QPE) for extreme lake-effect snowfall events over the United States lower Great Lakes region. GPM Microwave Imager (GMI) high frequency channels can clearly detect intense shallow convective snowfall events. However, GMI Goddard PROfiling (GPROF) QPE retrievals produce inconsistent results when compared against the Multi-Radar/Multi-Sensor (MRMS) ground-based radar reference dataset. While GPROF retrievals adequately capture intense snowfall rates and spatial patterns of one event, GPROF systematically underestimates intense snowfall rates in another event. Furthermore, GPROF produces abundant light snowfall rates that do not conform with MRMS observations. Ad-hoc precipitation rate thresholds are suggested to partially mitigate GPROF’s overproduction of light snowfall rates. The sensitivity and retrieval efficiency of GPROF to key parameters (2-meter temperature, total precipitable water, and background surface type) used to constrain the GPROF a-priori retrieval database are investigated. Results demonstrate that typical lake-effect snow environmental and surface conditions, especially coastal surfaces, are underpopulated in the database and adversely affect GPROF retrievals. For the two presented case studies, using snow cover a-priori database in the locations of originally deemed as coastline improves retrieval. This study suggests that it is particularly important to have more accurate GPROF surface classifications and better representativeness of the a-priori databases to improve intense lake-effect snow detection and retrieval performance.

Lisa Milani↗

Exploring drought-responsive crucial genes in Sorghum

Drought severely affects global food production. Sorghum is a typical drought-resistant model crop. Based on RNA-seq data for Sorghum with multiple time points and the gray correlation coefficient, this paper firstly selects candidate genes via mean variance test and constructs weighted gene differential co-expression networks (WGDCNs); then, based on guilt-by-rewiring principle, the WGDCNs and the hidden Markov random field model, drought-responsive crucial genes are identified for five developmental stages respectively. Enrichment and sequence alignment analysis reveal that the screened genes may play critical functional roles in drought responsiveness. A multilayer differential co-expression network for the screened genes reveals that Sorghum is very sensitive to pre-flowering drought. Furthermore, a crucial gene regulatory module is established, which regulates drought responsiveness via plant hormone signal transduction, MAPK cascades, and transcriptional regulations. The proposed method can well excavate crucial genes through RNA-seq data, which have implications in breeding of new varieties with improved drought tolerance.

60 APPLIED LIFE SCIENCES↗

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear global stability analysis of compressor stall phenomena

Compressor stall phenomena are analyzed from the point of view of nonlinear control theory, based on bifurcation-catastrophe techniques. This new approach appears promising and offers insight into such well-known compressor instability problems as surge and rotating stall and suggests strategies for recovery. Three interlocking dynamic nonlinear state space models are developed. It is shown that the problem of rotating stall can be viewed as an induced bifurcation of solution of the unstalled model. Hysteresis effects are shown to exist in the stall/recovery process. Surge cycles are observed for some critical parameter values. The oscillatory behavior is seen to be due to development of limit cycles, generated by Hopf bifurcation of solutions. More specifically, it is observed that at certain critical values of parameters, a family of stable limit cycles with growning and then diminishing amplitudes is generated, then giving rise to an unstable family of limit cycles. This unstable family in turn bifurcates into other unstable families. To further illustrate the utility of the methodology, some partial computation of domains is carried out, and parameter sensitivity analysis is performed.

Razavi, H.↗