Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate gradient”

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 199 records · Page 11

Satellite Observations of Desert Dust-induced Himalayan Snow Darkening

The optically thick aerosol layer along the southern edge of the Himalaya has been subject of several recent investigations relating to its radiative impacts on the South Asian summer monsoon and regional climate forcing. Prior to the onset of summer monsoon, mineral dust from southwest Asian deserts is transported over the Himalayan foothills on an annual basis. Episodic dust plumes are also advected over the Himalaya, visible as dust-laden snow surface in satellite imagery, particularly in western Himalaya. We examined spectral surface reflectance retrieved from spaceborne MODIS observations that show characteristic reduction in the visible wavelengths (0.47 nm) over western Himalaya, associated with dust-induced solar absorption. Case studies as well as seasonal variations of reflectance indicate a significant gradient across the visible (0.47 nm) to near-infrared (0.86 nm) spectrum (VIS-NIR), during premonsoon period. Enhanced absorption at shorter visible wavelengths and the resulting VIS-NIR gradient is consistent with model calculations of snow reflectance with dust impurity. While the role of black carbon in snow cannot be ruled out, our satellite-based analysis suggests the observed spectral reflectance gradient dominated by dust-induced solar absorption during premonsoon season. From an observational viewpoint, this study underscores the importance of mineral dust deposition toward darkening of the western Himalayan snow cover, with potential implications to accelerated seasonal snowmelt and regional snow albedo feedbacks.

Remote Sensing↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Evaluation of water vapor distribution in general circulation models using satellite observations

This paper presents a comparison of the water vapor distribution obtained from two general circulation models, the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM), with satellite observations of total precipitable water (TPW) from Special Sensor Microwave/Imager (SSM/I) and upper tropospheric relative humidity (UTH) from GOES. Overall, both models are successful in capturing the primary features of the observed water vapor distribution and its seasonal variation. For the ECMWF model, however, a systematic moist bias in TPW is noted over well-known stratocumulus regions in the eastern subtropical oceans. Comparison with radiosonde profiles suggests that this problem is attributable to difficulties in modeling the shallowness of the boundary layer and large vertical water vapor gradients which characterize these regions. In comparison, the CCM is more successful in capturing the low values of TPW in the stratocumulus regions, although it tends to exhibit a dry bias over the eastern half of the subtropical oceans and a corresponding moist bias in the western half. The CCM also significantly overestimates the daily variability of the moisture fields in convective regions, suggesting a problem in simulating the temporal nature of moisture transport by deep convection. Comparison of the monthly mean UTH distribution indicates generally larger discrepancies than were noted for TPW owing to the greater influence of large-scale dynamical processes in determining the distribution of UTH. In particular, the ECMWF model exhibits a distinct dry bias along the Intertropical Convergence Zone (ITCZ) and a moist bias over the subtropical descending branches of the Hadley cell, suggesting an underprediction in the strength of the Hadley circulation. The CCM, on the other hand, demonstrates greater discrepancies in UTH than are observed for the ECMWF model, but none that are as clearly correlated with well-known features of the large-scale circulation.

Soden, Brian J.↗

Mangrove Canopy Height Globally Related to Precipitation and Cyclone Frequency

Mangrove wetlands are among the most productive and carbon-dense ecosystems in the world. Their structural attributes vary considerably across spatial scales, yielding large uncertainties in regional and global estimates of carbon stocks. Here, we present a global analysis of mangrove canopy height gradients and aboveground carbon stocks based on remotely sensed measurements and field data. Our study highlights that precipitation, temperature and cyclone frequency explain 74% of the global trends in maximum canopy height, with other geophysical factors influencing the observed variability at local and regional scales. We find the tallest mangrove forests in Gabon, equatorial Africa, where stands attain 62.8 m. The total global mangrove carbon stock (above- and belowground biomass, and soil) is estimated at 5.03 Pg, with a quarter of this value stored in Indonesia. Our analysis implies sensitivity of mangrove structure to climate change, and offers a baseline to monitor national and regional trends in mangrove carbon stocks.

Carbon cycle↗

A numerical simulation of seasonal stratospheric climate. I - Zonal temperatures and winds

