Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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 1,243 records · Page 69

Recent Updates in the SMAP Level-4 Soil Moisture Algorithm

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level-4 Soil Moisture (L4_SM) product. The L4_SM data are published with a mean latency of ~2.5 days from the time of observation and provide global, three-hourly, 9-km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. The L4_SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF). In 2018, the L4_SM algorithm was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high-quality, space-borne remote sensing observations. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of surface soil moisture from below under non-equilibrium conditions, which brings the model's surface soil moisture more in line with the SMAP Level-2 and in situ soil moisture. Moreover, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4, and the "catchment deficit" model variable was removed from the EnKF state vector to avoid degrading the model's groundwater estimates.Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation-minus-forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4_SM product and the lessons learned from the effort.

Reichle, Rolf H.↗

Real-Time Assimilation of Goes-Derived Products into A Mesoscale Model and It's Impact on Short-Term (06-36hr) Forecasts from 17 October 1998 through the Present

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. Spatial heterogeneity of these parameters makes correct specification problematic since measurements are not routinely available. A technique has been developed for assimilating GOES-IR skin temperature tendencies, solar insolation, and surface albedo into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. The technique has been successfully employed in a number of mesoscale models in case-study mode. We have taken the next step and developed a study to determine if assimilating these types of data into mesoscale models in real-time can improve short-term (648h) forecasts of temperature, relative humidity, and QPF on a daily basis over relatively large regions. Therefore, an operational modeling/assimilation system has been developed at the GHCC during the past summer that allows us to produce simulations out to 48 hours in a timely manor. The PSU/NCAR MM5 is used in a nested configuration with a 25 km grid covering the southeastern third of the US. The model has been on-line since 1 July 1998 and forecast products are posted on our web site. The satellite algorithms that generate data to be assimilated came on-line 17 October 1998. Quantitative assessment of the forecast quality is performed via traditional verification statistics. In addition, invaluable qualitative information is obtained through close collaboration with several NWSFO's who are using the MM5 products in real-time on a daily basis. The assimilation technique has been applied in an off-line mode since 17 October. Results based on bulk statistical verification of surface meteorology over the entire Southeastern US show that assimilating the GOES-derived land surface tendencies and solar radiation results in a significant reduction of the shelter air temperature and RH bias on a daily basis. In fact, the assimilation technique has produced improved temperature and RH forecasts for 97% of the 100 simulations performed to date. Work is currently underway to determine the sensitivity of the assimilation procedure to the availability of satellite data, length of assimilation period, model initialization, and synoptic-scale meteorological conditions. In addition, results from a detailed energy budget analysis using the Early Eta, our operational MM5, and the assimilation runs will help us to better understand the satellite assimilation the land-surface energy budge. Research during the spring-summer of 1999 will focus on the impact of the assimilation technique during the warm season where it is hypothesized that it can have a positive impact on QPF during conditions of weak synoptic-scale forcing.

Lapenta, William M.↗

Methods for Validation and Intercomparison of Remote Sensing and In situ Ice Water Measurements: Case Studies from CRYSTAL-FACE and Model Results

