Vertical Profiles Research Articles

Cadmium/zinc stresses and plant cultivation influenced soil microflora: a pot experiment conducted in field

It’s of great challenge to address for heavy metal-contaminated soil. Once the farmland is contaminated with heavy metals, the microbial ecology of the plant rhizosphere will change, which in turn impacts crop productivity and quality. However, few studies have explored the effects of heavy metals on plant rhizosphere microbes in farmland and the role that plant cultivation plays in such a phytoremediation practice. In this study, the impacts of comfrey (Symphytum officinale L.) cultivation and the stresses of cadmium/zinc (Cd/Zn) on rhizosphere soil microflora were examined. Microbial DNA was collected from soils to evaluate the prevalence of bacteria and fungi communities in rhizosphere soils. High-throughput 16 S rRNA sequencing was used to determine the diversity of the bacterial and fungal communities. The results showed that growing comfrey on polluted soils reduced the levels of Cd and Zn from the vertical profile. Both the comfrey growth and Cd/Zn stresses affected the community of rhizosphere microorganisms (bacteria or fungi). Additionally, the analysis of PCoA and NMDS indicated that the cultivation of comfrey significantly changed the bacterial composition and structure of unpolluted soil. Comfrey cultivation in polluted and unpolluted soils did not result in much variance in the fungi’s species composition, but the fungal compositions of the two-type soils were noticeably different. This work provided a better understanding of the impacts of Cd/Zn stresses and comfrey cultivation on rhizosphere microbial community, as well as new insight into phytoremediation of heavy metal-contaminated soils.

Large inter-individual variability in force-velocity profile changes in response to acute high-load resistance training.

While the acute effects of high-load resistance training on the force generating capacity of muscles have been widely examined, limited data exist on the relationship with the force-velocity profile (FV). Evidence suggests high sensitivity of the vertical FV profile to monitor changes in the muscle's mechanical properties according to the type of the exercise protocol. However, the interpretation of the findings seems not as straightforward. Therefore, the purpose of this study was to examine the effects of a high-load resistance training protocol on the muscle's mechanical properties during loaded jumps and on the vertical force-velocity profile (FV) in relation to maximal strength. 29 resistance-trained male (mean age±SD: 35.4 ± 7.8 years) and 29 female athletes (mean age±SD: 32.5 ± 7.0 years) participated in the study. Five-repetition maximum (5RM) in back squat, unloaded countermovement jump (CMJ) and FV profile were assessed. Loaded jumps were performed against 25, 50, 75, and 100 percent of body mass. Participants performed exercise protocols corresponding to their 5RM. Immediately after, unloaded CMJ and FV profile measurements were repeated. A significant decrease in CMJ height (∼5-6%) and in average power (∼4%) was recorded for both men and women. The FV profile did not change after the exercise protocol; however, there was a significant decrease in theoretical maximal power (from 4 to 5%) and in theoretical maximal velocity (∼3%). Maximal strength was not associated with the changes in FV profile. Findings suggest that an acute high-load exercise decreased vertical jump performance and maximal power output, but without a concomitant change in FV profile. The large interindividual variability in FV measures indicates a less straightforward connection of the applied exercise with the acute response in the FV profile, highlighting the complexity of the FV profile to monitor changes in response to an acute training load.

A self-consistent model for dust settling and the vertical shear instability in protoplanetary disks

