Measurements of soil moisture are important for studies of climate and weather forecasting, flood prediction, and aquifer
recharge studies. Although soil moisture measurement networks exist, most are sparsely distributed and lack standardized instrumentation.
Measurements of soil moisture from satellites have extremely large spatial footprints (40–60 km). A methodology is described
here that uses existing networks of continuously-operating GPS receivers to measure soil moisture fluctuations. In this technique,
incoming signals are reflected off and attenuated by the ground before reception by the GPS receiver. These multipath reflections
directly affect signal-to-noise ratio (SNR) data routinely collected by GPS receivers, creating amplitude variations that
are a function of ground reflectivity and therefore soil moisture content. After describing this technique, multipath reflection
amplitudes at a GPS site in Tashkent, Uzbekistan are compared to estimates of soil moisture from the Noah land surface model.
Although the GPS multipath amplitudes and the land surface model are uncalibrated, over the 70-day period studied, they both
rise sharply following each rainfall event and slowly decrease over a period of ∼10 days. 相似文献
Climatic Change - Salinity response to sea-level rise is evaluated for a low-gradient, tidally active estuary, the lower St. Johns River, Florida. A high-resolution numerical model is forced by... 相似文献
When formulating a hydrologic model, scientists rely on parameterizations of multiple processes based on field data, but literature review suggests that more frequently people select parameterizations that were included in pre-existing models rather than re-evaluating the underlying field experiments. Problems arise when limited field data exist, when “trusted” approaches do not get reevaluated, and when sensitivities fundamentally change in different environments. The physics and dynamics of snow interception by conifers is just such a case, and it is critical to simulation of the water budget and surface albedo. The most commonly used interception parameterization is based on data from four trees from one site, but results from this field study are not directly transferable to locations with relatively warmer winters, where the dominant processes differ dramatically. Here, we combine a literature review with model experiments to demonstrate needed improvements. Our results show that the choice of model form and parameters can vary the fraction of snow lost through interception by as much as 30%. In most simulations, the warming of mean winter temperatures from −7 to 0°C reduces the modelled fraction of snow under the canopy compared to the open, but the magnitude of simulated decrease varies from about 10% to 40%. The range of results is even larger when considering models that neglect the melting of in-canopy snow in higher-humidity environments where canopy sublimation plays less of a role. Thus, we recommend that all models represent canopy snowmelt and include representation of increased loading due to increased adhesion and cohesion when temperatures rise from −3 to 0°C. In addition to model improvements, field experiments across climates and forest types are needed to investigate how to best model the combination of dynamically changing forest cover and snow cover to better understand and predict changes to albedo and water supplies. 相似文献
A quantitative physical model is presented which includes the factors that control the presence, or absence, of internally derived excess 40Ar or excess 4He in geological systems. In particular, the model incorporates the transport and partitioning properties of the rock surrounding the mineral of thermochronologic interest and illuminates the related effects on the amount of excess 40Ar or 4He preserved in the system. Modeling of a simplified 1-D rock column bounded by an external sink for 40Ar or 4He shows that a steady-state excess 40Ar or 4He profile develops, the magnitude of which is determined by a system parameter called the ‘transmissive timescale’, τT. The characteristic time required to reach this steady state depends upon τT and the ‘total local sink capacity’, TLSC, wherein the important role of local matrix mineral and fluid phases is incorporated. Together, these two system parameters (τT and TLSC) determine the evolution of excess 40Ar or 4He buildup within a system above the closure temperatures of all minerals involved. An analytical expression for the 1-D system describing the evolution of excess 40Ar (or by analogy 4He) in a particular potassium-bearing (or U-Th-bearing) mineral located at a distance, L, from an external sink has been derived empirically from model results:
The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels. 相似文献
We present a compact, high-order Richards’ equation solver using a local discontinuous Galerkin finite element method in space and a dual-time stepping method in time. Dual-time stepping methods convert a transient problem to a steady state problem, enabling direct evaluation of residual terms and resolve implicit equations in a step-wise manner keeping the method compact and amenable to parallel computing. Verification of our solver against an analytical solution shows high-order error convergence and demonstrates the solvers ability to maintain high accuracy using low spatial resolution; the method is robust and accurately resolves numerical solutions with time steps that are much larger than what is normally required for lower-order implicit schemes. Resilience of our solver (in terms of nonlinear convergence) is demonstrated in ponded infiltration into homogeneous and layered soils, for which HYDRUS-1D solutions are used as qualitative references to gauge performance of two slope limiting schemes.
It has been established that idealized western boundary currents, which encounter a gap in their supporting boundary, will assume one of two dominant steady states: a loop current state and a gap leaping state, and that transitions between these states display hysteresis. However, a question of whether the idealized geometries considered to date apply to the Gulf of Mexico Loop Current (LC) remained. Here, the nonlinear potential vorticity advection-diffusions equations are solved, for Gulf of Mexico topography, using Newton’s method. We demonstrate that, in application to the LC in the Gulf of Mexico, the original conclusions do hold and additionally describe peculiarities of the more realistic steady states. The existence of our numerically calculated steady LC states in the actual Gulf of Mexico are supported through analysis of historical sea surface height data, and implications of our results for LC modeling and forecasting are discussed.