Two numerical experiments were designed to isolate the effects of seasonal and latitudinal variations in longwave radiative transfer processes. The first experiment aimed at computing the rate of temperature change in the stratosphere due to longwave radiative transfer by CO2, H2O and O3 from Ramanathan's (1976) radiative-convective model applied to the troposphere-stratosphere circulation. The second experiment employs the Newtonian cooling approximation in which the rate of temperature change in the stratosphere by longwave radiative transfer is set equal to the product of the Newtonian cooling coefficient 'h' and the departure of the local temperature from a reference temperature. It is shown that the latitudinal temperature distribution of the lower stratosphere is maintained by the combined effects of dynamics, O3 solar heating and the longwave radiative coupling between troposphere and stratosphere. The latitudinal gradient in the troposphere-stratosphere longwave radiative coupling is maximum during winter and spring. The radiative response time (1/h) of the middle and upper stratosphere undergoes significant latitudinal and seasonal variations, largely due to the temperature dependence of h.

Ramanathan, V.↗

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) Contribution to CMIP6: Investigation of Sea-Level and Ocean Climate Change in Response to CO2 Forcing

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sealevel rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic meridional overturning circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea-surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, while the weakening of the AMOC causes redistribution of heat towards lower latitudes. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models.

thermosteric sea-level change↗

Venus: The case for a wet origin and a runaway greenhouse

To one interested in atmospheric evolution, the most intriguing aspect of our neighboring planet Venus is its lack of water. Measurements made by Pioneer Venus and by Several Venera spacecraft indicate that the present water abundance in Venus' lower atmosphere is of the order of 20 to 200 ppmv, or 3 x 10( exp -6) to 3 x 10 (exp -5) of the amount of water in Earth's oceans. The exact depletion factor is uncertain, in part because of an unexplained vertical gradient in H2O concentration in the lowest 10 km of the venusian atmosphere, but the general scarcity of water is well established. The interesting question, then, is: Was venus deficient in water when it formed and, if not, where did its water go? The conclusion that Venus was originally wet is consistent with its large endowment of other volatiles and with the enhanced D/H ratio in the present atmosphere. The most likely mechanism by which Venus could have lost its water is by the development of a runaway or moist greenhouse atmosphere followed by photodissociation of water vapor and escape of hydrogen to space. Climate model calculations that neglect cloud albedo feedback predict the existence of two critical transitions in atmospheric behavior at high solar fluxes: (1) at a solar flux of approximately 1.1 times the value at Earth's orbit, S(o), the abundance of stratospheric water vapor increases dramatically, permitting rapid escape of hydrogen to space (termed a moist greenhouse) and (2) at a solar flux of approximately 1.4 S(o), the oceans vaporize entirely, creating a true runaway greenhouse. If cloudiness increases at high surface temperatures, as seems likely, and if the dominant effect of clouds is to cool the planet by reflecting incident solar radiation, the actual solar flux required to create moist or runaway conditions would be higher than the values quoted above. Early in solar system history, solar luminosity was about 25 percent to 30 percent less than today, putting the flux at Venus' orbit in the range of 1.34 S(o) to 1.43 S(o). Thus, it is possible that Venus had liquid water on its surface for several hundred million years following its formation. Paradoxically, this might have facilitated water loss by sequestering atmospheric CO2 in carbonate rocks and by providing an effective medium for surface oxidation.

Kasting, J. F.↗

Reevaluating Suitability Estimates Based on Dynamics of Cropland Expansion in the Brazilian Amazon

Agricultural suitability maps are a key input for land use zoning and projections of cropland expansion. Suitability assessments typically consider edaphic conditions, climate, crop characteristics, and sometimes incorporate accessibility to transportation and market infrastructure. However, correct weighting among these disparate factors is challenging, given rapid development of new crop varieties, irrigation, and road networks, as well as changing global demand for agricultural commodities. Here, we compared three independent assessments of cropland suitability to spatial and temporal dynamics of agricultural expansion in the Brazilian state of Mato Grosso during 2001 2012. We found that areas of recent cropland expansion identified using satellite data were generally designated as low to moderate suitability for rainfed crop production. Our analysis highlighted the abrupt nature of suitability boundaries, rather than smooth gradients of agricultural potential, with little additional cropland expansion beyond the extent of the flattest areas (0-2% slope). Satellite-based estimates of the interannual variability in the use of existing crop areas also provided an alternate means to assess suitability. On average, cropland areas in the Cerrado biome had higher utilization (84%) than croplands in the Amazon region of northern Mato Grosso (74%). Areas of more recent expansion had lower utilization than croplands established before 2002, providing empirical evidence for lower suitability or alternative management strategies (e.g., pasture soya rotations) for lands undergoing more recent land use transitions. This unplanted reserve constitutes a large area of potentially available cropland (PAC)without further expansion, within the management limits imposed for pest management and fallow cycles. Using two key constraints on future cropland expansion, slope and restrictions on further deforestation of Amazon or Cerrado vegetation, we found little available flat land for further legal expansion of crop production in Mato Grosso. Dynamics of cropland expansion from more than a decade of satellite observations indicated narrow ranges of suitability criteria, restricting PAC under current policy conditions, and emphasizing the advantages of field-scale information to assess suitability and utilization.

