Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Total variation 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.

A Lagrange multiplier based divide and conquer finite element algorithm

A novel domain decomposition method based on a hybrid variational principle is presented. Prior to any computation, a given finite element mesh is torn into a set of totally disconnected submeshes. First, an incomplete solution is computed in each subdomain. Next, the compatibility of the displacement field at the interface nodes is enforced via discrete, polynomial and/or piecewise polynomial Lagrange multipliers. In the static case, each floating subdomain induces a local singularity that is resolved very efficiently. The interface problem associated with this domain decomposition method is, in general, indefinite and of variable size. A dedicated conjugate projected gradient algorithm is developed for solving the latter problem when it is not feasible to explicitly assemble the interface operator. When implemented on local memory multiprocessors, the proposed methodology requires less interprocessor communication than the classical method of substructuring. It is also suitable for parallel/vector computers with shared memory and compares favorably with factorization based parallel direct methods.

Farhat, C.↗

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↗

The Budget and Partitioning of Stratospheric Chlorine During the 1997 Arctic Summer

Volume mixing ratio profiles of HCl, HOCl, ClNO3, CH3Cl, CFC-12, CFC-11, CCl4, HCFC-22, and CFC-113 were measured simultaneously from 9 to 38 km by the Jet Propulsion Laboratory MkIV Fourier Transform Infrared solar absorption spectrometer during two balloon flights from Fairbanks, Alaska (64.8 N), on May 8 and July 8, 1997. The altitude variation of total organic chlorine (CCly) total inorganic chlorine (Cly), and the nearly constant value (3.7 +/- 0.2 ppbv) of their sum (Cl(sub TOT)) demonstrates that the stratospheric chlorine species available to react with O3 are supplied by the decomposition of organic chlorinated compounds whose abundances are well quantified. Measured profiles of HCl and ClNO3 agree well with profiles found by photochemical model (differences < 10% for altitudes below 35 km) constrained by various other constituents measured by MkIV. The production of HCl by ClO + OH plays a relatively small role in the partitioning of HCl and ClNO3 for the sampled air masses. However, better agreement with the measured profiles of HCl and ClNO3 is obtained when this source of HCl is included in the model. Both the measured and calculated [ClNO3]/[HCl] ratios exhibit the expected near linear variation with [O3]2/[CH4] over a broad range of altitudes. MkIV measurements of HCl, ClNO3, and CC1y agree well with ER-2 in situ observations of these quantities for directly comparable air masses. These results demonstrate good understanding of the budget of stratospheric chlorine and that the partitioning of inorganic chlorine is accurately described by photochemical models that employ JPL97 reaction rates and production of HCl from ClO + OH for the environmental conditions encountered: relatively warm temperatures, long periods of solar illumination, and relatively low aerosol surface areas.

Sen, B.↗

The Budget and Partitioning of Stratospheric Chlorine During the 1997 Arctic Summer

Volume mixing ratio profiles of HCl, HOCl, ClNO3, CH3Cl, CFC-12, and CFC-11, CCl4, HCFC-22, and CFC-113 were measured simultaneously from 9 to 38 km by the Jet Propulsion Laboratory MkIV Fourier Transform Infrared solar absorption spectrometer during two balloon flights from Fairbanks, Alaska (64.8 N), on May 8 and July 8, 1997. The altitude variation of total organic chlorine (CCl(sub y)), total inorganic chlorine (Cl(sub y)), and the nearly constant value (3.7 +/- 0.2 ppbv) of their sum (Cl(sub TOT)) demonstrates that the stratospheric chlorine species available to react with O3 are supplied by the decomposition of organic chlorinated compounds whose abundances are well quantified. Measured profiles of HCl and ClNO3 agree well with profiles found by photochemical model (differences < 10% for altitudes below 35 km) constrained by various other constituents measured by MkIV. The production of HCl by ClO + OH plays a relatively small role in the partitioning of HCl and ClNO3 for the sampled air masses. However, better agreement with the measured profiles of HCl and ClNO3 is obtained when this source of HCl is included in the model. Both the measured and calculated [ClNO3]/[HCl] ratios exhibit the expected near linear variation with [O3](sup 2)/[CH4] over a broad range of altitudes. MkIV measurements of HCl, ClNO3, and CCl(sub y) agree well with ER-2 in situ observations of these quantities for directly comparable air masses. These results demonstrate good understanding of the budget of stratospheric chlorine and that the partitioning of inorganic chlorine is accurately described by photochemical models that employ JPL97 reaction rates and production of HCl from ClO + OH for the environmental conditions encountered: relatively warm temperatures, long periods of solar illumination, and relatively low aerosol surface areas.

