Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

Timing Studies of X Persei and the Discovery of Its Transient Quasi-periodic Oscillation Feature

We present a timing analysis of X Persei (X Per) using observations made between 1998 and 2010 with the Proportional Counter Array (PCA) onboard the Rossi X-ray Timing Explorer (RXTE) and with the INTEGRAL Soft Gamma-Ray Imager (ISGRI). All pulse arrival times obtained from the RXTE-PCA observations are phase-connected and a timing solution is obtained using these arrival times. We update the long-term pulse frequency history of the source by measuring its pulse frequencies using RXTE-PCA and ISGRI data. From the RXTEPCA data, the relation between the frequency derivative and X-ray flux suggests accretion via the companion's stellar wind. However, the detection of a transient quasi-periodic oscillation feature, peaking at approximately 0.2 Hz, suggests the existence of an accretion disc. We find that doublebreak models fit the average power spectra well, which suggests that the source has at least two different accretion flow components dominating the overall flow. From the power spectrum of frequency derivatives, we measure a power-law index of approximately − 1, which implies that, on short time-scales, disc accretion dominates over noise, while on time-scales longer than the viscous time-scales, the noise dominates. From pulse profiles, we find a correlation between the pulse fraction and the count rate of the source.

accretion↗

Application of Whole Genome Sequencing and MALDI-TOF to the Identifcation of Bacillus Species Isolated from Cleanrooms at NASA Johnson Space Center

Astromaterial cleanrooms at NASA Johnson Space Center are built environments that hold samples, such as lunar rocks, from different space exploration missions. Bacillus sp. are frequently detected in routine microbial monitoring of these facilities. Since this, and related genera, can form endospores that can withstand harsh conditions, they could contaminate astromaterials. This could confound searches for extraterrestrial life. Whole genome sequencing (WGS) is widely used for identifying bacterial strains and tracking their source; however, WGS is expensive and time consuming. Matrix-assisted laser desorption ionization– time of flight mass spectrometry (MALDI-TOF) shows promise as a low-cost, rapid method of identifying strains of bacteria, but few studies have compared this proteomics method to WGS. To evaluate a high throughput method of tracking the source of contamination of this built environment, WGS and MALDI-TOF was conducted on 18 bacterial strains isolated from surfaces in astromaterials cleanrooms. WGS identified 14 Bacillus, 2 Paenibacillus, 1 Solibacillus and 1 Alcaligenes strains. These isolates showed similarity to strains commonly observed in spacecraft assembly cleanrooms at other facilities. Cluster analysis of mass spectra generated by MALDI-TOF grouped strains together that were greater than 94% similar to each other in terms of amino acid sequences of single copy core genes, as assessed by WGS. This suggests that MALDI-TOF and WGS results are consistent with each other and MALDI-TOF can rapidly identify strains of Bacillus sp. isolated from cleanroom environments with a resolution comparable to WGS. Based on phylogenomic analysis, these results also suggest the presence of a cosmopolitan class of Bacillus sp. that are more likely to be found in cleanrooms and similar built environments than in natural systems.

Farnaz Mazhari↗

Analysis of shape of porous cooled medium for an imposed surface heat flux and temperature

The surface of a porous cooled medium is to be maintained at a specified design temperature while being subjected to uniform heating by an external source. An analytical method is given for determining the shape of the medium surface that will satisfy these boundary conditions. The analysis accounts for temperature dependent variations of fluid density and viscosity and for temperature dependent matrix thermal conductivity. The energy equation is combined with Darcy's law in such a way that a potential can be defined that satisfies Laplace's equation. All of the heat-transfer and flow quantities are expressed in terms of this potential. The determination of the shape of the porous cooled region is thereby reduced to a free-boundary problem such as in inviscid free jet theory. Two illustrative examples are carried out: a porous leading edge with coolant supplied through a slot and a porous cooled duct with a rectangular outer boundary.

Siegel, R.↗

Extension of Case's method to the two-region critical cylinder and its application as an analytic standard

The normal mode expansion technique is applied to the transformed monoenergetic integral transport equation to develop a solution for the rotationally invariant and axially infinite, critical, two-region cylinder with a finite outer reflector boundary. Isotropic scattering and identical neutron mean free paths in the core and reflector regions are assumed. The solution in terms of singular integral equations is obtained by applying a completeness theorem found for the singular eigenfunctions. Numerical results for a variety of core and reflector multiplying properties and reflector thickness are presented and compared with the results of other method. An example of this type of application is given in a study of approximations inherent in the neutronic design analysis of a small, fast-neutron-spectrum reactor concept proposed as a space power source.

Westfall, R. M.↗

Noncontact Determination of Antisymmetric Plate Wave Velocity in Ceramic Matrix Composites

