Zenith Research Articles

Evaluation of the SAIL Radiative Transfer Model for Simulating Canopy Reflectance of Row Crop Canopies

The widely used SAIL (Scattering by Arbitrarily Inclined Leaves) radiative transfer model (RTM) is designed for canopies that can be considered as homogeneous turbid media and thus should be inadequate for row canopies. However, numerous studies have employed the SAIL model for row crops (e.g., wheat and maize) to simulate canopy reflectance or retrieve vegetation properties with satisfactory accuracy. One crucial reason may be that under certain conditions, a row crop canopy can be considered as a turbid medium, fulfilling the assumption of the SAIL model. Yet, a comprehensive analysis about the performance of SAIL in row canopies under various conditions is currently absent. In this study, we employed field datasets of wheat canopies and synthetic datasets of wheat and maize canopies to explore the impacts of the vegetation cover fraction (fCover), solar angle and soil background on the performance of SAIL in row crops. In the numerical experiments, the LESS 3D RTM was used as a reference to evaluate the performance of SAIL for various scenarios. The results show that the fCover is the most significant factor, and the row canopy with a high fCover has a low soil background influence. For a non-black soil background, both the field measurement and simulation datasets showed that the SAIL model accuracy initially decreased, and then increased with an increasing fCover, with the most significant errors occurring when the fCover was between about 0.4 and 0.7. As for the solar angles, the accuracy of synthetic wheat canopy will be higher with a larger SZA (solar zenith angle), but that of a synthetic maize canopy is little affected by the SZA. The accuracy of the SAA (solar azimuth angle) in an across-row direction is always higher than that in an along-row direction. Additionally, when the SZA ranges from 65° to 75° and the fCover of wheat canopies are greater than 0.6, SAIL can simulate the canopy reflectance with satisfactory accuracy (rRMSE < 10%); the same accuracy can be achieved in maize canopies as long as the fCover is greater than 0.8. These findings provide insight into the applicability of SAIL in row crops and support the use of SAIL in row canopies under certain conditions (with rRMSE < 10%).

Evaluation of total ozone measurements from Geostationary Environmental Monitoring Spectrometer (GEMS)

Abstract. The continued interest in air pollution and stratospheric ozone variability has motivated the development of a Geostationary Environmental Monitoring Spectrometer (GEMS) for hourly ozone monitoring. This paper provides the atmospheric science community with the world's first assessment of GEMS total column ozone (TCO) retrieval performance and diurnal ozone variation. The algorithm used for GEMS is a more advanced version of its predecessor, the Total Ozone Mapping Spectrometer (TOMS) V8, that incorporates several improvements, including a new lookup table, a simple Lambertian-equivalent reflectivity model, and a spectral dependence correction. The GEMS algorithm also uses the optimal estimation method (OEM) to make error analysis more accessible and robust. The estimated retrieval errors range from 1.5 to 2 DU in September and 2 DU in December, with a constant degree of freedom of the signal (DFS) of 1 in September and a variable DFS of 1.25 to 1.4 in December throughout the day, depending on solar zenith angle (SZA). To assess the performance of the GEMS algorithm, the hourly GEMS total ozone was compared with ground-based measurements from Pandora instruments and other satellite platforms from TROPOMI (TROPOspheric Monitoring Instrument) and OMPS (Ozone Mapping and Profiler Suite Nadir Mapper). GEMS has a high correlation of 0.97 and small RMSE values compared to Pandora TCO at Busan and Seoul in South Korea. It is notable that despite exhibiting seasonal dependence in the mean bias of GEMS with Pandora, GEMS is capable of observing daily variations in ozone that are highly consistent with Pandora measurements, with a bias of approximately 1 %. The comparison of GEMS TCO data with TROPOMI and OMPS TCO data shows a high correlation of 0.99 and low RMSE compared to TROPOMI and OMPS TCO data, but the data have a negative bias of −2.38 % and −2.17 %, with standard deviations of 1.33 % and 1.57 %, respectively. Similar to OMPS, the influence of SO2 from volcanic eruptions is not properly removed in some regions, leading to GEMS overestimating TCO in those areas. The mean biases of GEMS TCO data with TROPOMI and OMPS TCO are within ±1 % at low latitudes but become negative at midlatitudes, with an increasingly negative dependence on latitude. Furthermore, this dependence becomes more prominent from summer to winter. The empirical correction applied to the GEMS irradiance data improves the dependence of the mean bias on season and latitude, but a consistent bias still remains, and a marginal positive trend was observed in December. Therefore, further investigation into correction methods is needed. The results are a meaningful scientific advance by providing the first validated, hourly UV ozone retrievals from a satellite in geostationary orbit. This experience can be used to advance research with future geostationary environmental satellite missions, including the incoming TEMPO (Tropospheric Emissions: Monitoring of Pollution) and Sentinel-4.

