Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance estimation”

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,207 records · Page 67

Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements

In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.

Lei, Jinlong↗

Little vertical and circumferential variations in stem xylem water Δ 2 H and Δ 18 O in three tree species

Vertical and circumferential patterns of xylem water isotope signatures in large trees are rarely characterized but may influence interpretation of water source investigations and soil/xylem water isotope offsets. Furthermore, to examine vertical and circumferential variation in xylem water δ 2 H and δ 18 O, we collected xylem tissue at 5 heights (0, 1, 4, 7 and 10 m) from three replicate trees of three species with contrasting xylem anatomy, two angiosperms (Liquidambar styraciflua and Quercus nigra) and one gymnosperm (Pinus taeda). Concurrently, we also determined δ 2 H and δ 18 O of groundwater, soil water, and recent precipitation. On a different day, we collected circumferential samples at 1 m from the same trees to test for sectorality effects. Water from stem and soil samples were cryogenically extracted and analyzed for δ 2 H and δ 18 O. Mean xylem water δ 2 H and δ 18 O were significantly different between species, both vertically and circumferentially. We did not find significant (p = 0.05) systematic variation in δ 2 H or δ 18 O with height. We found no significant evidence for sectorality effects on δ 2 H, δ 18 O. Variances of vertical synoptic xylem water δ 2 H or δ 18 O were similar between species, 7.2–10.4‰ for δ 2 H and 0.58–0.81‰ for δ 18 O. Circumferential variances were also similar between species, 4.0–6.0 for δ 2 H and 0.37–0.44 for δ 18 O. A mixing model showed that sweetgum, water oak and loblolly pine, were drawing most of their water from deep soil from 45 to 190 cm (84.7, 68.4 and 53.2%, respectively) however, soil water-excess values indicate δ 2 H fractionation effects on these estimates. Dual isotope mixing model evaluation with single and multiple sample configurations showed that source water estimates were not affected by within tree variability of xylem water signatures. Xylem water δ 2 H and δ 18 O variability with height or circumference, was 3.2 and 2.7 times less than between tree variation because of transient temporal and spatial processes and is, therefore, not likely to affect interpretations of water sourcing in these three species.

59 BASIC BIOLOGICAL SCIENCES↗

Exact Bayesian computation using H-functions

In the Bayesian Statistics, the use of conjugate prior distributions allows to obtain exact posterior quantities, avoiding the use of approximate methods. However, there are only a few conjugate structures available, which limits the use in practical problems. In this work we propose a general procedure to obtain the posterior distribution in an exact form in non-conjugate location parameter models with known (possibly, heteroskedastic) variance. The theory is based on special functions, specifically the H-function, which embraces most of the probability distributions. We express explicitly in a computable form all the Bayesian tools for posterior inferences, that is the posterior and the posterior predictive distributions and their moments, as well as the cumulative posterior distribution. As an illustration we consider a problem in astronomy of estimating the distance to the Large Magellanic Cloud, the largest galaxy now orbiting the Milky Way Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Online eco-routing for electric vehicles using combinatorial multi-armed bandit with estimated covariance

Identifying energy-efficient routes in real-time has significant implications for the energy-optimal operations of electric vehicles (EVs). Here, this study proposes a novel model for EV online eco-routing problem, which obtains the minimal expected energy consumption paths (MECPs) for multiple origin-destination (OD) pairs simultaneously. Specifically, we formulate the routing problem as a bandit problem and solve it with online algorithms. We extend the algorithms by implementing a path elimination mechanism to reduce the candidate path set and introducing the variance and covariance of the energy consumption to reduce the uncertainties. The numerical results show that the proposed algorithms can efficiently obtain near-optimal MECPs, and the solution is significantly better than the widely used shortest trip time path algorithm (STTP) and shortest trip distance path algorithm (SDP). The variation considering link energy covariance and path elimination generates paths that save 4.1% of energy compared to the SDP and 5.4% to the STTP.

33 ADVANCED PROPULSION SYSTEMS↗

BAYESIAN CONFIDENCE LIMITS OF ELECTRON SPECTRA OBTAINED THROUGH REGULARIZED INVERSION OF SOLAR HARD X-RAY SPECTRA