Remote sensing observations, such as those from AURA, are necessary to understand the role of cirrus in determining the radiative and humidity budgets of the upper troposphere. Using these measurements quantitatively requires comparisons with in situ measurements that have previously been validated. However, a direct comparison of remote and in situ measurements is difficult due to the requirement that the spatial and temporal overlap be sufficient in order to guarantee that both instruments are measuring the same air parcel. A difficult as this might be for gas phase intercomparisons, cloud inhomogeneities significantly exacerbate the problem for cloud ice water content measurements. The CRYSTAL-FACE mission provided an opportunity to assess how well such intercomparisons can be performed and to establish flight plans that will be necessary for validation of future satellite instruments. During CRYSTAL-FACE, remote and in situ instruments were placed on different aircraft (NASA's ER-2 and WB-59, and the two planes flew in tandem so that the in situ payload flew in the field of view of the remote instruments. We show here that, even with this type of careful flight planning, it is not always possible to guarantee that remote and in situ instruments are viewing the same air parcel. We use ice water data derived from the in situ Harvard Total Water (HV-TW) instrument, and the remote Goddard Cloud Radar System (CRS) and show that agreement between HV-TW and CRS is a strong function of the horizontal separation and the time delay between the aircraft transects. We also use a cloud model to simulate possible trajectories through a cloud and evaluate the use of statistical analysis in determining the agreement between the two instruments. This type of analysis should guide flight planning for future intercomparison efforts, whether for aircraft or satellite-borne instrumentation.

Sayres, D.S.↗

Using Multi-Scale Modeling Systems and Satellite Data to Study the Precipitation Processes

In recent years, exponentially increasing computer power has extended Cloud Resolving Model (CRM) integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 km2 in three-dimensions. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique. Recently, a multi-scale modeling system with unified physics was developed at NASA Goddard. It consists of (l) a cloud-resolving model (Goddard Cumulus Ensemble model, GCE model), (2) a regional scale model (a NASA unified weather research and forecast, WRF), (3) a coupled CRM and global model (Goddard Multi-scale Modeling Framework, MMF), and (4) a land modeling system. The same microphysical processes, long and short wave radiative transfer and land processes and the explicit cloud-radiation, and cloud-land surface interactive processes are applied in this multi-scale modeling system. This modeling system has been coupled with a multi-satellite simulator to use NASA high-resolution satellite data to identify the strengths and weaknesses of cloud and precipitation processes simulated by the model. In this talk, the recent developments and applications of the multi-scale modeling system will be presented. In particular, the results from using multi-scale modeling system to study the precipitating systems and hurricanes/typhoons will be presented. The high-resolution spatial and temporal visualization will be utilized to show the evolution of precipitation processes. Also how to use of the multi-satellite simulator tqimproy~ precipitation processes will be discussed.

Tao, Wei-Kuo↗

Vibration Response Predictions for Heavy Panel Mounted Components from Panel Acreage Environment Specifications

The development of new launch vehicles in the Aerospace industry often relies on response measurements taken from previously developed vehicles during various stages of liftoff and ascent, and from wind tunnel models. These measurements include sound pressure levels, dynamic pressures in turbulent boundary layers and accelerations. Rigorous statistical scaling methods are applied to the data to derive new environments and estimate the performance of new skin panel structures. Scaling methods have proven to be reliable, particularly for designs similar to the vehicles used as the basis for scaling, and especially in regions of smooth acreage without exterior protuberances or heavy components mounted to the panel. To account for response attenuation of a panel-mounted component due to its apparent mass at higher frequencies, the vibroacoustics engineer often reduces the acreage vibration according to a weight ratio first suggested by Barrett. The accuracy of the reduction is reduced with increased weight of the panel-mounted component, and does not account for low-frequency amplification of the component/panel response as a system. A method is proposed that combines acreage vibration from scaling methods with finite element analysis to account for the frequency-dependent dynamics of heavy panel-mounted components. Since the acreage and mass-loaded skins respond to the same dynamic input pressure, such pressure may be eliminated in favor of a frequency-dependent scaling function applied to the acreage vibration to predict the mass-loaded panel response. The scaling function replaces the Barrett weight ratio, and contains all of the dynamic character of the loaded and unloaded skin panels. The solution simplifies for spatially uncorrelated and fully correlated input pressures. Since the prediction uses finite element models of the loaded and unloaded skins, a rich suite of response data are available to the design engineer, including interface forces, stress and strain, as well as acceleration and displacement. An extension of the method is also developed to incorporate the effect of a local protuberance near a heavy component. Acreage environments from traditional scaling methods with and without protuberance effects serve as the basis for the extension. Authors:

Harrison, Phillip↗

A Novel Machine Learning Method for Surface PM2.5 Estimations from Geostationary Satellites

Particulate matter (PM) with a diameter of less or equal to 2.5 μm, known as PM , affects human health as it penetrates the respiratory system. The Environmental Protection Agency (EPA) measures the atmospheric concentration of PM using air quality monitors stationed throughout the Continental United States (CONUS). Such measurements are points on a spatial domain and therefore, might not be representative of the air quality at nearby areas considering that the composition of the atmosphere is highly variable from place to place. Satellite based AOD permits a spatially uniform means of estimating PM and new geostationary satellites provide high temporal and spatial resolution estimation of AOD. However, the concentration of PM is non-linearly dependent on other atmospheric parameters that include relative humidity, temperature, and height of the planetary boundary layer. This information may be estimated at similar spatial and temporal resolutions as AOD from numerical modeling such as from the National Oceanic and Atmospheric Administration’s (NOAA) High Resolution Rapid Refresh (HRRR) model which resolves near real-time atmospheric conditions over the CONUS. The estimation of PM concentration is a multi-parametric problem that considers the effect of temporal dependencies among the different parameters. Deep learning approaches are appropriate for such complex estimation problems as they intrinsically capture relations among multiple non-linear parameters. This study compares deep-learning methods to traditional regression analysis to demonstrate the capabilities of these methods in predicting PM2.5 concentrations. Additionally, a novel ensemble learning approach is employed to identify scientific processes that could further improve the estimation of PM concentration. Utilizing Long Short-Term Memory (LSTM) neural networks, which are suitable for multivariate time series estimation problems as they are capable of learning long-term dependencies, individual models are created for each EPA station and trained on the aforementioned dataset collocated over each station. Individual station models are merged if the model's performance is improved by reducing the root mean squared error (RMSE) metric. This ensemble training method ultimately reduces the RMSE value. Evaluation of these results provide insights into physical processes and related observable parameters that may contribute to PM concentrations. Identified parameters evaluated to be statistically different between the merged and unmerged models are expected to improve overall performance. These new parameters are then utilized for reevaluation of the deep learning methods with an extreme gradient boosting model with an RMSE of 5.5 providing the best results.

George Priftis↗

Analysis of fractal dimensions of rat bones from film and digital images

OBJECTIVES: (1) To compare the effect of two different intra-oral image receptors on estimates of fractal dimension; and (2) to determine the variations in fractal dimensions between the femur, tibia and humerus of the rat and between their proximal, middle and distal regions. METHODS: The left femur, tibia and humerus from 24 4-6-month-old Sprague-Dawley rats were radiographed using intra-oral film and a charge-coupled device (CCD). Films were digitized at a pixel density comparable to the CCD using a flat-bed scanner. Square regions of interest were selected from proximal, middle, and distal regions of each bone. Fractal dimensions were estimated from the slope of regression lines fitted to plots of log power against log spatial frequency. RESULTS: The fractal dimensions estimates from digitized films were significantly greater than those produced from the CCD (P=0.0008). Estimated fractal dimensions of three types of bone were not significantly different (P=0.0544); however, the three regions of bones were significantly different (P=0.0239). The fractal dimensions estimated from radiographs of the proximal and distal regions of the bones were lower than comparable estimates obtained from the middle region. CONCLUSIONS: Different types of image receptors significantly affect estimates of fractal dimension. There was no difference in the fractal dimensions of the different bones but the three regions differed significantly.

NASA Discipline Musculoskeletal↗

The Transition of High-Resolution NASA MODIS Sea Surface Temperatures into the WRF Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composite at 2-km resolution that has been implemented in version 3 of the National Weather Service (NWS) Weather Research and Forecasting (WRF) Environmental Modeling System (EMS). The WRF EMS is a complete, full physics numerical weather prediction package that incorporates dynamical cores from both the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). The installation, configuration, and execution of either the ARW or NMM models is greatly simplified by the WRF EMS to encourage its use by NWS Weather Forecast Offices (WFOs) and the university community. The WRF EMS is easy to run on most Linux workstations and clusters without the need for compilers. Version 3 of the WRF EMS contains the most recent public release of the WRF-NMM and ARW modeling system (version 3 of the ARW is described in Skamarock et al. 2008), the WRF Pre-processing System (WPS) utilities, and the WRF Post-Processing program. The system is developed and maintained by the NWS National Science Operations Officer Science and Training Resource Coordinator. To initialize the WRF EMS with high-resolution MODIS SSTs, SPoRT developed the composite product consisting of MODIS SSTs over oceans and large lakes with the NCEP Real-Time Global (RTG) filling data over land points. Filling the land points is required due to minor inconsistencies between the WRF land-sea mask and that used to generate the MODIS SST composites. This methodology ensures a continuous field that adequately initializes all appropriate arrays in WRF. MODIS composites covering the Gulf of Mexico, western Atlantic Ocean and the Caribbean are generated daily at 0400, 0700, 1600, and 1900 UTC corresponding to overpass times of the NASA Aqua and Terra polar orbiting satellites. The MODIS SST product is output in gridded binary-1 (GRIB-1) data format for a seamless incorporation into WRF via the WPS utilities. The full-resolution, 1-km MODIS product is sub-sampled to 2-km grid spacing due to limitations in handling very large dimensions in the GRIB-1 data format. The GRIB-1 files are posted online at ftp://ftp.nsstc.org/sstcomp/WRF/, which is directly accessed by the WRF EMS scripts. The MODIS SST composites are also downloaded to the EMS data server, which is accessible by the WRF EMS users and NWS WFOs. The SPoRT MODIS SST composite provides the model with superior detail of the ocean gradients around Florida and surrounding waters, whereas the operational RTG SST typically depicts a relatively smooth field and is not able to capture sharp horizontal gradients in SST. Differences of 2-3 C are common over small horizontal distances, leading to enhanced SST gradients on either side of the Gulf Stream and along the edges of the cooler shelf waters. These sharper gradients can in turn produce atmospheric responses in simulated temperature and wind fields as depicted in LaCasse et al. Differences in atmospheric verification statistics over a several month study were generally small in the vicinity of south Florida; however, the validation of SSTs at specific buoy locations revealed important improvements in the biases and RMS errors, especially in the vicinity of the cooler shelf waters off the east-central Florida coast. A current weakness in the MODIS SST product is the occurrence of occasional discontinuities caused by high latency in SST coverage due to persistent cloud cover. An enhanced method developed by Jedlovec et al. (2009, GHRSST User Symposium) reduces the occurrence of these problems by adding Advanced Microwave Scanning Radiometer -- EOS (AMSR-E) SST data to the compositing process. Enhanced SST composites are produced over the ocean regions surrounding the Continental U.S. at four times each day corresponding to Terra and Aqua equator crossing times. For a given day and overpass time, both MODInd AMSR-E data from the previous seven days form a collection used in the compositing. At each MODIS pixel, cloud-free SST values from the collection are used to form a weighted average based on their latency (number of days from the current day). In this way, recent SST data are given more weight than older data. One of the primary issues involved in incorporating the AMSR-E microwave data in the composites is the tradeoff between the decreased spatial resolution of the AMSR-E data (25 km) and the increased coverage due to its near all-weather capability. Currently, the AMSR-E is given a weight of 20% compared to MODIS data, thereby preserving the spatial structure observed in the MODIS data. Day-time (night-time) AMSR-E SST data from Aqua are used with both Terra and Aqua MODIS day-time (night-time) SST data sets.

