Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “threshold”

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,045 records · Page 58

Radiometry Measurements of Mars at 1064 nm Using the Mars Orbiter Laser Altimeter

Measurements by the Mars Orbiter Laser Altimeter (MOLA) on board the Mars Global Surveyor (MGS) may be used to provides a radiometric measurement of Mars in addition to the topographic measurement. We will describe the principle of operation, a mathematical model, and the receiver calibration in this presentation. MOLA was designed primarily to measure Mars topography, surface roughness end the bidirectional reflectance to the laser beam. To achieve the highest sensitivity the receiver detection threshold is dynamically adjusted to be as low as possible while keeping a predetermined false alarm rate. The average false alarm rate 29 monitored in real time on board MOLA via a noise counter, whose output is fed to the threshold control loop. The false alarm rate at a given threshold is a function of the detector output noise which is the sum of the photo detector, shot noise due to the background light seen by the detector and the dark noise. A mathematical model has been developed that can be used to numerically solve for the optical background power given the MOLA threshold setting and the average noise count. The radiance of Mars can then be determined by dividing the optical power by the solid angle subtended by the MOLA receiver, the receiver optical band-width, end the Mars surface area within the receiver field of view. The phase angle which is the sun-Mars-MOLA angle is available from the MGS database. MOLA also measures simultaneously the bidirectional reflectance of Mars vie its 106-lum loser beam at nadir with nearly zero phase angle. The optical bandwidth of the MOLA receiver is 2um full width at half maximum (FWHM) and centered at 106-lum. The receiver field of view is 0.95mrad FWHM. The nominated spacecraft altitude is 100km and the ground track speed is about 3km/s. Under normal operation, the noise counter are read and the threshold levels are updated at 1Hz. The receiver sensitivity is limited by the detector dark noise to about 0.1nW, which corresponds to less than 2% the maximum radiance during daytime from the brightest area on Mars. The results from the mathematical model agree well with the prelaunch measurements at several calibrated optical power levels. The radiance of sunlit Mars estimated with this technique correlates well with the measurement from the MGS. Thermal Emission Spectrometer (TES) and the Hubble Space Telescope at similar wavelength.

Sun, Xiao-Li↗

Visibility of Wavelet Quantization Noise

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp)-L , where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We describe a mathematical model to predict DWT noise detection thresholds as a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗

Parafoveal Target Detectability Reversal Predicted by Local Luminance and Contrast Gain Control

This project is part of a program to develop image discrimination models for the prediction of the detectability of objects in a range of backgrounds. We wanted to see if the models could predict parafoveal object detection as well as they predict detection in foveal vision. We also wanted to make our simplified models more general by local computation of luminance and contrast gain control. A signal image (0.78 x 0.17 deg) was made by subtracting a simulated airport runway scene background image (2.7 deg square) from the same scene containing an obstructing aircraft. Signal visibility contrast thresholds were measured in a fully crossed factorial design with three factors: eccentricity (0 deg or 4 deg), background (uniform or runway scene background), and fixed-pattern white noise contrast (0%, 5%, or 10%). Three experienced observers responded to three repetitions of 60 2IFC trials in each condition and thresholds were estimated by maximum likelihood probit analysis. In the fovea the average detection contrast threshold was 4 dB lower for the runway background than for the uniform background, but in the parafovea, the average threshold was 6 dB higher for the runway background than for the uniform background. This interaction was similar across the different noise levels and for all three observers. A likely reason for the runway background giving a lower threshold in the fovea is the low luminance near the signal in that scene. In our model, the local luminance computation is controlled by a spatial spread parameter. When this parameter and a corresponding parameter for the spatial spread of contrast gain were increased for the parafoveal predictions, the model predicts the interaction of background with eccentricity.

Ahumada, Albert J., Jr.↗

Visibility of wavelet quantization noise

The discrete wavelet transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that we call DWT uniform quantization noise; it is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2-lambda, where r is display visual resolution in pixels/degree, and lambda is the wavelet level. Thresholds increase rapidly with wavelet spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from lowpass to horizontal/vertical to diagonal. We construct a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

NASA Center ARC↗

Binary-Signal Recovery

Binary communication through long cables, opto-isolators, isolating transformers, or repeaters can become distorted in characteristic ways. The usual solution is to slow the communication rate, change to a different method, or improve the communication media. It would help if the characteristic distortions could be accommodated at the receiving end to ease the communication problem. The distortions come from loss of the high-frequency content, which adds slopes to the transitions from ones to zeroes and zeroes to ones. This weakens the definition of the ones and zeroes in the time domain. The other major distortion is the reduction of low frequency, which causes the voltage that defines the ones or zeroes to drift out of recognizable range. This development describes a method for recovering a binary data stream from a signal that has been subjected to a loss of both higher-frequency content and low-frequency content that is essential to define the difference between ones and zeroes. The method makes use of the frequency structure of the waveform created by the data stream, and then enhances the characteristics related to the data to reconstruct the binary switching pattern. A major issue is simplicity. The approach taken here is to take the first derivative of the signal and then feed it to a hysteresis switch. This is equivalent in practice to using a non-resonant band pass filter feeding a Schmitt trigger. Obviously, the derivative signal needs to be offset to halfway between the thresholds of the hysteresis switch, and amplified so that the derivatives reliably exceed the thresholds. A transition from a zero to a one is the most substantial, fastest plus movement of voltage, and therefore will create the largest plus first derivative pulse. Since the quiet state of the derivative is sitting between the hysteresis thresholds, the plus pulse exceeds the plus threshold, switching the hysteresis switch plus, which re-establishes the data zero to one transition, except now at the logic levels of the receiving circuit. Similarly, a transition from a one to a zero will be the most substantial and fastest minus movement of voltage and therefore will create the largest minus first derivative pulse. The minus pulse exceeds the minus threshold, switching the hysteresis switch minus, which re-establishes the data one to zero transition. This innovation has a large increase in tolerance for the degradation of the binary pattern of ones and zeroes, and can reject the introduction of noise in the form of low frequencies that can cause the voltage pattern to drift up or down, and also higher frequencies that are beyond the recognizable content in the binary transitions.

Griebeler, Elmer L.↗

Tropical Cyclone Diurnal Cycle as Observed by TRMM

Using infrared satellite data, previous work has shown a consistent diurnal cycle in the pattern of cold cloud tops around mature tropical cyclones. In particular, an increase in the coverage by cold cloud tops often occurs in the inner core of the storm around the time of sunset and subsequently propagates outward to several hundred kilometers over the course of the following day. This consistent cycle may have important implications for structure and intensity changes of tropical cyclones and the forecasting of such changes. Because infrared satellite measurements are primarily sensitive to cloud top, the goal of this study is to use passive and active microwave measurements from the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) and Precipitation Radar (PR), respectively, to examine and better understand the tropical cyclone diurnal cycle throughout a larger depth of the storm's clouds. The National Hurricane Center's best track dataset was used to extract all PR and TMI pixels within 1000 km of each tropical cyclone that occurred in the Atlantic basin between 1998-2011. Then the data was composited according to radius (100-km bins from 0-1000 km) and local standard time (LST; 3-hr bins). Specifically, PR composites involved finding the percentage of pixels with reflectivity greater than or equal to 20 dBZ at various heights (i.e., 2-14 km in increments of 2 km) as a function of radius and time. The 37- and 85- GHz TMI channels are especially sensitive to scattering by precipitation-sized ice in the mid to upper portions of clouds. Hence, the percentage of 37- and 85-GHz polarization corrected temperatures less than various thresholds were calculated using data from all storms as a function of radius and time. For 37 GHz, thresholds of 260 K, 265 K, 270 K, and 275 K were used, and for 85 GHz, thresholds of 200-270 K in increments of 10 K were utilized. Note that convection forced by the interactions of a tropical cyclone with land (e.g., due to frictional convergence) may disrupt the natural convective cycle of a cyclone. Hence, only data pertaining to storms whose centers were greater than 300 km from land were included in the composites. Early results suggest the presence of a diurnal cycle in the PR composites of all Atlantic basin tropical cyclones from a height of 2-12 km from approximately 0-400 km radius, but the cycle is most apparent above 6 km. At a height of 8 km, there is a peak (minimum) in the percentage of PR pixels greater than or equal to 20 dBZ near 0 (21) LST in the inner core with some indication that this signal propagates outward with time. In contrast, the 37- and 85-GHz composites show little indication of a diurnal cycle at any radii, regardless of the threshold used. Ongoing work with this project will involve sub-setting the composites according to storm intensity to see if the diurnal cycle varies with storm strength. Moderate to strong vertical wind shear often leads to asymmetries in tropical cyclone convection and may disrupt the cyclone's natural diurnal cycle. Therefore, wind shear thresholds will be applied to the composites to determine if the diurnal cycle becomes more apparent in a low shear environment. Finally, other work to be completed will involve developing composites for other tropical cyclone basins, including the East Pacific, Northwest Pacific, South Pacific, and Indian Ocean.

