SEARCH · Engineering Papers

Results for “source term estimation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Efficient method for approximating nonlinear dynamics: applications to uncertainty propagation and estimation

High-order Taylor series expansions can be used to model nonlinear dynamics at the cost of integrating a large set of variational equations to obtain high-order state-transition tensors (STTs). This paper presents an innovative technique for approximating the high-order STTs that reduces significantly the computational cost by retaining only the dominant secular terms. We propagate the low-order partial derivatives of Kepler’s equation, which only requires the integration of six additional equations to extend an n-th order approximation to order (n + 1). The approximation stems from the Lindstedt-Poincare procedure and exploits the stability properties of orbital motion. Since the method makes no dynamical assumptions, it can accommodate any source of orbital perturbations. We show how the approximation of the second-order STT significantly increases the accuracy of the linear method for uncertainty propagation with only a small computational overhead. Finally, we derive a high-order approximate extended Kalman filter that implements the proposed approximation of the STT and improves the performance of linear filters. Examples of application with different perturbation sources include the heliocentric orbit of an asteroid, an orbiter around Europa, and an Earth-orbiting satellite.

Park, Ryan S↗

DOE OSTI2022

Photometric redshift-aided classification using ensemble learning

We present SHEEP, a new machine learning approach to the classic problem of astronomical source classification, which combines the outputs from the XGBoost, LightGBM, and CatBoost learning algorithms to create stronger classifiers. A novel step in our pipeline is that prior to performing the classification, SHEEP first estimates photometric redshifts, which are then placed into the data set as an additional feature for classification model training; this results in significant improvements in the subsequent classification performance. SHEEP contains two distinct classification methodologies: (i) Multi-class and (ii) one versus all with correction by a meta-learner. We demonstrate the performance of SHEEP for the classification of stars, galaxies, and quasars using a data set composed of SDSS and WISE photometry of 3.5 million astronomical sources. The resulting F1 -scores are as follows: 0.992 for galaxies; 0.967 for quasars; and 0.985 for stars. In terms of the F1-scores for the three classes, SHEEP is found to outperform a recent RandomForest-based classification approach using an essentially identical data set. Our methodology also facilitates model and data set explainability via feature importances; it also allows the selection of sources whose uncertain classifications may make them interesting sources for follow-up observations.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1978

Mercury's helium exosphere after Mariner 10's third encounter

From a comparison of the Mariner 10 third encounter UV spectrometer data with intensities generated from a newly constructed model exosphere, a new value of 4.5 x 10 to the -4th power for the fraction of the solar wind He(++) flux to be intercepted and captured by Mercury's magnetosphere, if the observed He atmosphere is maintained by the solar wind, has been derived. If an internal source for He prevails, the corresponding upper bound for the global outgassing rate is estimated to be 4.5 x 10 to the 22nd power/s. These values differ from those given earlier because of the present use of a surface temperature distribution satisfying the heat equation over Mercury's entire surface which employs Mariner-10-determined mean surface thermal characteristics. The mean standoff distance of Mercury's magnetopause averaged over Mercury's orbit is also used. Agreement between the observed and calculated intensities is found to be good. Minor discrepancies on the nightside of the terminator are explicable in terms of differences between actual and computed temperatures and scale height structure changes.

Curtis, S. A.↗

NASA NTRS2012

Magnetization of the Lunar Crust

Magnetic fields measured by the satellite Lunar Prospector show large scale features resulting from remanently magnetized crust. Vector data synthesized at satellite altitude from a spherical harmonic model of the lunar crustal field, and the radial component of the magnetometer data, have been used to produce spatially continuous global magnetization models for the lunar crust. The magnetization is expressed in terms of localized basis functions, with a magnetization solution selected having the smallest root-mean square magnetization for a given fit to the data, controlled by a damping parameter. Suites of magnetization models for layers with thicknesses between 10 and 50 km are able to reproduce much of the input data, with global misfits of less than 0.5 nT (within the uncertainties of the data), and some surface field estimates. The magnetization distributions show robust magnitudes for a range of model thicknesses and damping parameters, however the magnetization direction is unconstrained. These global models suggest that magnetized sources of the lunar crust can be represented by a 30 km thick magnetized layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the measured magnetizations of the Apollo samples and significantly weaker than crustal magnetizations for Mars and the Earth. These are the first global magnetization models for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in the absence of multiple sample returns, and can be used to predict the crustal contribution to the lunar magnetic field at a particular location.