Abstract The spatial distribution of dust particles in protoplanetary disks affects dust evolution and planetesimal formation processes. The vertical shear instability (VSI) is one of the candidate hydrodynamic mechanisms that can generate turbulence in the outer disk region and affect dust diffusion. Turbulence driven by the VSI has a predominant vertical motion that can prevent dust settling. On the other hand, the dust distribution controls the spatial distribution of the gas cooling rate, thereby affecting the strength of VSI-driven turbulence. Here, we present a semi-analytic model that determines the vertical dust distribution and the strength of VSI-driven turbulence in a self-consistent manner. The model uses an empirical formula for the vertical diffusion coefficient in VSI-driven turbulence obtained from our recent hydrodynamical simulations. The formula returns the vertical diffusion coefficient as a function of the vertical profile of the cooling rate, which is determined by the vertical dust distribution. We use this model to search for an equilibrium vertical dust profile where settling balances with turbulent diffusion for a given maximum grain size. We find that if the grains are sufficiently small, there exists a stable equilibrium dust distribution where VSI-driven turbulence is sustained at a level of αz ∼ 10−3, where αz is the dimensionless vertical diffusion coefficient. However, as the maximum grain size increases, the equilibrium solution vanishes because the VSI can no longer stop the settling of the grains. This runaway settling may explain highly settled dust rings found in the outer part of some protoplanetary disks.

The Recovery and Re-Calibration of a 13-Month Aerosol Extinction Profiles Dataset from Searchlight Observations from New Mexico, after the 1963 Agung Eruption

The recovery and re-calibration of a dataset of vertical aerosol extinction profiles of the 1963/64 stratospheric aerosol layer measured by a searchlight at 32°N in New Mexico, US, is reported. The recovered dataset consists of 105 aerosol extinction profiles at 550 nm that cover the period from December 1963 to December 1964. It is a unique record of the portion of the aerosol cloud from the March 1963 Agung volcanic eruption that was transported into the Northern Hemisphere subtropics. The data-recovery methodology involved re-digitizing the 105 original aerosol extinction profiles from individual Figures within a research report, followed by the re-calibration. It involves inverting the original equation used to compute the aerosol extinction profile to retrieve the corresponding normalized detector response profile. The re-calibration of the original aerosol extinction profiles used Rayleigh extinction profiles calculated from local soundings. Rayleigh and aerosol slant transmission corrections are applied using the MODTRAN code in transmission mode. Also, a best-estimate aerosol phase function was calculated from observations and applied to the entire column. The tropospheric aerosol phase function from an AERONET station in the vicinity of the searchlight location was applied between 2.76 to 11.7 km. The stratospheric phase function, applied for a 12.2 to 35.2 km altitude range, is calculated from particle-size distributions measured by a high-altitude aircraft in the vicinity of the searchlight in early 1964. The original error estimate was updated considering unaccounted errors. Both the re-calibrated aerosol extinction profiles and the re-calibrated stratospheric aerosol optical depth magnitudes showed higher magnitudes than the original aerosol extinction profiles and the original stratospheric aerosol optical depth, respectively. However, the magnitudes of the re-calibrated variables show a reasonable agreement with other contemporary observations. The re-calibrated stratospheric aerosol optical depth demonstrated its consistency with the tropics-to-pole decreasing trend, associated with the major volcanic eruption stratospheric aerosol pattern when compared to the time-coincident stratospheric aerosol optical depth lidar observations at Lexington at 42° N.

An automated geometry generation procedure for CWE applications in OpenFOAM: the case study of a thunderstorm downburst moving over Genoa city

AbstractSince a decade, Computational Fluid Dynamics (CFD) simulations for Wind Engineering (WE) on large-scale environments became a norm. Meso-$$\gamma$$ γ scales and microscale can be characterized by a large number of structures as well as complex terrain topographies. CFD simulations over such an area would commonly imply the involvement of a time-consuming geometry preparation in the pre-processing stage. This is also the innovative goal of the present paper for which an automated geometry generation procedure for Computational Wind Engineering (CWE) applications is developed in open-source environments. This paper further expands on the unsteady Reynolds-averaged Navier–Stokes (URANS) simulation of a thunderstorm downburst immersed in a background Atmospheric Boundary Layer (ABL) wind moving over Genoa city, Italy. Recommendations for geometry and grid generation of the present case study are also proposed. The study considers two simultaneous approaching wind conditions, a background ABL wind and a moving downburst (DB), the latter modeled as an impinging jet. The simulated results are compared with available full-scale LiDAR profiler data in terms of the vertical profiles of outflow velocity, showing relatively good agreement in terms of magnitude and profile shape.