Case, Jonathan L.↗

Recent Spatial and Temporal Anomalies and Trends of OLR as Observed by CERES and Computed Based on AIRS Retrievals

We show that a recent CERES-observed negative trend in OLR of approx.-0.1 W/sq m/yr averaged over the globe, for the time period of September 2002 through February 2010 used in this study, is found in the AIRS OLR data as well. Most importantly, even minute details (down to 1 x 1 Degree GCM-scale resolution) of spatial and temporal anomalies and trends of OLR as observed by CERES and computed based on AIRS-retrieved surface and atmospheric geophysical parameters over this time period are essentially the same. We see this correspondence even in the very large spatial variations of these trends with local values ranging from -2.6 W/sq m/yr to +3.0 W/sq m/yr in the tropics. This essentially perfect agreement of OLR anomalies and even local trends derived from observations by two different instruments, in totally independent and different manners, implies that both sets of results must be highly accurate; and indirectly validates the anomalies and trends of other AIRS derived products as well. These products show that global and regional anomalies and trends of OLR, water vapor and cloud cover over the last 7+ years are strongly influenced by El-Nino-La Nina cycles . We use the anomalies and trends of AIRS derived products to explain why the global OLR has a large negative trend over this time period; Global and tropical OLR began to decrease significantly at the onset of a strong La Nina in mid-2007. AIRS products show that cloudiness and mid-tropospheric water vapor began to increase in the tropics at roughly the same time, especially in the region 5degN - 20degS latitude extending eastward from 150degW to 30degE longitude, with a corresponding very large drop in OLR in this region. Late 2009 is characterized by a strong El-Nino, with a corresponding change in sign of observed tropical water vapor, cloud cover, and OLR anomalies. If one excludes the area 5degN - 20degS, 150degW - 30degE from the statistics, area mean OLR trends over the rest of the globe are substantially reduced over the time period under study.