The Impact and Correction of Sensitive Environmental Factors on Spectral Reflectance Measured In Situ

The spectral reflectance measured in situ is often regarded as the “truth” of objects, which plays an important role in Earth observation applications. However, in situ measurements are influenced by several factors such as atmospheric conditions, illumination and view geometry (I&VG), cloud coverage, and adjacency effects. In order to avoid the influence of these factors, in situ measurements are usually carried out under sunny days and close to noon. However, the impact of I&VG is still present in most cases. At present, people still know little about the influence mechanism of I&VG. Moreover, correcting the impact of I&VG is also a problem that needs to be urgently solved in reflectance spectroscopy. In this work, experiments are carried out using the multi-directional hyperspectral remote sensing simulation facility (MHSRS2F), which allows adjustment and control of the I&VG parameters. This paper proposes an uncertainty evaluation model for I&VG and quantifies the uncertainty caused by different I&VG parameters. Then, the sensitivity of reflectance to I&VG at different wavelengths is explored based on uncertainty models. Finally, a correction model for reflectance under different I&VG conditions is proposed. The results reveal that the uncertainty and sensitivity caused by observation height are relatively high, regardless of the surface heterogeneity. It directly affects the size of the field of view and the physicochemical characteristics of the object. For objects that approximate the Lambertian surface, more attention should be paid to the selection and variation of solar and view zenith angles and view azimuth angles. For objects with surface heterogeneity, the selection and variation of solar azimuth angle, view azimuth angle, and solar zenith angle are more crucial. The correction model proposed in this paper has a 41.25% correction effect on different view zenith angles, but the correction effect on other environmental factors is not significant.

Enhancing Leaf Area Index Estimation for Maize with Tower-Based Multi-Angular Spectral Observations.

The leaf area index (LAI) played a crucial role in ecological, hydrological, and climate models. The normalized difference vegetation index (NDVI) has been a widely used tool for LAI estimation. However, the NDVI quickly saturates in dense vegetation and is susceptible to soil background interference in sparse vegetation. We proposed a multi-angular NDVI (MAVI) to enhance LAI estimation using tower-based multi-angular observations, aiming to minimize the interference of soil background and saturation effects. Our methodology involved collecting continuous tower-based multi-angular reflectance and the LAI over a three-year period in maize cropland. Then we proposed the MAVI based on an analysis of how canopy reflectance varies with solar zenith angle (SZA). Finally, we quantitatively evaluated the MAVI's performance in LAI retrieval by comparing it to eight other vegetation indices (VIs). Statistical tests revealed that the MAVI exhibited an improved curvilinear relationship with the LAI when the NDVI is corrected using multi-angular observations (R2 = 0.945, RMSE = 0.345, rRMSE = 0.147). Furthermore, the MAVI-based model effectively mitigated soil background effects in sparse vegetation (R2 = 0.934, RMSE = 0.155, rRMSE = 0.157). Our findings demonstrated the utility of tower-based multi-angular spectral observations in LAI retrieval, having the potential to provide continuous data for validating space-borne LAI products. This research significantly expanded the potential applications of multi-angular observations.

