Spatial interpolation of snow water equivalency using surface observations and remotely sensed images of snow‐covered area

As demand for water continues to escalate in the western Unites States, so does the need for accurate monitoring of the snowpack in mountainous areas. In this study, we describe a simple methodology for generating gridded‐estimates of snow water equivalency (SWE) using both surface observations of SWE and remotely sensed estimates of snow‐covered area (SCA). Multiple regression was used to quantify the relationship between physiographic variables (elevation, slope, aspect, clear‐sky solar radiation, etc.) and SWE as measured at a number of sites in a mountainous basin in south‐central Idaho (Big Wood River Basin). The elevation of the snowline, obtained from the SCA estimates, was used to constrain the predicted SWE values. The results from the analysis are encouraging and compare well to those found in previous studies, which often utilized more sophisticated spatial interpolation techniques. Cross‐validation results indicate that the spatial interpolation method produces accurate SWE estimates [mean R² = 0·82, mean mean absolute error (MAE) = 4·34 cm, mean root mean squared error (RMSE) = 5·29 cm]. The basin examined in this study is typical of many mid‐elevation mountainous basins throughout the western United States, in terms of the distribution of topographic variables, as well as the number and characteristics of sites at which the necessary ground data are available. Thus, there is high potential for this methodology to be successfully applied to other mountainous basins. Copyright © 2010 John Wiley & Sons, Ltd.     

Thresholds for post-wildfire debris flows: Insights from the Pinal Fire,Arizona, USA

Wildfire significantly alters the hydrologic properties of a burned area, leading to increases in overland flow, erosion, and the potential for runoff-generated debris flows. The initiation of debris flows in recently burned areas is well characterized by rainfall intensity-duration (ID) thresholds. However, there is currently a paucity of data quantifying the rainfall intensities required to trigger post-wildfire debris flows, which limits our understanding of how and why rainfall ID thresholds vary in different climatic and geologic settings. In this study, we monitored debris-flow activity following the Pinal Fire in central Arizona, which differs from both a climatic and hydrogeomorphic perspective from other regions in the western United States where ID thresholds for post-wildfire debris flows are well established, namely the Transverse Ranges of southern California. Since the peak rainfall intensity within a rainstorm may exceed the rainfall intensity required to trigger a debris flow, the development of robust rainfall ID thresholds requires knowledge of the timing of debris flows within rainstorms. Existing post-wildfire debris-flow studies in Arizona only constrain the peak rainfall intensity within debris-flow-producing storms, which may far exceed the intensity that actually triggered the observed debris flow. In this study, we used pressure transducers within five burned drainage basins to constrain the timing of debris flows within rainstorms. Rainfall ID thresholds derived here from triggering rainfall intensities are, on average, 22 mm h⁻¹ lower than ID thresholds derived under the assumption that the triggering intensity is equal to the maximum rainfall intensity recorded during a rainstorm. We then use a hydrologic model to demonstrate that the magnitude of the 15-min rainfall ID threshold at the Pinal Fire site is associated with the rainfall intensity required to exceed a recently proposed dimensionless discharge threshold for debris-flow initiation. Model results further suggest that previously observed differences in regional ID thresholds between Arizona and the San Gabriel Mountains of southern California may be attributed, in large part, to differences in the hydraulic properties of burned soils. © 2019 John Wiley & Sons, Ltd.     

Using analytic element models to delineate drinking water source protection areas

Since 1999, Ohio EPA hydrogeologists have used two analytic element models (AEMs), the proprietary software GFLOW and U.S. EPA's WhAEM, to delineate protection areas for 535 public water systems. Both models now use the GFLOW2001 solution engine, integrate well with Geographic Information System (GIS) technology, have a user-friendly graphical interface, are capable of simulating a variety of complex hydrogeologic settings, and do not rely upon a model grid. These features simplify the modeling process and enable AEMs to bridge the gap between existing simplistic delineation methods and more complex numerical models. Ohio EPA hydrogeologists demonstrated that WhAEM2000 and GFLOW2000 were capable of producing capture zones similar to more widely accepted models by applying the AEMs to eight sites that had been previously delineated using other methods. After the Ohio EPA delineated protection areas using AEMs, more simplistic delineation methods used by other states (volumetric equation and arbitrary fixed radii) were applied to the same water systems to compare the differences between various methods. GIS software and two-tailed paired t-tests were used to quantify the differences in protection areas and analyze the data. The results of this analysis demonstrate that AEMs typically produce significantly different protection areas than the most simplistic delineation methods, in terms of total area and shape. If the volumetric equation had been used instead of AEMs, Ohio would not have protected 265 km2 of critical upgradient area and would have overprotected 269 km2 of primarily downgradient land. Since an increasing number of land-use restrictions are being tied to drinking water protection areas, this analysis has broad policy implications.     