Sen, B.↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

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 the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗

Ecosystem Impacts of Woody Encroachment In Texas: A Spatial Analysis Using AVIRIS

Woody encroachment, the increase of woody plant density relative to herbaceous vegetation, has been documented in drylands of Texas as well as worldwide (Archer 1994, Harrington and Harman 1995, Moleele et al. 2002). Over-grazing, fire suppression and climate change are implicated in the shift from open grasslands to ecosystems now populated by trees and shrubs (Scholes and Archer 1997, Archer et al. 2001), such as Prosopis glandulosa var. glandulosa (honey mesquite) in north Texas (Teague et al. 1997, Ansley et al. 2001, Asner et al. 2003a). Several studies have examined changes in ecosystem properties accompanying woody vegetation encroachment in the Southwest U.S., with research focused on increases in plant and soil carbon (C) and nitrogen (N) stores (Hoffman and Jackson 2000, Asner et al. 2003a), isotopic shifts in these pools (Boutton 1999, Archer et al. 2001), and increases in N cycling rates (Rundel et al. 1982, Hibbard et al. 2001). However, little is known regarding the impact of woody encroachment on N trace gas emissions from dryland regions such as Texas. NOx is produced in the soil during the processes of nitrification and denitrification (Firestone and Davidson 1989). The total N efflux from soils is most directly influenced by the internal cycling of N, which at a regionalscale, is controlled by the inputs and availability of N from vegetation via litterfall and subsequent decomposition (Robertson et al. 1989). Although plot-scale studies are critical to understanding controls over N oxide emissions, regionalization of the measurements is impeded by spatial variation in the factors contributing most to N cycling processes: soil properties (affecting soil moisture regimes and N stocks) and vegetation cover (affecting litter inputs and N uptake). While broad patterns in ecosystem structure and vegetation composition co-vary with general patterns of trace gas emissions (Matson 1997), there is no easily measured index of N availability that can be applied for regional-scale studies of N oxide fluxes. Remote sensing is arguably the only approach available to develop a spatially-explicit understanding of ecosystem processes. More specifically, remotely detectable spatial patterns in the distal controls over soil N properties, such as vegetation cover, land use and soil type (Robertson et al. 1989), should be exploited for regional studies of N oxide emissions. The woody encroachment phenomenon provides an opportunity to test the strength of the relationship between N oxide emissions and those factors controlling the fluxes that can be remotely measured. If such linkages can be firmly established, and if the spatial pattern of distal controls is relevant, then the combination of field measurements and remote sensing offers to improve regional-scale N oxide estimates. The paper presents the utility of linking field based sampling of soil NOx emissions with very high resolution remote sensing estimates of woody vegetation cover from the NASA AVIRIS, Airborne Visible-Infrared Imaging Spectrometer (Green et al. 1998, Asner and Green 2001) and automated spectral mixture analysis (Asner and Lobell 2000, Asner and Heidebrecht 2002) that provide a means to spatially extrapolate soil NOx emissions to the regional scale.

Martin, Roberta E.↗

Using Enhanced Grace Water Storage Data to Improve Drought Detection by the U.S. and North American Drought Monitors

NASA's Gravity Recovery and Climate Experiment (GRACE) satellites measure time variations of the Earth's gravity field enabling reliable detection of spatio-temporal variations in total terrestrial water storage (TWS), including groundwater. The U.S. and North American Drought Monitors rely heavily on precipitation indices and do not currently incorporate systematic observations of deep soil moisture and groundwater storage conditions. Thus GRACE has great potential to improve the Drought Monitors by filling this observational gap. GRACE TWS data were assimilating into the Catchment Land Surface Model using an ensemble Kalman smoother enabling spatial and temporal downscaling and vertical decomposition into soil moisture and groundwater components. The Drought Monitors combine several short- and long-term drought indicators expressed in percentiles as a reference to their historical frequency of occurrence. To be consistent, we generated a climatology of estimated soil moisture and ground water based on a 60-year Catchment model simulation, which was used to convert seven years of GRACE assimilated fields into drought indicator percentiles. At this stage we provide a preliminary evaluation of the GRACE assimilated moisture and indicator fields.

Houborg, Rasmus↗

Satellite Sampling and Retrieval Errors in Regional Monthly Rain Estimates from TMI AMSR-E, SSM/I, AMSU-B and the TRMM PR