High-temperature materials are of increasing importance in the development of more efficient engines and components for the aeronautics industry. In particular, ceramic matrix composite (CMC) and metal matrix composite (MMC) structures are under active development for these applications. The acousto-ultrasonic (AU) method has been shown to be useful for assessing mechanical properties in composite structures. In particular, plate wave analysis can characterize composites in terms of their stiffness moduli. It is desirable to monitor changes in mechanical properties that occur during thermomechanical testing and to monitor the health of components whose geometry or position make them hard to reach with conventional ultrasonic probes. In such applications, it would be useful to apply AU without coupling directly to the test surface. For a number of years, lasers have been under investigation as remote ultrasonic input sources and ultrasound detectors. The use of an ultrasonic transducer coupled through an air gap has also been under study. So far at the NASA Lewis Research Center, we have been more successful in using lasers as ultrasonic sources than as output devices. On the other hand, we have been more successful in using an air-coupled piezoelectric transducer as an output device than as an input device. For this reason, we studied the laser in/air-coupled-transducer out combination-using a pulsed NdYAG laser as the ultrasonic source and an air-coupled-transducer as the detector. The present work is focused on one of the AU parameters of interest, the ultrasonic velocity of the antisymmetric plate-wave mode. This easily identified antisymmetric pulse can be used to determine shear and flexure modulus. It was chosen for this initial work because the pulse arrival times are likely to be the most precise. The following schematic illustrates our experimental arrangement for using laser in/air-transducer out on SiC/SiC composite tensile specimens. The NdYAG pulse was directed downward by a 90 infrared prism to the top of the specimen, but at the edge of one end. An energy sensor measured a single pulse at 13 millijoules (mJ) before it passed through the prism, which attenuated 15 percent of its energy. It also provided an output trigger for the waveform time-delay synthesizer.

Kautz, Harold E.↗

Sparse Solutions for Single Class SVMs: A Bi-Criterion Approach

In this paper we propose an innovative learning algorithm - a variation of One-class nu Support Vector Machines (SVMs) learning algorithm to produce sparser solutions with much reduced computational complexities. The proposed technique returns an approximate solution, nearly as good as the solution set obtained by the classical approach, by minimizing the original risk function along with a regularization term. We introduce a bi-criterion optimization that helps guide the search towards the optimal set in much reduced time. The outcome of the proposed learning technique was compared with the benchmark one-class Support Vector machines algorithm which more often leads to solutions with redundant support vectors. Through out the analysis, the problem size for both optimization routines was kept consistent. We have tested the proposed algorithm on a variety of data sources under different conditions to demonstrate the effectiveness. In all cases the proposed algorithm closely preserves the accuracy of standard one-class nu SVMs while reducing both training time and test time by several factors.

Das, Santanu↗

Remote sensing of sea state by analysis of backscattered microwave phase fluctuations

The relative phase of a normally incident microwave signal scattered off the sea surface is a random function of space and time. The statistics of these random phase fluctuations contain information about the wave-height statistics. This study demonstrates that the wave-height spectra can be deduced directly in terms of the spectra of the phase fluctuations without recourse to inversion techniques even when surface roughness exceeds many wavelengths of the incident signal. In the particular case of a nadir-directed satellite-mounted microwave source operating in the backscatter mode, the wave-height and phase spatial spectra are proportional with a constant proportionality of (2k0) squared, where k0 is the wavenumber of the scattered signal.

Clifford, S. F.↗

Delta V errors due to attitude correction maneuvers of a spinning spacecraft

A statistical error analysis has been performed to evaluate the change in velocity of a spin stabilized spacecraft due to thruster induced attitude correction maneuvers. Analytical expressions were developed to determine the mean and variance of the velocity in terms of the nominal values and uncertainties of key spacecraft parameters. These parameters include thrust, thruster timing, thruster alignment, plume impingement force and estimated spacecraft mass. Typical values of these parameters and their uncertainties for NASA's Galileo spacecraft were used to evaluate these expressions. The results of the analysis allow the effect of attitude correction maneuvers to be evaluated with respect to the functional requirements on nongravitational acceleration. Other results include the identification of major sources of delta V error and studies of the sensitivity of delta V to variations in the error sources.

Zipursky, H. J.↗

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗

Ecological Impacts of the Space Shuttle Program at John F. Kennedy Space Center, Florida