Morton, Douglas C.↗

On the Importance of Neutral Composition and Temperature Measurements in the 100–200 km Altitude Region

Currently, thermospheric species densities and temperatures between ∼100 and 200 km are not known to the accuracy needed to fully characterize how the thermosphere transitions from a well-mixed atmosphere to a diffusively separated atmosphere with zero temperature gradient. This greatly inhibits scientific discovery attainable from either models or observations in this region, especially the understanding of mechanisms that drive thermosphere and ionospheric variability from space weather to climatological time scales. The purpose of this paper is to highlight the importance and critical need for new, global, height-resolved neutral composition (O, O2, N2) and temperature measurements in the new ignorosphere: the 100–200 km region of the thermosphere. We conclude with observation recommendations and requirements for new comprehensive composition and temperature measurements in the 100–200 km altitude region that would lead to significant advances in thermosphere-ionosphere science, space weather, and space climate.

thermosphere↗

An Earth longwave radiation climate model

An Earth outgoing longwave radiation (OLWR) climate model was constructed for radiation budget study. Required information is provided by on empirical 100mb water vapor mixing ratio equation of the mixing ratio interpolation scheme. Cloud top temperature is adjusted so that the calculation would agree with NOAA scanning radiometer measurements. Both clear sky and cloudy sky cases are calculated and discussed for global average, zonal average and world-wide distributed cases. The results agree well with the satellite observations. The clear sky case shows that the OLWR field is highly modulated by water vapor, especially in the tropics. The strongest longitudinal variation occurs in the tropics. This variation can be mostly explained by the strong water vapor gradient. Although in the zonal average case the tropics have a minimum in OLWR, the minimum is essentially contributed by a few very low flux regions, such as the Amazon, Indonesian and the Congo.

Yang, S. K.↗

Climate change and the middle atmosphere. II - The impact of volcanic aerosols

The response of the middle atmosphere to an increase in stratospheric aerosols, normally associated with increased volcanic activity, is investigated. The aerosols are found to induce a direct stratospheric response, with warming in the tropical lower stratosphere, and cooling at higher latitudes. On the shorter time scales, this radiative effect increases tropospheric static stability at low- to midlatitudes, which reduces the intensity of the Hadley cell and Ferrel cell. There is an associated increase in tropospheric standing wave energy and a decrease in midlatitude west winds, which result in additional wave energy propagation into the stratosphere at lower midlatitudes in both hemispheres. On the longer time scale, a strong hemispheric asymmetry arises. In the Northern Hemisphere eddy energy decreases, as does the middle-atmosphere residual circulation, and widespread stratospheric cooling results. In the Southern Hemisphere, the large increase in sea ice increases the tropospheric latitudinal temperature gradient, leading to increased eddy energy, an increased middle-atmosphere residual circulation, and some high-latitude stratospheric warming.

Rind, D.↗

Photosynthetic Control of Atmospheric Carbonyl Sulfide during the Growing Season

Climate models incorporate photosynthesis-climate feedbacks, yet we lack robust tools for large-scale assessments of these processes. Recent work suggests that carbonyl sulfide (COS), a trace gas consumed by plants, could provide a valuable constraint on photosynthesis. Here we analyze airborne observations of COS and carbon dioxide concentrations during the growing season over North America with a three-dimensional atmospheric transport model. We successfully modeled the persistent vertical drawdown of atmospheric COS using the quantitative relation between COS and photosynthesis that has been measured in plant chamber experiments. Furthermore, this drawdown is driven by plant uptake rather than other continental and oceanic fluxes in the model. These results provide quantitative evidence that COS gradients in the continental growing season may have broad use as a measurement-based photosynthesis tracer.

Campbell, J. Elliott↗

Reconstruction of Snow on Arctic Sea Ice

