Three-dimensional Radiative Transfer Model Research Articles

Simulating High-Resolution Sun-Induced Chlorophyll Fluorescence Image of Three-Dimensional Canopy Based on Photon Mapping

The remote sensing of sun-induced chlorophyll fluorescence (SIF) is an emerging technique with immense potential for terrestrial vegetation sciences. However, the interpretation of fluorescence data is often hindered by the complexity of observed land surfaces. Therefore, advanced remote sensing models, particularly physically based simulations, are critical to accurately interpret SIF data. In this work, we propose a three-dimensional (3D) radiative transfer model that employs the Monte Carlo ray-tracing technique to simulate the excitation and transport of SIF within plant canopies. This physically based approach can quantify the various radiative processes contributing to the observed SIF signal with high fidelity. The model’s performance is rigorously evaluated by comparing the simulated SIF spectra and angular distributions to field measurements, as well as conducting systematic comparisons with an established radiative transfer model. The results demonstrate the proposed model’s ability to reliably reproduce the key spectral and angular characteristics of SIF, with the coefficient of determination (R2) exceeding 0.98 and root mean square error (RMSE) being less than 0.08 mW m−2 sr−1 nm−1 for both the red and far-red fluorescence peaks. Furthermore, the model’s versatile representation of canopy structures, enabled by the decoupling of radiation and geometry, is applied to study the impact of 3D structure on SIF patterns. This capability makes the proposed model a highly attractive tool for investigating SIF distributions in realistic, heterogeneous canopy environments.

A dataset of topographic correction coefficients for shortwave downward radiation over the Pan-Third Pole

Prevalent Shortwave downward radiation (SWDR) estimates assume a flat surface, neglecting topographic effects and leading to significant errors in mountainous regions. We introduce SWDR topography correction coefficients (TCCs), based on the mountain radiative transfer model tailored for the Pan-Third Pole region. This dataset effectively bridges the disparities between flat-surface SWDR and rugged-surface SWDR, forming part of the Long-term Earth System spatiotemporally Seamless Radiation budget dataset (LessRad). Validation results using a three-dimensional radiative transfer model demonstrate the efficacy of this method in correcting solar direct radiation, sky diffuse radiation, and SWDR under diverse conditions. At a spatial resolution of 2.5 arc-minutes, the correction accuracy for solar direct radiation is characterized by a coefficient of determination (R²) of 0.998, a relative root mean square error (rRMSE) of 2.4%, and a relative bias (rbias) of 0.8%. For sky diffused radiation, an R² of 0.965, a rRMSE of 1.2%, and a rbias of −0.8%. SWDR corrections under clear and cloudy skies also show high accuracy, demonstrating the robustness of the TCCs approach.

Numerical investigation of optical characterization of polycarbonate panels

Polycarbonate panels (PC panels) are state-of-the-art transparent insulating materials widely used in the construction industry due to their cavity structure, which provides exceptional thermal insulation and optimal optical performance. However, the inherent anisotropy of the three-dimensional cavity structure complicates radiative transfer and requires consideration of both azimuth and zenith angles in optical performance evaluation. This aspect has received limited attention in existing research. This study aims to accurately characterize the optical performance of PC panels through numerical simulations. A three-dimensional radiative transfer model based on the discrete ordinate radiation model is developed to solve the radiation transfer equation. The model's independency regarding mesh division, angular discretization, and accuracy is validated. The effects of incidence angle, geometric parameters, and optical properties of PC panels on optical performance are analyzed. The findings reveal a strong correlation between transmittance and absorption with variations in incident zenith and azimuth angles. The transmittance exhibits a consistent monotonic variation expressible as a rational bifunction. Notably, absorption peaks occur within specific solid angle ranges, with increased structural complexity resulting in heightened absorption and greater uncertainty. For conventional PC materials, maximum transmittance ranges from 46.9 % to 73 %, while maximum absorption ranges from 2.3 % to 13.5 %. Increasing absorption coefficients, refractive index, and surface scattering coefficients nonlinearly decrease transmittance while increasing absorption. Additionally, deviations in transmittance and absorption with azimuth angle amplify with an increase in non-horizontal structures. Sensitivity analysis indicates a significant influence of zenith angle on transmittance, and absorption coefficient predominantly affects absorption.

