Broad scale remote sensing promises to build forest inventories at unprecedented scales. A crucial step in this process is to associate sensor data into individual crowns. While dozens of crown ...detection algorithms have been proposed, their performance is typically not compared based on standard data or evaluation metrics. There is a need for a benchmark dataset to minimize differences in reported results as well as support evaluation of algorithms across a broad range of forest types. Combining RGB, LiDAR and hyperspectral sensor data from the USA National Ecological Observatory Network’s Airborne Observation Platform with multiple types of evaluation data, we created a benchmark dataset to assess crown detection and delineation methods for canopy trees covering dominant forest types in the United States. This benchmark dataset includes an R package to standardize evaluation metrics and simplify comparisons between methods. The benchmark dataset contains over 6,000 image-annotated crowns, 400 field-annotated crowns, and 3,000 canopy stem points from a wide range of forest types. In addition, we include over 10,000 training crowns for optional use. We discuss the different evaluation data sources and assess the accuracy of the image-annotated crowns by comparing annotations among multiple annotators as well as overlapping field-annotated crowns. We provide an example submission and score for an open-source algorithm that can serve as a baseline for future methods.
We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk ...Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of mG, where mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of J m−3 in OMC 1, comparable both to the gravitational potential energy density of OMC 1 (∼10−7 J m−3) and to the energy density in the Orion BN/KL outflow (∼10−7 J m−3). We find that neither the Alfvén velocity in OMC 1 nor the velocity of the super-Alfvénic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the ∼500 yr lifetime of the outflow. Hence, we propose that the hourglass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.
Functional ecology has increasingly focused on describing ecological communities based on their traits (measurable features affecting individuals’ fitness and performance). Analyzing trait ...distributions within and among forests could significantly improve understanding of community composition and ecosystem function. Historically, data on trait distributions are generated by (1) collecting a small number of leaves from a small number of trees, which suffers from limited sampling but produces information at the fundamental ecological unit (the individual), or (2) using remote-sensing images to infer traits, producing information continuously across large regions, but as plots (containing multiple trees of different species) or pixels, not individuals. Remote-sensing methods that identify individual trees and estimate their traits would provide the benefits of both approaches, producing continuous large-scale data linked to biological individuals. We used data from the National Ecological Observatory Network (NEON) to develop a method to scale up functional traits from 160 trees to the millions of trees within the spatial extent of two NEON sites. The pipeline consists of three stages: (1) image segmentation, to identify individual trees and estimate structural traits; (2) an ensemble of models to infer leaf mass area (LMA), nitrogen, carbon, and phosphorus content using hyperspectral signatures, and DBH from allometry; and (3) predictions for segmented crowns for the full remote-sensing footprint at the NEON sites. The R² values on held-out test data ranged from 0.41 to 0.75 on held-out test data. The ensemble approach performed better than single partial least-squares models. Carbon performed poorly compared to other traits (R² of 0.41). The crown segmentation step contributed the most uncertainty in the pipeline, due to over-segmentation. The pipeline produced good estimates of DBH (R² of 0.62 on held-out data). Trait predictions for crowns performed significantly better than comparable predictions on pixels, resulting in improvement of R² on test data of between 0.07 and 0.26. We used the pipeline to produce individual-level trait data for ~5 million individual crowns, covering a total extent of ~360 km². This large data set allows testing ecological questions on landscape scales, revealing that foliar traits are correlated with structural traits and environmental conditions.
Forests provide biodiversity, ecosystem, and economic services. Information on individual trees is important for understanding forest ecosystems but obtaining individual-level data at broad scales is ...challenging due to the costs and logistics of data collection. While advances in remote sensing techniques allow surveys of individual trees at unprecedented extents, there remain technical challenges in turning sensor data into tangible information. Using deep learning methods, we produced an open-source data set of individual-level crown estimates for 100 million trees at 37 sites across the United States surveyed by the National Ecological Observatory Network's Airborne Observation Platform. Each canopy tree crown is represented by a rectangular bounding box and includes information on the height, crown area, and spatial location of the tree. These data have the potential to drive significant expansion of individual-level research on trees by facilitating both regional analyses and cross-region comparisons encompassing forest types from most of the United States.
Agricultural land now exceeds forests as the dominant global biome. Because of their global dominance, and potential expansion or loss, methods to estimate biomass and carbon in agricultural areas ...are necessary for monitoring global terrestrial carbon stocks and predicting carbon dynamics. Agricultural areas in the tropics have substantial tree cover and associated above ground biomass (AGB) and carbon. Active remote sensing data, such as airborne LiDAR (light detection and ranging), can provide accurate estimates of biomass stocks, but common plot-based methods may not be suitable for agricultural areas with dispersed and heterogeneous tree cover. The objectives of this research are to quantify AGB of a tropical agricultural landscape using a tree-based method that directly incorporates the size of individual trees, and to understand how landscape estimates of AGB from a tree-based method compare to estimates from a plot-based method. We use high-resolution (1.12 m) airborne LiDAR data collected on a 9280-ha region of the Azuero Peninsula of Panama. We model individual tree AGB with canopy dimensions from the LiDAR data. We apply the model to individual tree crown polygons and aggregate AGB estimates to compare with previously developed plot-based estimates. We find that agricultural trees are a distinct and dominant part of our study site. The tree-based approach estimates greater AGB in pixels with low forest cover than the plot-based approach, resulting a 2-fold difference in landscape AGB estimates between the methods for non-forested areas. Additionally, one third of the total landscape AGB exists in areas having <10% cover, based on a global tree cover product. Our study supports the continued use and development of allometric models to predict individual tree biomass from LiDAR-derived canopy dimensions and demonstrates the potential for spatial information from high-resolution data, such as relative isolation of canopies, to improve allometric models.
Display omitted
•Substantial tree biomass exists in areas with zero to little forest cover.•Landscape estimates from a tree-based method capture biomass of trees outside forests.•The plot-based method underestimates biomass in low forest areas relative to tree-based estimates.•The degree of tree canopy isolation influences the allometry of tree size with biomass.
The ecology of forest ecosystems depends on the composition of trees. Capturing fine-grained information on individual trees at broad scales provides a unique perspective on forest ecosystems, forest ...restoration, and responses to disturbance. Individual tree data at wide extents promises to increase the scale of forest analysis, biogeographic research, and ecosystem monitoring without losing details on individual species composition and abundance. Computer vision using deep neural networks can convert raw sensor data into predictions of individual canopy tree species through labeled data collected by field researchers. Using over 40,000 individual tree stems as training data, we create landscape-level species predictions for over 100 million individual trees across 24 sites in the National Ecological Observatory Network (NEON). Using hierarchical multi-temporal models fine-tuned for each geographic area, we produce open-source data available as 1 km 2 shapefiles with individual tree species prediction, as well as crown location, crown area, and height of 81 canopy tree species. Site-specific models had an average performance of 79% accuracy covering an average of 6 species per site, ranging from 3 to 15 species per site. All predictions are openly archived and have been uploaded to Google Earth Engine to benefit the ecology community and overlay with other remote sensing assets. We outline the potential utility and limitations of these data in ecology and computer vision research, as well as strategies for improving predictions using targeted data sampling.
• Premise of the study: In ecosystems maintained by low-intensity surface fires, tree bark thickness is a determinant of fire-survival because it protects underlying tissues from heat damage. ...However, it has been unclear whether relatively thick bark is maintained at all heights or only near the ground where damage is most likely.• Methods: We studied six Quercus species from the red and white clades, with three species characteristic of fire-maintained savannas and three species characteristic of forests with infrequent fire. Inner and outer bark (secondary phloem and rhytidome, respectively) thicknesses were measured at intervals from 10 to 300 cm above the ground. We used linear mixed-effects models to test for relationships among height, habitat, and clade on relative thickness (stem proportion) of total, inner, and outer bark. Bark moisture and tissue density were measured for each species at 10 cm.• Key results: Absolute and relative total bark thickness declined with height, with no difference in height-related changes between habitat groups. Relative outer bark thickness showed a height-by-habitat interaction. There was a clade effect on relative thickness, but no interaction with height. Moisture contents were higher in inner than outer bark, and red oaks had denser bark than white oaks, but neither trait differed by habitat.• Conclusions: Quercus species characteristic of fire-prone habitats invest more in outer bark near the ground where heat damage to outer tissues is most likely. Future investigations of bark should consider the height at which measurements are made and distinguish between inner and outer bark.
Mapping species through classification of imaging spectroscopy data is facilitating research to understand tree species distributions at increasingly greater spatial scales. Classification requires a ...dataset of field observations matched to the image, which will often reflect natural species distributions, resulting in an imbalanced dataset with many samples for common species and few samples for less common species. Despite the high prevalence of imbalanced datasets in multiclass species predictions, the effect on species prediction accuracy and landscape species abundance has not yet been quantified. First, we trained and assessed the accuracy of a support vector machine (SVM) model with a highly imbalanced dataset of 20 tropical species and one mixed-species class of 24 species identified in a hyperspectral image mosaic (350–2500 nm) of Panamanian farmland and secondary forest fragments. The model, with an overall accuracy of 62% ± 2.3% and F-score of 59% ± 2.7%, was applied to the full image mosaic (23,000 ha at a 2-m resolution) to produce a species prediction map, which suggested that this tropical agricultural landscape is more diverse than what has been presented in field-based studies. Second, we quantified the effect of class imbalance on model accuracy. Model assessment showed a trend where species with more samples were consistently over predicted while species with fewer samples were under predicted. Standardizing sample size reduced model accuracy, but also reduced the level of species over- and under-prediction. This study advances operational species mapping of diverse tropical landscapes by detailing the effect of imbalanced data on classification accuracy and providing estimates of tree species abundance in an agricultural landscape. Species maps using data and methods presented here can be used in landscape analyses of species distributions to understand human or environmental effects, in addition to focusing conservation efforts in areas with high tree cover and diversity.
Though there has been a significant amount of work investigating the early stages of low-mass star formation in recent years, the evolution of the mass assembly rate onto the central protostar ...remains largely unconstrained. Examining in depth the variation in this rate is critical to understanding the physics of star formation. Instabilities in the outer and inner circumstellar disk can lead to episodic outbursts. Observing these brightness variations at infrared or submillimeter wavelengths constrains the current accretion models. The JCMT Transient Survey is a three-year project dedicated to studying the continuum variability of deeply embedded protostars in eight nearby star-forming regions at a one-month cadence. We use the SCUBA-2 instrument to simultaneously observe these regions at wavelengths of 450 and 850 m. In this paper, we present the data reduction techniques, image alignment procedures, and relative flux calibration methods for 850 m data. We compare the properties and locations of bright, compact emission sources fitted with Gaussians over time. Doing so, we achieve a spatial alignment of better than 1″ between the repeated observations and an uncertainty of 2%-3% in the relative peak brightness of significant, localized emission. This combination of imaging performance is unprecedented in ground-based, single-dish submillimeter observations. Finally, we identify a few sources that show possible and confirmed brightness variations. These sources will be closely monitored and presented in further detail in additional studies throughout the duration of the survey.
We present the B-fields mapped in IRDC G34.43+0.24 using 850 m polarized dust emission observed with the POL-2 instrument at the James Clerk Maxwell telescope. We examine the magnetic field ...geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14″ and SHARP at 10″) to small scales (TADPOL at 2 5 and SMA at 1 5) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350 m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470 190 G, 100 40 G, and 60 34 G in the central, northern, and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter λ values 0.8 0.4, 1.1 0.8, and 0.9 0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfvénic with Alfvénic Mach numbers of 0.34 0.13, 0.53 0.30, and 0.49 0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity.