Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical mode decomposition”

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 73 records · Page 4

An Iterated Global Mascon Solution with Focus on Land Ice Mass Evolution

Land ice mass evolution is determined from a new GRACE global mascon solution. The solution is estimated directly from the reduction of the inter-satellite K-band range rate observations taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons are estimated with 10-day and 1-arc-degree equal area sampling, applying anisotropic constraints for enhanced temporal and spatial resolution of the recovered land ice signal. The details of the solution are presented including error and resolution analysis. An Ensemble Empirical Mode Decomposition (EEMD) adaptive filter is applied to the mascon solution time series to compute timing of balance seasons and annual mass balances. The details and causes of the spatial and temporal variability of the land ice regions studied are discussed.

Luthcke, S. B.↗

Antarctica, Greenland and Gulf of Alaska Land-ice Evolution from an Iterated GRACE Global Mascon Solution

We have determined the ice mass evolution of the Antarctica and Greenland ice sheets (AIS and GIS) and Gulf of Alaska (GOA) glaciers from a new GRACE global solution of equal-area surface mass concentration parcels (mascons) in equivalent height of water. The mascons were estimated directly from the reduction of the inter-satellite K-band range-rate (KBRR) observations, taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons were estimated with 10 day and 1 arc degree equal-area sampling, applying anisotropic constraints. An ensemble empirical mode decomposition adaptive filter was applied to the mascon time series to compute annual mass balances. The details and causes of the spatial and temporal variability of the land-ice regions studied are discussed. The estimated mass trend over the total GIS, AIS and GOA glaciers for the time period 1 December 2003 to 1 December 2010 is -380 plus or minus 31 Gt a(exp -1), equivalent to -1.05 plus or minus 0.09 mma(exp -1) sea-level rise. Over the same time period we estimate the mass acceleration to be -41 plus or minus 27 Gt a(exp -2), equivalent to a 0.11 plus or minus 0.08 mm a(exp -2) rate of change in sea level. The trends and accelerations are dependent on significant seasonal and annual balance anomalies.

MASCON↗

On Hilbert-Huang Transform Based Synthesis of a Signal Contaminated by Radio Frequency Interference or Fringes

