Surface freshwater storage and variability in the Amazon basin from multi-satellite observations, 1993-2007

The amount of water stored and moving through the surface water bodies of large river basins (river, floodplains, wetlands) plays a major role in the global water and biochemical cycles and is a critical parameter for water resources management. However, the spatiotemporal variations of these freshwater reservoirs are still widely unknown at the global scale. Here, we propose a hypsographic curve approach to estimate surface freshwater storage variations over the Amazon basin combining surface water extent from a multi-satellite-technique with topographic data from the Global Digital Elevation Model (GDEM) from Advance Spaceborne Thermal Emission and Reflection Radiometer (ASTER). Monthly surface water storage variations for 1993-2007 are presented, showing a strong seasonal and interannual variability, and are evaluated against in situ river discharge and precipitation. The basin-scale mean annual amplitude of similar to 1200 km(3) is in the range of previous estimates and contributes to about half of the Gravity Recovery And Climate Experiment (GRACE) total water storage variations. For the first time, we map the surface water volume anomaly during the extreme droughts of 1997 (October-November) and 2005 (September-October) and found that during these dry events the water stored in the river and floodplains of the Amazon basin was, respectively, similar to 230 (similar to 40%) and 210 (similar to 50%) km(3) below the 1993-2007 average. This new 15 year data set of surface water volume represents an unprecedented source of information for future hydrological or climate modeling of the Amazon. It is also a first step toward the development of such database at the global scale.

Copula-Based Modeling of TMI Brightness Temperature With Rainfall Type

Overland rain retrieval using spaceborne microwave radiometer offers a myriad of complications as land presents itself as a radiometrically warm and highly variable background. Hence, land rainfall algorithms of the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) have traditionally incorporated empirical relations of microwave brightness temperature (Tb) with rain rate, rather than relying on physically based radiative transfer modeling of rainfall (as implemented in the TMI ocean algorithm). In this paper, sensitivity analysis is conducted using the Spearman rank correlation coefficient as benchmark, to estimate the best combination of TMI low-frequency channels that are highly sensitive to the near surface rainfall rate from the TRMM Precipitation Radar (PR). Results indicate that the TMI channel combinations not only contain information about rainfall wherein liquid water drops are the dominant hydrometeors but also aid in surface noise reduction over a predominantly vegetative land surface background. Furthermore, the variations of rainfall signature in these channel combinations are not understood properly due to their inherent uncertainties and highly nonlinear relationship with rainfall. Copula theory is a powerful tool to characterize the dependence between complex hydrological variables as well as aid in uncertainty modeling by ensemble generation. Hence, this paper proposes a regional model using Archimedean copulas, to study the dependence of TMI channel combinations with respect to precipitation, over the land regions of Mahanadi basin, India, using version 7 orbital data from the passive and active sensors on board TRMM, namely, TMI and PR. Studies conducted for different rainfall regimes over the study area show the suitability of Clayton and Gumbel copulas for modeling convective and stratiform rainfall types for the majority of the intraseasonal months. Furthermore, large ensembles of TMI Tb (from the most sensitive TMI channel combination) were generated conditional on various quantiles (25th, 50th, 75th, and 95th) of the convective and the stratiform rainfall. Comparatively greater ambiguity was observed to model extreme values of the convective rain type. Finally, the efficiency of the proposed model was tested by comparing the results with traditionally employed linear and quadratic models. Results reveal the superior performance of the proposed copula-based technique.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Modeling of solar radiation management: a comparison of simulations using reduced solar constant and stratospheric sulphate aerosols

The climatic effects of Solar Radiation Management (SRM) geoengineering have been often modeled by simply reducing the solar constant. This is most likely valid only for space sunshades and not for atmosphere and surface based SRM methods. In this study, a global climate model is used to evaluate the differences in the climate response to SRM by uniform solar constant reduction and stratospheric aerosols. Our analysis shows that when global mean warming from a doubling of CO2 is nearly cancelled by both these methods, they are similar when important surface and tropospheric climate variables are considered. However, a difference of 1 K in the global mean stratospheric (61-9.8 hPa) temperature is simulated between the two SRM methods. Further, while the global mean surface diffuse radiation increases by similar to 23 % and direct radiation decreases by about 9 % in the case of sulphate aerosol SRM method, both direct and diffuse radiation decrease by similar fractional amounts (similar to 1.0 %) when solar constant is reduced. When CO2 fertilization effects from elevated CO2 concentration levels are removed, the contribution from shaded leaves to gross primary productivity (GPP) increases by 1.8 % in aerosol SRM because of increased diffuse light. However, this increase is almost offset by a 15.2 % decline in sunlit contribution due to reduced direct light. Overall both the SRM simulations show similar decrease in GPP (similar to 8 %) and net primary productivity (similar to 3 %). Based on our results we conclude that the climate states produced by a reduction in solar constant and addition of aerosols into the stratosphere can be considered almost similar except for two important aspects: stratospheric temperature change and the consequent implications for the dynamics and the chemistry of the stratosphere and the partitioning of direct versus diffuse radiation reaching the surface. Further, the likely dependence of global hydrological cycle response on aerosol particle size and the latitudinal and height distribution of aerosols is discussed.

