Indigenous forests cover 23.9% of New Zealand’s land area and provide highly valued ecosystem services, including climate regulation, habitat for native biota, regulation of soil erosion and recreation. Despite their importance, information on the number of tall trees and the tree height distribution across different forest classes is scarce. We present the first region-wide spatial inventory of tall trees (>30 m) based on airborne LiDAR (Light Detection and Ranging) measurements in New Zealand—covering the Greater Wellington region. This region has 159,000 ha of indigenous forest, primarily on steep mountainous land. We implement a high-performance tree mapping algorithm that uses local maxima in a canopy height model (CHM) as initial tree locations and accurately identifies the tree top positions by combining a raster-based tree crown delineation approach with information from the digital surface and terrain models. Our algorithm includes a check and correction for over-estimated heights of trees on very steep terrain such as on cliff edges. The number of tall trees (>30 m) occurring in indigenous forest in the Wellington Region is estimated to be 286,041 (±1%) and the number of giant trees (>40 m tall) is estimated to be 7340 (±1%). Stereo-analysis of aerial photographs was used to determine the accuracy of the automated tree mapping. The giant trees are mainly in the beech-broadleaved-podocarp and broadleaved-podocarp forests, with density being 0.04 and 0.12 (trees per hectare) respectively. The inventory of tall trees in the Wellington Region established here improves the characterization of indigenous forests for management and provides a useful baseline for long-term monitoring of forest conditions. Our tree top detection scheme provides a simple and fast method to accurately map overstory trees in flat as well as mountainous areas and can be directly applied to improve existing and build new tree inventories in regions where LiDAR data is available.

We use the 2B-GEOPROF-LIDAR R04 (2BGL4) and R05 (2BGL5) products and the 2B-CLDCLASS-LIDAR R04 (2BCL4) product, all generated by combining CloudSat radar and CALIPSO lidar satellite measurements with auxiliary data, to examine the vertical distribution of cloud occurrence around the Ross Ice Shelf (RIS) and Ross Sea region. We find that the 2BGL4 product, used in previous studies in this region, displays a discontinuity at 8.2 km which is not observable in the other products. This artefact appears to correspond to a change in the horizontal and vertical resolution of the CALIPSO dataset used above this level. We then use the 2BCL4 product to examine the vertical distribution of cloud occurrence, phase, and type over the RIS and Ross Sea. In particular we examine how synoptic conditions in the region, derived using a previously developed synoptic classification, impact the cloud environment and the contrasting response in the two regions. We observe large differences between the cloud occurrence as a function of altitude for synoptic regimes relative to those for seasonal variations. A stronger variation in the occurrence of clear skies and multi-layer cloud and in all cloud type occurrences over both the Ross Sea and RIS is associated more with synoptic type than seasonal composites. In addition, anomalies from the mean joint histogram of cloud top height against thickness display significant differences over the Ross Sea and RIS sectors as a function of synoptic regime, but are near identical over these two regions when a seasonal analysis is completed. However, the frequency of particular phases of cloud, notably mixed phase and water, is much more strongly modulated by seasonal than synoptic regime compositing, which suggests that temperature is still the most important control on cloud phase in the region.

Episodic deposition of light absorbing impurities on glaciers reduces albedo and exacerbates snow melt. In 2019/2020 a devastating Australian bushfire and desert dust event combined with favorable meteorological conditions transported an unprecedented mass of impurities across the Tasman Sea turning the Southern Alps of Aotearoa New Zealand red. Here we use time lapse cameras, airmass back trajectories, snow impurity geochemistry, and remote sensing to quantify the timing, provenance, and mass deposition of the event. Deposited in late November 2019, the impurities were dominated by mineral dust with a distinct southeastern Australian geochemical fingerprint. The event deposited ∼4,500 ± 500 tons of red dust to Southern Alps permanent snow and ice with a mean dust mass concentration of 6.5 ± 0.7 g m ⁻² . A southeast Australian desert dust storm generated by the same type of meteorological conditions as the 2020 New Year bushfires was the main driver of the glacier discoloration.

