Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Snow model”

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 1,513 records · Page 84

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗

Sensitivity of Climate to Changes in NDVI

The sensitivity of global and regional climate to changes in vegetation density is investigated using a coupled biosphere-atmosphere model. The magnitude of the vegetation changes and their spatial distribution are based on natural decadal variability of the normalized difference vegetation index (ndvi). Different scenarios using maximum and minimum vegetation cover were derived from satellite records spanning the period 1982-1990. Albedo decreased in the northern latitudes and increased in the tropics with increased ndvi. The increase in vegetation density revealed that the vegetation's physiological response was constrained by the limits of the available water resources. The difference between the maximum and minimum vegetation scenarios resulted in a 46% increase in absorbed visible solar radiation and a similar increase in gross photosynthetic C02 uptake on a global annual basis. This caused the canopy transpiration and interception fluxes to increase, and reduced those from the soil. The redistribution of the surface energy fluxes substantially reduced the Bowen ratio during the growing season, resulting in cooler and moister near-surface climate, except when soil moisture was limiting. Important effects of increased vegetation on climate are : (1) A cooling of about 1.8 K in the northern latitudes during the growing season and a slight warming during the winter, which is primarily due to the masking of high albedo of snow by a denser canopy. and (2) A year round cooling of 0.8 K in the tropics. These results suggest that increases in vegetation density could partially compensate for parallel increases in greenhouse warming . Increasing vegetation density globally caused both evapotranspiration and precipitation to increase. Evapotranspiration, however increased more than precipitation resulting in a global soil-water deficit of about 15 %. A spectral analysis on the simulated results showed that changes in the state of vegetation could affect the low-frequency modes of the precipitation distribution and might reduce its low frequency variability in the tropics while increasing it in northern latitudes.

Bounoua, L.↗

Greenland Ice Sheet Surface Temperature, Melt, and Mass Loss: 2000-2006

Extensive melt on the Greenland Ice Sheet has been documented by a variety of ground and satellite measurements in recent years. If the well-documented warming continues in the Arctic, melting of the Greenland Ice Sheet will likely accelerate, contributing to sea-level rise. Modeling studies indicate that an annual or summer temperature rise of 1 C on the ice sheet will increase melt by 20-50% therefore, surface temperature is one of the most important ice-sheet parameters to study for analysis of changes in the mass balance of the ice-sheet. The Greenland Ice Sheet contains enough water to produce a rise in eustatic sea level of up to 7.0 m if the ice were to melt completely. However, even small changes (centimeters) in sea level would cause important economic and societal consequences in the world's major coastal cities thus it is extremely important to monitor changes in the ice-sheet surface temperature and to ultimately quantify these changes in terms of amount of sea-level rise. We have compiled a high-resolution, daily time series of surface temperature of the Greenland Ice Sheet, using the I-km resolution, clear-sky land-surface temperature (LST) standard product from the Moderate-Resolution Imaging Spectroradiometer (MODIS), from 2000 - 2006. We also use Gravity Recovery and Climate Experiment (GRACE) data, averaged over 10-day periods, to measure change in mass of the ice sheet as it melt and snow accumulates. Surface temperature can be used to determine frequency of surface melt, timing of the start and the end of the melt season, and duration of melt. In conjunction with GRACE data, it can also be used to analyze timing of ice-sheet mass loss and gain.

Hall, Dorothy K.↗

Low-Level Liquid-Bearing Clouds Contribute to Seasonal Lower Atmosphere Stability and Surface Energy Forcing over a High-Mountain Watershed Environment

Measurements of atmospheric structure and surface energy budgets distributed along a high-altitude moun tain watershed environment near Crested Butte, Colorado, from two separate, but coordinated, field campaigns, Surface Atmosphere Integrated field Laboratory (SAIL) and Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH), are analyzed. This study identifies similarities and differences in how clouds influence the radiative budget over one snow-free summer season (2022) and two snow-covered seasons (2021/22; 2022/23) for this alpine location. A relationship between lower-tropospheric stability stratification and longwave radiative flux from the presence or absence of clouds is identified. When low clouds persisted, often with signatures of supercooled liquid in win ter, the lower troposphere experienced weaker stability, while radiatively clear skies that are less likely to be influenced by liquid droplets were associated with appreciably stronger lower-tropospheric stratification. Corresponding surface turbu lent heat fluxes partitioned differently based upon the cloud–stability stratification regime derived from early morning radiosounding profiles. Combined with the differences in the radiative budget largely resulting from dramatic seasonal dif ferences in surface albedo, the lower atmosphere stratification, surface energy budget, and near-surface thermodynamics are shown to be modified by the effective longwave radiative forcing of clouds. The diurnal evolution of thermodynamics and surface energy components varied depending on the early morning stratification state. Thus, the importance of quies cent versus synoptically active large-scale meteorology is hypothesized as a critical forcing for cloud properties and associ ated surface energy budget variations. The physical relationships between clouds, radiation, and stratification can provide a useful suite of metrics for process understanding and to evaluate numerical models in such an undersampled, highly com plex terrain environment.

54 ENVIRONMENTAL SCIENCES↗

TRMM Precipitation Radar and Microwave Imager Observations of Convective and Stratiform Rain and Their Theoretical Implications

Theoretical models of radiative transfer are developed to simulate the 85 GHz brightness temperature, T85, observed by the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer as a function of rain rate. These simulations are performed separately over regions of convective and stratiform rain. For the purpose of constructing vertical profiles of hydrometeors in these regions, guidance is taken from TRMM Precipitation Radar (PR) observations, and from recent investigations of LDR and Z(sub DR) measurements made by land-based polarimetric radars. We find that for a given rain rate, the extinction in 85 GHz due to hydrometeors above the freezing level is relatively weak in convective regions compared to that in stratiform regions. The hydrometeor profile above the freezing level responsible for this weak extinction in convective regions is inferred, from LDR and Z(sub DR) observations and theoretical considerations, to contain two layers: 1) a layer of two km thickness with mixed phase particles, liquid drops and graupel, directly above the freezing level, and 2) a layer of graupel extending from the top of the mixed-phase layer to the cloud top. Strong extinction in the stratiform regions is inferred to result from slowly-falling, low-density ice aggregates (snow) above the freezing level. These theoretical results are consistent with T85 measured by TMI and the rain rate deduced from PR for convective and stratiform rain regions. On the basis of this study we infer that the accuracy of the rain rate sensed by TMI depends critically on the specification of the convective or stratiform nature of the rain.

Prabhakara, C.↗

Applications of ArcticDEM for measuring volcanic dynamics, landslides, retrogressive thaw slumps, snowdrifts, and vegetation heights

Topographical changes are of fundamental interest to a wide range of Arctic science disciplines faced with the need to anticipate, monitor, and respond to the effects of climate change, including geohazard management, glaciology, hydrology, permafrost, and ecology. This study demonstrates several geomorphological, cryospheric, and biophysical applications of ArcticDEM – a large collection of publicly available, time-dependent digital elevation models (DEMs) of the Arctic. Our study illustrates ArcticDEM's applicability across different disciplines and five orders of magnitude of elevation derivatives, including measuring volcanic lava flows, ice cauldrons, post-failure landslides, retrogressive thaw slumps, snowdrifts, and tundra vegetation heights. We quantified surface elevation changes in different geological settings and conditions using the time series of ArcticDEM. Following the 2014–2015 Bárðarbunga eruption in Iceland, ArcticDEM analysis mapped the lava flow field, and revealed the post-eruptive ice flows and ice cauldron dynamics. The total dense-rock equivalent (DRE) volume of lava flows is estimated to be (1431 ± 2) million m 3 . Then, we present the aftermath of a landslide in Kinnikinnick, Alaska, yielding a total landslide volume of (400 ± 8) × 103 m 3 and a total area of 0.025 km 2 . ArcticDEM is further proven useful for studying retrogressive thaw slumps (RTS). The ArcticDEM-mapped RTS profile is validated by ICESat-2 and drone photogrammetry resulting in a standard deviation of 0.5 m. Volume estimates for lake-side and hillslope RTSs range between 40,000 ± 9000 m 3 and 1,160,000 ± 85,000 m 3 , highlighting applicability across a range of RTS magnitudes. A case study for mapping tundra snow demonstrates ArcticDEM's potential for identifying high-accumulation, late-lying snow areas. The approach proves effective in quantifying relative snow accumulation rather than absolute values (standard deviation of 0.25 m, bias of −0.41 m, and a correlation coefficient of 0.69 with snow depth estimated by unmanned aerial systems photogrammetry). Furthermore, ArcticDEM data show its feasibility for estimating tundra vegetation heights with a standard deviation of 0.3 m (no bias) and a correlation up to 0.8 compared to the light detection and ranging (LiDAR). The demonstrated capabilities of ArcticDEM will pave the way for the broad and pan-Arctic use of this new data source for many disciplines, especially when combined with other imagery products. The wide range of signals embedded in ArcticDEM underscores the potential challenges in deciphering signals in regions affected by various geological processes and environmental influences.

