Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observation Data”

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,369 records · Page 76

A Study of the Kepler K2 Variable EPIC 211957146 Exhibiting a Variable O’ Connell Effect

We present the multi-band photometric and spectroscopic study of an over-contact binary system, EPIC 211957146. The light curves exhibit a variable O’ Connell effect, confirmed from our observational data and the Kepler K2 data. The best photometric solution incorporating a dark spot over the primary component unveils that the system has a low-mass ratio ( q ∼ 0.17) and a high inclination ( i ∼ 85°). To confirm the solution and constrain the uncertainty, Monte-Carlo simulations are performed and the results are reported. Based on the O–C diagram analysis, we see that the variable shows a period increase at the rate of dP / dt ∼ 1.06 × 10{sup −6} days yr{sup −1}, which is higher than the theoretically predicted value. Presence of a third body having a period of ∼16.23 years is evident from the O–C diagram. No filled-in effect is observed in the H α line, while the effect is vividly present in the Na line. From the Kepler K2 data, we found that the primary and secondary minima exhibit an anti-correlated O–C variation followed by an erratic behavior. This is possibly caused by the longitudinal motion of the spot, and hence, we set a lower limit of ∼40 days for the spot modulation. We also observe a possibly associated photometric difference in the primary depth by comparing our light curves with Kepler K2 normalized light curves. This system has a low-mass ratio and a high fill-out factor, and, theoretically, such a physical configuration would lead to a merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Heterogeneous local dynamics revealed by classification analysis of spatially disaggregated time series data

Time series data provide a crucial window into infectious disease dynamics, yet their utility is often limited by the spatially aggregated form in which they are presented. When working with time series data, violating the implicit assumption of homogeneous dynamics below the scale of spatial aggregation could bias inferences about underlying processes. We tested this assumption in the context of the 2015–2016 Zika epidemic in Colombia, where time series of weekly case reports were available at national, departmental, and municipal scales. First, we performed a descriptive analysis, which showed that the timing of departmental-level epidemic peaks varied by three months and that departmental-level estimates of the time-varying reproduction number, R(t), showed patterns that were distinct from a national-level estimate. Second, we applied a classification algorithm to six features of proportional cumulative incidence curves, which showed that variability in epidemic duration, the length of the epidemic tail, and consistency with a cumulative normal density curve made the greatest contributions to distinguishing groups. Third, we applied this classification algorithm to data simulated with a stochastic transmission model, which showed that group assignments were consistent with simulated differences in the basic reproduction number, R0. This result, along with associations between spatial drivers of transmission and group assignments based on observed data, suggests that the classification algorithm is capable of detecting differences in temporal patterns that are associated with differences in underlying drivers of incidence patterns. Overall, this diversity of temporal patterns at local scales underscores the value of spatially disaggregated time series data.

60 APPLIED LIFE SCIENCES↗

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.

79 ASTRONOMY AND ASTROPHYSICS↗

A dynamo model of magnetic activity in solar-like stars with different rotational velocities

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1 M{sub ☉} stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f{sub m} of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f{sub m}. Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

79 ASTRONOMY AND ASTROPHYSICS↗

Supermassive dark-matter Q-balls in galactic centers?

Though widely accepted, it is not proven that supermassive compact objects (SMCOs) residing in galactic centers are black holes. In particular, the Milky Way’s SMCO can be a giant nontopological soliton, Q-ball, made of a scalar field: this fits perfectly all observational data. Similar but tiny Q-balls produced in the early Universe may constitute, partly or fully, the dark matter. This picture explains in a natural way, why our SMCO has very low accretion rate and why the observed angular size of the corresponding radio source is much smaller than expected. Interactions between dark-matter Q-balls may explain how SMCOs were seeded in galaxies and resolve well-known problems of standard (non-interacting) dark matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

West Flank Coso FORGE: 3D Temperature Model