Molnar, Gyula↗

Using NASA Environmental Data to Enhance Public Health Decision Making

The Universities Space Research Association at the NASA Marshall Space Flight Center is collaborating with the University of Alabama at Birmingham (UAB) School of Public Health and the Centers for Disease Control and Prevention (CDC) to address issues of environmental health and enhance public health decision making by utilizing NASA remotely sensed data and products. The objectives of this collaboration are to develop high-quality spatial data sets of environmental variables, and deliver the data sets and associated analyses to local, state and federal end-user groups. These data can be linked spatially and temporally to public health data, such as mortality and disease morbidity, for further analysis and decision making. Three daily environmental data sets have been developed for the conterminous U.S. on different spatial resolutions for the time period 2003-2008: (1) spatial surfaces of estimated fine particulate matter (PM2.5) exposures on a 10-km grid utilizing the US Environmental Protection Agency (EPA) ground observations and NASA s MODerate-resolution Imaging Spectroradiometer (MODIS) data; (2) a 1-km grid of Land Surface Temperature (LST) using MODIS data; and (3) a 12-km grid of daily Solar Insolation (SI) and maximum and minimum air temperature using the North American Land Data Assimilation System (NLDAS) forcing data. These environmental data sets will be linked with public health data from the UAB REasons for Geographic And Racial Differences in Stroke (REGARDS) national cohort study to determine whether exposures to these environmental risk factors are related to cognitive decline and other health outcomes. These environmental datasets and public health linkage analyses will be made available to public health professionals, researchers and the general public through the CDC Wide-ranging Online Data for Epidemiologic Research (WONDER) system and through peer reviewed publications. To date, two of the data sets have been released to the public in CDC WONDER, Daily Air Temperature and Heat Index for years 1979-2010, and Daily Fine Particulate Matter (PM2.5) air quality measures for years 2003-2008. These data in CDC WONDER can be aggregated to the county-level, state-level, or regional-level as per users need and downloaded in tabular, graphical, and map formats. The summary statistical output are available to web and app developers via the WONDER Application Programming Interface (API). The linkage of these data with the CDC WONDER system provides a significant addition to CDC WONDER, allowing public health researchers and policy makers to better include environmental exposure data in the context of other health data available in CDC WONDER online system. It also substantially expands public access to NASA environmental data, making their use by a wide range of decision makers feasible.