Chunli Dai↗

Ecological and genomic variation in ectomycorrhizal fungal exploration types

Ectomycorrhizal fungi (EMF) produce mycelia with variable extension and complexity, which can be classified according to soil ‘exploration types’ (ETs). ETs have received attention as one of the few mycorrhizal trait frameworks, but without an empirical classification of ET functional diversity and environmental preferences, understanding and interpreting EMF biogeographic patterns has been difficult. We conducted a synthesis combining: comparative EMF genomics to describe functional divergence in decomposition and nutrient cycling genes across ETs; and EMF trait distribution modeling across continental Europe, pairing soil and root EMF surveys to establish biogeographic ET niche profiles. We demonstrate a signature of ETs encoded in EMF genomes, which is independent from phylogeny and linked to biomass production strategies. EMF ET relative abundances were separated by soil, root, and dominant tree leaf type habitats and exhibited unique correlations with forest biotic (e.g. plant productivity and plant pathogen densities) and abiotic (e.g. nitrogen deposition and soil pH) conditions. These findings support a theory that EMF niche partitioning can be partially explained by extraradical mycelial traits, with underlying variation in ET biogeography likely arising from distinct decomposition and nutrient cycling potentials. We also identify important limitations to this trait framework and provide a guided outlook for future research.

biogeography↗

Contributions of GRACE to Climate Monitoring

The NASA/German Gravity Recovery and Climate Experiment (GRACE) was launched in March 2002. Rather than looking downward, GRACE continuously monitors the locations of and precise distance between twin satellites which orbit in tandem about 200 km apart. Variations in mass near Earth's surface cause heterogeneities in its gravity field, which in turn affect the orbits of satellites. Thus scientists can use GRACE data to map Earth's gravity field with enough accuracy to discern month to month changes caused by ocean circulation and redistribution of water stored on and in the land. Other gravitational influences, such as atmospheric circulation, post-glacial rebound, and solid earth movements are either independently determined and removed or are negligible on a monthly to sub-decadal timescale. Despite its coarse spatial (>150,000 sq km at mid-latitudes) and temporal (approx monthly) resolutions, GRACE has enabled significant advancements in the oceanic, hydrologic, and cryospheric science, and has great potential for climate monitoring, because it is the only global observing system able to measure ocean bottom pressures, total terrestrial water storage, and ice mass changes. The best known GRACE results are estimates of Greenland and Antarctic ice sheet loss rates. Previously, scientists had estimated ice mass losses using ground and satellite based altimetry and surface mass balance estimates based on snowfall accumulation and glacier discharge. While such measurements are still very useful for their spatial detail, they are imperfectly correlated with large-scale ice mass changes, due to snow and ice compaction and incomplete spatial coverage. GRACE enables scientists to generate monthly time series of Greenland and Antarctic ice mass, which have confirmed the shrinking of the polar ice sheets, one of the most obvious and indisputable manifestations of climate change. Further, GRACE has located and quantified hot spots of ice loss in southeastern Greenland and western Antarctica. For 2002 to present, the rate of ice mass loss has been 200 to 300 GT/yr in Greenland and 70 to 210 GT/yr in Antarctica, and some scientists are suggesting that the rates are accelerating. Similarly, GRACE has been used to monitor mass changes in alpine glaciers. Tamisiea et al. first characterized glacier melt along the southern coast of Alaska, more recently estimated to be occurring at a rate of 84 GT/yr. Chen et al. estimated that Patagonian glaciers are melting at a rate of 28 GT/yr, and estimated that the high mountains of central Asia lose ice at a rate of 47 GT/yr. Tapley et al. and Wahr et al. presented the first GRACE based estimates of changes in column-integrated terrestrial water storage (TWS; the sum of ground-water, soil moisture, surface waters, snow, ice, and water stored in vegetation) at continental scales. Since then, dozens of studies have shown that GRACE based estimates of regional to continental scale TWS variations agree with independent information, and some innovative uses of GRACE data have been developed. Rodell et al. (2004) and Swenson and Wahr (2006) demonstrated that by combining GRACE derived terrestrial water storage changes with observations of precipitation and runoff in a river basin scale water budget, it was possible to produce new estimates of evapotranspiration and atmospheric moisture convergence, essential climate variables that are difficult to estimate accurately. Similarly, GRACE has been used to constrain estimates of global river discharge and the contribution of changes in TWS to sea level rise. Crowley et al. observed a negative correlation between interannual TWS anomalies in the Amazon and the Congo River basin. Yeh et al. and Rodell et al. estimated regionally averaged groundwater storage variations based on GRACE and auxiliary observations. Rodell et al. and Tiwari et al. applied that method to quantify massive groundwater depletion in northern India caused by over reliance on aquifers for irration, and Famiglietti et al. found a similar situation in California's Central Valley. Zaitchik et al. and Lo et al. described approaches to use GRACE to constrain hydrological models, enabling integration of GRACE data with other observations and achieving much higher spatial and temporal resolutions than GRACE alone. Such approaches are now supporting applications including drought and water resources monitoring. Oceanography has likewise benefitted from the independent nature of GRACE observations. One application is measurement of the mass component of sea level rise, which complements radar altimetry and in situ measurements. GRACE also measures ocean bottom pressures (OBP), which help to refine understanding and modeling of ocean circulation and the ocean's fresh water budget, among other things. For example, Hayakawa et al. showed that GRACE observes OBP patterns absent from the background models of oceanic variability. Morison et al. used GRACE to describe important decadal scale shifts in circulation and an ongoing trend of freshening of the western Arctic, important indicators of climate variability. The research of Song and Zlotnicki and Chambers and Willis on GRACE-derived ocean bottom pressures in the sub-polar gyre led to the discovery of an ENSO teleconnection and a long-term change in OBP in the North Pacific sub-polar gyre that was not predicted by an ocean model. Further, Chambers and Willis were able to identify an internal redistribution of mass between Atlantic and Pacific Oceans lasting at least six years, which was not predicted by ocean models and was the first direct evidence of sustained mass transport from one ocean basin to another on periods longer than a year. Boening et al. observed a record increase in OBP over part of the southeastern Pacific in late 2009 and early 2010, primarily caused by wind stress curl associated with a strong and persistent anticyclone and likely related to the concurrent Central Pacific El Nino. GRACE has far surpassed its 5-year design lifetime, but it will likely succumb to the aging of batteries and instrument systems sometime in the next few years. NASA has begun initial development of a follow-on to GRACE with very similar design, which could launch as soon as 2016 and would provide continuity in the data record while improving resolution slightly. Higher resolution time variable gravity missions are also on the drawing board .

Rodell, Matthew↗

Multiscale Characterization of Additive Manufacturing Components with Computed Tomography, 3D X-ray Microscopy, and Deep Learning

Additive manufacturing (AM) facilitates the creation of complex-geometry parts, driving advancements in lightweight aerospace components, high-efficiency engine cooling channels, and customized medical implants. However, ensuring the quality and reliability of AM parts remains challenging due to internal defects, surface irregularities, porosity, and residual trapped powder, which are often inaccessible to traditional inspection methods. Recent developments in X-ray computed tomography (XCT) and 3D X-ray microscopy (XRM), particularly systems equipped with resolution-at-a-distance (RaaD™) capabilities, enable high-resolution, non-destructive evaluation of AM components across multiple scales, from sub-micrometer to macroscopic levels. This paper explores modern XCT and XRM techniques for multiscale characterization of AM parts, focusing on their ability to detect and analyze defects such as porosity, cracks, inclusions, and surface roughness, while offering insights into defect formation mechanisms, material properties, and process-induced variations. The integration of deep learning (DL) frameworks, including Simurgh, DeepRecon, and DeepScout, enhances XCT/XRM workflows by reducing scan times, improving resolution recovery, and enabling accurate defect detection even with limited projection data. These DL-based methods overcome limitations of traditional reconstruction techniques, enabling faster, more reliable characterization of dense materials like Inconel 718 and novel alloys such as AlCe. Applications include process parameter optimization, high-throughput quality control, and multistage AM process evaluation, with DL-enhanced workflows accelerating analysis times from weeks to days. Correlative imaging approaches further validate XCT and XRM data against scanning electron microscopy (SEM) images of physically sectioned samples, confirming the accuracy of DL-based reconstructions and enabling comprehensive defect analysis. While challenges remain in generalizing DL models to diverse materials and imaging conditions, improvements in resolution, noise reduction, and defect detection highlight the transformative potential of these methods. This multiscale and correlative approach enables precise identification and correlation of microstructural features with the overall performance of AM components. By integrating advanced XCT, XRM, and DL techniques, this paper demonstrates a significant leap forward in AM characterization, offering valuable insights into the relationships between processing parameters, microstructure, and part performance, and driving innovations that enhance the quality and reliability of AM products for demanding industrial applications.