Sensitivity of cirrus and contrail radiative effect on cloud microphysical and environmental parameters

Abstract. Natural cirrus clouds and contrails cover about 30 % of the Earth's mid-latitudes and up to 70 % of the tropics. Due to their widespread occurrence, cirrus clouds have a considerable impact on the Earth energy budget, which, on average, leads to a warming net radiative effect (solar + thermal infrared). However, whether the instantaneous radiative effect (RE), which in some cases corresponds to a radiative forcing, of natural cirrus or contrails is positive or negative depends on their microphysical, macrophysical, and optical properties, as well as the radiative properties of the environment. This is further complicated by the fact that the actual ice crystal shape is often unknown, and thus, ice clouds remain one of the components that are least understood in the Earth's radiative budget. The present study aims to investigate the dependency of the effect on cirrus RE on eight parameters, namely solar zenith angle, ice water content, ice crystal effective radius, cirrus temperature, surface albedo, surface temperature, cloud optical thickness of an underlying liquid water cloud, and three ice crystal shapes. In total, 283 500 plane-parallel radiative transfer simulations have been performed, not including three-dimensional scattering effects. Parameter ranges are selected that are typically associated with natural cirrus and contrails. In addition, the effect of variations in the relative humidity profile and the ice cloud geometric thickness have been investigated for a sub-set of the simulations. The multi-dimensionality and complexity of the eight-dimensional parameter space makes it impractical to discuss all potential configurations in detail. Therefore, specific cases are selected and discussed. For a given parameter combination, the largest impact on solar, thermal-infrared (TIR), and net RE is related to the ice crystal effective radius. The second most important parameter is ice water content, which equally impacts the solar and terrestrial RE. The solar RE of cirrus is also determined by solar zenith angle, surface albedo, liquid cloud optical thickness, and ice crystal shape (in descending priority). RE in the TIR spectrum is dominated by surface temperature, ice cloud temperature, liquid water cloud optical thickness, and ice crystal shape. Net RE is controlled by surface albedo, solar zenith angle, and surface temperature in decreasing importance. The relative importance of the studied parameters differs, depending on the ambient conditions. Furthermore, and during nighttime the net RE is equal to the TIR RE. The data set generated in this work is publicly available. It can be used as a lookup table to extract the RE of cirrus clouds, contrails, and contrail cirrus instead of full radiative transfer calculations.

Martian M2 peak behavior in the dayside near-terminator ionosphere during interplanetary coronal mass ejections

The interplanetary coronal mass ejections (ICMEs) can pose significant impacts on the Martian ionosphere, resulting in plasma depletion, variability, and escape to space. However, the connections between the ICMEs and the associated responses of the dayside near-terminator Martian ionospheric primary peak (M2) are not well understood. The present study primarily investigates the behavior of the ionospheric peak density (Nm) and height (hm) during the passage of ICMEs using observations from the Radio Occultation Science Experiment (ROSE) aboard MAVEN spacecraft. We have selected 8 such ICMEs (during 2017–2022) at Mars from the existing catalogs and studied the ROSE electron density profiles during quiet and disturbed time (ICMEs) for identical solar zenith angle range. We observed the elevation of the M2 peak (hm ∼ 4–16 km) during disturbed time (ICMEs) with a decrease in Nm (0.41–2.8 × 1010 m−3) in comparison to the quiet time. The present study, for the first time, addressed the influence of ICMEs on the M2 peak parameters (Nm and hm). We have proposed that the development of large vertical pressure gradient and electron temperature enhancement are plausible causes for ionospheric variability. Therefore, the present study provides new insights to understand peak plasma behavior in the dayside near-terminator ionosphere during ICMEs.

Variations of the magnetic field and plasma parameters near the Martian induced magnetosphere boundary