Transpacific Transport of Asian Peroxyacetyl Nitrate (PAN) Observed from Satellite: Implications for Ozone.

Peroxyacetyl nitrate (PAN) is produced in the atmosphere by photochemical oxidation of non-methane volatile organic compounds in the presence of nitrogen oxides (NOx), and it can be transported over long distances at cold temperatures before decomposing thermally to release NOx in the remote troposphere. It is both a tracer and a precursor for transpacific ozone pollution transported from East Asia to North America. Here, we directly demonstrate this transport with PAN satellite observations from the infrared atmospheric sounding interferometer (IASI). We reprocess the IASI PAN retrievals by replacing the constant prior vertical profile with vertical shape factors from the GEOS-Chem model that capture the contrasting shapes observed from aircraft over South Korea (KORUS-AQ) and the North Pacific (ATom). The reprocessed IASI PAN observations show maximum transpacific transport of East Asian pollution in spring, with events over the Northeast Pacific offshore from the Western US associated in GEOS-Chem with elevated ozone in the lower free troposphere. However, these events increase surface ozone in the US by less than 1 ppbv because the East Asian pollution mainly remains offshore as it circulates the Pacific High.

A Lidar Biomass Index of Tidal Marshes from Drone Lidar Point Cloud

Accompanying climate change and sea level rise, tidal marsh mortality in coastal wetlands has been globally observed that urges the documentation of high-resolution, 3D marsh inventory to assist resilience planning. Drone Lidar has proven useful in extracting the fine-scale bare earth terrain and canopy height. Beyond that, this study performed marsh biomass mapping from drone Lidar point cloud in a S. alterniflora-dominated estuary on the Southeast U.S. coast. Three point classes (ground, low-veg, and high-veg) were classified via point cloud deep learning. Considering only vegetation points in the vertical profile, a profile area-weighted height (HPA) was extracted at a grid size of 50 cm × 50 cm. Vegetation point densities were also extracted at each grid. Adopting the plant-level allometric equations of stem biomass from long-term S. alterniflora surveys, a Lidar biomass index (Lidar_BI) was built to represent the relative quantity of marsh biomass in a range of [0, 1] across the estuary. Compared with the clipped dry biomass samples, it achieved a comparable and slightly better performance (R2 = 0.5) than the commonly applied spectral index approaches (R2 = 0.4) in the same marsh field. This study indicates the feasibility of the drone Lidar point cloud for marsh biomass mapping. More advantageously, the drone Lidar approach yields information on plant community architecture, such as canopy height and plant density distributions, which are key factors in evaluating marsh habitat and its ecological services.

Evolution of 3-D chlorophyll in the northwestern Pacific Ocean using a Gaussian-activation deep neural network model

Insufficient studies in characterizing vertical structure of Chlorophyll-a (Chl-a) in the ocean critically limit better understanding about marine ecosystem based on global climate change. In this study, we developed a Gaussian-activation deep neural network (Gaussian-DNN) model to assess vertical Chl-a structure in the upper ocean at high spatial resolution. Our Gaussian-DNN model used the input variables including satellite data of sea surface Chl-a and in-situ vertical physics profiles (temperature and salinity) in the northwestern Pacific Ocean (NWPO). After validation test based on two independent datasets of BGC-Argo and ship measurement, we applied the Gaussian-DNN model to reconstruct temporal evolution of 3-D Chl-a structure in the NWPO. Our modelling results successfully explain over 80% of the Chl-a vertical profiles in the NWPO at a horizontal resolution of 1° × 1° and 1 m vertical resolution within upper 300 meters during 2004 to 2022. Moreover, according to our modelling results, the Subsurface Chlorophyll Maxima (SCMs) and total Chl-a within 0-300 m depths were extracted and presented seasonal variability overlapping longer-time trends of spatial discrepancies all over the NWPO. In addition, our sensitivity testing suggested that sea-water temperatures predominantly control 3-D structures of the Chl-a in the tropical NWPO, while salinity played a key role in the temperate gyre of the NWPO. Here, our development of the Gaussian-DNN model may also be applied to craft long term, 3-D Chl-a products in the global ocean.