Additive manufacturing↗

Analysis of Satellite-Derived Arctic Tropospheric BrO Columns in Conjunction with Aircraft Measurements During ARCTAS and ARCPAC

We derive tropospheric column BrO during the ARCTAS and ARCPAC field campaigns in spring 2008 using retrievals of total column BrO from the satellite UV nadir sensors OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a photochemical simulation. We conduct a comprehensive comparison of satellite-derived tropospheric BrO column to aircraft in-situ observations ofBrO and related species. The aircraft profiles reveal that tropospheric BrO, when present during April 2008, was distributed over a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). Perturbations to the total column resulting from tropospheric BrO are the same magnitude as perturbations due to longitudinal variations in the stratospheric component, so proper accounting of the stratospheric signal is essential for accurate determination of satellite-derived tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radiances were obtained over bright surfaces (albedo> 0.7), for solar zenith angle < 80 and clear sky conditions. The rapid activation of BrO due to surface processes (the bromine explosion) is apparent in both the OMI and GOME-2 based tropospheric columns. The wide orbital swath of OMI allows examination of the evolution of tropospheric BrO on about hourly time intervals near the pole. Low surface pressure, strong wind, and high PBL height are associated with an observed BrO activation event, supporting the notion of bromine activation by high winds over snow.

Choi, S.↗

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗

Hydrological projections over the Upper Indus Basin at 1.5 °C and 2.0 °C temperature increase

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.

54 ENVIRONMENTAL SCIENCES↗

Modelling Temporal Variability in the Carbon Balance of a Spruce/Moss Boreal Forest

A model of the daily carbon balance of a black spruce/feathermoss boreal forest ecosystem was developed and results compared to preliminary data from the 1994 BOREAS field campaign in northern Manitoba, Canada. The model, driven by daily weather conditions, simulated daily soil climate status (temperature and moisture profiles), spruce photosynthesis and respiration, moss photosynthesis and respiration, and litter decomposition. Model agreement with preliminary field data was good for net ecosystem exchange (NEE), capturing both the asymmetrical seasonality and short-term variability. During the growing season simulated daily NEE ranged from -4 g C m(exp -2) d(exp -1) (carbon uptake by ecosystem) to + 2 g C m(exp -2) d(exp -1) (carbon flux to atmosphere), with fluctuations from day to day. In the early winter simulated NEE values were + 0.5 g C m(exp -2) d(exp -1), dropping to + 0.2 g C m(exp -2) d(exp -1) in mid-winter. Simulated soil respiration during the growing season (+ 1 to + 5 g C m(exp -2) d(exp -1)) was dominated by metabolic respiration of the live moss, with litter decomposition usually contributing less than 30% and live spruce root respiration less than 10% of the total. Both spruce and moss net primary productivity (NPP) rates were higher in early summer than late summer. Simulated annual NEE for 1994 was -51 g C m(exp -2) y(exp -1), with 83% going into tree growth and 17% into the soil carbon accumulation. Moss NPP (58 g C m(exp -2) d(exp -1)) was considered to be litter (i.e. soil carbon input; no net increase in live moss biomass). Ecosystem respiration during the snow-covered season (84 g Cm(exp -2)) was 58% of the growing season net carbon uptake. A simulation of the same site for 1968-1989 showed about 10-20% year-to-year variability in heterotrophic respiration (mean of + 113 g C m-2 y@1). Moss NPP ranged from 19 to 114 g C m(exp -2) y(exp -1); spruce NPP from 81 to 150 g C nt-2 y,@l; spruce growth (NPP minus litterfall) from 34 to 103 g C m(exp -2) y(exp -1); NEE ranged from +37 to -142 g C m(exp -2) y(exp -1). Values for these carbon balance terms in 1994 were slightly smaller than the 1969 - 89 means. Higher ecosystem productivity years (more negative NEE) generally had early springs and relatively wet summers; lower productivity years had late springs and relatively dry summers.