Modeling the co-transport of viruses and colloids in unsaturated porous media

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations,are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Cluysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWL A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AIM, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWL The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SW!, maximum adsorption capacity of the AWI for viruses and inlet virus concentration. (C) 2015 Elsevier B.V. All rights reserved.

Modeling of the non-isothermal liquid droplet impact on a heated solid substrate with heterogeneous wettability

A comprehensive numerical investigation on the impingement and spreading of a non-isothermal liquid droplet on a solid substrate with heterogeneous wettability is presented in this work. The time-dependent incompressible Navier-Stokes equations are used to describe the fluid flow in the liquid droplet, whereas the heat transfer in the moving droplet and in the solid substrate is described by the energy equation. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the time-dependent incompressible Navier-Stokes equation and the energy equation in the time-dependent moving domain. Moreover, the Marangoni convection is included in the variational form of the Navier-Stokes equations without calculating the partial derivatives of the temperature on the free surface. The heterogeneous wettability is incorporated into the numerical model by defining a space-dependent contact angle. An array of simulations for droplet impingement on a heated solid substrate with circular patterned heterogeneous wettability are presented. The numerical study includes the influence of wettability contrast, pattern diameter, Reynolds number and Weber number on the confinement of the spreading droplet within the inner region, which is more wettable than the outer region. Also, the influence of these parameters on the total heat transfer from the solid substrate to the liquid droplet is examined. We observe that the equilibrium position depends on the wettability contrast and the diameter of the inner surface. Consequently. the heat transfer is more when the wettability contrast is small and/or the diameter of inner region is large. The influence of the Weber number on the total heat transfer is more compared to the Reynolds number, and the total heat transfer increases when the Weber number increases. (C) 2015 Elsevier Ltd. All rights reserved.

Microkinetic Modeling of CO Oxidation on Ionic Palladium-Substituted Ceria

This study presents a plausible dual-site mechanism and microkinetic model for CO oxidation over palladium-substituted ceria incorporating the theoretical oxygen storage capacity of different-catalysts into the kinetic model. A rate expression without prior assumption of rate-determining steps has been developed for the proposed microkinetic model using reaction route analysis. Experiments were conducted using various percentages of palladium in ceria that were synthesized by solution combustion. Obtained catalysts were characterized by X-ray diffraction, X-ray photoelectron spectra, and Brunauer-Emmett-Teller surface area measurements. A detailed mechanism was, developed, and the kinetic parameters and rate expression were validated with the conversion data obtained in the presence of the catalysts. Furthermore, a reduced rate expression based on rate-determining step and most abundant reactive intermediate approximation was obtained and tested against the original rate expression for different experimental conditions. From the results obtained it was concluded that the simulated rate predictions matched the experimental trend with reasonable accuracy, validating the kinetic parameters proposed it this study.

Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

COUPLED OPTICAL-THERMAL-FLUID MODELING OF A DIRECTLY HEATED TUBULAR SOLAR RECEIVER FOR SUPERCRITICAL CO2 BRAYTON CYCLE

Recent studies have evaluated closed-loop supercritical carbon dioxide (s-CO2) Brayton cycles to be a higher energy density system in comparison to conventional superheated steam Rankine systems. At turbine inlet conditions of 923K and 25 MPa, high thermal efficiency (similar to 50%) can be achieved. Achieving these high efficiencies will make concentrating solar power (CSP) technologies a competitive alternative to current power generation methods. To incorporate a s-CO2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 923 K (at 25 MPa) is necessary. The s-CO2 will need to increase in temperature by similar to 200 K as it passes through the solar receiver to satisfy the temperature requirements of a s-CO2 Brayton cycle with recuperation and recompression. In this study, an optical-thermal-fluid model was developed to design and evaluate a tubular receiver that will receive a heat input similar to 2 MWth from a heliostat field. The ray-tracing tool SolTrace was used to obtain the heat-flux distribution on the surfaces of the receiver. Computational fluid dynamics (CFD) modeling using the Discrete Ordinates (DO) radiation model was used to predict the temperature distribution and the resulting receiver efficiency. The effect of flow parameters, receiver geometry and radiation absorption by s-CO2 were studied. The receiver surface temperatures were found to be within the safe operational limit while exhibiting a receiver efficiency of similar to 85%.