Recherches complémentaires sur le bases du calculateur météorologique «Temp»

Résumé La formule de base, traduisant une propriété analytique d'une classe très générale de fonctions, est un corollaire du théorème fondamental démontré dans un mémoire précédent, d'après lequel, étant donnés une fonction continue,p(, ,t) des points (, ) d'une surface régulière fermée et du temps et le champ d'un vecteur vitesse de transfert ou d'advection tangent à et ayant des lignes de flux fermées et régulières, il existe un opérateur spatial, linéaire, non singulierA tel que la fonctionA(p+Const.) soit purement advective par rapport a (sans creusement ni comblement). Ce théorème peut être exprimé par l'équation , où est un opérateur spatial, linéaire et non singulier, fonction deA.La détermination de peut être faite, soit en comparant deux formes différentes de la solution générale de l'équation en , soit en utilisant un raisonnement a priori très simple. On arrive ainsi au résultat pour un certain scalaireu(, ).Dans le cas oùp(, ,t) est la perturbation de la pression sur la surface du géoïde l'équation résulte aussi, comme nous l'avons montré dans le mémoire précédent, de notre théorie hydrodynamique des perturbations. On montre ici que la même équation peut encore être déduite de l'équation de continuité associée à la condition d'équilibre quasi statique selon la verticale.Comme applications de la formule de base (solution générale de l'équation enM), on étudie les problèmes suivants: 1^o creusement et comblement en général; 2^o creusement et comblement des centres et des cols; 3^o mouvement des centres et des cols; 4^o instabilité d'un champ moyen; 5^o propriétés spatiales des champsp(, ,t) et des vecteurs d'advection analytiques.Après une discussion des erreurs de la prévision d'un champp(, ,t) par la formule de base, du fait des erreurs des observations et du fonctionnement du calculateur, on examine quelques particularités du transfert ou advection d'un champf ₀(, ) par le vecteur . Enfin, le dernier chapitre du mémoire donne des éclaircissements complémentaires sur la structure du calculateur électronique «Temp» (qui effectue automatiquement les opérations mathématiques de la formule de base) et expose l'état actuel de sa construction.

---

Summary The basic formula, expressing an analytical property of a very general class of functions, is a corollary of the fundamental theorem, proved in a previous paper, according to which, given a functionp(, ,t) of the points (, ) of a closed regular surface and of the time, and a transfer or advection velocity vector tangent to and having regular closed streamlines, there is a spatial, linear, non singular operatorA such thatA(p+const.) is a purely advective function in respect to (no deepening). This theorem can be expressed by the equation where is a spatial, linear, non singular operator depending onA.The determination of can be attained, either by the comparison of two different forms of the general solution of the -equation, or by a simple a priori reasonning. The conclusion is thus reached that for a certain scalaru(, ).Whenp(, ,t) is the pressure perturbation at sea level, it was shown, in the preceding paper, that the equation can also be derived from our hydrodynamical perturbation theory. We now show that for this particular case, the same equation is also a consequence of the equation of continuity together with the condition of quasi statical vertical equilibrium.The following problems are then analysed by means of the basic formula: 1^o deepening and filling in general; 2^o deepening and filling of the centres and cols; 3^o motion of the centres and cols; 4^o instability of a mean field; 5^o spatial properties of the analytical fields and advection vectors .The errors in the forecast of a field,p(, ,t) by means of the basic formula, due to the observational and computational errors, are discussed, and some peculiarities of the transfer or advection of a fieldf ₀(, ) by are examined. Finally, complementary points are disclosed on the structure of the electronic computer «Temp» which performs automatically the mathematical operations of the basic formula, and a brief report is given of the present state of its construction.

     

Nitrogen uptake and denitrification in restored and unrestored streams in urban Maryland, USA