Frolking, S.↗

Mesoscale simulations of the November 25-26 and December 5-6 cirrus cases using the RAMS model

The Regional Atmospheric Modeling System (RAMS), developed at Colorado State University, was used during the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) 2 (13 Nov. through 6 Dec. 1991) to provide real time forecasts of cirrus clouds. Forecasts were run once a day, initializing with the 0000 UTC dataset provided by NOAA (Forecast Systems Laboratory (FSL) Mesoscale Analysis and Prediction System (MAPS)). In order to obtain better agreement with observations, a second set of simulations were done for the FIRE 2 cases that occurred on 25-26 Nov. and 5-6 Dec. In this set of simulations, a more complex radiation scheme was used, the Chen/Cotton radiation scheme, along with the nucleation of ice occurring at ice supersaturations as opposed to nucleation occurring at water supersaturations that was done in the actual forecast version. The runs using these more complex schemes took longer wall clock time (7-9 hours for the actual forecasts as compared to 12-14 hrs for the runs using the more complex schemes) however, the final results of the simulations were definitely improved upon. Comparisons between these two sets of simulations are given. Now underway are simulations of these cases using a closed analytical solution for the auto-conversion of ice from a pristine ice class (sizes less than about 50 microns in effective diameter) to a snow class (effective diameters on the order of several hundred microns). This solution is employed along with a new scheme for the nucleation of ice crystals due to Meyers et al and Demott et al. The scheme is derived assuming complete gamma distributions for both the pristine and snow classes. The time rate of change of the number concentration and mass mixing-ratio of each distribution is found by calculating either the flux of crystals that grow beyond a certain critical diameter by vapor deposition in an ice supersaturated regime or by calculating the flux of crystals that evaporate to sizes below that same critical effective diameter.

Harrington, J. L.↗

Introduction to this Special Issue on Geostatistics and Scaling of Remote Sensing

The germination of this special PE&RS issue began at the Royal Geographical Society (with the Institute of British Geographers)(RCS-IBC) annual meeting in January, 1997 held at the University of Exeter in Exeter, England. The cold and snow of an England winter were greatly tempered by the friendly and cordial discussions that ensued at the meeting on possible ways to foster both dialog and research across "the Big Pond" between geographers in the US and the UK on the use of geostatistics and geospatial techniques for remote sensing of land surface processes. It was decided that one way to stimulate and enhance cooperation on the application of geostatistics and geospatial methods in remote sensing was to hold parallel sessions on these topics at appropriate meeting venues in 1998 in both the US and the UK Selected papers given at these sessions would be published as a special issue of PE&RS on the US side, and as a special issue of Computers and Geosciences (C&G) on the UK side, to highlight the commonality in research on geostatistics and geospatial methods in remote sensing and spatial data analysis on both sides of the Atlantic Ocean. As a consequence, a session on "Ceostatistics and Geospatial Techniques for Remote Sensing of Land Surface Processes" was held at the Association of American Geographers (AAG) annual meeting in Boston, Massachusetts in March, 1998, sponsored by the AAG's Remote Sensing Specialty Group (RSSG). A similar session was held at the RGS-IBG annual meeting in Guildford, Surrey, England in January 1998, organized by the Modeling and Advanced Techniques Special Interest Group (MAT SIG) of the Remote Sensing Society (RSS). The six papers that in part, comprise this issue of PE&RS, are the US complement to such a dual journal publication effort. Both of us are co-editors of each of the journal special issues, with the lead editor of each journal being from their respective side of the Atlantic where the journals are published. The special issue of C&G that constitutes the other half of this co-edited journal series will be published in early 1999, with 3 papers by US authors being published along with 6 papers authored by individuals from the UK and other places in Europe.