Initial evaluation of humidity profiles retrieved by the EOS-07 millimetre-wave humidity sounder

ABSTRACT The vertical profiles of specific humidity in the troposphere measured by the millimetre-wave humidity sounder (MHS) onboard of the Earth Observation Satellite (EOS)-07 (Microsoft 2B) are validated using ground-based radiosonde observations during July-August 2023. MHS onboard EOS-07 (hereafter referred as EOS-07/MHS) is a 6-channel cross-track scanning radiometer operating in the 183.31 ± 16.25 GHz band. The EOS-07 observations are used to retrieve three-dimensional specific humidity profiling from 1000 hPa (surface) to 100 hPa (16 km) atmospheric pressure levels with a spatial resolution of 10 km at nadir. The EOS-07/MHS measurements over India have a high correlation coefficient of 0.98 with the radiosonde observations. The mean relative bias is found to be 0.44 ± 0.55 g kg−1 below 350hPa and 2.8 ± 2.3 g kg−1 above 350 hPa. It is also found that above the 350 hPa level, MHS seems to be systematically overestimating the radiosonde measurements. After that, the EOS-07/MHS measurements are used to construct the diurnal variation of humidity during the active phases of the Indian summer monsoon. The observations are consistent with the present understanding of the Indian summer monsoon system. The present results are very encouraging and demonstrate the great potential of EOS-07/MHS observations of humidity for meteorological application, especially in understanding the hydrological cycle.

An Asymmetric Galactic Stellar Disk Traced by OB-type Stars from LAMOST DR7

Using 9943 OB-type stars from LAMOST DR7 in the solar neighborhood, we fit the vertical stellar density profile with the model including a single exponential distribution at different positions (R, Φ). The distributions of the scale heights and scale length show that the young disk traced by the OB-type stars is not axisymmetric. The scale length decreases versus the azimuthal angle Φ, i.e., from 1.14−0.19+0.30 kpc with Φ = − 3° to 0.77−0.06+0.08 kpc with Φ = 9°. Meanwhile we find signal of non-symmetry in the distribution of the scale height of the north and south of the disk plane. The scale height in the north side shows signal of flaring of the disk, while that of the south disk stays almost constant around h s = 130 pc. The distribution of the displacement of the disk plane Z 0 also shows variance versus the azimuthal angle Φ, which displays significant differences with the warp model constrained by the Cepheid stars. We also test different values for the position of the Sun, and the distance between the Sun and the Galactic center affects the scale heights and the displacement of the disk significantly, but that does not change our conclusion that the disk is not axisymmetric.

A Comparison of WaveNet and XGBoost for traditional direct wave propagation and seismic inversion using horizontal layer models

The application of machine learning in geophysics has steeply increased in the last decade, with the quality of its results varying according to the type of seismic problem in focus and the employed computational method. Deep Learning methods are achieving impressive results in this area, but we note that there is still a lack of certainty on whether classical machine learning methods can provide similar results. In the present paper, the objective was to attempt to fill part of that gap, by comparing a well-known non-DL machine-learning method with a DL method for the direct wave propagation and the seismic inversion problems for 2D horizontally-layered models. Both methods are evaluated in different scenarios, but under similar conditions, so that it is possible to understand the effect of parameter configuration on their final results. The dataset has 20,000 samples, each consisting of three vectors: a velocity vector with 236 values (representing a vertical profile of a randomly generated 2D layered model), a reflective vector with 600 values obtained directly from the velocity vector, and the associated seismogram vector with 11 traces, each containing 600 values. The overall results show that the WaveNet produces a lower Mean Squared Error between predicted and correct outputs than that of XGBoost. One challenge yet not dealt with is that the WaveNet can train well in GPU, but we did not succeed in doing the same with the XGBoost, due to the amount of data to be processed.