Nuclear spin ratios of deuterated ammonia in prestellar cores

Context. Molecules containing two or more hydrogen or deuterium atoms have different nuclear spin states which behave as separate chemical species. The relative abundances of these species can give clues to their origin. Formation on grains is believed to yield statistical spin ratios whereas gas-phase reactions are predicted to result in clear deviations from them. This is also true for ammonia and its deuterated forms NH2D, NHD2, and ND3. Aims. Here we aim to determine the ortho/para ratios of NH2D and NHD2 in dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions. Methods. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores H-MM1 and Oph D in Ophiuchus in the ground-state lines of ortho and para NH2D and NHD2. The fractional abundances of these molecules were derived employing three-dimensional radiative transfer modelling, using different assumptions about the abundance profiles as functions of density. We also ran gas-grain chemistry models with different scenarios concerning proton or deuteron exchanges and chemical desorption from grains to find out if one of these models can reproduce the observed spin ratios. Results. The observationally deduced ortho/para ratios of NH2D and NHD2 are in both cores within 10% of their statistical values 3 and 2, respectively, and taking 3 σ limits, deviations from these of about 20% are allowed. Of the chemistry models tested here, the model that assumes proton hop (as opposed to full scrambling) in reactions contributing to ammonia formation, and a constant efficiency of chemical desorption, comes nearest to the observed abundances and spin ratios. Conclusions. The nuclear spin ratios derived here are in contrast with spin-state chemistry models that assume full scrambling in proton donation and hydrogen abstraction reactions leading to deuterated ammonia. The efficiency of chemical desorption strongly influences the predicted abundances of NH3, NH2D, and NHD2, but has a lesser effect on their ortho/para ratios. For these the proton exchange scenario in the gas is decisive. We suggest that this is because of rapid re-processing of ammonia and related cations by gas-phase ion-molecule reactions.

Calibration of DART 3D model with UAV and Sentinel-2 for studying the radiative budget of conventional and agro-ecological maize fields

Quantifying the radiative budget (RB) of maize, the world's predominant cereal, is crucial to managing its growth in a context of global warming. Its microclimate including temperature, water, and energy flux depend on crop management practices. This paper compares the influence of agroecological practices (AE) and conventional agriculture (CA) on the RB of maize using three-dimensional (3D) radiative transfer modelling and remote sensing (RS) observations. We studied two neighbour maize fields in south-west France, respectively with AE and CA practices. Time series of photosynthetically active radiation absorbed (APAR) by plants (APARplant) and soil (APARsoil) were simulated with the DART radiative transfer model. First, realistic 3D models of the CA and AE fields were created and validated with a new method based on DART and Sentinel 2 and UAV reflectance and vegetation indices acquired in July 2019. At that date, around midday, APARplant was larger by 21.5 W/m2 in the AE field and APARsoil was larger by 20.1 W/m2 in the CA field. We explained these differences by differences in geometric and optical parameters (OP) between the two fields. OPsoil causes ≈45% of the differences in APARplant and APARsoil. The shape of plant causes ≈14% of the difference in APARplant and ≈17% in APARsoil. OPplant causes ≈14% of the difference in APARplant and ≈1% in APARsoil. Field geometry (row orientation, inter-row and plant distance) causes ≈7% of the difference in APARplant and 2% in APARsoil. It stresses the great role of OPsoil associated to AE practices. Because it integrates the non-linear effects of all above parameters, the LAI (Leaf Area Index) accounts for ≈40% difference in APARplant and APARsoil. The sensitivity of APAR to LAI is twice as low in the AE field than in the CA field. We showed that APAR differences between the AE and CA fields can be even greater with other maize row orientations. We also highlighted the role of the field model on simulated APAR by comparing the use of APAR simulated by a 3D field model and its corresponding 1D turbid field: 1D turbid models resulted in a 16% increase in APARplant and in similar APARsoil, compared to the 3D CA and AE models.

Possible episodic infall towards a compact disk in B335