The induced magnetosphere boundary (IMB) forms between the Martian magnetosheath and the induced magnetosphere/ionosphere and is important in energy and momentum transport and particle escape at Mars. Based on MAVEN multi-instrument observations in the year 2015, we find the IMBs usually appear with distinct variations of magnetic field and plasma, though they could be disturbed occasionally, especially at large solar zenith angles (SZAs). This study focuses on the properties of the distinct IMBs. The magnetic field strength is enhanced and the magnetic fluctuations are suppressed in IMB. The solar wind particles decrease and the planetary heavy ions start to be dominant in IMB. The variation of plasma is more effective in identifying the IMB's location since it is consistent at all SZAs, while the enhancement of the magnetic field strength or the suppression of the magnetic fluctuations would become ambiguous at large SZAs. The currents in IMB exhibit the signature of field-aligned-like and have larger magnitudes at the subsolar region. The kinetic effect of multiple ion species could be important near IMB. The properties of the Martian IMB investigated here are helpful in identifying this region and understanding the Mars-solar wind interaction and the Martian ionosphere/atm evolution.

A portable reflected-sunlight spectrometer for CO2 and CH4

Abstract. Mapping the greenhouse gases (GHGs) carbon dioxide (CO2) and methane (CH4) above source regions such as urban areas can deliver insights into the distribution and dynamics of local emission patterns. Here, we present the prototype development and an initial performance evaluation of a portable spectrometer that allows for measuring CO2 and CH4 concentrations integrated along a long (>10 km) horizontal path component through the atmospheric boundary layer above a target region. To this end, the spectrometer is positioned at an elevated site from which it points downward at reflection targets in the region, collecting the reflected sunlight at shallow viewing angles. The path-integrated CO2 and CH4 concentrations are inferred from the absorption fingerprint in the shortwave–infrared (SWIR) spectral range. While mimicking the concept of the stationary California Laboratory for Atmospheric Remote Sensing – Fourier Transform Spectrometer (CLARS-FTS) in Los Angeles, our portable setup requires minimal infrastructure and is straightforward to duplicate and to operate in various locations. For performance evaluation, we deployed the instrument, termed EM27/SCA, side by side with the CLARS-FTS at the Mt. Wilson Observatory (1670 m a.s.l.) above Los Angeles for a 1-month period in April/May 2022. We determined the relative precision of the retrieved slant column densities (SCDs) for urban reflection targets to be 0.36 %–0.55 % for O2, CO2 and CH4, where O2 is relevant for light path estimation. For the partial vertical column (VCD) below instrument level, which is the quantity carrying emission information, the propagated precision errors amount to 0.75 %–2 % for the three gases depending on the distance to the reflection target and solar zenith angle. The comparison to simultaneous CLARS-FTS measurements shows good consistency, but the observed diurnal patterns highlight the need to take light scattering into account to enable detection of emission patterns.

Combined modelling of annual and diurnal land surface temperature cycles

The land surface's thermal dynamics follows annual and diurnal cycles that are to a large extent controlled by solar geometry. Therefore, annual and diurnal variations of land surface temperature (LST) can be modelled with relatively simple functions controlled by a small number of parameters, typically from three to six. The parameter values of the models can be determined by fitting the respective functions to time series of LST observations. Commonly either annual or diurnal LST variations are modelled or they are modelled sequentially in a two-step process. Here, we combine an annual temperature cycle (ATC) model controlled by the solar zenith angle (ATCsza) with a four-parameter version of the diurnal temperature cycle (DTC) model GOT09: this yields a new annual-diurnal temperature cycle (ADTC) model that simultaneously describes the annual and diurnal surface temperature dynamics. The proposed ADTC model is controlled by physically meaningful parameters: annual minimum temperature, annual temperature amplitude, annual maximum of daily temperature amplitude, mean time of thermal noon and time lag of maximum temperature with respect to summer solstice. Thus, the entire annual-diurnal LST dynamics is described with only five parameters. The new model was tested by fitting it to one year of LST observations obtained for five globally representative tiles of the Moderate Resolution Imaging Spectroradiometer (MODIS), onboard EOS – TERRA and EOS – AQUA satellites. For these tiles, the mean of the root mean square error (RMSE) was 4.2 K. ADTC modelled LSTs were also compared against those obtained with the standard ATC model for the four MODIS overpass times at five representative sites: for these, an overall RMSE of 1.2 K between the two models was obtained. The ADTC derived LSTs were validated against in-situ measurements from three different sites, which yielded an overall RMSE of 3.4 K. Additional investigations over five areas with different land covers (i.e. urban, lake, forest, mountain area and desert) revealed the potential of the ADTC parameters to describe the corresponding surface and climate properties. Since it is driven by solar geometry, the ADTC model reproduces double LST peaks in the tropics naturally. Furthermore, all available observations are modelled simultaneously, which means that a single set of parameters is obtained for each pixel and year.