Lunar↗

DOE OSTI2026

Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES)

The C2SAGES project evaluated the feasibility of a community geothermal system for the planned Windy Ridge affordable housing development in Hinesburg, Vermont. Led by GTI Energy with Vermont Gas Systems, LN Consulting, NREL, and Frontier Energy, the work assessed technical design, energy performance, costs, business models, community engagement, maintenance, workforce development, and permitting. The proposed system was designed to serve 100% of the development’s heating, cooling, and domestic hot water loads. Compared with a baseline using air-source heat pumps and natural gas water heating, the geothermal system was estimated to reduce HVAC and domestic hot water energy use by about 45% to 48%, lower operating and maintenance costs, and reduce 30-year life-cycle costs by 37% for Phase 1 and 10% for Phase 2. Technical testing and modeling indicated that the Windy Ridge site is suitable for a community-scale geothermal system. The project also developed borehole field layouts, piping concepts, pump house designs, controls, maintenance plans, and supporting engineering drawings. The business model analysis found that first cost, ownership structure, and customer affordability remain major deployment challenges. Utility-led maintenance and operation were viewed favorably, but traditional utility cost-recovery models may require subsidy or revised financing structures to be practical for affordable housing. Community engagement highlighted the need for clear public education, transparent financing, reliable long-term maintenance, trained technicians, and the potential to pair geothermal systems with weatherization. Overall, the report concludes that community geothermal is technically feasible and offers meaningful energy, emissions, and life-cycle cost benefits, but broader deployment will depend on workable financing models and workforce readiness.

15 GEOTHERMAL ENERGY↗

NASA NTRS2019

Historical (1700–2012) Global Multi-Model Estimates of the Fire Emissions from the Fire Modeling Intercomparison Project (FireMIP)

Fire emissions are critical for carbon and nutrient cycles, climate, and air quality. Dynamic Global Vegetation Models (DGVMs) with interactive fire modeling provide important estimates for long-term and large-scale changes of fire emissions. Here we present the first multi-model estimates of global gridded historical fire emissions for 1700-2012, including carbon and 33 species of trace gases and aerosols. The dataset is based on simulations of nine DGVMs with different state-of-the-art global fire models that participated in the Fire Modeling Intercomparison Project (FireMIP), using the same and standardized protocols and forcing data, and the most up-to-date fire emission factor table from field and laboratory studies over various land cover types. We evaluate the simulations of present-day fire emissions by comparing them with satellite-based products. Evaluation results show that most DGVMs simulate present-day global fire emission totals within the range of satellite-based products, and can capture the high emissions over the tropical savannas, low emissions over the arid and sparsely vegetated regions, and the main features of seasonality. However, most of the models fail to simulate the interannual variability, partly due to a lack of modeling peat fires and tropical deforestation fires. Historically, all models show only a weak trend in global fire emissions before ~1850s, consistent with multi-source merged historical reconstructions. The long-term trends among DGVMs are quite different for the 20th century, with some models showing an increase and others a decrease in fire emissions, mainly as a result of the discrepancy in their simulated responses to human population density change and land-use and land-cover change (LULCC). Our study provides a basic dataset for developing regional and global multi-source merged historical reconstructions and merging methods, and analyzing historical changes of fire emissions and their uncertainties as well as their role in the Earth system. It also highlights the importance of accurately modeling the responses of fire emissions to LULCC and population density change in reducing uncertainties in historical reconstructions of fire emissions and providing more reliable future projections.

Li, Fang↗

NASA NTRS2012

SBUV version 8.6 Retrieval Algorithm: Error Analysis and Validation Technique