Context. Previous observations of the isolated Class 0 source B335 have presented evidence of ongoing infall in various molecular lines, such as HCO+, HCN, and CO. There have been no confirmed observations of a rotationally supported disk on scales greater than ~12 au. Aims. The presence of an outflow in B335 suggests that a disk is also expected to be present or undergoing formation. To constrain the earliest stages of protostellar evolution and disk formation, we aim to map the region where gas falls inwards and observationally constrain its kinematics. Furthermore, we aim to put strong limits on the size and orientation of any disk-like structure in B335. Methods. We used high angular resolution 13CO data from the Atacama Large Millimeter/submillimeter Array (ALMA) and combined it with shorter-baseline archival data to produce a high-fidelity image of the infall in B335. We also revisited the imaging of high-angular resolution Band 6 continuum data to study the dust distribution in the immediate vicinity of B335. Results. Continuum emission shows an elliptical structure (10 by 7 au) with a position angle 5 degrees east of north, consistent with the expectation for a forming disk in B335. A map of the infall velocity (as estimated from the 13CO emission), shows evidence of asymmetric infall, predominantly from the north and south. Close to the protostar, infall velocities appear to exceed free-fall velocities. Three-dimensional (3D) radiative transfer models, where the infall velocity is allowed to vary within the infall region, may explain the observed kinematics. Conclusions. The data suggest that a disk has started to form in B335 and that gas is falling towards that disk. However, kinematically-resolved line data towards the disk itself is needed to confirm the presence of a rotationally supported disk around this young protostar. The high infall velocities we measured are not easily reconcilable with a magnetic braking scenario, suggesting that there is a pressure gradient that allows the infall velocity to vary in the region.

Evaluating the saturation effect of vegetation indices in forests using 3D radiative transfer simulations and satellite observations

Vegetation indices (VIs) have been used extensively for qualitative and quantitative remote sensing monitoring of vegetation vigor and growth dynamics. However, the saturation phenomenon of VIs (i.e., insignificant change at moderate to high vegetation densities) poses a known limitation to their ability to characterize surface vegetation over the dense canopy. Although the mechanisms underlying saturation are relatively straightforward and several VIs have been proposed to mitigate the saturation effect, the assessment of the saturation effect of VIs remains insufficient. Notably, no unified metric has been proposed to quantify the VI saturation phenomenon, limiting VI selection in practical applications. In this study, we proposed two indicators to describe the saturation phenomenon and utilized a well-validated three-dimensional (3D) canopy radiative transfer (RT) model large-scale remote sensing data and image simulation framework (LESS) to simulate the bidirectional reflectance factor (BRF) of six forests scenes and assessed the variations in VIs in relation to leaf area index (LAI) values over different backgrounds, sun-sensor geometries, and spatial distribution types. The saturation characteristics of 36 VIs were evaluated in combination with simulation results and satellite observations from multiple sensors. The ranking of VI saturation from simulated and satellite results revealed a good agreement. Our results indicated that the simple ratio vegetation index (SR) performed best with the highest saturation point and can well characterize the surface vegetation condition until LAI reaches 4. Besides, we found that the saturation effect of VIs was influenced by soil brightness, sun-sensor geometry, and canopy structure. SR, modified simple ratio (MSR) and normalized green red difference index (NGRDI) were the most susceptible to these disturbing factors, although they had higher resistance to saturation. Modified triangular vegetation index 1 (MTVI1), modified non-linear vegetation index (MNLI), triangular greenness index (TGI), and triangular vegetation index (TriVI) performed well overall, combining the ability to resist saturation and disturbance factors. Appropriate application of VIs can help better understand vegetation responses to climate change and accurately assess ecosystem status. Our results contribute to the understanding of the VI saturation effect and provide a combined model and satellite data experimental workflow in appropriate VI selection to accurately characterize vegetation.

Models of Rotating Infall for the B335 Protostar

Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show redshifted absorption against a central continuum source indicative of infall in the HCO+ and HCN J = 4 → 3 transitions. The data are combined with a new estimate of the distance to provide strong constraints to three-dimensional radiative transfer models including a rotating, infalling envelope, outflow cavities, and a very small disk. The models favor ages since the initiation of collapse between 3 × 104 and 4 × 104 yr for both the continuum and the lines, resolving a conflict found in one-dimensional models. The models underpredict the continuum emission seen by ALMA, suggesting an additional component such as a pseudo-disk. The best-fitting model is used to convert variations in the 4.5 μm flux in recent years into a model for a variation of a factor of 5–7 in luminosity over the last 8 yr.