Complex determination of automatic robotic total stations' measurements' accuracy in underground spaces and comparison with results on the surface

The total station is the most basic geodetic measuring instrument, locally the most accurate and versatile. Its accuracy is the cornerstone of its use and is defined by the standard deviations of horizontal direction, zenith angle and slope distance measurements. These accuracy parameters are given by the manufacturer, but these are only valid under optimum measurement conditions. To ensure the credibility and reliability of the measurements, these values must be periodically ascertained or determined for atypical measurement configurations or measurement conditions. Standardised procedures are used for this purpose, but in our opinion, they do not reflect the full influence of the measuring conditions and other measuring aids. A comprehensive determination of the accuracies (variation components) from the alignment, where all possible influences in a given situation are applied, may be considered the most appropriate for determining the angular accuracy of measurements. Such atypical conditions are certainly represented by geodetic measurements in the confined spaces of an underground mine. Thus, an experimental determination of the accuracy of four different robotic total stations was carried out at the Center of experimental geotechnics in a mine Josef teaching centre (CTU in Prague), and the Fӧrstner method was used to determine the variation components. A network of 6 stations and 8 target points was designed. The grid size was approximately 32x21 m with 4 - 31 m lengths. A network with the same configuration was also duplicated at the surface to assess whether the accuracy is different in underground and how the results will correspond to the accuracy claimed by the manufacturers. The result of the testing is that the accuracy claimed by the manufacturers is maintained even under such difficult measurement conditions in narrow corridors and with short sights. The overall evaluation also found that the accuracy achieved underground and on the surface is identical, although it varies from instrument to instrument.

Impact of the cloud zenith angle of the pickup coil upon the sensitivity of a flexible rosette eddy current sensor

Bolt links and riveted joints, although widely used in aerospace and other mechanical systems, are susceptible to stress concentrations and fatigue cracks, which are undesirable. The flexible rosette eddy current sensor is typically used for the structural monitoring of bolt connections and riveted joints. This study proposes a method to improve the sensitivity of a flexible rosette eddy current sensor by reducing the cloud zenith angle of the pickup coil. First, the effect of the cloud zenith angle upon the sensitivity of the flexible rosette eddy current sensor is analyzed theoretically; the results show that a larger cloud zenith angle corresponds to a larger background signal. This diminishes the relative rate of change of the crack signal when the sensitivity of the sensor is evaluated based upon the maximum rate of change of the output. Therefore, the sensitivity is optimized by reducing the cloud zenith angle. Subsequently, through simulation and experimental methods, the relationship between sensor sensitivity and cloud zenith angle was studied. The results indicate that the cloud zenith angle significantly affects the sensitivity and that reducing the cloud zenith angle effectively reduces the background and improves the sensitivity.

3-D muographic inversion in the exploration of cavities and low-density fractured zones