Passive and active microwave rain sensors onboard earth-orbiting satellites estimate monthly rainfall from the instantaneous rain statistics collected during satellite overpasses. It is well known that climate-scale rain estimates from meteorological satellites incur sampling errors resulting from the process of discrete temporal sampling and statistical averaging. Sampling and retrieval errors ultimately become entangled in the estimation of the mean monthly rain rate. The sampling component of the error budget effectively introduces statistical noise into climate-scale rain estimates that obscure the error component associated with the instantaneous rain retrieval. Estimating the accuracy of the retrievals on monthly scales therefore necessitates a decomposition of the total error budget into sampling and retrieval error quantities. This paper presents results from a statistical evaluation of the sampling and retrieval errors for five different space-borne rain sensors on board nine orbiting satellites. Using an error decomposition methodology developed by one of the authors, sampling and retrieval errors were estimated at 0.25 resolution within 150 km of ground-based weather radars located at Kwajalein, Marshall Islands and Melbourne, Florida. Error and bias statistics were calculated according to the land, ocean and coast classifications of the surface terrain mask developed for the Goddard Profiling (GPROF) rain algorithm. Variations in the comparative error statistics are attributed to various factors related to differences in the swath geometry of each rain sensor, the orbital and instrument characteristics of the satellite and the regional climatology. The most significant result from this study found that each of the satellites incurred negative longterm oceanic retrieval biases of 10 to 30%.

Fisher, Brad↗

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↗

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗

Satellite Monitoring of Global Surface Soil Organic Carbon Dynamics Using the SMAP Level 4 Carbon Product

Soil organic carbon (SOC) is an important metric of soil health and the terrestrial carbon balance. Short‐term climate variations affect SOC through changes in temperature and moisture, which control vegetation growth and soil decomposition. We evaluated a satellite data‐driven carbon model, operating under the NASA Soil Moisture Active‐Passive (SMAP) mission, as a means of monitoring global surface SOC dynamics. The SMAP Level 4 Carbon (L4C) product estimates a daily global carbon budget including surface (0‐ to 5‐cm depth) SOC. We found that the L4C mean latitudinal SOC distribution is generally consistent with alternative assessments from static soil inventory records and dynamic global vegetation models (r ≥ 0.89). Within forest systems, based on inventory data, L4C SOC is most similar in magnitude to litterfall but is correlated with coarse woody debris ( urn:x-wiley:jgrg:media:jgrg21790:jgrg21790-math-0001) and total SOC ( urn:x-wiley:jgrg:media:jgrg21790:jgrg21790-math-0002). L4C SOC is sensitive to seasonal and annual climate variability, with mean residence times that range from 1.5 years in the wet tropics to 17 years in the cold tundra. Incorporating soil moisture retrievals from the SMAP L‐band (1.4 GHz) microwave radiometer within the L4C algorithm provides enhanced soil moisture sensitivity under low‐to‐moderate vegetation cover (<5 kg/sq.m vegetation water content). The L‐band soil moisture had the greatest impact on the L4C carbon budget in semiarid regions, which span almost 60% of the globe and account for substantial variability in the terrestrial carbon sink. The L4C operational product enables prognostic investigations into effects of recent climate trends and anomalies (e.g., droughts and pluvials) on shallow soil carbon dynamics.

K. Arthur Endsley↗

Quantifying Errors in TRMM-Based Multi-Sensor QPE Products Over Land in Preparation for GPM

Determining uncertainties in satellite-based multi-sensor quantitative precipitation estimates over land of fundamental importance to both data producers and hydro climatological applications. ,Evaluating TRMM-era products also lays the groundwork and sets the direction for algorithm and applications development for future missions including GPM. QPE uncertainties result mostly from the interplay of systematic errors and random errors. In this work, we will synthesize our recent results quantifying the error characteristics of satellite-based precipitation estimates. Both systematic errors and total uncertainties have been analyzed for six different TRMM-era precipitation products (3B42, 3B42RT, CMORPH, PERSIANN, NRL and GSMap). For systematic errors, we devised an error decomposition scheme to separate errors in precipitation estimates into three independent components, hit biases, missed precipitation and false precipitation. This decomposition scheme reveals hydroclimatologically-relevant error features and provides a better link to the error sources than conventional analysis, because in the latter these error components tend to cancel one another when aggregated or averaged in space or time. For the random errors, we calculated the measurement spread from the ensemble of these six quasi-independent products, and thus produced a global map of measurement uncertainties. The map yields a global view of the error characteristics and their regional and seasonal variations, reveals many undocumented error features over areas with no validation data available, and provides better guidance to global assimilation of satellite-based precipitation data. Insights gained from these results and how they could help with GPM will be highlighted.