Questions: (1) What can be learned by extending a national classification into unsampled forest types? (2) Are both remotely sensed and environmental predictors needed to model and map associations? (3) For mapping, are LiDAR‐generated canopy structure parameters or reflectance from spectral imagery more useful? (4) How can we assess uncertainty of a final map?Location: Warawara Forest, Northland, New Zealand.Methods: We sampled 205 vegetation plots and assigned them to an existing national classification using the fuzzy classification framework of noise clustering. Plots too distinct to be assigned were used to define new associations. We produced spatial models of each association using boosted regression trees. Predictors included 11 environmental, 11 canopy reflectance, and 17 canopy structure variables. We created a composite map by assigning each map pixel to the association with the highest occurrence probability. We evaluated uncertainty by examining locations where no class was predicted with probability above 0.2 and by creating a confusion map based on entropy.Results: Forty‐five plots were assigned to six of 79 existing national associations and 147 plots were used to define two new forest associations. Three shrubland types are widespread nationally, whereas two young forest types are northern. Three mature forest types are narrowly distributed nationally, with the new “High‐elevation hardwood forest” largely restricted to Warawara Forest. Three associations were mapped using remotely sensed predictors alone, whereas two also required environmental predictors. Overall, canopy reflectance predictors explained more deviance than canopy structure. Examining locations where no association was predicted well and where multiple associations were predicted equally showed areas mapped as younger forests to have greatest uncertainty.Conclusions: In answering our questions, we present a vegetation classification and map for Warawara Forest that provides a framework to guide the indigenous people's management responses to threats to valued communities and their species.

A polynya is an area of open water or reduced concentration of sea ice surrounded by either concentrated sea ice or land ice. They are often seen as sites of intense ocean–atmosphere heat exchange and as ice production factories. Given their importance, it is crucial to quantify the accuracy of satellite-derived polynya information. Polynyas in their early evolution phase are generally narrow and occur at scales likely too fine to be detected by widely used passive microwave (PMW) radiometric sensors. We derived 40 m scale polynya information over the western Ross Sea from high-resolution Synthetic Aperture Radar (SAR) Sentinel-1 C-band data and examined discrepancies with larger-scale estimates. We utilized two automated algorithms, supervised (a rule-based approach) and unsupervised (a combination of texture analysis with k-means clustering), to accurately identify the polynya areas. We generated data for validation using Sentinel-1 data at instances where polynyas can be visually delineated. Results from PMW sensors (NSIDC and AMSR2) and SAR-based algorithms (rule-based and texture-based) are compared with manually delineated polynya areas obtained through Sentinel-1. Analysis using PMW sensors revealed that NSIDC overestimates larger polynyas and underestimates smaller polynyas compared to AMSR2. We were more accurately able to identify polynya presence and area using Sentinel-1 SAR observations, especially in clear cases and cases when PMW data miscalculates the polynya’s presence. Of our SAR-based algorithms, the rule-based approach was more accurate than the texture-based approach at identifying clear polynyas when validated against manually delineated regions. Altogether, we emphasize the need for finer spatio-temporal resolution data for polynya studies.

This study compares high-resolution output (1.1-km horizontal grid length) from twice-daily forecasts produced by the Antarctic Mesoscale Prediction System (AMPS) with a dense observational network east of Ross Island. Covering 10 000 km², 15 SNOWWEB stations significantly increased the number of observation stations in the area to 19 during the 2014–15 austral summer. Collocated “virtual stations” created from AMPS output are combined with observations, producing a single dataset of zonal and meridional wind components used to train a self-organizing map (SOM). The resulting SOM is used to individually classify the observational and AMPS datasets, producing a time series of classifications for each dataset directly comparable to the other. Analysis of class composites shows two dominant weather patterns: low but directionally variable winds and high but directionally constant winds linked to the Ross Ice Shelf airstream (RAS). During RAS events the AMPS and SNOWWEB data displayed good temporal class alignment with good surface wind correlation. SOM analysis shows that AMPS did not accurately forecast surface-level wind speed or direction during light wind conditions where synoptic forcing was weak; however, it was able to forecast the low wind period occurrence accurately. Coggins’s regimes provide synoptic-scale context and show a reduced synoptic pressure gradient during these classes, increasing reliance on the ability of Polar WRF to resolve mesoscale dynamics. Available initialization data have insufficient resolution for the region’s complex topography, which likely impacts performance. The SOM analysis methods used are shown to be effective for model validation and are widely applicable.