tropical cyclones↗

Relationships Between Vestibular Measures as Potential Predictors for Spaceflight Sensorimotor Adaptation

Introduction: During space exploration missions astronauts are exposed to a series of novel sensorimotor environments, requiring sensorimotor adaptation. Until adaptation is complete, sensorimotor decrements occur, affecting critical tasks such as piloted landing or docking. Of particularly interest are locomotion tasks such as emergency vehicle egress or extra-vehicular activity. While nearly all astronauts eventually adapt sufficiently, it appears there are substantial individual differences in how quickly and effectively this adaptation occurs. These individual differences in capacity for sensorimotor adaptation are poorly understood. Broadly, we aim to identify measures that may serve as pre-flight predictors of and individual's adaptation capacity to spaceflight-induced sensorimotor changes. As a first step, since spaceflight is thought to involve a reinterpretation of graviceptor cues (e.g. otolith cues from the vestibular system) we investigate the relationships between various measures of vestibular function in humans. Methods: In a set of 15 ground-based control subjects, we quantified individual differences in vestibular function using three measures: 1) ocular vestibular evoked myogenic potential (oVEMP), 2) computerized dynamic posturography and 3) vestibular perceptual thresholds. oVEMP responses are elicited using a mechanical stimuli approach. Computerized dynamic posturography was used to quantify Sensory Organization Tests (SOTs), including SOT5M which involved performing pitching head movements while balancing on a sway-reference support surface with eyes closed. We implemented a vestibular perceptual threshold task using the tilt capabilities of the Tilt-Translation Sled (TTS) at JSC. On each trial, the subject was passively roll-tilted left ear down or right ear down in the dark and verbally provided a forced-choice response regarding which direction they felt tilted. The motion profile was a single-cycle sinusoid of angular acceleration with a duration of 5 seconds (frequency of 0.2 Hz), which was selected as it requires sensory integration of otolith and semicircular canal cues. Stimuli direction was randomized and magnitude was determined using an adaptive sampling procedure. One hundred trials were provided and each subject's responses were fit with a psychometric curve to estimate the subject's threshold. Results: Roll tilt perceptual thresholds at 0.2 Hz ranged from 0.5 degrees to 1.82 degrees across the 15 subjects (geometric mean of 1.04 degrees), consistent with previous studies. The inter-individual variability in thresholds may be able to help explain individual differences observed in sensorimotor adaptation to spaceflight. Analysis is ongoing for the oVEMPS and computerized dynamic posturography to identify relationships between the various vestibular measures. Discussion: Predicting individual differences in sensorimotor adaptation is critical both for the development of personalized countermeasures and mission planning. Here we aim to develop a basis of vestibular tests and parameters which may serve as predictors of individual differences in sensorimotor adaptability through studying the relationship between these measures.

Clark, T. K.↗

Detection and Use of Load and Gage Output Repeats of Wind Tunnel Strain-Gage Balance Data

