FITTING BOLE-VOLUME EQUATIONS TO SPATIALLY CORRELATED WITHIN-TREE DATA

Equations to predict the volume of an individual tree bole between stump height and the height at which its diameter has tapered to a specified minimum are common in forestry. When fitting such a regression equation, a sample of trees which span the range of sizes needed for eventual application of the equation is selected. Bole diameter is measured at ascending heights on the bole. Each tree, therefore, contributes multiple measurements to the data fitted to the equation. In contrast to past practice, we model these data in a manner which accounts for the likely spatial correlation among measurements within a tree. The resulting mixed-effects nonlinear model is fitted by REML and also by generalized estimating equations (GEE). Results from the two approaches are nearly identical, which suggests that the computationally less demanding GEE may be acceptable as a routine alternative to a fully parameterized approach

Spatiotemporal Calibration of Atmospheric Nitrogen Dioxide Concentration Estimates From an Air Quality Model for Connecticut

A spatiotemporal calibration and resolution refinement model was fitted to calibrate nitrogen dioxide (NO$_2$) concentration estimates from the Community Multiscale Air Quality (CMAQ) model, using two sources of observed data on NO$_2$ that differed in their spatial and temporal resolutions. To refine the spatial resolution of the CMAQ model estimates, we leveraged information using additional local covariates including total traffic volume within 2 km, population density, elevation, and land use characteristics. Predictions from this model greatly improved the bias in the CMAQ estimates, as observed by the much lower mean squared error (MSE) at the NO$_2$ monitor sites. The final model was used to predict the daily concentration of ambient NO$_2$ over the entire state of Connecticut on a grid with pixels of size 300 x 300 m. A comparison of the prediction map with a similar map for the CMAQ estimates showed marked improvement in the spatial resolution. The effect of local covariates was evident in the finer spatial resolution map, where the contribution of traffic on major highways to ambient NO$_2$ concentration stands out. An animation was also provided to show the change in the concentration of ambient NO$_2$ over space and time for 1994 and 1995.Comment: 23 pages, 8 figures, supplementary materia

ESTIMATING VARIANCE FUNCTIONS FOR WEIGHTED LINEAR REGRESSION

For linear models with heterogeneous error structure, four variance function models are examined for predicting the error structure in two loblolly pine data sets and one white oak data set. An index of fit and a simulation study were used to determine which models were best. The size of coefficients for linear and higher order terms varied drastically across different data sets, thus it is not desirable to recommend a general model containing both linear and higher order terms. The unspecified exponent model σ2vi = σ2(Di2 Hi)k 1 is recommended for all data sets considered. The k1 values ranged from 1.8 to 2.1. We recommend k1 = 2.0 for simplicity

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI

Cut-off importance sampling of bole volume.

Cut-off importance sampling (CIS) is introduced as a means of sampling individual trees for the purpose of estimating bole volume. The novel feature of this variant of importance sampling is the establishment on the bole of a cut-off height, HC, above which sampling is precluded. An estimator of bole volume between predetermined heights HL and HU > HC is proposed, and its design-based bias and mean square error are derived. In an application of CIS as the second stage of a two-stage sample to estimate aggregate bole volume, the gain in precision realized from CIS more than offset its bias when compared to the precision of importance sampling when HC = HU

Use of models in large-area forest surveys: comparing model-assisted, model-based and hybrid estimation

This paper focuses on the use of models for increasing the precision of estimators in large-area forest surveys. It is motivated by the increasing availability of remotely sensed data, which facilitates the development of models predicting the variables of interest in forest surveys. We present, review and compare three different estimation frameworks where models play a core role: model-assisted, model-based, and hybrid estimation. The first two are well known, whereas the third has only recently been introduced in forest surveys. Hybrid inference mixes designbased and model-based inference, since it relies on a probability sample of auxiliary data and a model predicting the target variable from the auxiliary data..We review studies on large-area forest surveys based on model-assisted, modelbased, and hybrid estimation, and discuss advantages and disadvantages of the approaches. We conclude that no general recommendations can be made about whether model-assisted, model-based, or hybrid estimation should be preferred. The choice depends on the objective of the survey and the possibilities to acquire appropriate field and remotely sensed data. We also conclude that modelling approaches can only be successfully applied for estimating target variables such as growing stock volume or biomass, which are adequately related to commonly available remotely sensed data, and thus purely field based surveys remain important for several important forest parameters. Keywords: Design-based inference, Model-assisted estimation, Model-based inference, Hybrid inference, National forest inventory, Remote sensing, Samplin

Unconstrained three-dimensional reaching in Rhesus monkeys

To better understand normative behavior for quantitative evaluation of motor recovery after injury, we studied arm movements by non-injured Rhesus monkeys during a food-retrieval task. While seated, monkeys reached, grasped, and retrieved food items. We recorded three-dimensional kinematics and muscle activity, and used inverse dynamics to calculate joint moments due to gravity, segmental interactions, and to the muscles and tissues of the arm. Endpoint paths showed curvature in three dimensions, suggesting that maintaining straight paths was not an important constraint. Joint moments were dominated by gravity. Generalized muscle and interaction moments were less than half of the gravitational moments. The relationships between shoulder and elbow resultant moments were linear during both reach and retrieval. Although both reach and retrieval required elbow flexor moments, an elbow extensor (triceps brachii) was active during both phases. Antagonistic muscles of both the elbow and hand were co-activated during reach and retrieval. Joint behavior could be described by lumped-parameter models analogous to torsional springs at the joints. Minor alterations to joint quasi-stiffness properties, aided by interaction moments, result in reciprocal movements that evolve under the influence of gravity. The strategies identified in monkeys to reach, grasp, and retrieve items will allow the quantification of prehension during recovery after a spinal cord injury and the effectiveness of therapeutic interventions