Peters-Lidard, Christa D.↗

Interannual Atmospheric Oscillations and Their Gravitational Effects

Past studies have examined the effects of the interannual atmospheric oscillations (IAO), or often called teleconnection patterns, on the Earth's rotation, such as El Nino/Southern Oscillation (ENSO), the quasi-biennial oscillation, and to a lesser extent the North Atlantic Oscillation (NAO). The present study focuses on IAO's effects on the gravitational field and geocenter; the purpose is to be able to identify and isolate the contributions of each prominent IAO in relation to the total contribution of the atmosphere, in terms of their magnitudes, their geographical patterns, and their interannual time history. We use the 40-year NCEP reanalysis of the monthly, global atmospheric surface pressure field as our basic data set. The method we apply to isolate the IAOs is the empirical orthogonal function (EOF) decomposition which is widely used in meteorological investigations. We do the EOF analysis first on IAOs' seasonal signals (by "collapsing" the 40-year series into 12 mean-months for each grid point) and obtain estimates for their respective contributions. Then we remove these seasonal signals from the data to focus on the (broad-band) interannual EOFs. We examine ENSO, NAO, North Pacific Oscillation, and other less prominent IAOs that can be identified from our data set, and compute their respective contributions to the variation of global gravitational field and geocenter motion. Finally, we compare the results with the available observational data, and discuss the implications w.r.t. the upcoming space missions such as GRACE.

Chao, Benjamin F.↗

Flux to the atmosphere of CH4 and CO2 from wetland ponds on the Hudson Bay lowlands (HBLs)

Ponds on peatlands of the Hudson Bay lowlands (HBLs) are complex ecosystems in which the fluxes to the atmosphere of CH4 and CO2 were controlled by interacting physical and biological factors. This resulted in strong diel variations of both dissolved gas concentrations and gas fluxes to the atmosphere, necessitating frequent sampling on a 24-hour schedule to enable accurate estimates of daily fluxes. Ponds at three sites on the HBL were constant net sources of CH4 and CO2 to the atmosphere at mean rates of 110-180 mg CH4 m(exp -2)/d and 3700-11,000 mg CO2 m(exp -2)/d. Rates peaked in August and September. For CH4 the pond fluxes were 3-30 times higher than adjacent vegetated surfaces. For CO2 the net pond fluxes were similar in magnitude to the vegetated fluxes but the direction of the flux was opposite, toward atmosphere. Even though ponds cover only 8-12% of the HBL area, they accounted for 30% of its total CH4 flux to the atmosphere. There is some circumstantial evidence that the ponds are being formed by decomposition of the underlying peat and that this decomposition is being stimulated by the activity of N2 fixing cyanobacteria that grow in mats at the peat-water interface. The fact that the gas fluxes from the ponds were so different from the surrounding vegetated surfaces means that any change in the ratio of pond to vegetated area, as may occur in response to climate change, would affect the total HBL fluxes.

Hamilton, J. David↗

A Data-Driven Reduced Order Model of an Isolated Rotor

There are numerous conceptual design stage rotorcraft analysis tasks which demand a high-fidelity and low cost method for rotor load distribution predictions. Considering Urban Air Mobility (UAM) vehicles aim to operate in close proximity to buildings and with unique rotor configurations, there is a significant challenge in quickly and accurately modeling rotors operating in complex, turbulent flow fields. One potential path for deriving a high-fidelity, low cost rotor model is with data-driven surrogate modeling. In this study, an initial investigation is taken to apply a proper orthogonal decomposition (POD) based reduced order model (ROM) for the purpose of pressure distribution prediction. In this study, a POD ROM was derived to produce distributed pressure predictions on rotor blades subjected to topology change due to variation in twist and taper ratio. Rotor twist was varied between 0◦, 10◦, 20◦, and 30◦ while taper ratio was varied between 1.0, 0.9, 0.8, and 0.7. All rotors consisted of a single blade. The POD ROM was validated for three demonstration cases; a high thrust rotor in hover, a low thrust rotor in hover, and a rotor in forward flight with a flight speed of M = 0.1. Results showed highly accurate distributed load predictions could be achieved at minimal computational cost. Computational cost for hovering blade surface pressure modeling was reduced from 12 hours on 440 cores to 10−5 seconds on a single core. For blade in forward flight cost was reduced from 20 hours on 440 cores to 0.6 seconds on a single core. For cases of high thrust and low thrust rotors, POD ROM was used to undergo a design optimization of the rotor such that figure of merit was maximized. Total optimization time for each case was 1 minute.

Data-Driven↗