SUMMARY Muography is an imaging tool based on the attenuation of cosmic muons for observing density anomalies associated with large objects, such as underground caves or fractured zones. Tomography based on muography measurements, that is, 3-D reconstruction of density distribution from 2-D muon flux maps, brings along special challenges. The detector field of view covering must be as balanced as possible, considering the muon flux drop at high zenith angles and the detector placement possibilities. The inversion from directional muon fluxes to a 3-D density map is usually underdetermined (more voxels than measurements). Therefore, the solution of the inversion can be unstable due to partial coverage. The instability can be solved by geologically relevant Bayesian constraints. However, the Bayesian principle results in parameter bias and artefacts. In this work, linearized (density-length based) inversion is applied by formulating the constraints associated with inversion to ensure the stability of parameter fitting. After testing the procedure on synthetic examples, an actual high-quality muography measurement data set from seven positions is used as input for the inversion. The resulting tomographic imaging provides details on the complicated internal structures of karstic fracture zone. The existence of low density zones in the imaged space was verified by samples from core drills, which consist of altered dolomite powder within the intact high density dolomite.

Simulating the 3D Distribution of Absorbed Shortwave Radiation in a Tree Crown: Comparison of Simplified and Monte Carlo Models

AbstractUrban trees are expected to alleviate severe thermal environments in summer. The absorbed shortwave radiation is a key variable to estimate the microclimatic effects of trees. In numerical simulation tools for the urban and built thermal environment, the difficulty lies in how to calculate the absorbed radiation in terms of accuracy and computational load. Here we investigated the difference between a radiative transfer model calculating scattering by the Monte Carlo method and a simpler model that accounts for scattering by modifying the optical depth. The calculations of absorbed radiation were performed under various conditions of leaf normal distribution and zenith angle of incident radiation, using tree models represented by voxels with different leaf area density values. For the photosynthetically active radiation, the relative difference between models was less than 10% for most cases. For the near‐infrared radiation, the difference was also small for a large canopy consisting of several trees. The difference was 30%–60% and 0%–20% for zenith angles of incident radiation of 0° and 60°, respectively, for small isolated trees with a crown height and width of less than 5 m, which was a significantly different condition from a uniform canopy. This large difference was related to the uniformity of the canopy and was explained by the transmittance from each voxel to the outside of the canopy. The present study enables us to evaluate the applicability of the simple model and correct the difference from the Monte Carlo model.

Experimental Investigation on Fragmentation Identification in Loose Slope Landslides by Infrared Emissivity Variability Features

Infrared radiation (IR) features that are influenced by infrared emissivity ε and physical temperature Td have been successfully applied to the early-warning of landslides. Although the infrared emissivity of a rock is a key parameter to determine its thermal radiation properties, the effect of particle size on the infrared emissivity of rock fragments is unknown. So in this paper, granite, marble, and sandstone were used as examples to conduct infrared imaging experiments on rock fragments. Their equivalent emissivity was used to interpret the detected infrared emission, including that from indoor backgrounds. In addition, the characteristics of changes in equivalent emissivity were discussed with reference to changes in observation direction and zenith angle. Then, a computation model of equivalent emissivity based on multiple observation directions and zenith angles was built to reveal the change in equivalent emissivity with particle sizes. The result indicates that the indoor background radiation has a predominant direction just above the rock fragments. The maximum deviation of infrared brightness temperature (IBT) was 0.260 K, and the maximum deviation of equivalent emissivity among different observation directions and zenith angles was 0.0065. After eliminating the influence of directional and angle effects with the operation of normalization, the general law of equivalent emissivity for all rock fragments that change with particle size is consistent. The maximum equivalent emissivity occurs at particle size 5 mm in the condition of particle size larger than 1 mm, while the equivalent emissivity changes inversely with particle size in the condition of particle size smaller than 1 mm. Above all, this study contributes new cognitions to Remote Sensing Rock Mechanics, and provides valuable evidence for better thermal infrared remote sensing monitoring on loose slope landslides.

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