Estimation of the suspended sediment concentration in ice-covered channels based on the gravitational theory

Rivers in cold high latitude northern regions usually freeze during winter, thereby forming an ice cover on the water surface. The additional ice cover significantly changes the flow characteristics of the rivers and hence affects the process of water and sediment transport. Understanding the vertical profile of suspended sediment concentration (SSC) in ice-covered channel flows is helpful for investigating the sediment deposition and the bed morphology of rivers located in cold regions. To predict the SSC distribution in ice-covered channels, this study is based on the gravitational theory and puts forward a new SSC model by dividing the flow underneath the ice cover into the upper ice layer and the lower bed layer with the maximum streamwise velocity and the zero Reynolds shear stress as the interface. The calculated results of the new SSC model agree well with the collected experimental data, and the accuracy of the prediction results is evaluated by error analysis. Besides, the differences in the governing equations and the prediction accuracy between the gravitational theory and the diffusion theory are discussed, and thus the pros and cons of the proposed SSC model are illustrated. Although the gravitational theory has the controversial point, the validation results and discussions demonstrate that adopting the gravitational theory to estimate the SSC distribution in ice-covered channels is effective. The theoretical approach used herein can provide some guidance for studying the sediment-laden flow of freezing rivers in high-latitude northern regions.

A radiative–convective model computing precipitation with the maximum entropy production hypothesis

Abstract. All climate models use parameterizations and tuning in order to be accurate. The different parameterizations and tuning processes are the primary source of difference between models. Because models are tuned with present observations of Earth, they may not accurately simulate climates of other planets or palaeoclimate. A model with no adjustable parameter that happens to fit today's observations is probably more universal and should be more appropriate to model palaeoclimate. However, to our knowledge, such a model does not exist or is yet to be developed. This paper aims to improve a parameter-free radiative–convective model that computes a realistic temperature vertical profile to compute the water cycle, giving a value on average tropical precipitation. Although it is known that the radiative transfer constrains the order of magnitude of precipitation, no parameter-free model has yet been able to compute precipitation. Our model finds a precipitation value closer to observations than similar radiative–convective models or some general circulation models (GCMs).

Fine Shear-Wave Velocity Structures of Subsurface beneath the Guangdong–Hong Kong–Macao Greater Bay Area with Dense Seismic Array and SPAC Method

Abstract We apply the spatial autocorrelation (SPAC) method to construct the 3D subsurface shear-wave velocity structure model using the short-period dense seismic array (containing 725 nodal geophones) located at the Guangdong–Hong Kong–Macao Greater Bay area (GBA). We first divided the dense array into numerous subarrays, with each subarray consisting of nine nodal geophones, and obtained 562 subarrays that can provide 1D VS profiles of the same quantity. Then, the SPAC method and genetic algorithm are utilized to extract the dispersion curve of the Rayleigh wave from the raw microtremor data and invert VS structure, respectively. Finally, a 3D VS structure model from the surface to 3.3 km depth is derived by combining all 1D VS structures. Relatively low-velocity anomalies above 700 m are considered unconsolidated shallow sediments as well as relatively high-velocity anomalies beneath 1100 m are attributed to consolidated granite bedrock. Meanwhile, low-velocity anomalies that are identified through the vertical VS profile at a depth of about 900–3000 m can be contributed to the fractured zone, and striped low-velocity anomalies in the horizontal VS maps reveal the location of the deeply buried faults in the study area. The results also mean that the SPAC method combined with the records of short-period dense seismic array can be effectively applied to image subsurface structures in high-populated urban area. The development of this noise-resistance and environment-friendly geophysical technique provides a reliable and effective way to explore the complicated subsurface geological structures, which is of great significance to urban engineering construction and earthquake disaster reduction work in densely populated urban agglomerations.

Balloon drift estimation and improved position estimates for radiosondes