x,y,z data of the 3D temperature model for the West Flank Coso FORGE site. Model grid spacing is 250m. The temperature model for the Coso geothermal field used over 100 geothermal production sized wells and intermediate-depth temperature holes. At the near surface of this model, two boundary temperatures were assumed: (1) areas with surface manifestations, including fumaroles along the northeast striking normal faults and northwest striking dextral faults with the hydrothermal field, a temperature of ~104 deg C was applied to datum at +1066 meters above sea level elevation, and (2) a near-surface temperature at about 10 meters depth, of 20 deg C was applied below the diurnal and annual conductive temperature perturbations. These assumptions were based on heat flow studies conducted at the CVF and for the Mojave Desert. On the edges of the hydrothermal system, a 73 deg C/km temperature gradient contour was established using conductive gradient data from shallow and intermediate-depth temperature holes. This contour was continued to all elevation datums between the 20 deg C surface and -1520 meters below mean sea level. Because the West Flank is outside of the geothermal field footprint, during Phase 1, the three wells inside the FORGE site were incorporated into the preexisting temperature model. To ensure a complete model was built based on all the available data sets, measured bottom-hole temperature gradients in certain wells were downward extrapolated to the next deepest elevation datum (or a maximum of about 25% of the well depth where conductive gradients are evident in the lower portions of the wells). After assuring that the margins of the geothermal field were going to be adequately modelled, the data was contoured using the Kriging method algorithm. Although the extrapolated temperatures and boundary conditions are not rigorous, the calculated temperatures are anticipated to be within ~6 deg C (20 deg F), or one contour interval, of the observed data within the Coso geothermal field. Based on a lack of temperature data west of 74-2TCH, the edges of this model still seem to have an effect on West Flank modeled temperatures.

15 GEOTHERMAL ENERGY↗

Surface temperature comparison of the Arctic winter MOSAiC observations, ERA5 reanalysis, and MODIS satellite retrieval

Atmospheric model systems, such as those used for weather forecast and reanalysis production, often have significant and systematic errors in their representation of the Arctic surface energy budget and its components. The newly available observation data of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (2019/2020) enable a range of model analyses and validation in order to advance our understanding of potential model deficiencies. In the present study, we analyze deficiencies in the surface radiative energy budget over Arctic sea ice in the ERA5 global atmospheric reanalysis by comparing against the winter MOSAiC campaign data, as well as, a pan-Arctic level-2 MODIS ice surface temperature remote sensing product. We find that ERA5 can simulate the timing of radiatively clear periods, though it is not able to distinguish the two observed radiative Arctic winter states, radiatively clear and opaquely cloudy, in the distribution of the net surface radiative budget. The ERA5 surface temperature over Arctic sea ice has a conditional error with a positive bias in radiatively clear conditions and a negative bias in opaquely cloudy conditions. The mean surface temperature error is 4°C for radiatively clear situations at MOSAiC and up to 15°C in some parts of the Arctic. The spatial variability of the surface temperature, given by 4 observation sites at MOSAiC, is not captured by ERA5 due to its spatial resolution but represented in the level-2 satellite product. The sensitivity analysis of possible error sources, using satellite products of snow depth and sea ice thickness, shows that the positive surface temperature errors during radiatively clear events are, to a large extent, caused by insufficient sea ice thickness and snow depth representation in the reanalysis system. A positive bias characterizes regions with ice thickness greater than 1.5 m, while the negative bias for thinner ice is partly compensated by the effect of snow.

54 ENVIRONMENTAL SCIENCES↗

Development and Validation of a Lifecycle-based Prognostics Architecture with Test Bed Validation