There is growing interest in rates of nitrate uptake and denitrification in restored streams to better understand the effects of restoration on nitrogen processing. This study quantified nitrate uptake in two restored and two unrestored streams in Baltimore, Maryland, USA using nitrate additions, denitrification enzyme assays, and a ¹⁵N isotope tracer addition in one of the urban restored streams, Minebank Run. Restoration included either incorporation of stormwater ponds below a storm drain and catch basins to attenuate flow or hydrologic “reconnection” of a stream channel to its floodplain. Stream restoration was conducted for restoring aging sanitary and bridge infrastructure and introducing some stormwater management in watersheds developed prior to current regulations. Denitrification potential in sediments was variable across streams, whereas nitrate uptake length appeared to be significantly correlated to surface water velocity, which was low in the restored streams during summer baseflow conditions. Uptake length of NO₃ ⁻–N in Minebank Run estimated by ¹⁵N tracer addition was 556 m. Whole stream denitrification rates in Minebank Run were 153 mg NO₃ ⁻–N m⁻² day⁻¹, and approximately 40% of the daily load of nitrate was estimated to be removed via denitrification over a distance of 220.5 m in a stream reach designed to be hydrologically “connected” to its floodplain. Increased hydrologic residence time in Minebank Run during baseflow likely influenced rates of whole stream denitrification, suggesting that hydrologic residence time may be a key factor influencing N uptake and denitrification. Restoration approaches that increase hydrologic “connectivity” with hyporheic sediments and increase hydrologic residence time may be useful for stimulating denitrification. More work is necessary, however, to examine changes in denitrification rates in restored streams across different seasons, variable N loads, and in response to the “flashy” hydrologic flow conditions during storms common in urban streams.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Ceres's global and localized mineralogical composition determined by Dawn's Visible and Infrared Spectrometer (VIR)

The Visible and Infrared Spectrometer (VIR) instrument on the Dawn mission observed Ceres’s surface at different spatial resolutions, revealing a nearly uniform global distribution of surface mineralogy. Clearly, Ceres experienced extensive water‐related processes and chemical differentiation. The surface is mainly composed of a dark component (carbon, magnetite?), Mg‐phyllosilicates, ammoniated clays, carbonates, and salts. The observed species suggest endogenous, global‐scale aqueous alteration. While mostly uniform at regional scale, Ceres’s surface shows small localized areas with different species and/or variations in abundances. Few local exposures of water ice are seen, especially at higher latitudes. Sodium carbonates have been identified in several areas on the surface, notably in Occator bright faculae. Organic matter has also been discovered in several places, most conspicuously in a large area close to the Ernutet crater. The observed mineralogies, with the presence of ammoniated species and sodium salts, have a strong resemblance to materials found on other bodies of the outer solar system, such as Enceladus. This poses some questions about the original material from which Ceres accreted, suggesting a colder environment for such material with respect to Ceres’s present position.     

Multi-wavelength observations of the 1998 September 27 flare spray

We report on observations of a large eruptive event associated with a flare that occurred on 27 September 1998 made with the Richard B. Dunn Solar Telescope at Sacramento Peak Observatory (several wave bands including off-line-center H), in soft and hard X-rays (GOES and BATSE), and in several TRACE wave bands (including Feix/x 171 Å, Fexii 195 Å, and Civ 1550 Å). The flare initiation is signaled by two H foot-point brightenings which are closely followed by a hard X-ray burst and a subsequent gradual increase in other wavelengths. The flare light curves show a complicated, three-component structure which includes two minor maxima before the main GOES class C5.2 peak after which there is a characteristic exponential decline. During the initial stages, a large spray event is observed within seconds of the hard X-ray burst which can be directly associated with a two-ribbon flare in H. The emission returns to pre-flare levels after about 35 min, by which time a set of bright post-flare loops have begun to form at temperatures of about 1.0–1.5 MK. Part of the flare plasma also intrudes into the penumbra of a large sunspot, generally a characteristic of very powerful flares, but the flare importance in GOES soft X-rays is in fact relatively modest. Much of the energy appears to be in the form of a second ejection which is observed in optical and ultraviolet bands, traveling out via several magnetic flux tubes from the main flare site (about 60° from Sun center) to beyond the limb.     

Geometric consistency for regression model estimation and testing in climatology and meteorology

Abstract

The effects of outliers on linear regression are examined. The sensitivity of classical least‐squares (LS) procedures to outliers is shown to be associated with the geometric inconsistency between the data space and the analysis space. This is illustrated for both estimation and inference. A geometrically consistent procedure based on the Euclidean distance is proposed. This procedure involves the least absolute deviation (LAD) regression and a new permutation test for matched pairs (PTMP). Comparisons made with LS techniques demonstrate that the proposed procedure is more resistant to the existence of outliers in the data set and leads to more intuitive results. Applications and illustrations using meteorological and climatological data are also discussed.