14 results
Search Results
Now showing 1 - 10 of 14
- ItemGrounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams(Katlenburg-Lindau : Copernicus, 2018) Reese, Ronja; Winkelmann, Ricarda; Gudmundsson, G. HilmarCurrently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.
- ItemPhotocurrent modelling and experimental confirmation for meteoric smoke particle detectors on board atmospheric sounding rockets(Katlenburg-Lindau : Copernicus, 2018-9-20) Giono, Gabriel; Strelnikov, Boris; Asmus, Heiner; Staszak, Tristan; Ivchenko, Nickolay; Lübken, Franz-JosefCharacterising the photoelectron current induced by the Sun's UV radiation is crucial to ensure accurate daylight measurements from particle detectors. This article lays out the methodology used to address this problem in the case of the meteoric smoke particle detectors (MSPDs), developed by the Leibniz Institute of Atmospheric Physics in Kühlungsborn (IAP) and flown on board the PMWEs (Polar Mesosphere Winter Echoes) sounding rockets in April 2018. The methodology focuses on two complementary aspects: modelling and experimental measurements. A detailed model of the MSPD photocurrent was created based on the expected solar UV flux, the atmospheric UV absorption as a function of height by molecular oxygen and ozone, the photoelectric yield of the material coating the MSPD as a function of wavelength, the index of refraction of these materials as a function of wavelength and the angle of incidence of the illumination onto the MSPD. Due to its complex structure, composed of a central electrode shielded by two concentric grids, extensive ray-tracing calculations were conducted to obtain the incidence angles of the illumination on the central electrode, and this was done for various orientations of the MSPD in respect to the Sun. Results of the modelled photocurrent at different heights and for different materials, as well as for different orientations of the detector, are presented. As a pre-flight confirmation, the model was used to reproduce the experimental measurements conducted by Robertson et al. (2014) and agrees within an order of magnitude. An experimental setup for the calibration of the MSPD photocurrent is also presented. The photocurrent induced by the Lyman-alpha line from a deuterium lamp was recorded inside a vacuum chamber using a narrowband filter, while a UV-sensitive photodiode was used to monitor the UV flux. These measurements were compared with the model prediction, and also matched within an order of magnitude. Although precisely modelling the photocurrent is a challenging task, this article quantitatively improved the understanding of the photocurrent on the MSPD and discusses possible strategies to untangle the meteoric smoke particles (MSPs) current from the photocurrent recorded in-flight.
- ItemVertical profiles of aerosol mass concentration derived by unmanned airborne in situ and remote sensing instruments during dust events(Katlenburg-Lindau : Copernicus, 2018) Mamali, Dimitra; Marinou, Eleni; Sciare, Jean; Pikridas, Michael; Kokkalis, Panagiotis; Kottas, Michael; Binietoglou, Ioannis; Tsekeri, Alexandra; Keleshis, Christos; Engelmann, Ronny; Baars, Holger; Ansmann, Albert; Amiridis, Vassilis; Russchenberg, Herman; Biskos, GeorgeIn situ measurements using unmanned aerial vehicles (UAVs) and remote sensing observations can independently provide dense vertically resolved measurements of atmospheric aerosols, information which is strongly required in climate models. In both cases, inverting the recorded signals to useful information requires assumptions and constraints, and this can make the comparison of the results difficult. Here we compare, for the first time, vertical profiles of the aerosol mass concentration derived from light detection and ranging (lidar) observations and in situ measurements using an optical particle counter on board a UAV during moderate and weak Saharan dust episodes. Agreement between the two measurement methods was within experimental uncertainty for the coarse mode (i.e. particles having radii > 0.5 μm), where the properties of dust particles can be assumed with good accuracy. This result proves that the two techniques can be used interchangeably for determining the vertical profiles of aerosol concentrations, bringing them a step closer towards their systematic exploitation in climate models.
- ItemSeparation of the optical and mass features of particle components in different aerosol mixtures by using POLIPHON retrievals in synergy with continuous polarized Micro-Pulse Lidar (P-MPL) measurements(Katlenburg-Lindau : Copernicus, 2018) Córdoba-Jabonero, Carmen; Sicard, Michaël; Ansmann, Albert; del Águila, Ana; Baars, HolgerThe application of the POLIPHON (POlarization-LIdar PHOtometer Networking) method is presented for the first time in synergy with continuous 24/7 polarized Micro-Pulse Lidar (P-MPL) measurements to derive the vertical separation of two or three particle components in different aerosol mixtures, and the retrieval of their particular optical properties. The procedure of extinction-to-mass conversion, together with an analysis of the mass extinction efficiency (MEE) parameter, is described, and the relative mass contribution of each aerosol component is also derived in a further step. The general POLIPHON algorithm is based on the specific particle linear depolarization ratio given for different types of aerosols and can be run in either 1-step (POL-1) or 2 steps (POL-2) versions with dependence on either the 2- or 3-component separation. In order to illustrate this procedure, aerosol mixing cases observed over Barcelona (NE Spain) are selected: a dust event on 5 July 2016, smoke plumes detected on 23 May 2016 and a pollination episode observed on 23 March 2016. In particular, the 3-component separation is just applied for the dust case: a combined POL-1 with POL-2 procedure (POL-1/2) is used, and additionally the fine-dust contribution to the total fine mode (fine dust plus non-dust aerosols) is estimated. The high dust impact before 12:00 UTC yields a mean mass loading of 0.6±0.1 g m'2 due to the prevalence of Saharan coarse-dust particles. After that time, the mean mass loading is reduced by two-thirds, showing a rather weak dust incidence. In the smoke case, the arrival of fine biomass-burning particles is detected at altitudes as high as 7 km. The smoke particles, probably mixed with less depolarizing non-smoke aerosols, are observed in air masses, having their origin from either North American fires or the Arctic area, as reported by HYSPLIT back-trajectory analysis. The particle linear depolarization ratio for smoke shows values in the 0.10-0.15 range and even higher at given times, and the daily mean smoke mass loading is 0.017±0.008 g m'2, around 3 % of that found for the dust event. Pollen particles are detected up to 1.5 km in height from 10:00 UTC during an intense pollination event with a particle linear depolarization ratio ranging between 0.10 and 0.15. The maximal mass loading of Platanus pollen particles is 0.011±0.003 g m'2, representing around 2 % of the dust loading during the higher dust incidence. Regarding the MEE derived for each aerosol component, their values are in agreement with others referenced in the literature for the specific aerosol types examined in this work: 0.5±0.1 and 1.7±0.2 m2 g'1 are found for coarse and fine dust particles, 4.5±1.4 m2 g'1 is derived for smoke and 2.4±0.5 m2 g'1 for non-smoke aerosols with Arctic origin, and a MEE of 2.4±0.8 m2 g'1 is obtained for pollen particles, though it can reach higher or lower values depending on predominantly smaller or larger pollen grain sizes. Results reveal the high potential of the P-MPL system, a simple polarization-sensitive elastic backscatter lidar working in a 24/7 operation mode, to retrieve the relative optical and mass contributions of each aerosol component throughout the day, reflecting the daily variability of their properties. In fact, this procedure can be simply implemented in other P-MPLs that also operate within the worldwide Micro-Pulse Lidar Network (MPLNET), thus extending the aerosol discrimination at a global scale. Moreover, the method has the advantage of also being relatively easily applicable to space-borne lidars with an equivalent configuration such as the ongoing Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) on board NASA CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) and the forthcoming Atmospheric Lidar (ATLID) on board the ESA EarthCARE mission.
- ItemIntercomparison of middle-atmospheric wind in observations and models(Katlenburg-Lindau : Copernicus, 2018-4-6) Rüfenacht, Rolf; Baumgarten, Gerd; Hildebrand, Jens; Schranz, Franziska; Matthias, Vivien; Stober, Gunter; Lübken, Franz-Josef; Kämpfer, NiklausWind profile information throughout the entire upper stratosphere and lower mesosphere (USLM) is important for the understanding of atmospheric dynamics but became available only recently, thanks to developments in remote sensing techniques and modelling approaches. However, as wind measurements from these altitudes are rare, such products have generally not yet been validated with (other) observations. This paper presents the first long-term intercomparison of wind observations in the USLM by co-located microwave radiometer and lidar instruments at Andenes, Norway (69.3∘ N, 16.0∘ E). Good correspondence has been found at all altitudes for both horizontal wind components for nighttime as well as daylight conditions. Biases are mostly within the random errors and do not exceed 5–10 m s−1, which is less than 10 % of the typically encountered wind speeds. Moreover, comparisons of the observations with the major reanalyses and models covering this altitude range are shown, in particular with the recently released ERA5, ECMWF's first reanalysis to cover the whole USLM region. The agreement between models and observations is very good in general, but temporally limited occurrences of pronounced discrepancies (up to 40 m s−1) exist. In the article's Appendix the possibility of obtaining nighttime wind information about the mesopause region by means of microwave radiometry is investigated.
- ItemThe Fifth International Workshop on Ice Nucleation phase 2 (FIN-02): Laboratory intercomparison of ice nucleation measurements(Katlenburg-Lindau : Copernicus, 2018) DeMott, Paul J.; Möhler, Ottmar; Cziczo, Daniel J.; Hiranuma, Naruki; Petters, Markus D.; Petters, Sarah S.; Belosi, Franco; Bingemer, Heinz G.; Brooks, Sarah D.; Budke, Carsten; Burkert-Kohn, Monika; Collier, Kristen N.; Danielczok, Anja; Eppers, Oliver; Felgitsch, Laura; Garimella, Sarvesh; Grothe, Hinrich; Herenz, Paul; Hill, Thomas C. J.; Höhler, Kristina; Kanji, Zamin A.; Kiselev, Alexei; Koop, Thomas; Kristensen, Thomas B.; Krüger, Konstantin; Kulkarni, Gourihar; Levin, Ezra J. T.; Murray, Benjamin J.; Nicosia, Alessia; O'Sullivan, Daniel; Peckhaus, Andreas; Polen, Michael J.; Price, Hannah C.; Reicher, Naama; Rothenberg, Daniel A.; Rudich, Yinon; Santachiara, Gianni; Schiebel, Thea; Schrod, Jann; Seifried, Teresa M.; Stratmann, Frank; Sullivan, Ryan C.; Suski, Kaitlyn J.; Szakáll, Miklós; Taylor, Hans P.; Ullrich, Romy; Vergara-Temprado, Jesus; Wagner, Robert; Whale, Thomas F.; Weber, Daniel; Welti, André; Wilson, Theodore W.; Wolf, Martin J.; Zenker, JakeThe second phase of the Fifth International Ice Nucleation Workshop (FIN-02) involved the gathering of a large number of researchers at the Karlsruhe Institute of Technology's Aerosol Interactions and Dynamics of the Atmosphere (AIDA) facility to promote characterization and understanding of ice nucleation measurements made by a variety of methods used worldwide. Compared to the previous workshop in 2007, participation was doubled, reflecting a vibrant research area. Experimental methods involved sampling of aerosol particles by direct processing ice nucleation measuring systems from the same volume of air in separate experiments using different ice nucleating particle (INP) types, and collections of aerosol particle samples onto filters or into liquid for sharing amongst measurement techniques that post-process these samples. In this manner, any errors introduced by differences in generation methods when samples are shared across laboratories were mitigated. Furthermore, as much as possible, aerosol particle size distribution was controlled so that the size limitations of different methods were minimized. The results presented here use data from the workshop to assess the comparability of immersion freezing measurement methods activating INPs in bulk suspensions, methods that activate INPs in condensation and/or immersion freezing modes as single particles on a substrate, continuous flow diffusion chambers (CFDCs) directly sampling and processing particles well above water saturation to maximize immersion and subsequent freezing of aerosol particles, and expansion cloud chamber simulations in which liquid cloud droplets were first activated on aerosol particles prior to freezing. The AIDA expansion chamber measurements are expected to be the closest representation to INP activation in atmospheric cloud parcels in these comparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX and potassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam (Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively, and a bacterial INP (Snomax®). Considered together, the agreement among post-processed immersion freezing measurements of the numbers and fractions of particles active at different temperatures following bulk collection of particles into liquid was excellent, with possible temperature uncertainties inferred to be a key factor in determining INP uncertainties. Collection onto filters for rinsing versus directly into liquid in impingers made little difference. For methods that activated collected single particles on a substrate at a controlled humidity at or above water saturation, agreement with immersion freezing methods was good in most cases, but was biased low in a few others for reasons that have not been resolved, but could relate to water vapor competition effects. Amongst CFDC-style instruments, various factors requiring (variable) higher supersaturations to achieve equivalent immersion freezing activation dominate the uncertainty between these measurements, and for comparison with bulk immersion freezing methods. When operated above water saturation to include assessment of immersion freezing, CFDC measurements often measured at or above the upper bound of immersion freezing device measurements, but often underestimated INP concentration in comparison to an immersion freezing method that first activates all particles into liquid droplets prior to cooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber-Portable Ice Nucleation Chamber), and typically slightly underestimated INP number concentrations in comparison to cloud parcel expansions in the AIDA chamber; this can be largely mitigated when it is possible to raise the relative humidity to sufficiently high values in the CFDCs, although this is not always possible operationally. Correspondence of measurements of INPs among direct sampling and post-processing systems varied depending on the INP type. Agreement was best for Snomax® particles in the temperature regime colder than -10°C, where their ice nucleation activity is nearly maximized and changes very little with temperature. At temperatures warmer than -10°C, Snomax® INP measurements (all via freezing of suspensions) demonstrated discrepancies consistent with previous reports of the instability of its protein aggregates that appear to make it less suitable as a calibration INP at these temperatures. For Argentinian soil dust particles, there was excellent agreement across all measurement methods; measures ranged within 1 order of magnitude for INP number concentrations, active fractions and calculated active site densities over a 25 to 30°C range and 5 to 8 orders of corresponding magnitude change in number concentrations. This was also the case for all temperatures warmer than -25°C in Tunisian dust experiments. In contrast, discrepancies in measurements of INP concentrations or active site densities that exceeded 2 orders of magnitude across a broad range of temperature measurements found at temperatures warmer than -25°C in a previous study were replicated for illite NX. Discrepancies also exceeded 2 orders of magnitude at temperatures of -20 to -25°C for potassium feldspar (K-feldspar), but these coincided with the range of temperatures at which INP concentrations increase rapidly at approximately an order of magnitude per 2°C cooling for K-feldspar. These few discrepancies did not outweigh the overall positive outcomes of the workshop activity, nor the future utility of this data set or future similar efforts for resolving remaining measurement issues. Measurements of the same materials were repeatable over the time of the workshop and demonstrated strong consistency with prior studies, as reflected by agreement of data broadly with parameterizations of different specific or general (e.g., soil dust) aerosol types. The divergent measurements of the INP activity of illite NX by direct versus post-processing methods were not repeated for other particle types, and the Snomax° data demonstrated that, at least for a biological INP type, there is no expected measurement bias between bulk collection and direct immediately processed freezing methods to as warm as -10°C. Since particle size ranges were limited for this workshop, it can be expected that for atmospheric populations of INPs, measurement discrepancies will appear due to the different capabilities of methods for sampling the full aerosol size distribution, or due to limitations on achieving sufficient water supersaturations to fully capture immersion freezing in direct processing instruments. Overall, this workshop presents an improved picture of present capabilities for measuring INPs than in past workshops, and provides direction toward addressing remaining measurement issues.
- ItemSimulation of the future sea level contribution of Greenland with a new glacial system model(Katlenburg-Lindau : Copernicus, 2018) Calov, Reinhard; Beyer, Sebastian; Greve, Ralf; Beckmann, Johanna; Willeit, Matteo; Kleiner, Thomas; Rückamp, Martin; Humbert, Angelika; Ganopolski, AndreyWe introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. The aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961-1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961-1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation-surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm until the year 2300, including our simulated additional sea level rise due to elevation-surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.
- ItemA confined-unconfined aquifer model for subglacial hydrology and its application to the Northeast Greenland Ice Stream(Katlenburg-Lindau : Copernicus, 2018) Beyer, Sebastian; Kleiner, Thomas; Aizinger, Vadym; Rückamp, Martin; Humbert, AngelikaSubglacial hydrology plays an important role in ice sheet dynamics as it determines the sliding velocity. It also drives freshwater into the ocean, leading to undercutting of calving fronts by plumes. Modeling subglacial water has been a challenge for decades. Only recently have new approaches been developed such as representing subglacial channels and thin water sheets by separate layers of variable hydraulic conductivity. We extend this concept by modeling a confined-unconfined aquifer system (CUAS) in a single layer of an equivalent porous medium (EPM). The advantage of this formulation is that it prevents unphysical values of pressure at reasonable computational cost. We performed sensitivity tests to investigate the effect of different model parameters. The strongest influence of model parameters was detected in terms of governing the opening and closure of the system. Furthermore, we applied the model to the Northeast Greenland Ice Stream, where an efficient system independent of seasonal input was identified about 500 km downstream from the ice divide. Using the effective pressure from the hydrology model, the Ice Sheet System Model (ISSM) showed considerable improvements in modeled velocities in the coastal region.
- ItemSALSA2.0: The sectional aerosol module of the aerosol-chemistry-climate model ECHAM6.3.0-HAM2.3-MOZ1.0(Katlenburg-Lindau : Copernicus, 2018) Kokkola, Harri; Kühn, Thomas; Laakso, Anton; Bergman, Tommi; Lehtinen, Kari E. J.; Mielonen, Tero; Arola, Antti; Stadtler, Scarlet; Korhonen, Hannele; Ferrachat, Sylvaine; Lohmann, Ulrike; Neubauer, David; Tegen, Ina; Siegenthaler-Le Drian, Colombe; Schultz, Martin G.; Bey, Isabelle; Stier, Philip; Daskalakis, Nikos; Heald, Colette L.; Romakkaniemi, SamiIn this paper, we present the implementation and evaluation of the aerosol microphysics module SALSA2.0 in the framework of the aerosol-chemistry-climate model ECHAM-HAMMOZ. It is an alternative microphysics module to the default modal microphysics scheme M7 in ECHAM-HAMMOZ. The SALSA2.0 implementation within ECHAM-HAMMOZ is evaluated against observations of aerosol optical properties, aerosol mass, and size distributions, comparing also to the skill of the M7 implementation. The largest differences between the implementation of SALSA2.0 and M7 are in the methods used for calculating microphysical processes, i.e., nucleation, condensation, coagulation, and hydration. These differences in the microphysics are reflected in the results so that the largest differences between SALSA2.0 and M7 are evident over regions where the aerosol size distribution is heavily modified by the microphysical processing of aerosol particles. Such regions are, for example, highly polluted regions and regions strongly affected by biomass burning. In addition, in a simulation of the 1991 Mt. Pinatubo eruption in which a stratospheric sulfate plume was formed, the global burden and the effective radii of the stratospheric aerosol are very different in SALSA2.0 and M7. While SALSA2.0 was able to reproduce the observed time evolution of the global burden of sulfate and the effective radii of stratospheric aerosol, M7 strongly overestimates the removal of coarse stratospheric particles and thus underestimates the effective radius of stratospheric aerosol. As the mode widths of M7 have been optimized for the troposphere and were not designed to represent stratospheric aerosol, the ability of M7 to simulate the volcano plume was improved by modifying the mode widths, decreasing the standard deviations of the accumulation and coarse modes from 1.59 and 2.0, respectively, to 1.2 similar to what was observed after the Mt. Pinatubo eruption. Overall, SALSA2.0 shows promise in improving the aerosol description of ECHAM-HAMMOZ and can be further improved by implementing methods for aerosol processes that are more suitable for the sectional method, e.g., size-dependent emissions for aerosol species and size-resolved wet deposition.
- ItemA generic pixel-to-point comparison for simulated large-scale ecosystem properties and ground-based observations: An example from the Amazon region(Katlenburg-Lindau : Copernicus, 2018) Rammig, Anja; Heinke, Jens; Hofhansl, Florian; Verbeeck, Hans; Baker, Timothy R.; Christoffersen, Bradley; Ciais, Philippe; De Deurwaerder, Hannes; Fleischer, Katrin; Galbraith, David; Guimberteau, Matthieu; Huth, Andreas; Johnson, Michelle; Krujit, Bart; Langerwisch, Fanny; Meir, Patrick; Papastefanou, Phillip; Sampaio, Gilvan; Thonicke, Kirsten; von Randow, Celso; Zang, Christian; Rödig, EdnaComparing model output and observed data is an important step for assessing model performance and quality of simulation results. However, such comparisons are often hampered by differences in spatial scales between local point observations and large-scale simulations of grid cells or pixels. In this study, we propose a generic approach for a pixel-to-point comparison and provide statistical measures accounting for the uncertainty resulting from landscape variability and measurement errors in ecosystem variables. The basic concept of our approach is to determine the statistical properties of small-scale (within-pixel) variability and observational errors, and to use this information to correct for their effect when large-scale area averages (pixel) are compared to small-scale point estimates. We demonstrate our approach by comparing simulated values of aboveground biomass, woody productivity (woody net primary productivity, NPP) and residence time of woody biomass from four dynamic global vegetation models (DGVMs) with measured inventory data from permanent plots in the Amazon rainforest, a region with the typical problem of low data availability, potential scale mismatch and thus high model uncertainty. We find that the DGVMs under- and overestimate aboveground biomass by 25% and up to 60%, respectively. Our comparison metrics provide a quantitative measure for model-data agreement and show moderate to good agreement with the region-wide spatial biomass pattern detected by plot observations. However, all four DGVMs overestimate woody productivity and underestimate residence time of woody biomass even when accounting for the large uncertainty range of the observational data. This is because DGVMs do not represent the relation between productivity and residence time of woody biomass correctly. Thus, the DGVMs may simulate the correct large-scale patterns of biomass but for the wrong reasons. We conclude that more information about the underlying processes driving biomass distribution are necessary to improve DGVMs. Our approach provides robust statistical measures for any pixel-to-point comparison, which is applicable for evaluation of models and remote-sensing products.