Norden E. Huang et al. had proposed and published the Hilbert-Huang Transform (HHT) concept correspondently in 1996, 1998. The HHT is a novel method for adaptive spectral analysis of non-linear and non-stationary signals. The HHT comprises two components: - the Huang Empirical Mode Decomposition (EMD), resulting in an adaptive data-derived basis of Intrinsic Mode functions (IMFs), and the Hilbert Spectral Analysis (HSA1) based on the Hilbert Transform for 1-dimension (1D) applied to the EMD IMF's outcome. Although paper describes the HHT concept in great depth, it does not contain all needed methodology to implement the HHT computer code. In 2004, Semion Kizhner and Karin Blank implemented the reference digital HHT real-time data processing system for 1D (HHT-DPS Version 1.4). The case for 2-Dimension (2D) (HHT2) proved to be difficult due to the computational complexity of EMD for 2D (EMD2) and absence of a suitable Hilbert Transform for 2D spectral analysis (HSA2). The real-time EMD2 and HSA2 comprise the real-time HHT2. Kizhner completed the real-time EMD2 and the HSA2 reference digital implementations respectively in 2013 & 2014. Still, the HHT2 outcome synthesis remains an active research area. This paper presents the initial concepts and preliminary results of HHT2-based synthesis and its application to processing of signals contaminated by Radio-Frequency Interference (RFI), as well as optical systems' fringe detection and mitigation at design stage. The Soil Moisture Active Passive (SMAP mission (SMAP) carries a radiometer instrument that measures Earth soil moisture at L1 frequency (1.4 GHz polarimetric - H, V, 3rd and 4th Stokes parameters). There is abundant RFI at L1 and because soil moisture is a strategic parameter, it is important to be able to recover the RFI-contaminated measurement samples (15% of telemetry). State-of-the-art only allows RFI detection and removes RFI-contaminated measurements. The HHT-based analysis and synthesis facilitates recovery of measurements contaminated by all kinds of RFI, including jamming [7-8]. The fringes are inherent in optical systems and multi-layer complex contour expensive coatings are employed to remove the unwanted fringes. HHT2-based analysis allows test image decomposition to analyze and detect fringes, and HHT2-based synthesis of useful image.

Spacecraft Electronics Equipment↗

The 27-Day Rotational Variations in Total Solar Irradiance Observations: from SORCE-TIM, ACRIMSAT-ACRIM III, and SOHO-VIRGO

During the last decade, observations from SORCE (Solar Radiation and Climate Experiment)/TIM (Total Irradiance Monitor), ACRIMSAT (Active Cavity Radiometer Irradiance Monitor Satellite)/ACRIM III, and SOHO (Solar and Heliospheric Observatory)VIRGO (Variability of IRradiance and Gravity Oscillations Sun PhotoMeter) provided the Total Solar Irradiance (TSI) measurements with unprecedented accuracy and stability to determine the amount of solar irradiance reaching the top of the atmosphere and how solar irradiance varies in different time scales. These three independent measurements are analyzed using the EEMD (Ensemble Empirical Mode Decomposition) method to characterize the phase and amplitude of 27-day solar rotational variation in TSI. The mode decomposition clearly identifies a 27-day solar rotational signature in TSI measurements. The rotational variations of TSI from the three independent observations are generally consistent with each other, despite different mean TSI values. During the declining phase of solar cycle 23, the amplitude of TSI 27-day variations is as high as 0.8 watts per square meter (approximately 0.05 percent), while during the rising phase of solar cycle 24, the amplitude is up to 0.4 watts per square meter (approximately 0.04 percent). During the minimum phase (2008-2009), the amplitude of the rotational mode is only 0.1 watts per square meter. The correlation of this rotational mode between TIM and ACRIM III is approximately 0.92 and the slope of the local peak values is approximately 0.98. The correlation between TIM and VIRGO is approximately 0.96 and the slope of the local peak values isapproximately 0.98, very similar to the slope with ACRIM III.

ACRIM III↗

Solar Rotational Modulations of Spectral Irradiance and Correlations with the Variability of Total Solar Irradiance

Aims: We characterize the solar rotational modulations of spectral solar irradiance (SSI) and compare them with the corresponding changes of total solar irradiance (TSI). Solar rotational modulations of TSI and SSI at wavelengths between 120 and 1600 nm are identified over one hundred Carrington rotational cycles during 2003-2013. Methods: The SORCE (Solar Radiation and Climate Experiment) and TIMED (Thermosphere Ionosphere Mesosphere Energetics and Dynamics)/SEE (Solar EUV Experiment) measured and SATIRE-S modeled solar irradiances are analyzed using the EEMD (Ensemble Empirical Mode Decomposition) method to determine the phase and amplitude of 27-day solar rotational variation in TSI and SSI. Results: The mode decomposition clearly identifies 27-day solar rotational variations in SSI between 120 and 1600 nm, and there is a robust wavelength dependence in the phase of the rotational mode relative to that of TSI. The rotational modes of visible (VIS) and near infrared (NIR) are in phase with the mode of TSI, but the phase of the rotational mode of ultraviolet (UV) exhibits differences from that of TSI. While it is questionable that the VIS to NIR portion of the solar spectrum has yet been observed with sufficient accuracy and precision to determine the 11-year solar cycle variations, the temporal variations over one hundred cycles of 27-day solar rotation, independent of the two solar cycles in which they are embedded, show distinct solar rotational modulations at each wavelength.

Spectral Solar Irradiance↗

On Holo-Hilbert Spectral Analysis: A Full Informational Spectral Representation for Nonlinear and Non-Stationary Data

The Holo-Hilbert spectral analysis (HHSA) method is introduced to cure the deficiencies of traditional spectral analysis and to give a full informational representation of nonlinear and non-stationary data. It uses a nested empirical mode decomposition and Hilbert-Huang transform (HHT) approach to identify intrinsic amplitude and frequency modulations often present in nonlinear systems. Comparisons are first made with traditional spectrum analysis, which usually achieved its results through convolutional integral transforms based on additive expansions of an a priori determined basis, mostly under linear and stationary assumptions. Thus, for non-stationary processes, the best one could do historically was to use the time- frequency representations, in which the amplitude (or energy density) variation is still represented in terms of time. For nonlinear processes, the data can have both amplitude and frequency modulations (intra-mode and inter-mode) generated by two different mechanisms: linear additive or nonlinear multiplicative processes. As all existing spectral analysis methods are based on additive expansions, either a priori or adaptive, none of them could possibly represent the multiplicative processes. While the earlier adaptive HHT spectral analysis approach could accommodate the intra-wave nonlinearity quite remarkably, it remained that any inter-wave nonlinear multiplicative mechanisms that include cross-scale coupling and phase-lock modulations were left untreated. To resolve the multiplicative processes issue, additional dimensions in the spectrum result are needed to account for the variations in both the amplitude and frequency modulations simultaneously. HHSA accommodates all the processes: additive and multiplicative, intra-mode and inter-mode, stationary and nonstationary, linear and nonlinear interactions. The Holo prefix in HHSA denotes a multiple dimensional representation with both additive and multiplicative capabilities.

Huang, Norden E.↗

Future Nuisance Flooding in Norfolk, VA, From Astronomical Tides and Annual to Decadal Internal Climate Variability

Abstract Increasing sea level rise will lead to more instances of nuisance flooding along the Virginia coastline in the coming decades, causing road closures and deteriorating infrastructure. These minor flood events can be caused by astronomical tides alone, in addition to internal climate variability on annual to decadal timescales. An assessment of nuisance flooding from these two effects is presented up until the year2050 for Norfolk, Virginia. The analysis of water levels indicates that nuisance flooding from tides alone in conjunction with a medium-high sea level scenario will result in flooding beginning in 2030 with frequency increasing thereafter. The addition of climate variability, by use of an empirical mode decomposition, leads to a substantial increase of flooding relative to the tides-alone analysis and shows flood events beginning as soon as 2020. High tides in the future will produce nuisance flooding without the need of other drivers such as coastal storms.

Burgos, A. G.↗

Solar Migrating Diurnal Tide in the Upper Mesosphere and Lower Thermosphere from SD-WACCM-X, Aura/MLS, and TIMED/SABER

Solar heating on the rotating Earth is a fundamental forcing of atmospheric tidal waves. The migrating diurnal tide, propagating from the troposphere to the upper mesosphere and lower thermosphere, is a large global disturbance that drives most of the daily variations in dynamics, thermal structures, chemistry, as well as atmospheric compositions. In this study, we present a comprehensive analysis of the (1,1) propagating diurnal tides derived from multi-year MLS and SABER observations, to characterize and better understand interannual and long-term tidal variability in the upper atmosphere. Although the (1,1) tide is driven by the solar heating, its interannual variations are largely determined by the internal variability of Earth’s climate system. The tidal amplitudes derived from MLS and SABER data agree well with each other in terms of monthly climatology and interannual variations, showing a consistent seasonal cycle and modulations from the Quasi Biennial Oscillation (QBO). The Ensemble Empirical Mode Decomposition (EEMD) analysis, employed to extract the low frequency variations, also reveals an ENSO-like (42 months) influence on tidal amplitudes in the upper mesosphere. While their (1,1) tidal amplitudes are less affected by the solar cycle, the mean mixing ratio of carbon monoxide (CO) and ozone (O3) is significantly modulated by the 11-year solar cycle due to their UV-dependent photochemistry.

mesosphere↗

Developmental Flight Instrumentation: Review of Space Shuttle, Ares I-X, and Artemis I

Ascent vehicles in the developmental stages of the program are instrumented with Developmental Flight Instrumentation (DFI) sensors. These sensors establish a link between a vehicle and engineers on the ground to communicate conditions experienced during the ascent. These data are then compared to pre-flight predictions used in the design process. The aerodynamic, acoustic, thermal, and structural data are either telemetered to ground stations during the ascent or stored on the vehicle for post-flight recovery and archived at the Huntsville Operations Support Center (HOSC). Following NASA's Artemis I Space Launch System (SLS) launch on November 16, 2020, data from three separate programs are available at the HOSC: Space Shuttle Program (Space Transport System (STS)), Constellation Program (Ares I-X), and Artemis Program (SLS). Availability of these data presents a unique opportunity to examine DFI data from three distinct vehicles and analyze the broad impact of the DFI data on the understanding of transonic aerodynamics. Classical spectrogram and Empirical Mode Decomposition techniques were used to present data in aerodynamically analogous regions on each vehicle. On the SLS and Ares I-X, a region downstream of the Launch Abort System motors was chosen. Comparing SLS and the STS, a region downstream of booster Froward Attach Hardware was selected as analogous flow region. Some other regions of interest were also identified. Although similarities in flow features on three vehicles were identified, some challenges in the comparison were also encountered, especially due to poor temporal and spatial resolution of Shuttle measurements.

Space Shuttle↗

Developmental Flight Instrumentation: Review of Space Shuttle, Ares I-X, and Artemis I

Ascent vehicles in the developmental stages of the program are instrumented with Developmental Flight Instrumentation (DFI) sensors. These sensors establish a link between a vehicle and engineers on the ground to communicate conditions experienced during the ascent. These data are then compared to pre-flight predictions used in the design process. The aerodynamic, acoustic, thermal, and structural data are either telemetered to ground stations during the ascent or stored on the vehicle for post-flight recovery and archived at the Huntsville Operations Support Center (HOSC). Following NASA's Artemis I Space Launch System (SLS) launch on November 16, 2020, data from three separate programs are available at the HOSC: Space Shuttle Program (Space Transport System (STS)), Constellation Program (Ares I-X), and Artemis Program (SLS). Availability of these data presents a unique opportunity to examine DFI data from three distinct vehicles and analyze the broad impact of the DFI data on the understanding of transonic aerodynamics. Classical spectrogram and Empirical Mode Decomposition techniques were used to present data in aerodynamically analogous regions on each vehicle. On the SLS and Ares I-X, a region downstream of the Launch Abort System motors was chosen. Comparing SLS and the STS, a region downstream of booster Froward Attach Hardware was selected as analogous flow region. Some other regions of interest were also identified. Although similarities in flow features on three vehicles were identified, some challenges in the comparison were also encountered, especially due to poor temporal and spatial resolution of Shuttle measurements.

Space Shuttle↗

Non-Stationary Power System Forced Oscillation Analysis using Synchrosqueezing Transform

Non-stationary forced oscillations (FOs) have been observed in power system operations. However, most detection methods assume that the frequency of FOs is stationary. In this paper, we present a methodology for the analysis of nonstationary FOs. Firstly, Fourier synchrosqueezing transform (FSST) is used to provide a concentrated time-frequency representation of the signals that allows identification and retrieval of non-stationary signal components. To continue, the Dissipating Energy Flow (DEF) method is applied to the extracted components to locate the source of forced oscillations. The methodology is tested using simulated as well as real PMU data. In conclusion, the results show that the proposed FSST-based signal decomposition provides a systematic framework for the application of DEF Method to non-stationary FOs.

42 ENGINEERING↗

An Alternate Method for Estimating Dynamic Height from XBT Profiles Using Empirical Vertical Modes

A technique is presented that applies modal decomposition to estimate dynamic height (0-450 db) from Expendable BathyThermograph (XBT) temperature profiles. Salinity-Temperature-Depth (STD) data are used to establish empirical relationships between vertically integrated temperature profiles and empirical dynamic height modes. These are then applied to XBT data to estimate dynamic height. A standard error of 0.028 dynamic meters is obtained for the waters of the Gulf of Alaska- an ocean region subject to substantial freshwater buoyancy forcing and with a T-S relationship that has considerable scatter. The residual error is a substantial improvement relative to the conventional T-S correlation technique when applied to this region. Systematic errors between estimated and true dynamic height were evaluated. The 20-year-long time series at Ocean Station P (50 deg N, 145 deg W) indicated weak variations in the error interannually, but not seasonally. There were no evident systematic alongshore variations in the error in the ocean boundary current regime near the perimeter of the Alaska gyre. The results prove satisfactory for the purpose of this work, which is to generate dynamic height from XBT data for coanalysis with satellite altimeter data, given that the altimeter height precision is likewise on the order of 2-3 cm. While the technique has not been applied to other ocean regions where the T-S relation has less scatter, it is suggested that it could provide some improvement over previously applied methods, as well.

Lagerloef, Gary S. E.↗

Principal modes of atmospheric circulation anomalies associated with global angular momentum fluctuations

This paper provides a description of the variability of global atmospheric angular momentum (GAM) and its relationship with principal modes of three-dimensional atmospheric circulation anomalies. The data used are 5-day mean global wind fields from the European Centre for Medium-Range Weather Forecasts initialized dataset for 1980-1989. Significant seasonal variation of GAM is observed with maxima in April and November and a minimum during late July. The amplitude of the annual cycle is largest in the upper troposphere and decreases toward the surface. Although the lower tropospheric contribution to the total angular momentum is relatively small, its annual cycle is out of phase with those of the upper atmosphere and GAM. Also identified is a distinct semiannual component, with double peaks appearing in April and November. This signal is most noticeable in the upper troposphere above the 300-mb level. The principal modes of zonal-mean angular momentum and meridional circulation anomalies and their coupled modes are obtained by using empirical orthogonal function analysis and singular value decomposition. It is shown that the leading modes of the angular momentum and meridional circulation are coupled with each other and are responsible for much of the variability in GAM. The coupled modes represent fluctuations of upper-level subtropical zonal flow, which are linked to the modulation of Hadley circulation intensity in both hemispheres. It is found that GAM is highly correlated with the first eigenvector of upper-level streamfunction anomalies, which consists of a superrotational flow in the tropics and subtropics, except over the central Pacific where a 'blocked' flow with two subtropical anticyclonic circulation cells straddling the equator is found. Much of the blocked flow is due to the establishment of dipole anomalies in the velocity potential with centers over the central Pacific and the Maritime Continent on the interannual time scale. On the intraseasonal time scale, GAM fluctuation is dominated by superrotational flow in the tropics, with the blocked flow present to a much lesser extent. The associated velocity potential anomaly has a weak dipole structure with centers over the Indian Ocean and the eastern Pacific. The implications of the above results on the total angular momentum balance of the earth-atmosphere system are also discussed.

Kang, In-Sik↗

Non-seasonal Gravitational Effects of Interannual Meteorological Oscillations in Atmosphere and Ocean

Via mass transport, the interannual meteorological oscillations (ENSO, NAO, etc.) produce slight variations in the global gravitational field, which can be observed by satellite laser ranging (SLR) and anticipated space missions such as GRACE. The meteorological oscillation modes are identified using empirical orthogonal function and principal component decomposition using de-seasoned surface fields produced from the 40-years of NCEP reanalysis. This decompositional analysis requires an accounting for issues associated with area-weighting and non-zero mean. The atmospheric contributions to gravitational variation, mode by mode and in combinations, are then computed via regional integration. To address the corresponding oceanic contributions, we examine the output of the numerical POCM ocean circulation model (1992-1998) in two ways: ocean bottom pressure field and ocean surface topography minus steric effect. The resulting lowest-degree Stokes coefficients for both NCEP reanalysis and POCM are then compared with space geodetic SLR observations, to identify the importance of each meteorological oscillations in gravitational variation signals.

Chao, Benjamin F.↗

Gravitational Signal of Mass Redistribution Due to Interannual Meteorological Oscillations in Atmosphere and Ocean

Interannual meteorological oscillations (ENSO, QBO, NAO, etc.) have demonstrable influences on Earth's rotation. Here we study their effects on global gravitational field, whose temporal variations are being studied using SLR (satellite laser ranging) data and in anticipation of the new space mission GRACE. The meteorological oscillation modes are identified using the EOF (empirical orthogonal function)/PC (principal component) decomposition of surface fields (in which we take care of issues associated with the area-weighting and non-zero mean). We examine two fields, one for the global surface pressure field for the atmosphere obtained from the NCEP reanalysis (for the past 40 years), one for the surface topography field for the ocean from the Topex/Poseidon (T/P) data (for the past 8 years). We use monthly maps, and remove the mean-monthly ("climatology") values from each grid point, hence focusing only on non-seasonal signals. The T/P data were first subject to a steric correction where the steric contribution to the ocean surface topography was removed according to output of the numerical POCM model. The respective atmospheric and oceanic contributions to the gravitational variation, in terms of harmonic Stokes coefficients, are then combined mode-by-mode. Since the T/P data already contain the oceanic response to overlying atmospheric pressure, no regards to the inverted-barometer behavior for the ocean need be considered. Results for the lowest-degree Stokes coefficients can then be compared with space geodetic observations including the Earth's rotation and the SLR data mentioned above, to identify the importance of each meteorological oscillations in gravitational variation signals.

Chao, B. F.↗

Interannual Sea Level Variations in the Tropical Indian Ocean from Geosat and Shallow Water Simulations

Sea level variations of the Indian Ocean north of 20 deg S are analyzed from Geosat satellite altimeter data over April 1985-September 1989. These variations are compared and interpreted with numerical simulations derived from a reduced gravity model forced by FSU observed winds over the same period. After decomposition into complex empirical orthogonal functions, the low-frequency anomalies are described by the first two modes for observations as well as for simulations. The sums of the two modes contain 34% and 40% of the observed and simulated variances, respectively. Averaged over the basin, the observed and simulated sea level changes are correlated by 0.92 over 1985-1988. The strongest change happens during the El Ninio 1986-1987: between winter 1986 and summer 1987 the basin-averaged sea level rises by approx. 1 cm. These low-frequency variations can partly be explained by changes in the Sverdrup circulation. The southern tropical Indian Ocean between 1O deg and 20 deg S is the domain where those changes are strongest: the averaged sea level rises by approx. 4.5 cm between winter 1986 and winter 1987. There, the signal propagates southwestward across the basin at a speed similar to free Rossby waves. Sensitivity of observed anomalies is examined over 1987-1988, with different orbit ephemeris, tropospheric corrections, and error reduction processes. The uncertainty of the basin-averaged sea level estimates is mostly due to the way the orbit error is reduced and reaches approx. 1 cm. Nonetheless, spatial correlation is good between the various observations and better than between observations and simulations. Sensitivity of simulated anomalies to the wind uncertainty, examined with Former Soviet Union (FSU) and European Center for Medium-Range Weather Forecasting (ECMWF) forcings over 1985-1988, shows that the variance of the simulations driven by ECMWF is 52% smaller, as FSU winds are stronger than ECMWF. Results show that the wind strength also affects the dynamic response of the ocean: anomalies propagate westward across the basin more than twice as fast with FSU than with ECMWF. It is found that the discrepancy is larger between ECMWF and FSU simulations than between observations and FSU simulations.

Perigaud, Claire↗

Simulation of high-frequency dissolved oxygen dynamics in a shallow estuary, the Corsica River, Chesapeake Bay

Understanding shallow water biogeochemical dynamics is a challenge in coastal regions, due to the presence of highly variable land-water interface fluxes, tight coupling with sediment processes, tidal dynamics, and diurnal variability in biogeochemical processes. While the deployment of continuous monitoring devices has improved our understanding of high-frequency (12 - 24 hours) variability and spatial heterogeneity in shallow regions, mechanistic modeling of these dynamics has lagged behind conceptual and empirical models. The inherent complexity of shallow water systems is represented in the Corsica River estuary, a small basin within the Chesapeake Bay ecosystem, where abundant monitoring data have been collected from long-term monitoring stations, continuous monitoring sensors, synoptic sensor surveys, and measurements of sediment-water fluxes. A state-of-the-art modeling system, the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), was applied to the Corsica domain with a high-resolution grid and nutrient loads from the most recent version of the Chesapeake Bay watershed model. The Corsica SCHISM model reproduced observed high-frequency variability in dissolved oxygen, as well as seasonal variability in chlorophyll-a and sediment-water fluxes. Time-series signal analyses using Empirical Model Decomposition and spectral analysis revealed that the diurnal and M2 tide frequencies are the dominant high-frequency modes and physical transport contributes a larger share to dissolved oxygen budgets than biogeochemical processes on an hourly time scale. Heterogeneity and patchiness in dissolved oxygen resulting from phytoplankton distributions and geometry-driven eddies amplify the physical transport effect, and on longer time scales oxygen is controlled more by photosynthesis and respiration. Our simulation demonstrates that interactions among physical and biological dynamics generate complex high-frequency variability in water quality and non-linear reposes to nutrient loading and environmental forcing in shallow water systems.

54 ENVIRONMENTAL SCIENCES↗