Al-Hamdan, Mohammad↗

Advanced Techniques for Assessment of Postural and Locomotor Ataxia, Spatial Orientation, and Gaze Stability

In addition to adapting to microgravity, major neurovestibular problems of space flight include postflight difficulties with standing, walking, turning corners, and other activities that require stable upright posture and gaze stability. These difficulties inhibit astronauts' ability to stand or escape from their vehicle during emergencies. The long-ter7n goal of the NSBRI is the development of countermeasures to ameliorate the effects of long duration space flight. These countermeasures must be tested with valid and reliable tools. This project aims to develop quantitative, parametric approaches for assessing gaze stability and spatial orientation during normal gait and when gait is perturbed. Two of this year's most important findings concern head fixation distance and ideal trajectory analysis. During a normal cycle of walking the head moves up and down linearly. A simultaneous angular pitching motion of the head keeps it aligned toward an imaginary point in space at a distance of about one meter in front of a subject and along the line of march. This distance is called the head fixation distance. Head fixation distance provides the fundamental framework necessary for understanding the functional significance of the vestibular reflexes that couple head motion to eye motion. This framework facilitates the intelligent design of counter-measures for the effects of exposure to microgravity upon the vestibular ocular reflexes. Ideal trajectory analysis is a simple candidate countermeasure based upon quantifying body sway during repeated up and down stair stepping. It provides one number that estimates the body sway deviation from an ideal sinusoidal body sway trajectory normalized on the subject's height. This concept has been developed with NSBRI funding in less than one year. These findings are explained in more detail below. Compared to assessments of the vestibuo-ocular reflex, analysis of vestibular effects on locomotor function is relatively less well developed and quantified. We are improving this situation by applying methodologies such as nonlinear orbital stability to quantify responses and by using multivariate statistical approaches to link together the responses across separate tests. In this way we can exploit the information available and increase the ability to discriminate between normal and pathological responses. Measures of stability and orientation are compared to measures such as dynamic visual acuity and with balance function tests. The responses of normal human subjects and of patients having well documented pathophysiologies are being characterized. When these studies are completed, we should have a clearer idea about normal and abnormal patterns of eye, head, and body movements during locomotion and their stability in a wide range of environments. We plan eventually to use this information to validate the efficacy of candidate neurovestibular and neuromuscular rehabilitative techniques. Some representative studies made during this year are summarized.