Abstract. The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (Nd) comparisons are highly variable depending on the Nd retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and Nd values due to a lower cloud effective radius than observed. However, the simulated albedo response to Nd is suppressed due to a correlation between the solar zenith angle (SZA) and Nd created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln Nd distribution, which increases the cloud albedo sensitivity to ln Nd to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of Nd to CCN is greater than observed for TOA-retrieved Nd, and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved Nd. The LWP response to Nd in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–Nd relationship, observed clouds become much more adiabatic as Nd increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–Nd relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–Nd relationship non-causal. Nd retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

Average Intensity of a Laguerre—Gaussian Vector Vortex Beam through Inhomogeneous Atmospheric Turbulence

We investigate the evolution properties of a partially coherent Laguerre–Gaussian vector vortex (LGVV) beam through inhomogeneous atmospheric turbulence. Analytical formulae for the elements of a cross-spectral density matrix of a partially coherent LGVV beam propagating in turbulence are derived with the help of the extended Huygens–Fresnel principle. Our outcomes demonstrate that the normalized initial profile of a partially coherent beam with concentric dark rings gradually evolves into a Gaussian-like beam profile in turbulence. We also find that the beam is emitted at a large zenith angle and quickly converts to a Gaussian-like beam. Furthermore, it is also shown that a propagation beam with a large topological charge has a stronger ability to resist atmospheric turbulence. In order to confirm our numerical results, we combine the complex screen method and multi-phase screen method to simulate the propagation of a partially coherent LGVV beam in atmospheric turbulence. It is indicated that the simulation results are in good agreement with theoretical predictions. Our results will pave the way for the development of free-space optical communications and remote sensing.

Multi-longitudinal-mode Doppler lidar technology based on quadri-channel Mach–Zehnder interferometer

The detection principle of a multi-longitudinal-mode (MLM) Doppler lidar based on a quadri-channel Mach–Zehnder interferometer (QMZI) was analysed. The measurement error formulas for the radial wind speed and backscatter ratio of this type of Doppler lidar system were derived. The advantages of this detection technology were thoroughly analysed based on the principle. The optical path difference (OPD) ‘l’ between the two interference arms of the QMZI was optimised. For a detection height of 0–20 km, ‘l’ was selected as 50 cm to balance the detection of radial wind speed and backscatter ratio. The overall structure of the MLM Doppler lidar based on QMZI was designed, and detailed design parameters of the system were provided. Further simulation was conducted to evaluate the detection performance of the designed QMZI- based MLM Doppler lidar system. The simulation results showed that when using a pulse laser with a power of 17.5 W@355 nm and a telescope with an aperture of 25 cm, a detection zenith angle of 30°, a vertical resolution of 26 m@0–10 km and 52 m@10–20 km and a pulse accumulation time of 1 min, the radial wind speed measurement errors of the system during the day and night were less than 2.1 m/s and 2.5 m/s; the backscatter ratio relative measurement errors were less than 0.81% and 0.97%.

Assessing the Allometric Scaling of Vectorized Branch Lengths of Trees with Terrestrial Laser Scanning and Quantitative Structure Modeling: A Case Study in Guyana

