Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wind direction variability”

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

Dynamic Testing of the NASA Hypersonic Project Combined Cycle Engine Testbed for Mode Transition Experiments

NASA is interested in developing technology that leads to more routine, safe, and affordable access to space. Access to space using airbreathing propulsion systems has potential to meet these objectives based on Airbreathing Access to Space (AAS) system studies. To this end, the NASA Fundamental Aeronautics Program (FAP) Hypersonic Project is conducting fundamental research on a Turbine Based Combined Cycle (TBCC) propulsion system. The TBCC being studied considers a dual flow-path inlet system. One flow-path includes variable geometry to regulate airflow to a turbine engine cycle. The turbine cycle provides propulsion from take-off to supersonic flight. The second flow-path supports a dual-mode scramjet (DMSJ) cycle which would be initiated at supersonic speed to further accelerate the vehicle to hypersonic speed. For a TBCC propulsion system to accelerate a vehicle from supersonic to hypersonic speed, a critical enabling technology is the ability to safely and effectively transition from the turbine to the DMSJ-referred to as mode transition. To experimentally test methods of mode transition, a Combined Cycle Engine (CCE) Large-scale Inlet testbed was designed with two flow paths-a low speed flow-path sized for a turbine cycle and a high speed flow-path designed for a DMSJ. This testbed system is identified as the CCE Large-Scale Inlet for Mode Transition studies (CCE-LIMX). The test plan for the CCE-LIMX in the NASA Glenn Research Center (GRC) 10- by 10-ft Supersonic Wind Tunnel (10x10 SWT) is segmented into multiple phases. The first phase is a matrix of inlet characterization (IC) tests to evaluate the inlet performance and establish the mode transition schedule. The second phase is a matrix of dynamic system identification (SysID) experiments designed to support closed-loop control development at mode transition schedule operating points for the CCE-LIMX. The third phase includes a direct demonstration of controlled mode transition using a closed loop control system developed with the data obtained from the first two phases. Plans for a fourth phase include mode transition experiments with a turbine engine. This paper, focusing on the first two phases of experiments, presents developed operational and analysis tools for streamlined testing and data reduction procedures.

Source record↗

Reconstructed Historical Land Cover and Biophysical Parameters for Studies of Land-Atmosphere Interactions within the Eastern United States

The local environment where we live within the Earth's biosphere is often taken for granted. This environment can vary depending on whether the land cover is a forest, grassland, wetland, water body, bare soil, pastureland, agricultural field, village, residential suburb, or an urban complex with concrete, asphalt, and large buildings. In general, the type and characteristics of land cover influence surface temperatures, sunlight exposure and duration, relative humidity, wind speed and direction, soil moisture amount, plant life, birds, and other wildlife in our backyards. The physical and biological properties (biophysical characteristics) of land cover help to determine our surface environment because they directly affect surface radiation, heat, and soil moisture processes, and also feedback to regional weather and climate. Depending on the spatial scale and land use intensity, land cover changes can have profound impacts on our local and regional environment. Over the past 350 years, the eastern half of the United States, an area extending from the grassland prairies of the Great Plains to the Gulf and Atlantic coasts, has experienced extensive land cover and land use changes that began with land clearing in the 1600s, led to extensive deforestation and intensive land use practices by 1920, and then evolved to the present-day landscape. Determining the consequences of such land cover changes on regional and global climate is a major research issue. Such research requires detailed historical land cover data and modeling experiments simulating historical climates. Given the need to understand the effects of historical land cover changes in the eastern United States, some questions include: - What were the most important land cover transformations and how did they alter biophysical characteristics of the land cover at key points in time since the mid-1600s? - How have land cover and land use changes over the past 350 years affected the land surface environment including surface weather, hydrologic, and climatic variability? - How do the potential effects of regional human-induced land cover change on the environment compare to similar changes that are caused by the natural variations of the Earth's climate system? To help answer these questions, we reconstructed a fractional land cover and biophysical parameter dataset for the eastern United States at 1650, 1850, 1920, and 1992 time-slices. Each land cover fraction is associated with a biophysical parameter class, a suite of parameters defining the biophysical characteristics of that kind of land cover. This new dataset is designed for use in computer models of land-atmosphere interactions, to understand and quantify the effects of historical land cover changes on the water, energy, and carbon cycles