SBUV version 8.6 algorithm was used to reprocess data from the Back Scattered Ultra Violet (BUV), the Solar Back Scattered Ultra Violet (SBUV) and a number of SBUV/2 instruments, which 'span a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s)[see Bhartia et al, 2012]. In the new version Daumont et al. [1992] ozone cross section were used, and new ozone [McPeters et ai, 2007] and cloud climatologies Doiner and Bhartia, 1995] were implemented. The algorithm uses the Optimum Estimation technique [Rodgers, 2000] to retrieve ozone profiles as ozone layer (partial column, DU) on 21 pressure layers. The corresponding total ozone values are calculated by summing ozone columns at individual layers. The algorithm is optimized to accurately retrieve monthly zonal mean (mzm) profiles rather than an individual profile, since it uses monthly zonal mean ozone climatology as the A Priori. Thus, the SBUV version 8.6 ozone dataset is better suited for long-term trend analysis and monitoring ozone changes rather than for studying short-term ozone variability. Here we discuss some characteristics of the SBUV algorithm and sources of error in the SBUV profile and total ozone retrievals. For the first time the Averaging Kernels, smoothing errors and weighting functions (or Jacobians) are included in the SBUV metadata. The Averaging Kernels (AK) represent the sensitivity of the retrieved profile to the true state and contain valuable information about the retrieval algorithm, such as Vertical Resolution, Degrees of Freedom for Signals (DFS) and Retrieval Efficiency [Rodgers, 2000]. Analysis of AK for mzm ozone profiles shows that the total number of DFS for ozone profiles varies from 4.4 to 5.5 out of 6-9 wavelengths used for retrieval. The number of wavelengths in turn depends on solar zenith angles. Between 25 and 0.5 hPa, where SBUV vertical resolution is the highest, DFS for individual layers are about 0.5.

Kramarova, N. A.↗

NASA NTRS1994

Deep VLA imaging of twelve extended 3CR quasars

A representative sample of 12 extended quasars from the 3CR catalog has been imaged at 4.9 GHz using the Very Large Array (VLA). These full synthesis observations typically achieve an rms noise of 20 micro-Jy per beam. Jets are detected on at least one side of every source. The jets are well collimated compared with those in less powerfull sources, but spreading is detected in most of them. The opening angles of several jets are not constant, but show recollimation after an initial regime of rapid spreading. Many of the jets contain quasiperiodic strings of knots, of which the knot closest to the central feature is usually the brightest (until the jet nears its hot spot). The degrees of linear polarization at the jet knots range from less than 5% to approx. 50%, but show no common trend with distance along the jets. In knots that are elongated in directions close to that of the jet, the E vectors tend to be orthogonal to the jet axis. The prominence of the inner, straighter jet segments relative to the extended lobes correlates significantly with the prominence of the milliarcsecond-scale central features, but the prominence of the more bent jet segments does not. Candidates for counterjet emission are detected in seven sources, but there is no unambiguous, continuous counterjet in any of them. Estimates of the flux density ratios between the straighter jet segments and the counterjets based on these tentative detections range from 1.2:1 to greater than 175:1. There is no evidence in this sample that counterjet detectability correlates with such putative inclination indicators as central feature prominence or projected linear size. There is also no evidence that the prominence of the counterjets anticorrelates with that of the jets as predicted by simple relativistic-beaming models for the jet/counterjet asymmetry. There is, however, strong evidence that large bends in the main jet favor counterjet detection, and there are no counterjet candidates opposite long, uninterrupted straight segments of the main jets. The detectability of the counterjets in these quasars may therefore be strongly influenced by interactions between the underlying beams and inhomogeneities in the surrounding material. We offer a new empirical definition of the term 'hot spot' that is intended to improve the distinction between such features and 'jet knots.' Both the compactness of hot spots and their position in the lobe are affected by whether they are fed by a detectable jet. When the hot spots differ significantly in compactness, the more compact one is always on the jetted side. Jetted hot spots are also more likely to be recessed deeply from the outer edge of their lobes than are their counterjetted counterparts. The jetted hot spot is less prominent relative to other extended emission if the jet bends through a large angle, particularly if a large bend occurs abruptly. The counterjetted hot spot is also less well defined if the jet is more bent. The lobes of several sources show considerable inhomogeneity, including filamentation.

Bridle, Alan H.↗

NASA NTRS1975

The fluorescence of chlorophyll and yellow substances in natural waters: A note on the problems of measurement and the importance of their remote sensing

There are two chromophylls which, if sensed remotely from high altitude, would revolutionize the ability to survey large areas of the world's oceans. The chromophylls of importance are: the photosynthetic pigments of plankton algae and a group of organic materials frequently termed dissolved yellow substances. These are derived from plants and carried into the ocean by fresh water inflow. The attenuation of light by phytoplankton is characterized by two distinctive bands (450, 675 nm) which represent absorption by chloroplastic pigments. Yellow substances are characterized by a strong ultraviolet absorption which tails over into the visible region. It is emphasized that chlorophyll determination could be a unique technique for estimating the extent of eutrophication in coastal waters, and that a high altitude observer equipped with temperature, chlorophyll and yellow substance sensors has the possibility of detecting the magnitude of eutrophication and its sources by using laser induced fluorescent devices.

Yentsch, C. S.↗

NASA NTRS1997

Lateral Shock of the R Aquarii Jet

The R Aqr jet was observed with the VLA B-configuration at two epochs separated by approximately 13.2 yr. Comparison of the resulting 6 cm continuum images show that the radio jet has undergone a lateral counterclockwise rotation of approximately 6 deg-12 deg on the plane of the sky. The model of jet parcels on independent trajectories is difficult to reconcile with these observations and leads us to consider a path-oriented jet (i.e., younger parcels follow the same path as older parcels). Comparison of the most recent radio image with a nearly contemporaneous HST/FOC ultraviolet image at approximately 2330 Angstroms suggests that the ultraviolet emission lies along the leading side of the rotating radio jet. In conjunction with a proper motion analysis of the jet material that yields empirical space-velocity and resulting acceleration-magnitude relationships as a function of distance from the central source, we evaluate the observational results in terms of a schematic model in which the jet emission consists of plane-parallel isothermal shocks along the leading edge of rotation. In such a radiating shock, the ultraviolet-emitting region is consistent with the adiabatic region in the form of a high-temperature, low-density sheath that surrounds the cooled postshock radio-emitting region. Within the context of the schematic model, we obtain the temperatures, densities, and pressures within the preshock, adiabatic, and postshock regions as a function of distance from the central source; the physical parameters so derived compare favorably to previously published estimates. We obtain a total jet mass of 3.1 x 10(exp -5) solar mass and an age of approximately 115 yr. We evaluate the model in the context of its density-boundary condition, its applicability to an episodic or quasi-continuous jet, and angular momentum considerations.

Hollis, J. M.↗

NASA NTRS2014

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↗

NASA NTRS2014

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↗

DOE OSTI2023

Lagrangian displacement field estimators in cosmology

The late-time nonlinear Lagrangian displacement field is highly correlated with the initial field, so reconstructing it could enable us to extract primordial cosmological information. Our previous work carefully studied the displacement field reconstructed from the late-time density field using the iterative method proposed and found that it does not fully converge to the true, underlying displacement field (e.g., ~8 % offset at k ~ 0.2 h Mpc –1 at z = 0.6 ). We also constructed the Lagrangian perturbation theory model for the reconstructed field, but the model could not explain the discrepancy between the true and the reconstructed fields in the previous work. Here, the main sources of the discrepancy were speculated to be a numerical artifact in the displacement estimator due to the discreteness of the sample. In this paper, we develop two new estimators of the displacement fields to reduce such a numerical discreteness effect, the normalized momentum estimator and the rescaled resumed estimator. We show that the discrepancy reported is not due to the numerical artifacts. We conclude that the method cannot fully reconstruct the shape of the nonlinear displacement field at the redshift we studied, while it is still an efficient baryon acoustic oscillation reconstruction method. In parallel, by properly accounting for the UV-sensitive term in a reconstruction procedure with an effective field theory approach, we improve the theoretical model for the reconstructed displacement field, by almost 5 times, from ~15% to the level of a few percent at k ~ 0.2 h Mpc –1 at the redshift z = 0.6.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗

NASA NTRS1977

Solar seismology. I - The stability of the solar p-modes

The stability of the radial p-modes of the sun is investigated by computing nonadiabatic eigenvalues and eigenfunctions for a solar envelope model which extends from an inner radius of about 0.3 solar radius out to an optical depth of about 0.0003. The calculations take into account in a crude fashion the response of the convective flux to the oscillation. The dynamical effect of turbulence in the convection zone is parametrized in terms of a turbulent shear viscosity. The results show that if damping by turbulent viscosity is neglected, all modes with periods longer than 6 minutes are unstable. The familiar kappa-mechanism, which operates in the H ionization-H(-) opacity region, is the dominant source of driving of the oscillations. Modes with periods shorter than 6 minutes are stabilized by radiative damping in the solar atmosphere. When turbulent dissipation of pulsational energy is included, all modes are predicted to be stable. However, the margin of stability is very small. In view of the large uncertainty that must be assigned to the estimate of turbulent damping, it is concluded that theoretical calculations cannot unequivocally resolve the question of the stability of the solar p-modes.