Allometric scaling is closely related to the morphology, function and behavior of trees, which are of great significance to the study of ecology. However, most of the traditional allometric scaling studies used the scalar attributes of trees, without considering the 3D vector mode of tree growth. In order to investigate the allometric relationships between branch lengths in 3D vector mode, in this study, an accurate and detailed quantitative structure model was used to reconstruct tree architectures from 3D point cloud data collected by terrestrial laser scanning and extract the structural parameters of each branch (length, branching level and zenith angle). The standardized major axis was used to establish and analyze the scalar and vector allometric relationships between branch lengths. Our results show that at the same branching level and using the same allometric model, the scaling exponents between the lengths of branches and the lengths of cumulative child (descendant) branches (no matter whether the lengths are in scalar or vector form) are similar among trees, and there is no significant difference between the scaling exponents of most trees. And the scaling exponents between the lengths of the horizontal components of branches and the cumulative lengths of the horizontal components of the child (descendant) branches are much larger than those between the lengths of the vertical components of branches and the cumulative lengths of vertical components of the child (descendant) branches. At different branching levels, the scaling exponents between the lengths of branches and the cumulative lengths of descendant branches tend to decrease with the increase in the branching level. The allometric models in terms of the cumulative lengths of horizontal components of the child (descendant) branches and the allometric models in terms of the cumulative lengths of child (descendant) branches have similar model accuracy and scaling exponents. The study results of allometric relationships between tree branch lengths in 3D mode are of great importance for understanding the crown morphology and the branching rule, which is helpful to further understand the growth strategies and adaptation mechanisms of trees and explain the growth and development mechanisms of trees from a physiological perspective.

The Method of Multi-Angle Remote Sensing Observation Based on Unmanned Aerial Vehicles and the Validation of BRDF

The measurement of bidirectional reflectivity for ground-based objects is a highly intricate task, with significant limitations in the capabilities of both ground-based and satellite-based observations from multiple viewpoints. In recent years, unmanned aerial vehicles (UAVs) have emerged as a novel remote sensing method, offering convenience and cost-effectiveness while enabling multi-view observations. This study devised a polygonal flight path along the hemisphere to achieve bidirectional reflectance distribution function (BRDF) measurements for large zenith angles and all azimuth angles. By employing photogrammetry’s principle of aerial triangulation, accurate observation angles were restored, and the geometric structure of “sun-object-view” was constructed. Furthermore, three BRDF models (M_Walthall, RPV, RTLSR) were compared and evaluated at the UAV scale in terms of fitting quality, shape structure, and reflectance errors to assess their inversion performance. The results demonstrated that the RPV model exhibited superior inversion performance followed, by M_Walthall; however, RTLST performed comparatively poorly. Notably, the M_Walthall model excelled in capturing smooth terrain object characteristics while RPV proved applicable to various types of rough terrain objects with multi-scale applicability for both UAVs and satellites. These methods and findings are crucial for an extensive exploration into the bidirectional reflectivity properties of ground-based objects, and provide an essential technical procedure for studying various ground-based objects’ in-plane reflection properties.

Exploring sensitivity: Unveiling the impact of input parameters on Venus ionosphere [formula omitted] layer characteristics

This study investigates the impact of model input parameters on the characteristic features of the dayside V2 layer in the Venus ionosphere using an in-house developed one-dimensional photochemical model (1D-PCM) and observations from Venus Express radio science experiments (VeRa). 1D-PCM is simulated for different combinations of neutral density models (VTS3, a global empirical model of the Venus Thermosphere, and VenusGRAM, the Venus Global Reference Atmospheric Model), solar flux models (Solar 2000 (S2K), and observations by SEE (Solar EUV Experiment) onboard the NASA’s TIMED, the Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics, satellite). The hmV2 (height of the V2 layer plasma density) is found to be better reproduced in VTS3-SEE and VTS3-S2K pairs of combination, indicating that VTS3 is a more representative model for Venusian neutral density. The nmV2 (peak plasma density of the V2 layer) remains consistently higher when using SEE data for solar fluxes, but the model deviates significantly for Solar Zenith Angles (SZA) between 65∘ and 85∘, possibly due to differences in neutral density at high SZA where the temperature is lower than predicted by the neutral density model. Quantitative analysis reveals that errors in the 1D-PCM, with deviations of 1 to 2 km in hmV2 and 5%–10% difference in nmV2 between the model and observations, could be due to uncertainties in the neutral density and solar flux models. Comparison with the Transplanet model suggests that accounting for input uncertainties can improve the model’s ability to reproduce observations. This study provides insights into the influence of input parameters on the Venus ionosphere and highlights the importance of considering uncertainties in modeling studies.