GF-1 leaf area index product across China based on three-dimensional stochastic radiation transfer model

GF-1 leaf area index product across China based on three-dimensional stochastic radiation transfer model

CANOPY STRUCTURAL EFFECTS ON BIDIRECTIONAL REFLECTANCE SIMULATED BY THE LESS MODEL: A CASE STUDY OF PICEA CRASSIFOLIA FORESTS

Abstract. Surface albedo is a dominant factor affecting the earth energy balance. Bidirectional reflectance distribution function (BRDF) describes the anisotropic surface reflection of solar radiation. Bidirectional effect is especially noticeable in the forest canopies with three-dimensional canopy structure. Thus, it is challenging to evaluate the canopy structural effects on BRDF using traditional radiative transfer modelling methods, where leaf in canopies is generally assumed to be homogeneous. LESS, a ray-tracing-based three-dimensional (3D) radiative transfer model can effectively simulate bidirectional reflectance factor while fully considering multi-component spectral and structural characteristics of vegetation. In this study, we applied the LESS model to evaluate the impacts of changes in leaf area index (LAI) on canopy reflectance in a Picea crassifolia forest. First, canopy structural parameters obtained from forest inventory were used to construct forest scenes and drive the LESS model. Then, we validated the simulation accuracy of the LESS model against BRDF calculated from MODIS BRDF parameter product using kernel-driven Ross-Li bi-directional reflectance function. Finally, we constructed forest scenes with different LAI to investigate the canopy structural effects on BRDF. We found that the BRDF simulated by the LESS model is in good accordance with the MODIS BRDF data in the red (R2 = 0.97, RMSE = 0.03) and NIR (R2 = 0.94, RMSE = 0.05) bands. The simulated canopy reflectance in the red band displayed ‘dome’ shape and decreased with increasing LAI. In contrast, canopy reflectance in the NIR band displayed ‘bowl’ shape and increased with increasing LAI. Our study highlighted the sensitivity of canopy BRDF to changes in canopy structure and have implications for retrieving LAI from BRDF data.

Evaluation of the impact of crop residue on fractional vegetation cover estimation by vegetation indices over conservation tillage cropland: a simulation study

ABSTRACT Accurate estimation of fractional vegetation cover (FVC) is of great significance to agricultural production. Crop residue management affect crop residue cover (CRC) over croplands. Crop and crop residue on the soil surface both contribute to overall canopy reflectance. Few studies, however, have examined the effect of crop residue on vegetation indices (VIs) and estimated FVC. The present study evaluated the response of eight commonly used VIs to crop residues and FVC uncertainty caused by crop residue based on the dimidiate pixel model (DPM) by using simulated reflectance of low-tilled cropland via a three-dimensional radiative transfer model. The absolute difference (AD) was used to quantify the spectral difference between crop residues and soils in red and near infrared wavelengths. Increases in normalized difference VI (NDVI), ratio VI (RVI), transformed soil-adjusted VI (TSAVI), and normalized difference phenology index (NDPI) were observed when green crops were mixed with crop residue that had negative ADs with soils, but decreases in enhance VI (EVI), perpendicular VI (PVI), SAVI, and litter-soil-adjusted VI (L-SAVI) were observed when crop residue was present under medium and high vegetation cover. The presence of crop residue with a positive AD with soils reduced NDVI, RVI, TSAVI, and NDPI while increased the other VIs. Crop residue had the least impact on EVI- and SAVI-based DPMs, with FVC-estimated uncertainty less than 0.1, followed by the NDPI- and L-SAVI-based model, while DPMs based on NDVI- and RVI performed poorly. Each VI-based DPM’s estimated uncertainty was highly correlated with AD values. Furthermore, the majority of the VI-based models were sensitive to solar position except for the NDPI-based model. Our findings highlight the need of considering the impact of crop residue on FVC retrieval over low-tilled cropland in future research.

Research on Space-Based Visible Detection for Conical Space Targets

We conducted research into the visible light detection of conical space targets and proposed a method to determine the caliber of a space-based detection system. We constructed a three-dimensional conical space target radiation transfer model and analyzed the effect of the phase angle between the solar incident direction, the detector’s bore sight direction, and the conical target inclination on the radiation intensity of the space target received by the detector. Finally, we obtained the threshold caliber which can meet the detection requirements of the system with the system signal-to-noise ratio as the evaluation standard, and we took the SBV sensor as an example to verify that the caliber determination method in this paper is effective.