Quattrochi, Dale A.↗

The Use of Sun Elevation Angle for Stereogrammetric Boreal Forest Height in Open Canopies

Stereogrammetry applied to globally available high resolution spaceborne imagery (HRSI; less than 5 m spatial resolution) yields fine-scaled digital surface models (DSMs) of elevation. These DSMs may represent elevations that range from the ground to the vegetation canopy surface, are produced from stereoscopic image pairs (stereo pairs) that have a variety of acquisition characteristics, and have been coupled with lidar data of forest structure and ground surface elevation to examine forest height. This work explores surface elevations from HRSI DSMs derived from two types of acquisitions in open canopy forests. We (1) apply an automated mass-production stereogrammetry workflow to along-track HRSI stereo pairs, (2) identify multiple spatially coincident DSMs whose stereo pairs were acquired under different solar geometry, (3) vertically co-register these DSMs using coincident spaceborne lidar footprints (from ICESat-GLAS) as reference, and(4) examine differences in surface elevations between the reference lidar and the co-registered HRSI DSMs associated with two general types of acquisitions (DSM types) from different sun elevation angles. We find that these DSM types, distinguished by sun elevation angle at the time of stereo pair acquisition, are associated with different surface elevations estimated from automated stereogrammetry in open canopy forests. For DSM values with corresponding reference ground surface elevation from spaceborne lidar footprints in open canopy northern Siberian Larix forests with slopes less than10, our results show that HRSI DSM acquired with sun elevation angles greater than 35deg and less than 25deg (during snow-free conditions) produced characteristic and consistently distinct distributions of elevation differences from reference lidar. The former include DSMs of near-ground surfaces with root mean square errors less than 0.68 m relative to lidar. The latter, particularly those with angles less than 10deg, show distributions with larger differences from lidar that are associated with open canopy forests whose vegetation surface elevations are captured. Terrain aspect did not have a strong effect on the distribution of vegetation surfaces. Using the two DSM types together, the distribution of DSM-differenced heights in forests (6.0 m, sigma = 1.4 m) was consistent with the distribution of plot-level mean tree heights (6.5m, sigma = 1.2 m). We conclude that the variation in sun elevation angle at time of stereo pair acquisition can create illumination conditions conducive for capturing elevations of surfaces either near the ground or associated with vegetation canopy. Knowledge of HRSI acquisition solar geometry and snow cover can be used to understand and combine stereogrammetric surface elevation estimates to co-register rand difference overlapping DSMs, providing a means to map forest height at fine scales, resolving the vertical structure of groups of trees from spaceborne platforms in open canopy forests.

Paul M Montesano↗

Relative Influence of Initial Surface and Atmospheric Conditions on Seasonal Water and Energy Balances

We constructed and analyzed wet and dry soil moisture composites for the mid-latitude GCIP region of the central US using long climate model simulations made with the NCAR CCM3 and reanalysis products from NCEP. Using the diagnostic composites as a guide, we have completed a series of predictability experiments in which we imposed soil water initial conditions in CCM3 for the GCIP region for June 1 from anomalously wet and dry years, with atmospheric initial conditions taken from June 1 of a year with 'near-normal' soil water, and initial soil water from the near-normal year and atmospheric initial conditions from the wet and dry years. Preliminary results indicate that the initial state of the atmosphere is more important than the initial state of soil water determining the subsequent late spring and summer evolution of sod water over the GCIP region. Surprisingly, neither the composites or the predictability experiments yielded a strong influence of soil moisture on the atmosphere. To explore this further, we have made runs with extreme dry soil moisture initial anomalies imposed over the GCIP region (the soil close to being completely dry). These runs did yield a very strong effect on the atmosphere that persisted for at least three months. We conclude that the magnitude of the initial soil moisture anomaly is crucial, at least in CCM3, and are currently investigating whether a threshold exists, below which little impact is seen. In a complementary study, we compared the impact of the initial condition of snow cover versus the initial atmospheric state over the western US (corresponding to the westward extension of the GAPP program follow-on to GCIP). In this case, the initial prescription of snow cover is far more important than the initial atmospheric state in determining the subsequent evolution of snow cover. We are currently working to understand the very different soil water and snow cover results.

Oglesby, Robert J.↗

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