Snow on Arctic sea ice is a poorly observed variable that plays an important role in the Arctic climate system and impacts the remote sensing systems that monitor Arctic sea ice. We present and validate a reconstruction of Arctic snow depth based on observed sea ice motion and snowfall derived from reanalysis data. Overall, the reconstruction is in good agreement with direct measurements of snow depth from Operation IceBridge, slightly underestimating mean IceBridge snow depth. At the local scale (10 km), the reconstruction is more skilled than a weighted climatology over first year ice, but underestimates deeper snow over multiyear ice. Reconstructions of single buoy snow depths are unskilled, but the reconstruction shows skill in simulating the mean snow depth across all buoys. Spring snow depths show a lowtohigh crossArctic gradient and tend to be greatest in the Atlantic sector of the eastern Arctic. The relationship between ice type (multiyear or firstyear ice) and snow depth previously documented in the western Arctic is not evident in the eastern Arctic. Using ice type to weight snow depths for satellite ice thickness retrievals may not be justifiable in the eastern Arctic. Reconstructed snow depth across the Arctic shows significant interannual variability, suggesting that use of a fixed snow depth climatology may lead to biases in retrieved ice thickness and its variability. However, interannual variability in panArctic mean snow depth is comparable or smaller than the uncertainty in both the reconstruction and IceBridge, highlighting the need for high accuracy snow depth products and reconstructions.

sea ice↗

Impact of Meteorological Factors on the Mesoscale Morphology of Cloud Streets during a Cold-Air Outbreak over the Western North Atlantic

Postfrontal clouds (PFC) are ubiquitous in the marine boundary layer, and their morphology is essential to estimating the radiation budget in weather and climate models. Here we examine the roles of sea surface temperature (SST) and meteorological factors in controlling the mesoscale morphology and evolution of shallow clouds associated with a cold-air outbreak that occurred on 1 March 2020 during phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). Our results show that the simulated PFC structure and ambient conditions by the Weather Research and Forecasting (WRF) Model are generally consistent with observations from GOES-16 and dropsonde measurements. We also examine the thermodynamical and dynamical influences in the cloud mesoscale morphology using WRF sensitivity experiments driven by two meteorological forcing datasets with different domain-mean SST and spatial gradients, which lead to dissimilar values of hydrometeor water path and cloud core fraction. The SST from ERA5 leads to weaker stability and higher inversion height than the SST from FNL does. In addition, the use of large-scale meteorological forcings from ERA5 yields a distinctive time evolution of wind direction shear in the inner domain, which favors the formation and persistence of longer cloud rolls. Both factors contribute to a change in the time evolution of domain-mean water path and cloud core fraction of cloud streets. Our study takes advantage of the simulation driven by the differences between two large-scale forcing datasets to illustrate the importance of SST and wind direction shear in the cloud street morphology in a realistic scenario

Jingyi Chen↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Impacts of the Mount Pinatubo eruption on ENSO in the GEOS seasonal-to-subseasonal forecasting system

The eruption of Mount Pinatubo in June 1991 introduced a perturbation of the Earth's global energy budget by increasing the stratospheric aerosol loading by an order of magnitude, with effects on the global climate. In this presentation we analyze the effects of the Mt. Pinatubo eruption on the seasonal forecast performed with Goddard Earth Observing System Seasonal-to Subseasonal (GEOS-S2S) system, an Earth System Model that includes an interactive ocean and a bulk aerosol model coupled to radiation. We performed 10-member ensembles for the year after the eruption (June 1991-May 1992) at ~0.5 horizontal resolution, with and without the inclusion of the Mt. Pinatubo eruption. In GEOS-S2S, the eruption leads to ta strengthening of El Nino peaking in January 1992. The strengthening is mainly due to the weakening of the trade winds, which is caused by a attening of the temperature gradient across the Pacic due to a differential response to the volcanic forcing between the central and eastern Pacic (ocean-dynamical thermostat). This response largely depends on the assumed size for the volcanic aerosols. Indeed, we performed simulations assuming a volcanic aerosol effective radius of 0.35 m (similar to tropospheric aerosol, and the default in GEOS) and 0.6 m (closer to observations of volcanic aerosol from Pinatubo-sized eruptions). We nd that in the latter case the tropical radiative forcing is lower, since smaller aerosols scatter shortwave radiation more eciently than larger ones. Accordingly, the impact on ENSO is not statistically signicant when a larger and more realistic particle radius is assumed.

Aquila, Valentina↗

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