Abstract. When comparing model output with historical radiosonde observations, it is usually assumed that a radiosonde has risen exactly above its starting point and has not been displaced by wind. This changed only relatively recently with the availability of Global Navigation Satellite System (GNSS) receivers aboard radiosondes in the late 1990s, but even then the balloon trajectory data were often not transmitted, although this information was the basis for estimating the wind in the first place. Depending on the conditions and time of year, radiosondes can sometimes drift a few hundred kilometres, particularly at the middle latitudes during the winter months. The position errors can lead to non-negligible representation errors when the corresponding observations are assimilated. This paper presents a methodology to compute changes in the balloon position during its vertical ascent, using only limited information, such as the vertical profile of wind contained in the historical observation reports. The sensitivity of the method to various parameters is investigated, such as the vertical resolution of the input data, the assumption about the vertical ascent speed of the balloon, and the departure of the surface of Earth from a sphere. The paper considers modern GNSS sonde data reports for validation, for which the full trajectory of the balloon is available, alongside the reported wind. Evaluation is also conducted by comparison with ERA5 and by conducting low-resolution data assimilation experiments. Overall, the results indicate that the trajectory of the radiosondes can be accurately reconstructed from original data of varying vertical resolutions and that the more accurate balloon position reduces representation errors and, in some cases, systematic errors.

Investigating the impact of coupling HARMONIE-WINS50 (cy43) meteorology to LOTOS-EUROS (v2.2.002) on a simulation of NO2 concentrations over the Netherlands

Abstract. Meteorological fields calculated by numerical weather prediction (NWP) models drive offline chemical transport models (CTMs) to solve the transport, chemical reactions, and atmospheric interaction over the geographical domain of interest. HARMONIE (HIRLAM ALADIN Research on Mesoscale Operational NWP in Euromed) is a state-of-the-art non-hydrostatic NWP community model used at several European weather agencies to forecast weather at the local and/or regional scale. In this work, the HARMONIE WINS50 (cycle 43 cy43) reanalysis dataset at a resolution of 0.025° × 0.025° covering an area surrounding the North Sea for the years 2019–2021 was coupled offline to the LOTOS-EUROS (LOng-Term Ozone Simulation-EURopean Operational Smog model, v2.2.002) CTM. The impact of using either meteorological fields from HARMONIE or from ECMWF on LOTOS-EUROS simulations of NO2 has been evaluated against ground-level observations and TROPOMI tropospheric NO2 vertical columns. Furthermore, the difference between crucial meteorological input parameters such as the boundary layer height and the vertical diffusion coefficient between the hydrostatic ECMWF and non-hydrostatic HARMONIE data has been studied, and the vertical profiles of temperature, humidity, and wind are evaluated against meteorological observations at Cabauw in The Netherlands. The results of these first evaluations of the LOTOS-EUROS model performance in both configurations are used to investigate current uncertainties in air quality forecasting in relation to driving meteorological parameters and to assess the potential for improvements in forecasting pollution episodes at high resolutions based on the HARMONIE NWP model.

Observing low-altitude features in ozone concentrations in a shoreline environment via uncrewed aerial systems

Abstract. Ozone is a pollutant formed in the atmosphere by photochemical processes involving nitrogen oxides (NOx) and volatile organic compounds (VOCs) when exposed to sunlight. Tropospheric boundary layer ozone is regularly measured at ground stations and sampled infrequently through balloon, lidar, and crewed aircraft platforms, which have demonstrated characteristic patterns with altitude. Here, to better resolve vertical profiles of ozone within the atmospheric boundary layer, we developed and evaluated an uncrewed aircraft system (UAS) platform for measuring ozone and meteorological parameters of temperature, pressure, and humidity. To evaluate this approach, a UAS was flown with a portable ozone monitor and a meteorological temperature and humidity sensor to compare to tall tower measurements in northern Wisconsin. In June 2020, as a part of the WiscoDISCO20 campaign, a DJI M600 hexacopter UAS was flown with the same sensors to measure Lake Michigan shoreline ozone concentrations. This latter UAS experiment revealed a low-altitude structure in ozone concentrations in a shoreline environment showing the highest ozone at altitudes from 20–100 m a.g.l. These first such measurements of low-altitude ozone via a UAS in the Great Lakes region revealed a very shallow layer of ozone-rich air lying above the surface.