Wall, Conrad., III↗

A polymer, random walk model for the size-distribution of large DNA fragments after high linear energy transfer radiation

DNA double-strand breaks (DSBs) produced by densely ionizing radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. Stochastic DSB clustering at large scales, from > 100 Mbp down to < 0.01 Mbp, is modeled using computer simulations and analytic equations. A random-walk, coarse-grained polymer model for chromatin is combined with a simple track structure model in Monte Carlo software called DNAbreak and is applied to data on alpha-particle irradiation of V-79 cells. The chromatin model neglects molecular details but systematically incorporates an increase in average spatial separation between two DNA loci as the number of base-pairs between the loci increases. Fragment-size distributions obtained using DNAbreak match data on large fragments about as well as distributions previously obtained with a less mechanistic approach. Dose-response relations, linear at small doses of high linear energy transfer (LET) radiation, are obtained. They are found to be non-linear when the dose becomes so large that there is a significant probability of overlapping or close juxtaposition, along one chromosome, for different DSB clusters from different tracks. The non-linearity is more evident for large fragments than for small. The DNAbreak results furnish an example of the RLC (randomly located clusters) analytic formalism, which generalizes the broken-stick fragment-size distribution of the random-breakage model that is often applied to low-LET data.

DNA/radiation effects↗

Calibration of TOPEX/POSEIDON at Platform Harvest

We present estimates for the mean bias of the TOPEX/POSEIDON NASA altimeter (ALT) and the Centre National d'Etudes Spatiales altimeter (SSALT) using in-situ data gathered at Platform Harvest during the first 36 cycles of the mission. Data for 21 overflights of the ALT and six overflights of the SSALT have been analyzed. The analysis includes an independent assessment of in-situ measurements of sea level, the radial component of the orbit, wet tropospheric path delay, and ionospheric path delay. (The sign convention used is such that, to correct the geophysical data record values for sea level, add the bias algebraically. Unless otherwise stated, the uncertainty in a given parameter is depicted by +/- sigma(sub x), where sigma(sub x) is the sample standard deviation of x about the mean.) Tide gauges at Harvest provide estimates of sea level with an uncertainty of +/- 1.5 cm. The uncertainty in the radial component of the orbit is estimated to be +/- 1.3 cm. In-situ measurements of tropopsheric path delay at Harvest compare to within +/- 1.3 cm of the TOPEX/POSEIDON microwave radiometer, and in-situ measurements of the ionospheric path delay compare to within -0.4 +/- 0.7 cm of the dual-frequency ALT and 1.1 +/- 0.6 cm of Doppler orbitography and radiopositioning integrated by satellite. We obtain mean bias estimates of -14.5 +/- 2.9 cm for the ALT and +0.9 +/- 3.1 cm for the SSALT (where the uncertainties are based on the standard deviation of the estimated mean (sigma(sub bar x/y), which is derived from sample statistics and estimates for errors that cannot be observed). These results are consistent with independent estimates for the relative bias between the two altimeters. A linear regression applied to the complete set of data shows that there is a discernable secular trend in the time series for the ALT bias estimates. A preliminary analysis of data obtained through cycle 48 suggests that the apparent secular drift may be the result of a poorly sampled annual signal.

Christensen, E. J.↗

Bayesian Analysis of the Power Spectrum of the Cosmic Microwave Background