Goldreich, P.↗

NASA NTRS1984

Friction- and mountain-torque estimates from global atmospheric data

Seasonal, zonal surface torques between the atmosphere and the earth are estimated and compared, using data from a number of independent sources. The mountain torque is computed both from surface pressure data and from isobaric height data. The friction torque is estimated from the oceanic stress data of Hellerman and Rosenstein. Results for the total torque are inferred from atmospheric angular momentum data. Finally, the globally integrated total torque is compared with astronomical observations of the earth's rotation rate. These comparisons help us to assess the quality of the different results. Zonal torques are also computed using results from a GFDL general circulation model of the atmosphere. A comparison with the corresponding results inferred from real data is presented and interpreted in terms of model accuracy.

Wahr, J. M.↗

NASA NTRS1986

Basaltic fissure eruptions, plume heights, and atmospheric aerosols

Convective plumes that rise above Hawaiian-style fire fountains consist of volcanic gases, aerosols, fine ash, and entrained heated air. Plume theory has been applied to observational estimates of the rate of thermal energy release from large fire fountains. The theoretically predicted heights of maintained plumes agree very well with the heights found from actual observations. Predicted plume heights for both central-vent (point-source) and fissure (line-source) eruptions indicate a stratospheric penetration by plumes that form over vents with very high magma-production rates. Flood basalt fissure eruptions that produce individual lava flows with volumes greater than 100 cu km at very high mass eruption rates are capable of injecting large quantities of sulfate aerosols into the lower stratosphere, with potentially drastic short-term atmospheric consequences, like acid precipitation, darkening of the sky, and climatic cooling.

Stothers, R. B.↗

NASA NTRS2010

Tidal Marsh Outwelling of Dissolved Organic Matter and Resulting Temporal Variability in Coastal Water Optical and Biogeochemical Properties

Coastal wetlands are highly dynamic environments at the land-ocean interface where human activities, short-term physical forcings and intense episodic events result in high biological and chemical variability. Long being recognized as among the most productive ecosystems in the world, tidally-influenced coastal marshes are hot spots of biogeochemical transformation and exchange. High temporal resolution observations that we performed in several marsh-estuarine systems of the Chesapeake Bay revealed significant variability in water optical and biogeochemical characteristics at hourly time scales, associated with tidally-driven hydrology. Water in the tidal creek draining each marsh was sampled every hour during several semi-diurnal tidal cycles using ISCO automated samplers. Measurements showed that water leaving the marsh during ebbing tide was consistently enriched in dissolved organic carbon (DOC), frequently by more than a factor of two, compared to water entering the marsh during flooding tide. Estimates of DOC fluxes showed a net DOC export from the marsh to the estuary during seasons of both low and high biomass of marsh vegetation. Chlorophyll amounts were typically lower in the water draining the marsh, compared to that entering the marsh during flooding tide, suggesting that marshes act as transformers of particulate to dissolved organic matter. Moreover, detailed optical and compositional analyses demonstrated that marshes are important sources of optically and chemically distinctive, relatively complex, high molecular weight, aromatic-rich and highly colored dissolved organic compounds. Compared to adjacent estuarine waters, marsh-exported colored dissolved organic matter (CDOM) was characterized by considerably stronger absorption (more than a factor of three in some cases), larger DOC-specific absorption, lower exponential spectral slope, larger fluorescence signal, lower fluorescence per unit absorbance, and higher fluorescence at visible wavelengths. Observed patterns in water optical and biogeochemical variables were very consistent among different marsh systems and throughout the year, despite continued tidal exchange, implying rapid transformation of marsh DOM in the estuary through both photochemical and microbial processes. These findings illustrate the importance of tidal marsh ecosystems as sources, sinks and/or transformers of biologically important nutrients, carbon and colored dissolved organic compounds, and their influence on short-term biological, optical and biogeochemical variability in coastal waters.

Tzortziou, Maria↗