A Gaussian Function Model of Mesoscale Eddy Temperature Anomalies and Research of Spatial Distribution Characteristics

Mesoscale eddies are ubiquitous oceanic phenomena and play an important role in ocean circulation, ocean dynamics, and the transport of material energy. Temperature anomalies are a crucial parameter that reflects the state of mesoscale eddies. This study proposes a Gaussian function model to fit the vertical temperature anomaly (TA) profile to facilitate the analysis of variations, and the principle of the model is based on the fact that each TA profile tends to fluctuate around one or more peaks. The model is extracted and validated using Argo profiles within cyclonic and anticyclonic eddies in the Northwest Pacific Ocean spanning over the period from 2002 to 2021. The validation demonstrates that the model can accurately recover the vertical TA profiles with a limited number of parameters. This makes it suitable for analysing the spatial distribution patterns that require a large sample count. The analysis indicates that eddies with different TA profiles have a spatial aggregation effect in geographic distribution. Eddies with lower extreme temperature anomalies, at depths of 200–300 m, are mainly distributed along two bands on the north side of the Kuroshio Extension (KE) and the North Equatorial Current. Eddies with extreme TAs at the deepest depth (500–600 m) are distributed along the KE.

Ozone Monitoring Spectrometer-Limb observation (OMSL) on-orbit polarization correction for atmospheric radiation measurements

OMSL is a passive, space-based spectrometer designed to obtain the vertical profiles of trace gases such as ozone, by measuring atmospheric radiation in limb-viewing mode. Due to the presence of internal components such as gratings and dichroic mirror, OMSL is sensitive to polarization of incident light, which can lead to errors in the radiometric measurements, especially when the incident light is fully polarized. To eliminate the errors, an on-orbit polarization correction method is applied to OMSL. The key to this method is to calculate the polarization correction factor by combining the on-orbit polarization measurement and the pre-flight polarization calibration. In this paper, the on-orbit polarization correction method is introduced in detail. Furthermore, the first results of the on-orbit correction are analyzed and compared with simulation of the SCIATRAN.

Dissolved organic matter quality in thermokarst lake water and sediments across a permafrost gradient, Western Siberia

Thermokarst (thaw) lakes of permafrost peatlands are among the most important sentinels of climate change and sizable contributors of greenhouse gas emissions (GHG) in high latitudes. These lakes are humic, often acidic and exhibit fast growing/drainage depending on the local environmental and permafrost thaw. In contrast to good knowledge of the thermokarst lake water hydrochemistry and GHG fluxes, the sediments pore waters remain virtually unknown, despite the fact that these are hot spots of biogeochemical processes including GHG generation. Towards better understating of dissolved organic matter (DOM) quality at the lake water – sediment interface and in the sediments pore waters, here we studied concentration and optical (UV, visual) properties of DOM of 11 thermokarst lakes located in four permafrost zones of Western Siberia Lowland. We found systematic evaluation of DOM concentration, SUVA and various optical parameters along the vertical profile of lake sediments. The lake size and hence, the stage of lake development, had generally weak control on DOM quality.The permafrost zone exhibited clear impact on DOM porewater concentration, optical characteristics, aromaticity and weight average molecular weight (WAMW). The lowest quality of DOM, reflected in highest SUVA and WAMW, corresponding to the dominance of terrestrial sources, was observed at the southern boundary of the permafrost, in the sporadic/discontinuous zone. This suggests active mobilization of organic matter leachates from the interstitial peat and soil porewaters to the lake, presumably via subsurface or suprapermafrost influx. Applying a substitute space for time scenario for future evolution of OM characteristics in thermokarst lake sediments of Western Siberia, we foresee a decrease of DOM quality, molecular weight and potential bioavailability in lakes of continuous permafrost zone, and an increase in these parameters in the sporadic/discontinuous permafrost zone.