On-line monitoring and tracking of nuclear plant system and component degradation is being investigated as a method for improving the safety, reliability, and maintainability of aging nuclear power plants. Accurate prediction of the current degradation state of system components and structures is important for accurate estimates of their remaining useful life (RUL). The correct quantification and propagation of both the measurement uncertainty and model uncertainty is necessary for quantifying the uncertainty of the RUL prediction. This research project developed and validated methods to perform RUL estimation throughout the lifecycle of plant components. Prognostic methods should seamlessly operate from beginning of component life (BOL) to end of component life (EOL). We term this "Lifecycle Prognostics." When a component is put into use, the only information available may be past failure times of similar components used in similar conditions, and the predicted failure distribution can be estimated with reliability methods such as Weibull Analysis (Type I Prognostics). As the component operates, it begins to degrade and consume its available life. This life consumption may be a function of system stresses, and the failure distribution should be updated to account for the system operational stress levels (Type II Prognostics). When degradation becomes apparent, this information can be used to again improve the RUL estimate (Type III Prognostics). This research focused on developing prognostics algorithms for the three types of prognostics, developing uncertainty quantification methods for each of the algorithms, and, most importantly, developing a framework using Bayesian methods to transition between prognostic model types and update failure distribution estimates as new information becomes available. The developed methods were then validated on a range of accelerated degradation test beds. The ultimate goal of prognostics is to provide an accurate assessment for RUL predictions, with as little uncertainty as possible. From a reliability and maintenance standpoint, there would be improved safety by avoiding all failures. Calculated risk would decrease, saving money by avoiding unnecessary maintenance. One major bottleneck for data-driven prognostics is the availability of run-to-failure degradation data. Without enough degradation data leading to failure, prognostic models can yield RUL distributions with large uncertainty or mathematically unsound predictions. To address these issues a "Lifecycle Prognostics" method was developed to create RUL distributions from Beginning of Life (BOL) to End of Life (EOL). This employs established Type I, II, and III prognostic methods, and Bayesian transitioning between each Type. Bayesian methods, as opposed to classical frequency statistics, show how an expected value, a priori, changes with new data to form a posterior distribution. For example, when you purchase a component you have a prior belief, or estimation, of how long it will operate before failing. As you operate it, you may collect information related to its condition that will allow you to update your estimated failure time. Bayesian methods are best used when limited data are available. The use of a prior also means that information is conserved when new data are available. The weightings of the prior belief and information contained in the sampled data are dependent on the variance (uncertainty) of the prior, the variance (uncertainty) of the data, and the amount of measured data (number of samples). If the variance of the prior is small compared to the uncertainty of the data, the prior will be weighed more heavily. However, as more data are collected, the data will be weighted more heavily and will eventually swamp out the prior in calculating the posterior distribution of model parameters. Fundamentally Bayesian analysis updates a prior belief with new data to get a posterior belief. The general approach to applying the Bayesian method to lifecycle prognostics consisted of identifying the prior, which is the RUL estimate and uncertainty from the previous prognostics type, and combining it with observational data related to the newer prognostics type. The resulting lifecycle prognostics algorithm uses all available information throughout the component lifecycle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On the Cloud Observations in JAXA's Next Coming Satellite Missions

The use of JAXA's next generation satellites, the EarthCARE and the GCOM-C, for observing overall cloud systems on the Earth is discussed. The satellites will be launched in the middle of 2010-era and contribute for observing aerosols and clouds in terms of climate change, environment, weather forecasting, and cloud revolution process study. This paper describes the role of such satellites and how to use the observing data showing concepts and some sample viewgraphs. Synergistic use of sensors is a key of the study. Visible to infrared bands are used for cloudy and clear discriminating from passively obtained satellite images. Cloud properties such as the cloud optical thickness, the effective particle radii, and the cloud top temperature will be retrieved from visible to infrared wavelengths of imagers. Additionally, we are going to combine cloud properties obtained from passive imagers and radar reflectivities obtained from an active radar in order to improve our understanding of cloud evolution process. This is one of the new techniques of satellite data analysis in terms of cloud sciences in the next decade. Since the climate change and cloud process study have mutual beneficial relationship, a multispectral wide-swath imagers like the GCOM-C SGLI and a comprehensive observation package of cloud and aerosol like the EarthCARE are both necessary.

satellite remote sensing↗

InSight Robotic Arm Testing Activities for HP3 Mole Anomaly Recovery on Mars