The Space Shuttle Program was one of NASAs first major undertakings to fall under the environmental impact analysis and documentation requirements of the National Environmental Policy Act of 1969 (NEPA). Space Shuttle Program activities at John F. Kennedy Space Center (KSC) and the associated Merritt Island National Wildlife Refuge (MINWR) contributed directly and indirectly to both negative and positive ecological trends in the region through the long-term, stable expenditure of resources over the 40 year program life cycle. These expenditures provided support to regional growth and development in conjunction with other sources that altered land use patterns, eliminated and modified habitats, and contributed to cultural eutrophication of the Indian River Lagoon. At KSC, most Space Shuttle Program related actions were conducted in previously developed facilities and industrial areas with the exception of the construction of the shuttle landing facility (SLF) and the space station processing facility (SSPF). Launch and operations impacts were minimal as a result of the low annual launch rate. The majority of concerns identified during the NEPA process such as potential weather modification, acid rain off site, and local climate change did not occur. Launch impacts from deposition of HCl and particulates were assimilated as a result of the high buffering capacity of the system and low launch and loading rates. Metals deposition from exhaust deposition did not display acute impacts. Sub-lethal effects are being investigated as part of the Resource Conservation and Recovery Act (RCRA) regulatory process. Major positive Space Shuttle Program effects were derived from the adequate resources available at the Center to implement the numerous environmental laws and regulations designed to enhance the quality of the environment and minimize impacts from human activities. This included reduced discharges of domestic and industrial wastewater, creation of stormwater management systems, remediation of past contamination sites, implementation of hazardous waste management systems, and creation of a culture of sustainability. Working with partners such as the USFWS and the St Johns River Water Management District (SJRWMD), wetlands and scrub restoration and management initiatives were implemented to enhance fish and wildlife populations at the Center. KSC remains the single largest preserve on the east coast of Florida in part due to NASAs commitment to stewardship. Ongoing Ecological Program projects are directed at development of information and knowledge to address future KSC management questions including the transition to a joint government and commercial launch facility, enhanced habitat management requirements for wetlands and scrub, potential impacts of emerging contaminants, and adaptation to climate change including projected sea level rise over the next 50-75 years.

Space Shuttle Program↗

Accretion torques in X-ray pulsars

An analysis of the accretion process in an X-ray pulsar, whereby angular momentum is transferred to the star and its rotation period is changed, is presented, and an expression for the fractional rate of change of the pulse period in terms of X-ray luminosity and other star parameters is derived. It is shown that observed characteristic spin-up time scales for seven X-ray pulsars strongly support the view that in every source (1) the pulse period reflects the rotation period of a compact object, (2) the accretion is mediated by a disk surrounding the compact object and rotating in the same sense, and (3) the compact object is a neutron star rather than a white dwarf.

Rappaport, S.↗

Impact of Emissions and Long-Range Transport on Multi-Decadal Aerosol Trends: Implications for Air Quality and Climate

We present a global model analysis of the impact of long-range transport and anthropogenic emissions on the aerosol trends in the major pollution regions in the northern hemisphere and in the Arctic in the past three decades. We will use the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model to analyze the multi-spatial and temporal scale data, including observations from Terra, Aqua, and CALIPSO satellites and from the long-term surface monitoring stations. We will analyze the source attribution (SA) and source-receptor (SR) relationships in North America, Europe, East Asia, South Asia, and the Arctic at the surface and free troposphere and establish the quantitative linkages between emissions from different source regions. We will discuss the implications for regional air quality and climate change.

Chin, Mian↗

High charge state carbon and oxygen ions in Earth's equatorial quasi-trapping region

Observations of energetic (1.5 - 300 keV/e) medium-to-high charge state (+3 less than or equal to Q less than or equal to +7) solar wind origin C and O ions made in the quasi-trapping region (QTR) of Earth's magnetosphere are compared to ion trajectories calculated in model equatorial magnetospheric magnetic and electric fields. These comparisons indicate that solar wind ions entering the QTR on the nightside as an energetic component of the plasma sheet exit the region on the dayside, experiencing little or no charge exchange on the way. Measurements made by the CHarge Energy Mass (CHEM) ion spectrometer on board the Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft at 7 less than L less than 9 from September 1984 to January 1989 are the source of the new results contained herein: quantitative long-term determination of number densities, average energies, energy spectra, local time distributions, and their variation with geomagnetic disturbance level as indexed by Kp. Solar wind primaries (ions with charge states unchanged) and their secondaries (ions with generally lower charge states produced from primaries in the magnetosphere via charge exchange)are observed throughout the QTR and have distinctly different local time variations that persist over the entire 4-year analysis interval. During Kp larger than or equal to 3 deg intervals, primary ion (e.g., O(+6)) densities exhibit a pronounced predawn maximum with average energy minimum and a broad near-local-noon density minimum with average energy maximum. Secondary ion (e.g., O(+5)) densities do not have an identifiable predawn peak, rather they have a broad dayside maximum peaked in local morning and a nightside minimum. During Kp less than or equal to 2(-) intervals, primary ion density peaks are less intense, broader in local time extent, and centered near midnight, while secondary ion density local time variations diminish. The long-time-interval baseline helps to refine and extend previous observations; for example, we show that ionospheric contribution to O(+3)) is negligible. Through comparison with model ion trajectories, we interpret the lack of pronounced secondary ion density peaks colocated with the primary density peaks to indicate that: (1) negligible charge exchange occurs at L greater than 7, that is, solar wind secondaries are produced at L less than 7, and (2) solar wind secondaries do not form a significant portion of the plasma sheet population injected into the QTR. We conclude that little of the energetic solar wind secondary ion population is recirculated through the magnetosphere.