Many astrophysical observations are characterized by a single, non-repeatable measurement of a source brightness or intensity, from which we are to construct estimates for the true intensity and its uncertainty. For example, the hard X-ray count spectrum from transient events such as solar flares can be observed only once, and from this single spectrum one must determine the best estimate of the underlying source spectrum I({epsilon}), and hence the form of the responsible electron spectrum F(E). Including statistical uncertainties on the measured count spectrum yields a 'confidence strip' that delineates the boundaries of electron spectra that are consistent with the observed photon spectrum. In this short article, we point out that the expectation values of the source brightness and its variance in a given photon energy bin are in general not (as has been assumed in prior works) equal to n, the number of counts observed in that energy bin. Rather, they depend both on n and on prior knowledge of the overall photon spectrum. Using Bayesian statistics, we provide an explicit procedure and formulas for determining the 'confidence strip' (Bayesian credible region) for F(E), thus providing rigorous bounds on the intensity and shape of the accelerated electron spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

A study on the clustering properties of radio-selected sources in the Lockman Hole region at 325 MHz

Studying the spatial distribution of extragalactic source populations is vital in understanding the matter distribution in the Universe. It also enables understanding the cosmological evolution of dark matter density fields and the relationship between dark matter and luminous matter. Clustering studies are also required for EoR foreground studies since it affects the relevant angular scales. This paper investigates the angular and spatial clustering properties and the bias parameter of radio-selected sources in the Lockman Hole field at 325 MHz. The data probes sources with fluxes ≳0.3 mJy within a radius of 1.8° around the phase centre of a 6° × 6° mosaic. Based on their radio luminosity, the sources are classified into Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs). Clustering and bias parameters are determined for the combined populations and the classified sources. The spatial correlation length and the bias of AGNs are greater than SFGs- indicating that more massive haloes host the former. This study is the first reported estimate of the clustering property of sources at 325 MHz, intermediate between the pre-existing studies at high and low-frequency bands. It also probes a well-studied deep field at an unexplored frequency with moderate depth and area. Clustering studies require such observations along different lines of sight, with various fields and data sets across frequencies to avoid cosmic variance and systematics. Thus, an extragalactic deep field has been studied in this work to contribute to this knowledge.

79 ASTRONOMY AND ASTROPHYSICS↗

Patterns of patient specific dosimetry in total body irradiation

Purpose: Total body irradiation (TBI) has been used for bone marrow transplant for hematologic and immune deficiency conditions. The goal of TBI is to deliver a homogeneous dose to the entire body, with a generally accepted range of dose uniformity being within {+-}10% of the prescribed dose. The moving table technique for TBI could make dose uniform in whole body by adjusting couch speed. However, it is difficult to accurately estimate the actual dose by calculation and hence in vivo dosimetry (IVD) is routinely performed. Here, the authors present patterns of patient-specific IVD in 161 TBI patients treated at our institution. Methods: Cobalt-60 teletherapy unit (Model C9 Cobalt-60 teletherapy unit, Picker X-ray Corporation) with customized moving bed (SITI Industrial Products, Inc., Fishers, IN) were used for TBI treatment. During treatment, OneDose{sup TM} (Sicel Technology, NC) Metal Oxide-silicon Semiconductor Field Effect Transistor detectors were placed at patient body surface; both entrance and exit side of the beam at patient head, neck, mediastinum, umbilicus, and knee to estimate midplane dose. When large differences (>10%) between the prescribed and measured dose were observed, dose delivery was corrected for subsequent fractions by the adjustment of couch speed and/or bolus placement. Under IRB exempt status, the authors retrospectively analyzed the treatment records of 161 patients who received TBI treatment between 2006 and 2011. Results: Across the entire cohort, the median {+-} SD (range) percent variance between calculated and measured dose for head, neck, mediastinum, umbilicus, and knee was -2.3 {+-} 10.2% (-66.2 to +35.3), 1.1 {+-} 11.5% (-62.2 to +40.3), -1.9 {+-} 9.5% (-66.4 to +46.6), -1.1 {+-} 7.2% (-35.2 to +42.9), and 3.4 {+-} 12.2% (-47.9 to +108.5), respectively. More than half of treatments were within {+-}10% of the prescribed dose for all anatomical regions. For 80% of treatments (10%-90%), dose at the umbilicus was within {+-}10%. However, some large differences greater than 35% were also found at several points. For one case, the knee received double the prescribed dose. When the dose differences for multiple fractions were averaged, compliance ({+-}10%) between the prescription and measured dose was improved compared to the dose difference of the first single fraction, for example, as at umbilicus, which improved from 83.9% to 98.5%. Conclusions: Actual dose measurement analysis of TBI patients revealed a potentially wide variance from the calculated dose. Based from their IVD method for TBI using Cobalt-60 irradiator and moving table, {+-}10% over entire body is hard to achieve. However, it can be significantly improved with immediate feedback after the first fraction prior to subsequent treatments.