The InSight lander’s Heat Flow and Physical Properties Package (HP3) was deployed on Mars in February 2019and began attempting to penetrate to its target depth range of3-5 meters shortly thereafter. However, the mole’s downwardprogress stopped after only 35 cm of penetration. In response,the project convened an Anomaly Response Team (ART) andsince then has been attempting to diagnose the problem andassist the mole using the tools available on Mars. The key assetused in this effort has been the Instrument Deployment System(IDS), which includes two cameras and a robotic InstrumentDeployment Arm (IDA). Since the IDS was originally intendedonly to deploy InSight’s primary instruments to the Martiansurface, new testbed setups, experiments, and operational protocols (e.g., command sequences) were required and had to bedeveloped on a short timeline. The HP3 Mole ART first focusedon gathering all observable data on Mars about the state ofthe mole and Support Structure Assembly (SSA), as well asthe physical properties of the Martian regolith. This includedusing the robotic arm to point the IDC at the SSA duringdiagnostic hammering attempts to observe motion of the SSAand science tether. Images taken during these attempts revealedsome motion of the SSA, but no apparent change in mole depth.At JPL, the IDS and Testbed teams re-created the hardwareconfiguration on Mars based on limited knowledge of the mole’sstate. They devised and tested techniques to use the roboticarm and cameras to accomplish previously untested activitieson Mars, including imaging the HP3, using the IDA to interactwith the terrain, and using the IDA to move the SSA awayfrom the partially-embedded mole. The team executed the morepromising techniques on Mars. After diagnostic hammering onMars, the team decided to move the SSA to gain visibility ofthe mole’s configuration and access to the soil around the mole.After developing the technique and practicing the maneuver inthe InSight testbed, the team lifted the SSA on Mars and placedit behind the mole. This revealed a pit surrounding the nowexposed mole, observations of which provided essential cluesfor determining the root cause of the mole’s lack of progress.The IDS and Testbed teams altered the testbed to match thesituation on Mars. They devised IDA techniques to determinethe Martian soil properties and assist the mole’s descent. Theytested these techniques in the testbed and executed the more978-1-7281-2734-7/20/$31.00 c 2021 IEEE. Copyright 2020 CaliforniaInstitute of Technology. U.S. Government sponsorship acknowledged.promising ones on Mars. These include using the robotic arm toalter the regolith near the mole and to push on the mole while ithammers. This paper discusses the anomaly resolution testing inthe testbed at JPL, describes how the IDS team prepared for theanomaly recovery activities on Mars, and provides preliminaryresults of the efforts to assist the HP3 mole on Mars.

Kim, Junggon↗

A Spatiotemporal-Aware Weighting Scheme for Improving Climate Model Ensemble Predictions

Multimodel ensembling has been widely used to improve climate model predictions, and the improvement strongly depends on the ensembling scheme. In this work, we propose a Bayesian neural network (BNN) ensembling method, which combines climate models within a Bayesian model averaging framework, to improve the predictive capability of model ensembles. Our proposed BNN approach calculates spatiotemporally varying model weights and biases by leveraging individual models' simulation skill, calibrates the ensemble prediction against observations by considering observation data uncertainty, and quantifies epistemic uncertainty when extrapolating to new conditions. More importantly, the BNN method provides interpretability about which climate model contributes more to the ensemble prediction at which locations and times. Thus, beyond its predictive capability, the method also brings insights and understanding of the models to guide further model and data development. In this study, we design experiments using an ensemble of CMIP6 climate model simulations to illustrate the BNN ensembling method's capability with respect to prediction accuracy, interpretability, and uncertainty quantification (UQ). We demonstrate that BNN can correctly assign larger weights to the regions and seasons where the individual model fits the observation better. Moreover, its offered interpretability is consistent with our understanding of localized climate model performance. Additionally, BNN shows an increasing uncertainty when the prediction is farther away from the period with constrained data, which appropriately reflects our trustworthiness of the models in the changing climate.

54 ENVIRONMENTAL SCIENCES↗

Demography, dynamics and data: building confidence for simulating changes in the world's forests

Vegetation demographic models (VDMs) are advanced tools for simulating forest responses to climate and land-use changes, and are essential for projecting carbon cycling and large-scale forest management strategies. Despite their increasing incorporation into Earth System Models, VDMs differ in their demographic assumptions, with no prior quantitative comparison of their performance. We benchmarked nine VDMs against observational data from boreal, temperate and tropical sites, assessing their accuracy in predicting tree growth, carbon turnover, biomass stocks and size distributions. Models were simulated under consistent climate conditions with postdisturbance recovery monitored for at least 420 yr. Postdisturbance carbon recovery trajectories showed significant variability while remaining within observational ranges. Initial regrowth rates varied substantially (0.03-0.60, 0.18-0.70 and 0.35-1.10 kgCm-2 yr-1 for boreal, temperate and tropical sites, respectively), influenced by each model's initial forest state. Models captured mature forest carbon content but showed compensating effects between overestimated growth and underestimated mortality rates. This first multi-model benchmarking identifies growth and mortality rates as critical calibration targets and highlights the need to refine postdisturbance establishment conditions for model development. We outline specific benchmarking variables needed to improve predictions of forest responses to environmental change.