Effects of Mixture Mode on the Canopy Bidirectional Reflectance of Coniferous–Broadleaved Mixed Plantations

One of the main initiatives for China to achieve the goal of being carbon neutral before 2060 is transforming monocultures into mixed plantations in subtropical China, because mixed forests possess a higher quality than monocultures in various ways. Very high spatial resolution (VHR) satellite imagery is very promising to precisely monitor the transformation process under the premise of clarifying the canopy reflectance anisotropy of mixed plantations. However, it is almost impossible to understand the canopy reflectance anisotropy of mixed plantations with real satellite data due to the extreme lack of multiangular VHR satellite images. In this study, the effects of the mixture mode on the canopy bidirectional reflectance factor (BRF) were comprehensively analyzed with simulated VHR images. The three-dimensional (3D) Discrete Anisotropic Radiative Transfer model (DART) was used to construct a pure coniferous scene, a pure broadleaved scene, and 27 coniferous–broadleaved mixed plantation scenes containing 3 mixture patterns (i.e., mixed by single trees, mixed by stripes, and mixed by patches) and 9 mixing proportions (i.e., from 10% to 90% with the interval of 10%), and to simulate red (R) and near-infrared (NIR) VHR images for these 3D scenes at both the solar principal plane (SPP) and perpendicular plane (PP) under different solar-viewing geometries. Negative correlations were generally found between the canopy BRF and the ratio of conifers in a mixed stand. The anisotropy of conifer dominated plantations is more prominent than broadleaf dominated plantations, especially for the single tree mixture. Although the level of anisotropy is much lower for PP than SPP, it should not be ignored, especially for the R band. Observations under large viewing zenith angles at PP are more preferred to study the effect of mixing proportions, followed by forward observations at SPP. The R band image has higher potential to distinguish mixture patterns for broadleaf-dominated situations, while the NIR band image has a higher potential for conifer-dominated situations. Furthermore, the canopy BRF generally increases with the solar zenith angle, and one meter can be considered as the optimal spatial resolution for the optical monitoring of the mixture mode. The findings of the current study add some valuable theoretical knowledge for the accurate monitoring of coniferous–broadleaved mixed plantations with VHR imagery.

Self-consistent modelling of the dust component in protoplanetary and circumplanetary disks: the case of PDS 70

Context.Direct observations of young stellar objects are important to test established theories of planet formation. PDS 70 is one of the few cases where robust evidence favours the presence of two planetary mass companions inside the gap of the transition disk. Those planets are believed to be going through the last stages of accretion from the protoplanetary disk, a process likely mediated by a circumplanetary disk (CPD).Aims.We aim to develop a three-dimensional radiative transfer model for the dust component of the PDS 70 system which reproduces the system’s global features observed at two different wavelengths: 855 μm with the Atacama large millimeter/submillimeter array (ALMA) and 1.25 μm with the Spectro-polarimetric high-contrast exoplanet research instrument at the Very large telescope (VLT/SPHERE). We use this model to investigate the physical properties of the planetary companion PDS 70 c and its potential circumplanetary disk.Methods.We base our modelling process on published ALMA and VLT/SPHERE observations of the dust continuum emission at 855 and 1.25 μm, respectively. We selected initial values for the physical properties of the planet and CPD through appropriate assumptions about the nature and evolutionary stage of the object. We modified the properties of the protoplanetary disk iteratively until the predictions retrieved from the model were consistent with both data sets. Simulations were carried out with the three-dimensional radiative transfer codeMCMax3D.Results.We provide a model that jointly explains the global features of the PDS 70 system seen in sub-millimetre and polarised-scattered light. Our model suggests that spatial segregation of dust grains is present in the protoplanetary disk. The sub-millimetre modelling of the PDS 70 c source favours the presence of an optically thick CPD and places an upper limit on its dust mass of 0.7M⊕. Furthermore, analysis of the thermal structure of the CPD demonstrates that the planet luminosity is the dominant heating mechanism of dust grains inside 0.6 au from the planet, while heating by stellar photons dominates at larger planetocentric distances.Conclusions.A CPD surrounding the planetary-mass companion PDS 70 c is a plausible scenario to explain the substructure observed with ALMA. The heating feedback from the protoplanetary disk has an non-negligible effect on the equilibrium temperature of dust grains in the outskirts of the CPD. The connection between the CPD properties and the planet mass still depends on a series of key assumptions. Further observations with high spatial and spectral resolution also for the gas component of the CPD are required to break the degeneracy between the properties of the planet and the disk.