Steyaert, Louis T.↗

Effects of spatial resolution on WRF v3.8.1 simulated meteorology over the central Himalaya

The sensitive ecosystem of the central Himalayan (CH) region, which is experiencing enhanced stress from anthropogenic forcing, requires adequate atmospheric observations and an improved representation of the Himalaya in the models. However, the accuracy of atmospheric models remains limited in this region due to highly complex mountainous topography. This article delineates the effects of spatial resolution on the modeled meteorology and dynamics over the CH by utilizing the Weather Research and Forecasting (WRF) model extensively evaluated against the Ganges Valley Aerosol Experiment (GVAX) observations during the summer monsoon. The WRF simulation is performed over a domain (d01) encompassing northern India at 15 km x 15 km resolution and two nests (d02 at 5 km x 5 km and d03 at 1 km x 1 km) centered over the CH, with boundary conditions from the respective parent domains. WRF simulations reveal higher variability in meteorology, e.g., relative humidity (RH = 70.3%-96.1%) and wind speed (WS = 1.1-4.2m s(-1)), compared to the ERA-Interim reanalysis (RH D 80.0%-85.0%, WS D 1.2-2.3m s(-1)) over northern India owing to the higher resolution. WRF-simulated temporal evolution of meteorological variables is found to agree with balloon-borne measurements, with stronger correlations aloft (r = 0.44-0.92) than those in the lower troposphere (r = 0.18-0.48). The model overestimates temperature (warm bias by 2.8 degrees C) and underestimates RH (dry bias by 6.4 %) at the surface in d01. Model results show a significant improvement in d03 (P = 827.6 hPa, T = 19.8 degrees C, RH = 92.3 %), closer to the GVAX observations (P = 801.4 hPa, T = 19.5 degrees C, RH = 94.7 %). Interpolating the output from the coarser domains (d01, d02) to the altitude of the station reduces the biases in pressure and temperature; however, it suppresses the diurnal variations, highlighting the importance of well-resolved terrain (d03). Temporal variations in near-surface P, T, and RH are also reproduced by WRF in d03 to an extent (r>0:5). A sensitivity simulation incorporating the feedback from the nested domain demonstrates the improvement in simulated P, T, and RH over the CH. Our study shows that the WRF model setup at finer spatial resolution can significantly reduce the biases in simulated meteorology, and such an improved representation of the CH can be adopted through domain feedback into regional-scale simulations. Interestingly, WRF simulates a dominant easterly wind component at 1 km x 1 km resolution (d03), which is missing in the coarse simulations; however, the frequency of southeasterlies remains underestimated. The model simulation implementing a highresolution (3 s) topography input (SRTM) improved the prediction of wind directions; nevertheless, further improvements are required to better reproduce the observed local-scale dynamics over the CH.

Singh, Jaydeep↗

The Solar Connections Observatory for Planetary Environments