demographic vegetation model benchmarking↗

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗

Hemispheric Asymmetry of Phase Partition in Mixed‐Phase Clouds Based on Near Global‐Scale Airborne Observations

Mixed-phase clouds contribute to substantial uncertainties in global climate models due to their complex microphysical properties. Former model evaluations almost exclusively rely on satellite observations to assess cloud phase distributions globally. This study investigated mixed-phase cloud properties using near global-scale in situ observation data sets from 14 flight campaigns in combination with collocated output from a global climate model. The Southern Hemisphere (SH) shows significantly higher occurrence frequencies and higher mass fractions of supercooled liquid water than Northern Hemisphere (NH) based on observations at 0.2 and 100 km horizontal scales. Such hemispheric asymmetry is not captured by the model. The model also consistently overestimates liquid water content (LWC) in all cloud phases but shows ice water content (IWC) biases that vary with phase. Key processes contributing to model biases in phase partition can be identified through the combination of evaluation of phase frequency, liquid mass fraction, LWC and IWC.

Yang, Ching An [San Jose State University, CA (Uni↗

Increased drought and extreme events over continental United States under high emissions scenario

The frequency, severity, and extent of climate extremes in future will have an impact on human well-being, ecosystems, and the effectiveness of emissions mitigation and carbon sequestration strategies. The specific objectives of this study were to downscale climate data for US weather stations and analyze future trends in meteorological drought and temperature extremes over continental United States (CONUS). We used data from 4161 weather stations across the CONUS to downscale future precipitation projections from three Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase Six (CMIP6), specifically for the high emission scenario SSP5 8.5. Comparing historic observations with climate model projections revealed a significant bias in total annual precipitation days and total precipitation amounts. The average number of annual precipitation days across CONUS was projected to be 205 ± 26, 184 ± 33, and 181 ± 25 days in the BCC, CanESM, and UKESM models, respectively, compared to 91 ± 24 days in the observed data. Analyzing the duration of drought periods in different ecoregions of CONUS showed an increase in the number of drought months in the future (2023–2052) compared to the historical period (1989–2018). The analysis of precipitation and temperature changes in various ecoregions of CONUS revealed an increased frequency of droughts in the future, along with longer durations of warm spells. Eastern temperate forests and the Great Plains, which encompass the majority of CONUS agricultural lands, are projected to experience higher drought counts in the future. Drought projections show an increasing trend in future drought occurrences due to rising temperatures and changes in precipitation patterns. Our high-resolution climate projections can inform policy makers about the hotspots and their anticipated future trajectories.

54 ENVIRONMENTAL SCIENCES↗

On the role of steepened structures and drift waves in equatorial spread F

A description is presented of a unified model of equatorial spread F which organizes much of the observational data and theoretical results in a self-consistent way. It is suggested that an initial bottomside instability nonlinearly evolves to the point that very steep gradients develop and that these steepened structures are responsible for the power spectra observed by probe experiments. This primary process is very likely the Rayleigh-Taylor instability. Observations indicate that some of the bottomside structures rise into the topside. It is shown that the observed gradients on the edges of these rising bubbles should be unstable to drift waves with a linear growth rate of about 1 per second and a perpendicular wavelength, for maximum growth, near the ion gyroradius.

Costa, E.↗

Estimation of effective hydrologic properties of soils from observations of vegetation density

An existing one-dimensional model of the annual water balance is reviewed. Slight improvements are made in the method of calculating the bare soil component of evaporation, and in the way surface retention is handled. A natural selection hypothesis, which specifies the equilibrium vegetation density for a given, water limited, climate-soil system, is verified through comparisons with observed data and is employed in the annual water balance of watersheds in Clinton, Ma., and Santa Paula, Ca., to estimate effective areal average soil properties. Comparison of CDF's of annual basin yield derived using these soil properties with observed CDF's provides excellent verification of the soil-selection procedure. This method of parameterization of the land surface should be useful with present global circulation models, enabling them to account for both the non-linearity in the relationship between soil moisture flux and soil moisture concentration, and the variability of soil properties from place to place over the Earth's surface.

Tellers, T. E.↗