Application of a Hypergeometric Model in Simulating Canopy Gap Fraction and BRF for Forest Plantations on Sloping Terrains

The influence of tree distribution and slope on canopy gap fraction (GF) and bidirectional reflectance factor (BRF) is shown here to be non-negligible. Trees are often assumed to be randomly distributed in natural forests due to random distribution of natural resources, but this assumption is not valid for forest plantations. A geometric optical model for forest plantations (GOFP) is a geometric optical model for forest plantations on horizontal surfaces based on the theory of exclusion distance among crowns. Sloping terrains change the exclusion distance among crowns, and inevitably affect the canopy GF and BRF. In this article, GOFP with a hypergeometric model (distances among trees are considered) on horizontal surfaces is modified as GOFP-T to simulate BRF for forest plantations on sloping terrains under two scenarios: (horizontal distances among crowns remain unchanged with slope) and (sloping distances among crowns remain unchanged with slope). Two three-dimensional (3-D) radiative transfer models (DART and LESS) and field measurements are used to evaluate and validate GOFP-T simulations. The results show that 1) the canopy GF, four component area ratios, and canopy BRF simulated by GOFP-T show high consistency with results from the two 3-D model: root-mean-square errors in GF, sunlit foliage, and sunlit ground are less than 0.02, 0.06, and 0.03, respectively; 2) forest coverage and canopy reflectance in GOFP-T are compared well with point cloud results from an airboard LiDAR system and Landsat 8 OLI surface reflectance products, respectively, indicating that GOFP-T has ability in simulating canopy reflectance for forest plantations on sloping terrains. GOFP-T with the hypergeometric model in this article is the first simple model to simulate canopy GF and BRF for forest plantations on sloping terrains.

A GPU-Based Solution for Ray Tracing 3-D Radiative Transfer Model for Optical and Thermal Images

Three-dimensional (3D) radiative transfer (RT) models are frequently recognized as a prerequisite when using high spatial resolution remote sensing data in heterogeneous surfaces. However, most studies of 3D RT models have been restricted to limited applications due to the low computational efficiency. Therefore, this study proposed a graphic processing unit (GPU)-based solution for the ray tracing 3D RT model. A state-of-the-art graphics and compute application programming interface, Vulkan, was introduced to implement the RT process. A bounding box method was adopted for the computation acceleration. By comparison with a central processing unit (CPU)-based solution, the performance efficiency of the proposed solution is significantly better: the simulation time of a GPU model is significantly reduced by more than 99% when facing a large-scale simulation mission. The simulation accuracy of the two solutions is similar, with root mean squared errors (RMSEs) lower than 0.005, 0.032 and 0.31 K for the red, near-infrared (NIR) and brightness temperature images, respectively. An evaluation based on airborne multiangle measurements also indicated that the accuracy of the proposed solution was satisfactory for simulating the red and NIR bidirectional reflectance factor and brightness temperature directional anisotropies, with RMSEs lower than 0.003, 0.020 and 0.20 K, respectively, when treating the whole scene as a pixel. Considering the simulation accuracy and efficiency, a GPU-based model will be an important supplement to the CPU model.

Three-dimensional radiative transfer modeling of forest: recent progress, applications, and future opportunities

Three-dimensional radiative transfer modeling of forest: recent progress, applications, and future opportunities

Estimating Leaf Area Index with Dynamic Leaf Optical Properties