60 APPLIED LIFE SCIENCES↗

Radar and Lidar Radar DEM

Using radar and lidar data, the aim is to improve 3D rendering of terrain, including digital elevation models (DEM) and estimates of vegetation height and biomass in a variety of forest types and terrains. The 3D mapping of vegetation structure and the analysis are useful to determine the role of forest in climate change (carbon cycle), in providing habitat and as a provider of socio-economic services. This in turn will lead to potential for development of more effective land-use management. The first part of the project was to characterize the Shuttle Radar Topography Mission DEM error with respect to ICESat/GLAS point estimates of elevation. We investigated potential trends with latitude, canopy height, signal to noise ratio (SNR), number of LiDAR waveform peaks, and maximum peak width. Scatter plots were produced for each variable and were fitted with 1st and 2nd degree polynomials. Higher order trends were visually inspected through filtering with a mean and median filter. We also assessed trends in the DEM error variance. Finally, a map showing how DEM error was geographically distributed globally was created.

remote sensing↗

Uncertainty in land carbon budget simulated by terrestrial biosphere models: the role of atmospheric forcing

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Magnetic field directional discontinuities. 1: Minimum variance errors

Errors associated with the minimum variance analysis of directional discontinuity normal components were investigated using both computer simulation of discontinuities with controlled properties and the examination of current sheets observed by the Mariner 10 spacecraft. The simulated discontinuities were created by adding fluctuations, represented by isotropic noise, to exactly known but varying (in a plane) magnetic field components. An empirical expression for the magnitude of the error in an estimated discontinuity normal component, relative to the total field across the discontinuity, was derived, as well as other relevant statistical properties. Use of the empirical relation in the analysis of 644 discontinuities observed by Mariner 10 provides a more precise, but probably conservative, estimate of an upper bound on the relative normal component value for tangential discontinuities that can be used to separate rotational from tangential discontinuities in studies using only magnetic field data from a single spacecraft, at least for the interplanetary region of space considered.

Lepping, R. P.↗

Analysis of satellite-derived ozone and water vapor measurements

Research under this grant has involved analyses of satellite-derived measurements of water vapor and total ozone. Upper tropospheric water vapor data from the Nimbus 7 temperature-humidity infrared radiometer (THIR) 6.7 micron channel were analyzed by Fourier transformation to provide wavenumber spectra for a case of an upper level system over Europe. The power law spectrum for horizontal scales from 60 to several hundred km suggests that when convective cloud energy sources are not present, the enstrophy-cascading process (with variance proportional to the minus three power of wavenumber) may hold to much smaller scales than previously thought. Several investigations of the middle atmosphere (stratosphere and mesosphere) were also conducted. Slow atmospheric oscillations with periods of 1-2 months were investigated in total ozone measurements and in ionospheric data. Using one-point correlation maps and total ozone mapping spectrometer (TOMS) data, strong statistical evidence was found for the planetary-scale oscillations at these periods. A dipole-like pattern in the tropical equatorial Indian ocean-western Pacific region was documented, along with wavetrain-like patterns emanating into the extratropics. While attempting to follow the oscillations upward above stratosphere, oscillations were studied in ionospheric data (D-region radio wave absorption, in the upper mesosphere, around 80-90 km altitude). Cross-spectral analyses with solar flux data revealed that much of the observed 1-2 month variance in the ionospheric D-region was well correlated to solar variations. Further research has involved investigations of the photochemical lifetime of greenhouse gases methane and nitrous oxide. Using Nimbus 7 satellite data, a new analysis technique was employed to determine the lifetime of these gases in situ in the upper stratosphere. The lifetimes are found to be in good agreement with theoretical estimates.

Stanford, John L.↗

Causal correlation of foliar biochemical concentrations with AVIRIS spectra using forced entry linear regression

