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.

At least 19 records

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗

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.↗

Point-scale organic-matter decomposition in streambeds is weakly associated with reach-scale respiration

Stream and river ecosystems play a central role in the movement and decomposition of particulate organic matter, serving as a conduit between terrestrial hillslopes and coastal environments. Microbial-catalyzed decomposition generates simpler organic molecules that fuel respiration, often in the sediments of these ecosystems. However, the degree of connection between sediment-associated respiration (ER??d) and organic-matter decomposition remains poorly understood. How that relationship compares to decomposition’s relationship with whole ecosystem (ERtot) and water column (ERwc) respiration is also not clear. We examined the link between particulate organic matter decomposition—using cellulose-based cotton strips as a standardized substrate—and all three components of respiration across 48 sites in the environmentally diverse Yakima River Basin (Washington State, USA). We hypothesized that decomposition within sediments would be most strongly related to ERsed, but decomposition rates were more closely associated with ERtot, with little connection to ERsed or ERwc. This suggests that particulate organic matter decomposition within stream/river sediments reflects integrated system respiration rather than processes confined to sediments or the water column alone. Further, across the basin, decomposition rates nearly spanned the previously reported global range for streams and rivers and were best explained by total dissolved nitrogen (TDN), sediment grain size, and aridity of the upstream drainage area. These results highlight the strong influence of land cover and basin-scale biophysical variation on sediment-associated decomposition processes and indicate that mechanistic models of organic matter decomposition in streams/rivers should account for coupled sediment–water–land interactions.

Stegen, James C. (ORCID:0000000191357424)↗

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.↗

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Ecohydrological controls on root and microbial respiration in the East River watershed of Colorado

The main objective of this project was to conduct exploratory work to quantify how snow and rain water inputs influence the CO2 coming from the soil surface (soil CO2 flux), and its plant and microbial sources, in the East River watershed, near Crested Butte, Colorado. Knowledge gained from this effort laid the groundwork for a more comprehensive (ongoing) follow-on grant that is using a combination of experiments, field observations, machine learning and simulation modeling to fully disentangle these relationships. New field measurements were made at four locations along an elevational transect on Snodgrass Mountain at Rocky Mountain Biological Laboratory. These sites were chosen to differ in total snowpack and water table depth, and included both deciduous (aspen) and evergreen (spruce/fir) forest types. We used automated measurements of soil CO2 concentrations to quantify the total soil CO2 flux, and the vertical CO2 production within the soil profile at each site. Isotope (radiocarbon, 14C) measurements determined how much of the CO2 emitted from the soil surface came from plant respiration (root metabolism) versus microbial respiration (decomposition of soil organic matter) sources. Supporting data on plant phenology and microbial dynamics provided context for the observed variation in respiration sources. This work was motivated by our overarching hypothesis that quantifying belowground plant and microbial processes separately, and how they are influenced by snow and rain inputs, is necessary for understanding and predicting how the East River watershed ecosystems will respond to future environmental change.

54 ENVIRONMENTAL SCIENCES↗

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↗

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗

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↗