Christon, S. P.↗

Dissertation Defense Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional "validation by test only" mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

CFD↗

Dissertation Defense: Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional validation by test only mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions.Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations. This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System spacecraft system.Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Fluid Dynamics↗

NIAC Phase 1 Final Study Report on Titan Aerial Daughtercraft

Saturns giant moon Titan has become one of the most fascinating bodies in the Solar System. Even though it is a billion miles from Earth, data from the Cassini mission reveals that Titan has a very diverse, Earth-like surface, with mountains, fluvial channels, lakes, evaporite basins, plains, dunes, and seas [Lopes 2010] (Figure 1). But unlike Earth, Titans surface likely is composed of organic chemistry products derived from complex atmospheric photochemistry [Lorenz 2008]. In addition, Titan has an active meteorological system with observed storms and precipitation-induced surface darkening suggesting a hydrocarbon cycle analogous to Earths water cycle [Turtle 2011].Titan is the richest laboratory in the solar system for studying prebiotic chemistry, which makes studying its chemistry from the surface and in the atmosphere one of the most important objectives in planetary science [Decadal 2011]. The diversity of surface features on Titan related to organic solids and liquids makes long-range mobility with surface access important [Decadal 2011]. This has not been possible to date, because mission concepts have had either no mobility (landers), no surface access (balloons and airplanes), or low maturity, high risk, and/or high development costs for this environment (e,g. large, self-sufficient, long-duration helicopters). Enabling in situ mobility could revolutionize Titan exploration, similarly to the way rovers revolutionized Mars exploration. Recent progress on several fronts has suggested that small-scale rotorcraft deployed as daughtercraft from a lander or balloon mothercraft may be an effective, affordable approach to expanding Titan surface access. This includes rapid progress on autonomous navigation capabilities of such aircraft for terrestrial applications and on miniaturization, driven by the consumer mobile electronics market, of high performance of sensors, processors, and other avionics components needed for such aircraft. Chemical analysis, for example with a mass spectrometer, will be important to any Titan surface mission. Anticipating that it may be more practical to host chemical analysis instruments on a mothership than a daughtercraft, we defined system and mission concepts that deploy a small rotorcraft, termed a Titan Aerial Daughtercraft (TAD), from a lander or balloon to perform high-resolution imaging and mapping, potentially land to acquire microscopic images or other in situ measurements, and acquire samples to return to analytical instruments on the mothership. In principle, the ability to recharge batteries in TAD from a radioisotope or other long-lived power source on the mothership could enable multiple sorties. For a lander-based mission, a variety of landing sites is conceivable, including near lake margins, in dry lake beds, or in regions of plains, dunes, or putative cryovolanic or impact melt features. Such missions may require landing with greater precision than in previous missions (Huygens) and mission studies; this could also enhance the ability of TAD to reach interesting terrain from the landing site. Precision descent may also benefit balloon missions, with or without a daughtercraft, by increasing the probability that the balloon will drift over desired terrain early in its mission. Given these potential benefits, the overall concept studied here includes brief consideration of precision descent for landing or balloon deployment, followed by one or more sorties by a rotorcraft deployed from the mothership, with the ability to return to the mothership.

Saturn↗

Watershed Modeling with Remotely Sensed Big Data: MODIS Leaf Area Index Improves Hydrology and Water Quality Predictions

Traditional watershed modeling often overlooks the role of vegetation dynamics. There is also little quantitative evidence to suggest that increased physical realism of vegetation dynamics in process-based models improves hydrology and water quality predictions simultaneously. In this study, we applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershedin the midwestern United States. We modified the SWAT source code to automatically override the model’s built-in semiempirical LAI with spatially distributed and temporally continuous estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that the overestimation or underestimation of LAI imparted a proportional cascading effect on how the model partitions hydrologic fluxes and nutrient pools. As such, assimilation of MODIS LAI data corrected the model’sLAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study fills the existing knowledge gap regarding vegetation dynamics in watershed modeling and confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.

Adnan Rajib↗