A major goal of airborne imaging spectrometry is to estimate the biochemical composition of vegetation canopies from reflectance spectra. Remotely-sensed estimates of foliar biochemical concentrations of forests would provide valuable indicators of ecosystem function at regional and eventually global scales. Empirical research has shown a relationship exists between the amount of radiation reflected from absorption features and the concentration of given biochemicals in leaves and canopies (Matson et al., 1994, Johnson et al., 1994). A technique commonly used to determine which wavelengths have the strongest correlation with the biochemical of interest is unguided (stepwise) multiple regression. Wavelengths are entered into a multivariate regression equation, in their order of importance, each contributing to the reduction of the variance in the measured biochemical concentration. A significant problem with the use of stepwise regression for determining the correlation between biochemical concentration and spectra is that of 'overfitting' as there are significantly more wavebands than biochemical measurements. This could result in the selection of wavebands which may be more accurately attributable to noise or canopy effects. In addition, there is a real problem of collinearity in that the individual biochemical concentrations may covary. A strong correlation between the reflectance at a given wavelength and the concentration of a biochemical of interest, therefore, may be due to the effect of another biochemical which is closely related. Furthermore, it is not always possible to account for potentially suitable waveband omissions in the stepwise selection procedure. This concern about the suitability of stepwise regression has been identified and acknowledged in a number of recent studies (Wessman et al., 1988, Curran, 1989, Curran et al., 1992, Peterson and Hubbard, 1992, Martine and Aber, 1994, Kupiec, 1994). These studies have pointed to the lack of a physical link between wavelengths chosen by stepwise regression and the biochemical of interest, and this in turn has cast doubts on the use of imaging spectrometry for the estimation of foliar biochemical concentrations at sites distant from the training sites. To investigate this problem, an analysis was conducted on the variation in canopy biochemical concentrations and reflectance spectra using forced entry linear regression.

Dawson, Terence P.↗

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗

NOISE IN THE CROSS-POWER SPECTRUM OF THE VELA PULSAR

We compare the noise in interferometric measurements of the Vela pulsar from ground- and space-based antennas with theoretical predictions. The noise depends on both the flux density and the interferometric phase of the source. Because the Vela pulsar is bright and scintillating, these comparisons extend into both the low and high signal-to-noise regimes. Furthermore, our diversity of baselines explores the full range of variation in interferometric phase. We find excellent agreement between theoretical expectations and our estimates of noise among samples within the characteristic scintillation scales. Namely, the noise is drawn from an elliptical Gaussian distribution in the complex plane, centered on the signal. The major axis, aligned with the signal phase, varies quadratically with the signal, while the minor axis, at quadrature, varies with the same linear coefficients. For weak signal, the noise approaches a circular Gaussian distribution. Both the variance and covariance of the noise are also affected by artifacts of digitization and correlation. In particular, we show that gating introduces correlations between nearby spectral channels.

79 ASTRONOMY AND ASTROPHYSICS↗

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor)

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ~5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ~33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.

59 BASIC BIOLOGICAL SCIENCES↗

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Soil respiration across a permafrost transition zone: spatial structure and environmental correlates

Soil respiration is a key ecosystem function whereby shifts in respiration rates can shift systems from carbon sinks to sources. Soil respiration in permafrost-associated systems is particularly important given climate change driven permafrost thaw that leads to significant uncertainty in resulting ecosystem carbon dynamics. Here we characterize the spatial structure and environmental drivers of soil respiration across a permafrost transition zone. We find that soil respiration is characterized by a non-linear threshold that occurs at active-layer depths greater than 140 cm. We also find that within each season, tree basal area is a dominant driver of soil respiration regardless of spatial scale, but only in spatial domains with significant spatial variability in basal area. Our analyses further show that spatial variation (the coefficient of variation) and mean-variance power-law scaling of soil respiration in our boreal system are consistent with previous work in other ecosystems (e.g., tropical forests) and in population ecology, respectively. Comparing our results to those in other ecosystems suggests that temporally stable features such as tree-stand structure are often primary drivers of spatial variation in soil respiration. If so, this provides an opportunity to better estimate the magnitude and spatial variation in soil respiration through remote sensing. Finally, combining such an approach with broader knowledge of thresholding behavior – here related to active layer depth – would provide empirical constraints on models aimed at predicting ecosystem responses to ongoing permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