There is a wealth of cosmological information encoded in the spatial power spectrum of temperature anisotropies of the cosmic microwave background. The sky, when viewed in the microwave, is very uniform, with a nearly perfect blackbody spectrum at 2.7 degrees. Very small amplitude brightness fluctuations (to one part in a million!!) trace small density perturbations in the early universe (roughly 300,000 years after the Big Bang), which later grow through gravitational instability to the large-scale structure seen in redshift surveys... In this talk, I will discuss a Bayesian formulation of this problem; discuss a Gibbs sampling approach to numerically sampling from the Bayesian posterior, and the application of this approach to the first-year data from the Wilkinson Microwave Anisotropy Probe. I will also comment on recent algorithmic developments for this approach to be tractable for the even more massive data set to be returned from the Planck satellite.

Bayesian inference↗

The Next Generation X-Ray Galaxy Survey with eROSITA

eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e. not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50–200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect ≳15 000 galaxies (3σ significance) at 50–200 Mpc, which is ~ 100× more normal galaxies than detected in any X-ray survey to date.

Surveys↗

An Assessment of Differences Between Cloud Effective Particle Radius Retrievals for Marine Water Clouds from Three MODIS Spectral Bands

The Moderate Resolution Imaging Spectroradiometer (MODIS) cloud product provides three separate 1 km resolution retrievals of cloud particle effective radii (r (sub e)), derived from 1.6, 2.1 and 3.7 micron band observations. In this study, differences among the three size retrievals for maritime water clouds (designated as r (sub e), 1.6 r (sub e), 2.1 and r (sub e),3.7) were systematically investigated through a series of case studies and global analyses. Substantial differences are found between r (sub e),3.7 and r (sub e),2.1 retrievals (delta r (sub e),3.7-2.l), with a strong dependence on cloud regime. The differences are typically small, within +/- 2 micron, over relatively spatially homogeneous coastal stratocumulus cloud regions. However, for trade wind cumulus regimes, r (sub e),3.7 was found to be substantially smaller than r (sub e),2.1, sometimes by more than 10 micron. The correlation of delta r(sub e),3.7-2.1 with key cloud parameters, including the cloud optical thickness (tau), r (sub e) and a cloud horizontal heterogeneity index (H-sigma) derived from 250 m resolution MODIS 0.86 micron band observations, were investigated using one month of MODIS Terra data. It was found that differences among the three r (sub e) retrievals for optically thin clouds (tau <5) are highly variable, ranging from - 15 micron to 10 micron, likely due to the large MODIS retrieval uncertainties when the cloud is thin. The delta r (sub e),3.7-2.1 exhibited a threshold-like dependence on both r (sub e),2.l and H-sigma. The re,3.7 is found to agree reasonably well with re,2.! when re,2.l is smaller than about 15J-Lm, but becomes increasingly smaller than re,2.1 once re,2.! exceeds this size. All three re retrievals showed little dependence when H-sigma < 0.3 (defined as standard deviation divided by the mean for the 250 m pixels within a 1 km pixel retrieval). However, for H-=sigma >0.3, both r (sub e),1.6 and r (sub e),2.1 were seen to increase quickly with H-sigma. On the other hand, r (sub e),3.7 statistics showed little dependence on H-sigma and remained relatively stable over the whole range of H-sigma values. Potential contributing causes to the substantial r (sub e),3.7 and r (sub e),2.1 differences are discussed. In particular, based on both 1-D and 3-D radiative transfer simulations, we have elucidated mechanisms by which cloud heterogeneity and 3-D radiative effects can cause large differences between r (sub e),3.7 and r (sub e),2.l retrievals for highly inhomogeneous clouds. Our results suggest that the contrast in observed delta r (sub e)3.7-2.1 between cloud regimes is correlated with increases in both cloud r (sub e) and H-sigma. We also speculate that in some highly inhomogeneous drizzling clouds, vertical structure induced by drizzle and 3-D radiative effects might operate together to cause dramatic differences between r (sub e),3.7 and r (sub e),2.1 retrievals.

Platnick, Steven↗

A Downscaling Analysis of the Urban Influence on Rainfall: TRMM Satellite Component AMS Conference on Satellite Meteorology and Oceanography