Criteria are discussed that may be used for the detection of load and gage output repeats of wind tunnel strain-gage balance data. First, empirical thresholds are introduced that help determine if the loads or electrical outputs of a pair of balance calibration or check load data points match. A threshold of 0.01 percent of the load capacity is suggested for the identification of matching loads. Similarly, a threshold of 0.1 microV/V is recommended for the identification of matching electrical outputs. Two examples for the use of load and output repeats are discussed to illustrate benefits of the implementation of a repeat point detection algorithm in a balance data analysis software package. The first example uses the suggested load threshold to identify repeat data points that may be used to compute pure errors of the balance loads. This type of analysis may reveal hidden data quality issues that could potentially be avoided by making calibration process improvements. The second example uses the electrical output threshold for the identification of balance fouling. Data from the calibration of a six-component force balance is used to illustrate the calculation of the pure error of the balance loads.

strain-gage balance↗

Improved Detection of Low Stratus and Fog at Dawn from Dual Geostationary (COMS and FY-2D) Satellites

A novel method was proposed to detect low stratus and fog (LSF) at dawn during the summertime using near simultaneous observations from dual geostationary-orbit satellites (GEOs): the Korean Communication, Ocean and Meteorological Satellite (COMS; 128.2 degrees East) and the Chinese Feng-Yun-2D (FY-2D; 86.5 degrees East). The orbital positions of the GEOs provided a large difference (around 46.5 degrees) in the viewing zenith angle (VZA) in the study region (122-132 degrees East, 32.5-42.5 degrees North) and high contrast observations at dawn. Numerical simulations were carried out to derive the optical properties of the LSF and the radiative differences between the GEOs due to both the VZA and spectral response function (SRF). The conventional threshold method, which used the visible reflectance at 0.67 microns (R (sub 0.67)) and the brightness temperature difference between 3.7 microns and 11 microns (BTD (sub 3.7-11), had limitations in detecting LSF at dawn due to weak visible light and fast night-to-day transition. We utilized the observed stereo differences (delta R (sub 0.67), delta BTD (sub 3.7-11) between the two GEOs to improve LSF detection. The dual-satellite observations were verified with ground-based data from 45 stations in South Korea co-located with each GEO pixel. The delta R (sub 0.67) threshold value showed better accuracy (78 percent vs. 67 percent) than the conventional R (sub 0.67) threshold method. In addition, the delta BTD (sub 3.7-11) threshold was better (55 percent vs. 38 percent) than the BTD (sub 3.7-11) threshold. The dual-satellite method allowed more reliable LSF detection using the combination of delta R (sub 0.67) and R (sub 0.67), particularly for LSF without cumuliform or high clouds. Our method is applicable to multiple geostationary satellites for continuous LSF monitoring.

geostationary satellite↗

Multi-phase covetic and methods of synthesis thereof

There are provided methods and systems for creating multi-phase covetics. For example, there is provided a process for making a composite material. The process includes forming a multi-phase covetic. The forming includes heating a melt including a metal in a molten state and a carbon source to a first temperature threshold to form metal-carbon bonds. The forming further includes subsequently heating the melt to a second temperature threshold, the second temperature threshold being greater than the first temperature threshold. The second temperature threshold is a temperature at or above which ordered multi-phase covetics form in the melt.

36 MATERIALS SCIENCE↗

Multi-phase covetic and methods of synthesis thereof

There are provided methods and systems for creating multi-phase covetics. For example, there is provided a process for making a composite material. The process includes forming a multi-phase covetic. The forming includes heating a melt including a metal in a molten state and a carbon source to a first temperature threshold to form metal-carbon bonds. The forming further includes subsequently heating the melt to a second temperature threshold, the second temperature threshold being greater than the first temperature threshold. The second temperature threshold is a temperature at or above which ordered multi-phase covetics form in the melt.

36 MATERIALS SCIENCE↗

Methods and system for a battery

Systems and methods for operating a redox flow battery system may include switching the redox flow battery system to an idle mode, wherein the idle mode includes operation of the redox flow battery system outside of a charging mode and outside of a discharge mode; in response to switching to the idle mode, repeatedly cycling operation of an electrolyte pump between an idling threshold flow rate less than a charging threshold flow rate and a deactivation threshold flow rate; and in response to switching to the charging mode, maintaining operation of the electrolyte pump at the charging threshold flow rate greater than the idling threshold flow rate. In this way, a responsiveness of the redox flow battery system to charging and discharging commands can be maintained while in idle, while reducing parasitic pumping losses due to pumping and heating, and reducing shunt current losses.

Song, Yang↗

Thermal stress in degraded forests in the Brazilian Amazon Arc of Deforestation

Understanding thermal stress in tropical forests has taken on new urgency in light of accelerating climate change and expansion of deforestation and forest degradation. Degraded tropical forests in particular may be approaching critical temperature thresholds even more rapidly than intact forests, with implications for tree survival and ecosystem recovery. We investigate thermal stress in degraded tropical forests within the Brazilian Amazon Arc of Deforestation. Using land surface temperature data from the ECOsystem Spaceborne Thermal Radiometer Experiment on the international Space Station (ECOSTRESS), we compared canopy temperatures of intact, selectively logged, and burned forests in Feliz Natal, Mato Grosso, Brazil. Upper canopy temperatures in previously burned forests were 4.1% higher (mean = 36.5 °C) and 50.9% more variable compared to intact and logged forests, which showed remarkably similar temperature distributions (means of 34.9 °C and 35.1 °C, respectively). Modeled leaf temperature distributions based on canopy temperature measurements from one of the warmest days in a 2 year dry season record indicated that 87% of leaves in the warmest burned forest patches exceeded the temperature threshold where respiration surpasses photosynthesis, compared to approximately 72%–74% in intact and logged forests, respectively, in the same time period. After controlling for environmental factors, burned forests were predicted to be 2.6 °C warmer on average [1.39 °C–3.96 °C, 95% credible interval] than intact forests across a 5–40 m canopy height range. Burned forests showed modest thermal recovery over time, with temperatures decreasing by approximately 1.2 °C over a 30 year recovery period. In contrast, logged forests showed minimal thermal differences from intact forests (−0.353 on average, [−1.125–0.274, 95% credible interval]) and negligible change (0.15 °C) over the same timeframe. While the absolute probabilities of exceeding damaging thermal thresholds remain low across all forest types under current climate conditions, the probability of leaves exceeding temperatures that cause permanent leaf damage was ten times higher in burned forests, with implications for the future of burned forest regeneration in water-limited regions of the Amazon basin. In particular, the combination of higher mean temperatures, greater temperature variability, and more frequent exposure to damaging thermal thresholds implies that burned tropical forests will experience substantially higher mortality rates and slower biomass recovery compared to intact and selectively logged forests, especially in water-limited regions where trees cannot rely on evaporative cooling to moderate canopy temperatures.

Brazilian Amazon Arc of Deforestation↗

The adequate stimulus for avian short latency vestibular responses to linear translation

Transient linear acceleration stimuli have been shown to elicit eighth nerve vestibular compound action potentials in birds and mammals. The present study was undertaken to better define the nature of the adequate stimulus for neurons generating the response in the chicken (Gallus domesticus). In particular, the study evaluated the question of whether the neurons studied are most sensitive to the maximum level of linear acceleration achieved or to the rate of change in acceleration (da/dt, or jerk). To do this, vestibular response thresholds were measured as a function of stimulus onset slope. Traditional computer signal averaging was used to record responses to pulsed linear acceleration stimuli. Stimulus onset slope was systematically varied. Acceleration thresholds decreased with increasing stimulus onset slope (decreasing stimulus rise time). When stimuli were expressed in units of jerk (g/ms), thresholds were virtually constant for all stimulus rise times. Moreover, stimuli having identical jerk magnitudes but widely varying peak acceleration levels produced virtually identical responses. Vestibular response thresholds, latencies and amplitudes appear to be determined strictly by stimulus jerk magnitudes. Stimulus attributes such as peak acceleration or rise time alone do not provide sufficient information to predict response parameter quantities. Indeed, the major response parameters were shown to be virtually independent of peak acceleration levels or rise time when these stimulus features were isolated and considered separately. It is concluded that the neurons generating short latency vestibular evoked potentials do so as "jerk encoders" in the chicken. Primary afferents classified as "irregular", and which traditionally fall into the broad category of "dynamic" or "phasic" neurons, would seem to be the most likely candidates for the neural generators of short latency vestibular compound action potentials.

NASA Discipline Neuroscience↗

Evaluation of the rheological and electrical percolation of high‐density polyethylene/carbon black composites using mathematical models

Abstract In this work, conductive polymer composites (CPCs) of bio‐based polyethylene (BioPe) containing different concentrations of carbon black (CB) were developed. By using oscillatory rheology analysis, a Newtonian plateau was observed in BioPe, and all BioPe/CB composites had a behavior of a pseudo‐solid and that composites with volume fractions ranging from 0.24 to 0.56 presented higher viscosity, storage, and loss modulus. This suggests the formation of a percolated network and by using the power‐law models, it was observed that the electrical percolation threshold was higher than the rheological percolation threshold. The electrical conductivity was measured using the four‐point probe method and a sigmoid model was used to predict the CPCs' electrical conductivity percolation threshold. The results indicated that the four‐point probe method presented satisfactory results according to the calculated standard deviations and voltage–current characteristics for each round of measurements considering the same ranging as used in rheology analysis. The analytical model used showed a coefficient of determination ( R 2 ) higher than 95%, allowing the prediction of the electrical conductivity of the CPC and the percolation threshold as a function of the volumetric fraction of the CB.

da SiIva, Moacy P.↗

Opportunities for Earth Observation to Inform Risk Management for Ocean Tipping Points

Abstract As climate change continues, the likelihood of passing critical thresholds or tipping points increases. Hence, there is a need to advance the science for detecting such thresholds. In this paper, we assess the needs and opportunities for Earth Observation (EO, here understood to refer to satellite observations) to inform society in responding to the risks associated with ten potential large-scale ocean tipping elements: Atlantic Meridional Overturning Circulation; Atlantic Subpolar Gyre; Beaufort Gyre; Arctic halocline; Kuroshio Large Meander; deoxygenation; phytoplankton; zooplankton; higher level ecosystems (including fisheries); and marine biodiversity. We review current scientific understanding and identify specific EO and related modelling needs for each of these tipping elements. We draw out some generic points that apply across several of the elements. These common points include the importance of maintaining long-term, consistent time series; the need to combine EO data consistently with in situ data types (including subsurface), for example through data assimilation; and the need to reduce or work with current mismatches in resolution (in both directions) between climate models and EO datasets. Our analysis shows that developing EO, modelling and prediction systems together, with understanding of the strengths and limitations of each, provides many promising paths towards monitoring and early warning systems for tipping, and towards the development of the next generation of climate models.

Wood, Richard A. (ORCID:0000000239609513)↗

Trace elements and heavy metals in black vultures ( Coragyps atratus ) and turkey vultures ( Cathartes aura ) in the southeastern United States

Many vulture species worldwide are declining at alarming rates due to a variety of anthropogenic causes, including exposure to pollutants and pharmaceuticals through consumption of contaminated carrion. However, little is known about the extent to which vultures are exposed to various contaminants as well as toxicity thresholds for trace elements and heavy metals. Our objective was to quantify levels of trace elements and heavy metals within black vulture (Coragyps atratus) and turkey vulture (Cathartes aura) tissues to determine the extent to which populations in the southeastern United States are exposed to carrion that contains high levels of contaminants. We collected black vulture liver samples and examined differences in trace element and heavy metal concentrations between sexes and age classes (adult and juvenile). Further, we collected 81 blood and 42 feather samples from additional black and turkey vultures and compared differences between species and age classes. We found similar element concentrations between juvenile and adult black vultures with the exception of Cu, where levels 40 were higher in juveniles compared to adults. However, we did observe substantial differences in element concentrations between species for both blood and feather samples, with black vultures generally having higher concentrations of most elements. Our data revealed higher element levels in both species compared to toxicity thresholds found in other bird of prey species, such as blood and liver toxicity threshold suggestions for Pb poisoning in Falconiformes. Further, while average contaminant levels were generally low, extreme outliers were observed for some elements, including Pb, suggesting some individuals were exposed to high levels of potentially toxic elements. In conclusion, more research is needed to better understand contaminant exposure in black and turkey vultures across a broader geographic region, as well as elucidate toxicity thresholds and non-lethal impacts of contaminant exposure in these species.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