Leaf area index (LAI) plays an important role in models of climate, hydrology, and ecosystem productivity. The physical model-based inversion method is a practical approach for large-scale LAI inversion. However, the ill-posed inversion problem, due to the limited constraint of inaccurate input parameters, is the dominant source of inversion errors. For instance, variables related to leaf optical properties are always set as constants or have large ranges, instead of the actual leaf reflectance of pixel vegetation in the current model-based inversions. This paper proposes to estimate LAI with the actual leaf optical property of pixels, calculated from the leaf chlorophyll content (Chlleaf) product, using a three-dimensional stochastic radiative transfer model (3D-RTM)-based, look-up table method. The parameter characterizing leaf optical properties in the 3D-RTM-based LAI inversion algorithm, single scattering albedo (SSA), is calculated with the Chlleaf product, instead of setting fixed values across a growing season. An algorithm to invert LAI with the dynamic SSA of the red band (SSAred) is proposed. The retrieval index (RI) increases from less than 42% to 100%, and the RMSE decreases to less than 0.28 in the simulations. The validation results show that the RMSE of the dynamic SSA decreases from 1.338 to 0.511, compared with the existing 3D-RTM-based LUT algorithm. The overestimation problem under high LAI conditions is reduced.

Characterizing reflectance anisotropy of background soil in open-canopy plantations using UAV-based multiangular images

The canopy bidirectional reflectance distribution function (BRDF) plays a pivotal role in estimating the biophysical parameters of plants, whereas soil background anisotropy creates challenges for their retrieval. Soil optical properties affect canopy anisotropic characteristics, especially in open-canopy areas. However, the remote sensing of background anisotropy is challenging due to the difficulties of information extraction in complex forest ecosystems and varying illumination conditions. This study develops an efficient photogrammetric technique to extract the background soil bidirectional reflectance factor (BRF) from unmanned aerial vehicle (UAV)-based multiangular images and to verify the need for accurate soil anisotropy information in canopy radiative transfer modeling. Soil BRF profiles were measured over three open-canopy sample plots from multiangular remotely sensed multispectral images collected with a hexacopter. As validation, reference soil BRF profiles were synchronously acquired by a ground-based multiangular imaging system.A high level of consistency between the ground- and UAV-measured soil BRF was observed with an RMSE of less than 0.012. Uncertainty analysis of the measured soil BRF showed that multiple scattering between sunlit soil in large sunflecks and foliage elements contributed less than 5%. Both results demonstrated that soil anisotropy can be accurately extracted from UAV multiangular measurements. To explicitly demonstrate that the use of soil anisotropy can reduce uncertainties in canopy radiative transfer simulations, we simulated the canopy BRF with Lambertian soil and with anisotropic soil using a three-dimensional (3D) radiative transfer model under different soil moisture content (SMC) levels, canopy cover (CC) levels and solar zenith angles (SZAs) with simulated realistic forest scenes. We found that less CC, lower SZAs and less SMC lead to a more significant influence of soil anisotropy on canopy reflectance; e.g., the reflectance bias reaches up to 0.3 in the hotspot direction. This illustrates that neglecting soil anisotropy can cause considerable errors in the modeling of the canopy BRF of open forests (i.e., CC levels of less than 0.5). The proposed technique facilitates the characterization of anisotropic forest background soil, which is important for advancing canopy radiative transfer modeling and validation and for the retrieval of vegetation parameters.

Investigation of the Geometric Shape Effect on the Solar Energy Potential of Gymnasium Buildings

Gymnasium are typically large-span buildings with abundant solar energy resources due to their extensive roof surface. However, relevant research on this topic has not been thoroughly conducted to investigate the effect of the geometric shape of gymnasium buildings on their solar potential. In this paper, an investigation of the geometric shape effect on the solar potential of gymnasium buildings is presented. A three-dimensional radiation transfer model coupled with historical meteorological data was established to estimate the real-time solar potential of the roof of a gymnasium building. The rooftop solar potential of three typical building foundation shapes and different types of roof shapes that have evolved was systematically analyzed. An annual solar potential cloud map of each gymnasium building is generated. The monthly and annual average solar radiation intensities of the different types of roof shapes are investigated. Compared to the optimal tilt angle, the maximum decrease in the average radiation intensity reached −20.42%, while the minimum decline was −8.64% for all types of building shapes. The solar energy potential fluctuated by up to 11% across the various roof shapes, which indicate that shape selection is of vital importance for integrated photovoltaic gymnasium buildings. The results presented in this work are essential for clarifying the effects of the geometric shape of gymnasium buildings on the solar potential of their roofs, which provide an important reference for building design.