The NASA Sun-Earth Connection theme roadmap calls for comparative study of how the planets, comets, and local interstellar medium (LISM) interact with the Sun and respond to solar variability. Through such a study we advance our understanding of basic physical plasma and gas dynamic processes, thus increasing our predictive capabilities for the terrestrial, planetary, and interplanetary environments where future remote and human exploration will occur. Because the other planets have lacked study initiatives comparable to the terrestrial ITM, LWS, and EOS programs, our understanding of the upper atmospheres and near space environments on these worlds is far less detailed than our knowledge of the Earth. To close this gap we propose a mission to study {\it all) of the solar interacting bodies in our planetary system out to the heliopause with a single remote sensing space observatory, the Solar Connections Observatory for Planetary Environments (SCOPE). SCOPE consists of a binocular EUV/FUV telescope operating from a remote, driftaway orbit that provides sub-arcsecond imaging and broadband medium resolution spectro-imaging over the 55-290 nm bandpass, and high (R>10$^{5}$ resolution H Ly-$\alpha$ emission line profile measurements of small scale planetary and wide field diffuse solar system structures. A key to the SCOPE approach is to include Earth as a primary science target. From its remote vantage point SCOPE will be able to observe auroral emission to and beyond the rotational pole. The other planets and comets will be monitored in long duration campaigns centered when possible on solar opposition when interleaved terrestrial-planet observations can be used to directly compare the response of both worlds to the same solar wind stream and UV radiation field. Using a combination of observations and MHD models, SCOPE will isolate the different controlling parameters in each planet system and gain insight into the underlying physical processes that define the solar connection.

Oliversen, Ronald J.↗

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

Multi-Sensor Cloud and Aerosol Retrieval Simulator and Remote Sensing from Model Parameters : Aerosols - Part 2

The Multi-sensor Cloud Retrieval Simulator (MCRS) produces a simulated radiance product from any high-resolution general circulation model with interactive aerosol as if a specific sensor such as the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of the atmospheric column and land ocean surface at a specific location. Previously the MCRS code only included contributions from atmosphere and clouds in its radiance calculations and did not incorporate properties of aerosols. In this paper we added a new aerosol properties module to the MCRS code that allows users to insert a mixture of up to 15 different aerosol species in any of 36 vertical layers. This new MCRS code is now known as MCARS (Multi-sensor Cloud and Aerosol Retrieval Simulator). Inclusion of an aerosol module into MCARS not only allows for extensive, tightly controlled testing of various aspects of satellite operational cloud and aerosol properties retrieval algorithms, but also provides a platform for comparing cloud and aerosol models against satellite measurements. This kind of two-way platform can improve the efficacy of model parameterizations of measured satellite radiances, allowing the assessment of model skill consistently with the retrieval algorithm. The MCARS code provides dynamic controls for appearance of cloud and aerosol layers. Thereby detailed quantitative studies of the impacts of various atmospheric components can be controlled. In this paper we illustrate the operation of MCARS by deriving simulated radiances from various data field output by the Goddard Earth Observing System version 5 (GEOS-5) model. The model aerosol fields are prepared for translation to simulated radiance using the same model sub grid variability parameterizations as are used for cloud and atmospheric properties profiles, namely the ICA technique. After MCARS computes modeled sensor radiances equivalent to their observed counterparts, these radiances are presented as input to operational remote-sensing algorithms. Specifically, the MCARS-computed radiances are input into the processing chain used to produce the MODIS Data Collection 6 aerosol product (MOYD04). TheMOYD04 product is of course normally produced from MOYD021KM MODIS Level-1B radiance product directly acquired by the MODIS instrument. MCARS matches the format and metadata of a MOYD021KM product. The resulting MCARS output can be directly provided to MODAPS (MODIS Adaptive Processing System) as input to various operational atmospheric retrieval algorithms. Thus the operational algorithms can be tested directly without needing to make any software changes to accommodate an alternative input source. We show direct application of this synthetic product in analysis of the performance of the MOD04 operational algorithm. We use biomass-burning case studies over Amazonia employed in a recent Working Group on Numerical Experimentation (WGNE)-sponsored study of aerosol impacts on numerical weather prediction (Freitas et al., 2015). We demonstrate that a known low bias in retrieved MODIS aerosol optical depth appears to be due to a disconnect between actual column relative humidity and the value assumed by the MODIS aerosol product.

aerosol retrieval↗

Electric Vehicle Managed Charging: Forward-Looking Estimates of Bulk Power System Value

When and where electric vehicle charging occurs has significant implications for power systems supporting widespread electric vehicle deployment with high shares of wind and solar generation. Numerous studies have estimated the value of scheduling or otherwise managing electric vehicle charging in such power systems. This study improves on those earlier works by leveraging detailed simulation models for electric vehicle adoption, electric vehicle use, electric vehicle charging, and bulk power system operations; and linking them with methods for describing charging flexibility at both the individual vehicle and aggregate levels. This study closely analyzes electric vehicle managed charging (EVMC) performance along the dimensions of flexibility type (within-charging session or within-week scheduling), dispatch mechanism (direct load control or one of several price-based mechanisms), and participation rate, under the assumptions of ubiquitous chargers and all trips completed on time. The study is located in a passenger light-duty vehicle adoption scenario with 100% electric vehicle sales by 2035, and in an envisioned 2038 New England power system for which within-region generation is 84% clean.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

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↗

VLTI-AMBER Velocity-Resolved Aperture-Synthesis Imaging of Eta Carinae with a Spectral Resolution of 12 000: Studies of the Primary Star Wind and Innermost Wind-Wind Collision Zone

The mass loss from massive stars is not understood well. Eta Carinae is a unique object for studying the massive stellar wind during the luminous blue variable phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Aims. We want to investigate the structure and kinematics of Car's primary star wind and wind-wind collision zone with a high spatial resolution of approx.6 mas (approx.14 au) and high spectral resolution of R = 12 000. Methods. Observations of Car were carried out with the ESO Very Large Telescope Interferometer (VLTI) and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 di erent spectral channels distributed across the Br(gamma) 2.166 micron emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is 126. The fan-shaped structure extends approximately 8.0 mas (approx.18:8 au) to the southeast and 5.8 mas (approx.13:6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

Velocity-resolved aperture-synthesis images↗

Market optimization and technoeconomic analysis of hydrogen-electricity coproduction systems

Decarbonization efforts across North America, Europe, and beyond rely on variable renewable energy sources such as wind and solar, as well as alternative fuels, such as hydrogen, to support the sustainable energy transition. These advancements have prompted a need for more flexibility in the electric grid to complement non-dispatchable energy sources and increased demand from electrification. Integrated energy systems are well suited to provide this flexibility, but conventional technoeconomic modeling paradigms neglect the time-varying dynamic nature of the grid and thus undervalue resource flexibility. In this work, we develop a computational optimization framework for dynamic market-based technoeconomic comparison of integrated energy systems that coproduce low-carbon electricity and hydrogen (e.g., solid oxide fuel cells, solid oxide electrolysis) against technologies that only produce electricity (e.g., natural gas combined cycle with carbon capture) or only produce hydrogen. Our framework starts with rigorous physics-based process models, built in the open-source Institute for the Design of Advanced Energy Systems (IDAES) modeling and optimization platform, for six energy process concepts. Using these rigorous models and a workflow to optimally design each technology, the framework is shown to be capable of evaluating new and emerging technologies in varying energy markets under a plethora of future scenarios (i.e., renewables penetration, carbon tax, etc.). Ultimately, our framework finds that solid oxide fuel cell-based coproduction systems achieve positive profits for 85% of the analyzed market scenarios. From these market optimization results, we use multivariate linear regression (R 2 values up to 0.99) to determine which electricity price statistics are most significant to predict the optimized annual profit of each system. The proposed framework provides a powerful tool for directly comparing flexible, multi-product energy process concepts to help discern optimal technology and integration options.

08 HYDROGEN↗

Southern Ocean Climate and Sea Ice Anomalies Associated with the Southern Oscillation

The anomalies in the climate and sea ice cover of the Southern Ocean and their relationships with the Southern Oscillation (SO) are investigated using a 17-year of data set from 1982 through 1998. We correlate the polar climate anomalies with the Southern Oscillation index (SOI) and examine the composites of these anomalies under the positive (SOI > 0), neutral (0 > SOI > -1), and negative (SOI < -1) phases of SOL The climate data set consists of sea-level pressure, wind, surface air temperature, and sea surface temperature fields, while the sea ice data set describes its extent, concentration, motion, and surface temperature. The analysis depicts, for the first time, the spatial variability in the relationship of the above variables and the SOL The strongest correlation between the SOI and the polar climate anomalies are found in the Bellingshausen, Amundsen and Ross sea sectors. The composite fields reveal anomalies that are organized in distinct large-scale spatial patterns with opposing polarities at the two extremes of SOI, and suggest oscillating climate anomalies that are closely linked to the SO. Within these sectors, positive (negative) phases of the SOI are generally associated with lower (higher) sea-level pressure, cooler (warmer) surface air temperature, and cooler (warmer) sea surface temperature in these sectors. Associations between these climate anomalies and the behavior of the Antarctic sea ice cover are clearly evident. Recent anomalies in the sea ice cover that are apparently associated with the SOI include: the record decrease in the sea ice extent in the Bellingshausen Sea from mid- 1988 through early 199 1; the relationship between Ross Sea SST and ENSO signal, and reduced sea ice concentration in the Ross Sea; and, the shortening of the ice season in the eastern Ross Sea, Amundsen Sea, far western Weddell Sea, and the lengthening of the ice season in the western Ross Sea, Bellingshausen Sea and central Weddell Sea gyre over the period 1988-1994. Four ENSO episodes over the last 17 years contributed to a negative mean in the SOI (-0.5). In each of these episodes, significant retreats in the Bellingshausen/Amundsen Sea were observed providing direct confirmation of the impact of SO on the Antarctic sea ice cover.

Kwok, R.↗

Waterway Ice Thickness Measurements

The ship on the opposite page is a U. S. Steel Corporation tanker cruising through the ice-covered waters of the Great Lakes in the dead of winter. The ship's crew is able to navigate safely by plotting courses through open water or thin ice, a technique made possible by a multi-agency technology demonstration program in which NASA is a leading participant. Traditionally, the Great Lakes-St. Lawrence Seaway System is closed to shipping for more than three months of winter season because of ice blockage, particularly fluctuations in the thickness and location of ice cover due to storms, wind, currents and variable temperatures. Shippers have long sought a system of navigation that would allow year-round operation on the Lakes and produce enormous economic and fuel conservation benefits. Interrupted operations require that industrial firms stockpile materials to carry them through the impassable months, which is costly. Alternatively, they must haul cargos by more expensive overland transportation. Studies estimate the economic benefits of year-round Great Lakes shipping in the hundreds of millions of dollars annually and fuel consumption savings in the tens of millions of gallons. Under Project Icewarn, NASA, the U.S. Coast Guard and the National Oceanic Atmospheric Administration collaborated in development and demonstration of a system that permits safe year-round operations. It employs airborne radars, satellite communications relay and facsimile transmission to provide shippers and ships' masters up-to-date ice charts. Lewis Research Center contributed an accurate methods of measuring ice thickness by means of a special "short-pulse" type of radar. In a three-year demonstration program, Coast Guard aircraft equipped with Side-Looking Airborne Radar (SLAR) flew over the Great Lakes three or four times a week. The SLAR, which can penetrate clouds, provided large area readings of the type and distribution of ice cover. The information was supplemented by short-pulse radar measurements of ice thickness. The radar data was relayed by a NOAA satellite to a ground station where NOAA analyzed it and created picture maps, such as the one shown at lower left, showing where icebreakers can cut paths easily or where shipping can move through thin ice without the aid of icebreakers. The ice charts were then relayed directly to the wheelhouses of ships operating on the Lakes. Following up the success of the Great Lakes program, the icewarn team applied its system in another demonstration, this one a similarly successful application designed to aid Arctic coast shipping along the Alaskan North Slope. Further improvement of the ice-monitoring system is planned. Although aircraft-mounted radar is effective, satellites could provide more frequent data. After the launch this year of Seasat, an ocean-monitoring satellite, NASA will conduct tests to determine the ice-mapping capability and accuracy of satellite radar images.

Source record↗

Decoding Golden Eagle Movement Behavior from High-Resolution, Variable-Rate Telemetry Data Through Bayesian Filtering

The recent advances in animal tracking technology have enabled the collection of a vast amount of in situ data regarding the movement of wildlife at high spatiotemporal resolution. These data are usually available at variable time resolutions and contains noise (error) originating from GPS fixes. Decoding movement characteristics, particularly of flying animals, from telemetry data while handling these factors is a challenging yet important task for conservation purposes. Typically, this task is broken into two subtasks: resampling, and model calibration. The resampling subtask converts the variable rate positional data into a constant time interval data, while the model calibration subtask uses the resampled data to tune time-invariant parameters of the proposed models. For telemetry data at high temporal resolutions (order of 1 second), it is very challenging to decouple noise from actual movements using interpolation-based resampling techniques. Any errors introduced during resampling can significantly alter the the calibration and prediction attributes of the movement model. We address this problem through a unified Bayesian state-space framework that can handle both the resampling and calibration tasks in a single step. In addition, we use the speed and heading of the bird from telemetry data to regularize the position information of the bird. We use a Kalman filtering approach to include these nonlinearly related motion parameters within the state space framework. We cross-validated to quantify how this inclusion affects the model performance in estimating true bird movements. The relationship between the true state of the bird and environmental and topographical covariates is then represented parametrically. These parameters are then tuned using stochastic sampling strategies like Markov Chain Monte Carlo (MCMC). We use the telemetry data collected from golden eagles in the western USA to demonstrate the applicability of this approach to build a predictive, probabilistic movement model. Our preliminary results show that this approach provides improved predictive performance in terms of capturing higher-order motion parameters such as angular and horizontal accelerations, which may have simpler and more direct relationships with environmental covariates than corresponding speeds. In this talk, we will demonstrate how this state-space approach benefits the prediction capabilities of a movement model in simulating golden eagle paths through a wind power plant in Wyoming given certain atmospheric conditions. The model outcomes are aimed at informing mitigation strategies that can minimize the potential for collisions of golden eagles with wind turbines.

Bayesian methods↗

Terrestrial Planets Comparative Climatology (TPCC) mission concept

The authors and co-signers of the Terrestrial Planets Comparative Climatology (TPCC) mission concept white paper advocate that planetary science in the next decade would greatly benefit from comparatively studying the fundamental behavior of the atmospheres of Venus and Mars, contemporaneously and with the same instrumentation, to capture atmospheric response to the same solar forcing, and with a minimum of instrument-related variability. NASA’s vision is to understand the effects of the Sun on the solar system and the potential for life elsewhere (Obj,1.1 of [1]). This vision directly motivates the study of terrestrial climate histories and evolution, with the ultimate goal of understanding what factors allow the emergence and sustainability of life/biosphere. It is widely accepted that Venus, Mars, and Earth were formed at the same time and in the same neighborhood of the solar nebula [2,3], and had water—a necessary ingredient for life. In spite of strong overlap in the atmospheric composition of these planets, their current climatic states are vastly different. While Earth retains water, both Mars and Venus’ surfaces are now dry. The manner in which these terrestrial bodies evolved, and whether Mars or Venus currently hosts or previously hosted life, are fundamental mysteries in planetary science. As such, the TPCC mission concept has been developed.The TPCC guiding science goals are:(1) to understand Venus’ evolutionary history compared to Mars and Earth, and (2) to elucidate the driving forces and mechanisms of global circulation at Mars and Venus. These goals address the science question: Why have Mars’ and Venus’ climates evolved so differently from each other and from Earth’s? Furthermore, a greater depth of understanding our own solar system, enables us to understand the variability and likelihood of habitable conditions within extra-solar planetary systems.We suggest a mission architecture and payload that would determine the noble gas content and isotopic ratios, including the D/H ratio, at Venus with sufficient accuracy to identify the timeline of Venus’ climatic evolution for comparison with Mars. In addition, our concept would systematically “follow the energy” at Mars and Venus, in order to establish how solar forcing on the two atmospheric systems leads to the current balance of dynamical, chemical, and radiative processes responsible for each planet’s climatic state.This would include acquisition of the first global atmospheric wind profiles on either planet;global, continuous profiles of aerosols, gases, and water vapor at vertical resolutions higher than any previous mission;and regular UV and IR night glow observations tracking global circulation on both planets.To increase the fidelity of the planned planetary comparisons, the TPCC mission concept baseline (estimated to be New Frontiers Class)would send a single spacecraft to visit one planet and then the next, allowing measurements of each planet with the same instrumentation within similar parts of the solar cycle. The concept uses commercial solar-electric propulsion (SEP; e.g., [4])similar to the Psyche mission [5], plus a probesent deep into Venus’ atmosphere(e.g., LLISSE concept; [6]).Alternatively, “twin” spacecraft (with added probe to the Venus craft),developed at the same time,and launched close in time,would increase the overlap in operation periods and should also be studied. Either architecture is a timely alternative to independently developed and launched spacecraft payloads, which ultimately are more costly and difficult to inter-calibrate.

Leslie K. Tamppari↗

Impacts of High Variable Renewable Energy Futures on Electric-Sector Decision Making: Demand-Side Effects

Previous work by the Berkeley Lab describes how high shares of variable renewable energy (VRE) such as wind and solar power could change wholesale electricity price dynamics. These include the timing of when electricity is cheap or expensive, locational differences in the cost of electricity, and the degree of regularity or predictability in those costs. Many decentralized decision-makers on the demand-side may not yet have considered the implications of these possible future changes. In this report, we evaluate the sensitivity of a set of demand-side decisions to different levels of VRE penetration ranging from a low of 5-20% to a high of 40-50%. The analysis builds on hourly wholesale energy and capacity prices in different VRE scenarios for four wholesale markets in the United States for the year 2030 (CAISO, ERCOT, NYISO, and SPP). The principal question for this exploration is whether private and public electric-sector decisions that are made based on assumptions reflecting low VRE levels still achieve their intended objective in a high VRE scenario with 40-50% wind and solar? This scoping report evaluates the impacts of changing patterns of peak system needs on the benefits of demand reductions by examining the altered value of different energy efficiency (EE) measures. Similarly, we investigate new opportunities for large energy consumers that may arise from periods with very low wholesale electricity prices. We calculate the value of new process investments (e.g., hydrogen production and other generalized electro-commodities), showcase the varying value of new product storage investments (such as reservoir extensions at a desalination plant), and estimate the benefits of increased process flexibility that uses electricity as a process-input in addition to traditional fossil fuels (e.g., district energy systems). Finally, many decentralized decision-makers and end-use customers are not directly exposed to wholesale electricity prices but instead receive price signals from their retail electricity rates. As wind and solar shares increase, we compare the economic efficiency of flat retail rates relative to more dynamic time-of-use tariffs with and without critical peak-pricing events.

14 SOLAR ENERGY↗

Eta Carinae in the Context of the Most Massive Stars

Eta Car, with its historical outbursts, visible ejecta and massive, variable winds, continues to challenge both observers and modelers. In just the past five years over 100 papers have been published on this fascinating object. We now know it to be a massive binary system with a 5.54-year period. In January 2009, Car underwent one of its periodic low-states, associated with periastron passage of the two massive stars. This event was monitored by an intensive multi-wavelength campaign ranging from -rays to radio. A large amount of data was collected to test a number of evolving models including 3-D models of the massive interacting winds. August 2009 was an excellent time for observers and theorists to come together and review the accumulated studies, as have occurred in four meetings since 1998 devoted to Eta Car. Indeed, Car behaved both predictably and unpredictably during this most recent periastron, spurring timely discussions. Coincidently, WR140 also passed through periastron in early 2009. It, too, is a intensively studied massive interacting binary. Comparison of its properties, as well as the properties of other massive stars, with those of Eta Car is very instructive. These well-known examples of evolved massive binary systems provide many clues as to the fate of the most massive stars. What are the effects of the interacting winds, of individual stellar rotation, and of the circumstellar material on what we see as hypernovae/supernovae? We hope to learn. Topics discussed in this 1.5 day Joint Discussion were: Car: the 2009.0 event: Monitoring campaigns in X-rays, optical, radio, interferometry WR140 and HD5980: similarities and differences to Car LBVs and Eta Carinae: What is the relationship? Massive binary systems, wind interactions and 3-D modeling Shapes of the Homunculus & Little Homunculus: what do we learn about mass ejection? Massive stars: the connection to supernovae, hypernovae and gamma ray bursters Where do we go from here? (future directions) The Science Organizing Committee: Co-chairs: Augusto Damineli (Brazil) & Theodore R. Gull (USA). Members: D. John Hillier (USA), Gloria Koenigsberger (Mexico), Georges Meynet (Switzerland), Nidia Morrell (Chile), Atsuo T. Okazaki (Japan), Stanley P. Owocki (USA), Andy M.T. Pol- lock (Spain), Nathan Smith (USA), Christiaan L. Sterken (Belgium), Nicole St Louis (Canada), Karel A. van der Hucht (Netherlands), Roberto Viotti (Italy) and GerdWeigelt (Germany)

Gull, Theodore R.↗

Scanning Backscatter Lidar Observations for Characterizing 4-D Cloud and Aerosol Fields to Improve Radiative Transfer Parameterizations

Clouds have a powerful influence on atmospheric radiative transfer and hence are crucial to understanding and interpreting the exchange of radiation between the Earth's surface, the atmosphere, and space. Because clouds are highly variable in space, time and physical makeup, it is important to be able to observe them in three dimensions (3-D) with sufficient resolution that the data can be used to generate and validate parameterizations of cloud fields at the resolution scale of global climate models (GCMs). Simulation of photon transport in three dimensionally inhomogeneous cloud fields show that spatial inhomogeneities tend to decrease cloud reflection and absorption and increase direct and diffuse transmission, Therefore it is an important task to characterize cloud spatial structures in three dimensions on the scale of GCM grid elements. In order to validate cloud parameterizations that represent the ensemble, or mean and variance of cloud properties within a GCM grid element, measurements of the parameters must be obtained on a much finer scale so that the statistics on those measurements are truly representative. High spatial sampling resolution is required, on the order of 1 km or less. Since the radiation fields respond almost instantaneously to changes in the cloud field, and clouds changes occur on scales of seconds and less when viewed on scales of approximately 100m, the temporal resolution of cloud properties should be measured and characterized on second time scales. GCM time steps are typically on the order of an hour, but in order to obtain sufficient statistical representations of cloud properties in the parameterizations that are used as model inputs, averaged values of cloud properties should be calculated on time scales on the order of 10-100 s. The Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE) provides exceptional temporal (100 ms) and spatial (30 m) resolution measurements of aerosol and cloud backscatter in three dimensions. HARLIE was used in a ground-based configuration in several recent field campaigns. Principal data products include aerosol backscatter profiles, boundary layer heights, entrainment zone thickness, cloud fraction as a function of altitude and horizontal wind vector profiles based on correlating the motions of clouds and aerosol structures across portions of the scan. Comparisons will be made between various cloud detecting instruments to develop a baseline performance metric.

Schwemmer, Geary K.↗

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