A recent publication by Shepherd et al. (2002) demonstrated the feasibility of using TRMM precipitation radar (PR) estimates to identify precipitation anomalies caused by urbanization. The approach is particularly useful for investigating this global process because TRMM data span large portions of the globe and comprise an extended temporal dataset. Recent literature suggests that urbanized regions of Houston, Texas may be influencing lightning and precipitation formation over and downwind of the city. Possible mechanisms include: (1) enhanced convergence through interactions between the sea breeze, Galveston bay breeze, and urban heat island circulations, (2) enhanced convergence due to increased surface roughness over the city and/or destabilization of the boundary layer by the UHI, or (3) enhanced cloud condensation nuclei due to urban and industrial aerosol sources. In this study, a downscaling analysis of spatial and temporal trends in rainfall around the Houston Area is being conducted. The downscaling analysis concept involves identifying and quantifying urban rainfall anomalies at progressively smaller spatial and temporal scales using the TRMM satellite, ground-based radar, and a dense network of rain gauges. The goal is to test the hypothesis that the Houston urban district and regions in the climatological downwind region of the city exhibit enhanced rainfall amounts relative to the climatological upwind regions. TRMM was launched in 1997 and currently operates in a low-inclination (35 deg), non-sun-synchronous orbit at an altitude of 402 km (350 km prior to August 2001). The satellite analysis follows the methodologies described in Shepherd et al. (2002). Nearly five years of TRMM PR-derived mean monthly rainfall estimates are utilized to produce annual and warm season isohyetal analyses around Houston. Early results indicate that rainfall rates (mm/h) for the entire period are largest within 100 km northeast and east of Houston (e.g. the "hypothesized downwind region"). The mean rainfall rate over the Houston urban center is 30.5% larger than the upwind control region. The mean rainfall rate in the downwind region is 34.4% larger than the upwind region. An analysis of a parameter called the urban rainfall ratio (URR) illustrates that 65% (88%) of the satellite-derived rainfall rates in the downwind (upwind control) region are greater (less) than the mean background rainfall rate of the entire study region. When the data is stratified by summer months from 1998 to 2001 (June-August), even greater influence over and downwind of the urban area is observed in the statistics. This result is consistent with published reports of urban-generated rainfall being more prevalent in the warm season. The research demonstrates that the evolving TRMM satellite climatology is a credible way to detect mesoscale precipitation signatures that may be linked to urbanization. Early results also corroborate recent findings on Houston-induced convection/drainfall anomalies. Burian and Shepherd will report on other aspects of the downscaling analysis in future forums, but early rain gauge results are consistent with the satellite-based observations.

Shepherd, J. Marshall↗

Comparison of Continuous Wave CO2 Doppler Lidar Calibration Using Earth Surface Targets in Laboratory and Airborne Measurements

Routine backscatter, beta, measurements by an airborne or space-based lidar from designated earth surfaces with known and fairly uniform beta properties can potentially offer lidar calibration opportunities. This can in turn be used to obtain accurate atmospheric aerosol and cloud beta measurements on large spatial scales. This is important because achieving a precise calibration factor for large pulsed lidars then need not rest solely on using a standard hard target procedure. Furthermore, calibration from designated earth surfaces would provide an inflight performance evaluation of the lidar. Hence, with active remote sensing using lasers with high resolution data, calibration of a space-based lidar using earth's surfaces will be extremely useful. The calibration methodology using the earth's surface initially requires measuring beta of various earth surfaces simulated in the laboratory using a focused continuous wave (CW) CO2 Doppler lidar and then use these beta measurements as standards for the earth surface signal from airborne or space-based lidars. Since beta from the earth's surface may be retrieved at different angles of incidence, beta would also need to be measured at various angles of incidences of the different surfaces. In general, Earth-surface reflectance measurements have been made in the infrared, but the use of lidars to characterize them and in turn use of the Earth's surface to calibrate lidars has not been made. The feasibility of this calibration methodology is demonstrated through a comparison of these laboratory measurements with actual earth surface beta retrieved from the same lidar during the NASA/Multi-center Airborne Coherent Atmospheric Wind Sensor (MACAWS) mission on NASA's DC8 aircraft from 13 - 26 September, 1995. For the selected earth surface from the airborne lidar data, an average beta for the surface was established and the statistics of lidar efficiency was determined. This was compared with the actual lidar efficiency determined with the standard calibrating hard target.

Jarzembski, Maurice A.↗