Search

We discuss robust estimations for the variance of normally distributed random variables in the presence of interference. The robust estimators are based on either ranking or the geometric mean. For the interference models used, estimators based on the geometric mean outperform the rank-based ones in both mitigating the effect of interference and reducing the statistical error when there is no interference. One reason for this is that estimators using the geometric mean do not suffer from the “heavy tail” phenomenon like the rank-based estimators do. The ratio of the standard deviation over the mean of the power random variable is sensitive to interference. It can thus be used as a criterion to combine the sample mean with a robust estimator to form a hybrid estimator. We apply the estimators to the Arecibo incoherent scatter radar signals to determine the total power and Doppler velocities in the ionospheric E-region altitudes. Although all the robust estimators selected deal with light contamination well, the hybrid estimator is most effective in all circumstances. It performs well in suppressing heavy contamination and is as efficient as the sample mean in reducing the statistical error. Accurate incoherent scatter radar measurements, especially at nighttime and at E-region altitudes, can improve studies of ionospheric dynamics and composition.

High-resolution remote sensing images can distinguish objects of smaller size, so as to more clearly express the texture features and structural information of objects, and provide a data source for the development of large-scale mapping, high-precision stereo measurement and other fields. The purpose of this paper is to complete the height estimation of the buildings by analyzing the stereoscopic observation formed by the front and rear view images of the Gaofen-7 line array CCD. After determining the roof profile of the building on the rear-view image, assuming a series of object elevations of the building, that is, searching for elevations with a certain step distance within a certain elevation search range, adopt the object-based image matching VLL algorithm, Through the RFM imaging model of the Gaofen-7 sensor, the rear-view contour is projected to the front-view image, and then the PSNR is selected as the similarity measure of the window, and the similarity of the image block area corresponding to the front- and rear-view contour is calculated. Corresponds to the hypothetical object elevation as the estimated height of the building. Under the technical route of this paper, the height of buildings on high-resolution images can be estimated to a level within 3 meters of accuracy.

Indole (ind) is a nitrogen-containing heterocyclic volatile organic compound commonly emitted from animal husbandry and from different plants like maize with global emissions of 0.1 Tg yr−1. The chemical composition and optical properties of indole secondary organic aerosol (SOA) and brown carbon (BrC) are still not well understood. To address this, environmental chamber experiments were conducted to investigate the oxidation of indole at atmospherically relevant concentrations of selected oxidants (OH radicals and O3) with or without NO2. In the presence of NO2, the SOA yields decreased by more than a factor of 2, but the mass absorption coefficient at 365 nm (MAC365) of ind-SOA was 4.3 ± 0.4 m2 g−1, which was 5 times higher than that in experiments without NO2. In the presence of NO2, C8H6N2O2 (identified as 3-nitroindole) contributed 76 % to all organic compounds detected by a chemical ionization mass spectrometer, contributing ∼ 50 % of the light absorption at 365 nm (Abs365). In the absence of NO2, the dominating chromophore was C8H7O3N, contributing to 20 %–30 % of Abs365. Indole contributes substantially to the formation of secondary BrC and its potential impact on the atmospheric radiative transfer is further enhanced in the presence of NO2, as it significantly increases the specific light absorption of ind-SOA by facilitating the formation of 3-nitroindole. This work provides new insights into an important process of brown carbon formation by interaction of two pollutants, NO2 and indole, mainly emitted by anthropogenic activities.

Stabilized Criegee intermediates (SCIs) and organic peroxy radicals (RO2) are critical in atmospheric oxidation processes and secondary organic aerosol (SOA) formation. However, the influence of temperature on their corresponding reaction mechanisms in SOA formation is unclear. Through utilizing formic acid as a SCI scavenger and regulating the ratio of hydroperoxyl radials (HO2) to RO2 ([HO2]/[RO2]) from ∼ 0.3 to ∼ 1.9 using different concentrations of CO, the roles of RO2 and SCIs in SOA formation were investigated from 248 to 298 K, particularly for dimer formation in β-pinene ozonolysis. The SOA yield increased by 21 % from 298 to 273 K, while it decreased by 40 % from 273 to 248 K. Both changing [HO2]/[RO2] and scavenging SCIs significantly affect SOA yield and composition. SCI reactions accounted for more than 40 % of dimer and SOA mass formation for all temperatures. Increasing [HO2]/[RO2] inhibited dimer and SOA formation, and this inhibition became larger with decreasing temperature. Compared to low [HO2]/[RO2] (0.30–0.34), the dimer abundance at high [HO2]/[RO2] (1.53–1.88) decreased by about 31 % at 298 K and 70 % at 248 K. [HO2]/[RO2] has a specific impact on SCI-controlled dimers at lower temperatures by especially influencing the C9–SCI reactions with RO2. The dimers formed from C9–SCI reactions with RO2 were estimated to decrease by 61 % at high [HO2]/[RO2] compared to low [HO2]/[RO2] at 248 K. The high reactivity and substantial contribution to SOA of β-pinene-derived SCIs at lower temperatures observed in this study suggest that monoterpene-derived SCI reactions should be accounted for in describing colder regions of the atmosphere.

In order to study the unsteady flow and vortex characteristics of tongue cavitation, numerical calculation is carried out for the whole flow channel given different conditions. Then, the calculation results are verified experimentally. The results show that after the occurrence of tongue cavitation in the centrifugal pump, it deteriorates with the decrease of NPSHa. However, when NPSHa is reduced to 3.78 m, it does not change significantly anymore with the decrease of NPSHa. The extrusion of fluid by the vapor at the tongue promotes the formation of the separation vortex, and the re-jet flow caused by the separation vortex leads to vapor shedding. The frequency of cavitation shedding is consistent with the frequency of vortex shedding. In the vorticity transport equation, the relative vortex stretching term and the relative vortex dilatation term dominate the vapor shedding by controlling the change in vorticity. The baroclinic torque term mainly affects the change of vorticity at the vapor-liquid interface, but to a much lesser extent than the first two terms.

Brown carbon (BrC) aerosol is light-absorbing organic carbon that affects radiative forcing and atmospheric photochemistry. The BrC chromophoric composition and its linkage to optical properties at the molecular level, however, are still not well characterized. In this study, we investigate the day–night differences in the chromophoric composition (38 species) and optical properties of water-soluble and water-insoluble BrC fractions (WS-BrC and WIS-BrC) in aerosol samples collected in Shijiazhuang, one of the most polluted cities in China. We found that the light absorption contribution of WS-BrC to total BrC at 365 nm was higher during the day (62±8 %) than during the night (47±26 %), which is in line with the difference in chromophoric polarity between daytime (more polar nitrated aromatics) and nighttime (more less-polar polycyclic aromatic hydrocarbons, PAHs). The high polarity and water solubility of BrC in the daytime suggests the enhanced contribution of secondary formation to BrC during the day. There was a decrease in the mass absorption efficiency of BrC from nighttime to daytime (2.88±0.24 vs. 2.58±0.14 for WS-BrC and 1.43±0.83 vs. 1.02±0.49 m2 g C−1 for WIS-BrC, respectively). Large polycyclic aromatic hydrocarbons (PAHs) with four- to six-ring PAHs and nitrophenols contributed to 76.7 % of the total light absorption between 300–420 nm at nighttime, while nitrocatechols and two- to three-ring oxygenated PAHs accounted for 52.6 % of the total light absorption during the day. The total mass concentrations of the identified chromophores showed larger day–night difference during the low-pollution period (day-to-night ratio of 4.3) than during the high-pollution period (day-to-night ratio of 1.8). The large day–night difference in BrC composition and absorption, therefore, should be considered when estimating the sources, atmospheric processes, and impacts of BrC.

The burning of incense for home use is a widespread practice that has been shown to have significant negative impacts on human health and air quality. However, there is a lack of understanding regarding its emission profiles and associated health risks. To address this knowledge gap, we utilized a state-of-the-art thermal-desorption comprehensive two-dimensional gas chromatography–mass spectrometer (TD-GC × GC-MS) to (semi-)quantify the emission factors (EFs) of 317 volatile compounds and thoroughly investigate the organic profiles of smoke from incense burning across a full-volatility range. Results showed that toluene (70.8±35.7 µg g−1) is the most abundant compound in smoke from incensing burning, followed by benzene, furfural, and phenol. Phenol, toluene, furfural, 2-furanmethanol, benzene, and benzyl alcohol are the main contributors to ozone and secondary organic aerosol (SOA) estimation. Intermediate volatility organic compounds (IVOCs) accounted for 19.2 % of the total EFs but 40.0 % of the estimated SOA. Additionally, a novel pixel-based method, combined with aroma analysis, revealed that furfural can act as a key tracer of incense burning and is responsible for the distinctive aroma of incense smoke. High-bioaccumulation-potential (BAP) assessment using pixel-based partition coefficient estimation revealed that acenaphthylene, dibenzofuran, and phthalate esters (PAEs) are chemicals of high-risk concern and warrant further control. Our results highlight the critical importance of investigating home-use incense burning and provide new insights into the health impacts of smoke from incense burning using novel approaches.

Introduction: Luspatercept reduced the severity of anemia in transfusion-dependent lower risk MDS with ring sideroblasts (RS) in the MEDALIST trial. Here we present real-world data in a prospective multicenter registry study of MDS patients treated with luspatercept (ChiCTR2300071857). Patients who were not enrolled in the MEDALIST trial were also evaluated, including those of higher IPSS-R risk, without ring sideroblasts, or failed to prior hypomethylating agent (HMA) therapy. Methods: Eligible patients were ≥ 18 years old and were RBC transfusion-dependent (defined as ≥2 units of RBC transfusion per 8 weeks). Primary endpoint was RBC transfusion independence (RBC-TI) for ≥8 weeks. Secondary endpoints are hematological improvement-erythroid (HI-E, defined as >50% reduction in RBC transfusion), neutrophil (HI-N) and platelet (HI-P) per the IWG 2006 criteria. Results: From June 2022 to August 2023, 25 patients were treated with luspatercept for a median of 7 cycles (Table 1). The median follow-up time was 19.4 weeks. The rate of RBC-TI was 36%. The median duration of RBC-TI was 227 [IQR: 129-324] days. The rates of HI-E, HI-N, HI-P were 52% (13/25), 21% (3/14) and 33% (5/15), respectively. Multi-variant analysis revealed that lower transfusion burden (

Internal solitary waves (ISWs) are a ubiquitous phenomenon in the dynamic ocean system, which play a crucial role in driving transport through turbulent mixing. Over the past few decades, numerical modelling has become a vital approach to investigate the generation mechanism and spatial distribution of ISWs. The northern South China Sea (NSCS) has been treated as a physical oceanographic focus of ISWs in massive numerical studies since the last century. However, there has been no systematic evaluation of a reliable three-dimensional (3D) model about accurately reproducing ISW characteristics in the NSCS. In this study, we implement a 3D ISW forecasting model in the NSCS and quantitatively evaluate the requirements of factors (i.e. model resolution, tidal forcing, and stratification selection) in accurately depicting ISW properties by comparison with observational data at a mooring station in the vicinity of the Dongsha Atoll. Firstly, the 500 m resolution model can basically reproduce the principal ISW characteristics, while the 250 m resolution model would be a better solution to identify wave properties, specifically increasing 40 % accuracy of predicting characteristic half-widths. Nonetheless, a 250 m resolution model spends nearly 5-fold the computational resources of a 500 m resolution model in the same model domain. Compared with the former two, the model with a lower resolution of 1000 m severely underestimates the nonlinearity of ISWs, resulting in an incorrect ISW field in the NSCS. Secondly, the model with 8 (or 13) primary tidal constituents can accurately reproduce the real ISW field in the NSCS, while the one with four main harmonics (M2, S2, K1 and O1) would underestimate averaged wave-induced velocity for about 38 % and averaged mode-1 wave amplitude for about 15 %. Thirdly, the model with the initial condition of field-extracted stratification gives a better performance in predicting some wave properties than the model with climatological stratification, namely 13 % improvement of arrival time and 46 % improvement of characteristic half-width. Finally, background currents, spatially varying stratification and external (wind) forcing are discussed to reproduce a more realistic ISW field in the future numerical simulations.

Innovation leads the development of technology and society, so the position of innovation and entrepreneurship training to nurture innovative talents is very important. The paper explores the ‘trinity’ innovation and entrepreneurship talent training model of ‘promoting learning through competition, education through competition, and innovation through competition’, focus on high level innovation and entrepreneurship competitions, enhance the stickiness among colleges, educators, students, and social enterprises, and enable students to participate in extracurricular academic activities that are highly consistent with the needs of the industrial market. To achieve this goal, School of Remote Sensing and Information Engineering of Wuhan University puts forward a ‘remote sensing +’ training model. Through the two-year practice of this model, 136 students from different schools and majors participated in the training of the ‘Remote Sensing +’ Innovation and Entrepreneurship Center in Wuhan University, and won three gold awards in the China's largest innovation and entrepreneurship competition. Due to the characteristics of cross-integration of remote sensing science and technology itself, it has played a supporting role in innovation education, integrating other multi-specialized knowledge, and providing more entrepreneurial and employment opportunities for students of related majors.

Sugars and amino acids are major classes of organic components in atmospheric fine particles and play important roles in atmospheric processes. However, the identification of their sources in different regions is explored little. To characterize local primary sources (biomass burning, plant, and soil sources) and evaluate their contributions to the total sugar compound and amino acid (AA) pool in different regions, fine particulate matter samples were collected from the urban, rural, and forest areas in Nanchang, China. The concentrations and compositions of sugar compounds (anhydrosugars, primary sugars, and sugar alcohols), free amino acids (FAAs), and combined amino acids (CAAs) were analysed by gas chromatography–mass spectrometry (GC-MS) after silylation derivatization. Urban areas had significantly higher average Σ sugar concentration (317±139 ng m−3) than that of the rural (181±72 ng m−3) and forest (275±154 ng m−3) areas (p<0.05). Overall, the distribution pattern of sugar compounds and CAAs in PM2.5 was generally similar in three areas. Levoglucosan accounted for 24.4 %, 22.0 %, and 21.7 %, respectively, of the total sugar pool in the urban, rural, and forest areas. This suggests that plant and soil sources, as well as biomass burning (BB), provide important contributions to aerosol sugars and CAAs in three areas. In the urban area, the concentrations of anhydrosugars showed a positive correlation with combined Gly concentrations, but no correlation was found between these two compounds in the rural and forest areas, indicating that the urban area is mainly affected by local combustion sources. This conclusion was also supported by the positive correlation between levoglucosan and non-sea-salt potassium, only observed in the urban area (Lev=0.07 K++37.7,r=0.6,p<0.05). In addition, the average levoglucosan / mannosan (L/M) ratio in the urban area (59.9) was much higher than in the rural (6.9) and forest areas (7.2), implying BB aerosols collected in the urban area originated from lignite burning, while the type of biofuels used in the rural and forest areas is mainly softwood. The concentrations of sugar alcohols in the rural and forest areas were positively correlated with that of CAAs, which are abundant in the topsoil (r=0.53∼0.62,p<0.05), suggesting that the contribution of local topsoil sources is large in these two areas. In the rural and forest areas, the concentrations of primary sugars were positively correlated with those of combined CAA species abundant in local dominant vegetation. Our findings suggest that combining specific sugar tracers and chemical profiles of CAAs in local emission sources can provide insight into primary source characteristics, including the types of biofuels burned, the contribution of topsoil sources, and local vegetation types.

Abstract In clinical medicine, the segmentation of blood vessels in retinal images is essential for subsequent analysis in clinical diagnosis. However, retinal images are often noisy and their vascular structure is relatively tiny, which poses significant challenges for vessel segmentation. To improve the performance of vessel segmentation, an improved model IterNet++ based on the architecture of IterNet is proposed. First, curvelet signal analysis is applied to enhance retinal images. Second, residual convolution (ResConv) blocks and guided filters are introduced to utilise the encoder features of previous iterations in the model to reduce overfitting. Third, offline hard‐sample mining is used to improve segmentation performance by utilising training samples with low segmentation accuracy as many possible on a few‐sample training set. In addition, a test‐time augmentation method is applied to testing samples in test dataset during inference. Extensive experiments show that this model achieves Dice scores of 0.8313, 0.8277, and 0.8372 on DRIVE, CHASE‐DB1, and STARE datasets, respectively, demonstrating the best performance compared with IterNet and other baseline models.

Enhancements of stationary planetary waves (SPWs) and traveling planetary waves (TPWs) are commonly observed in the middle atmosphere during sudden stratospheric warming (SSW) events. Based on the least squares fitting method (Wu et al., 1995), numerous studies have used satellite measurements to investigate the characteristics of TPWs during SSWs, but they have ignored the effect of the SPWs. However, a rapid and large change in the SPWs during SSWs may lead to significant disturbances in the amplitude of derived TPWs. In this study, we present a new methodology for obtaining the amplitudes and wave numbers of traveling quasi-5-day oscillations (Q5DOs) in the middle atmosphere during major SSWs. Our new fitting method is developed by inhibiting the effect of a rapid and large change in SPWs during SSWs. We demonstrate the effectiveness of the new method using both synthetic data and satellite observations. The results of the simulations indicate that the new method can suppress the aliasing from SPWs and capture the real variations in TPWs during SSWs. Based on the geopotential height data measured by the Aura satellite from 2004 to 2021, the variations in traveling Q5DOs during eight midwinter major SSWs are reevaluated using the new method. The differences in the fitted amplitudes between the least squares fitting method and the new method are usually over 100 m during the SSW onsets. Our analysis indicates that previously reported Q5DOs during SSWs might be contaminated by SPWs, leading to both overestimation and underestimation of the amplitudes of the traveling Q5DOs.

To elucidate the molecular chemical compositions, volatility–polarity distributions, and influencing factors of Chinese cooking emissions, a comprehensive cooking emission experiment was conducted. Volatile organic compounds (VOCs), intermediate volatility, and semi-volatile organic compounds (I/SVOCs) from cooking fumes were analysed by a thermal desorption comprehensive two-dimensional gas chromatography coupled with quadrupole mass spectrometer (TD-GC × GC-qMS). Emissions from four typical Chinese dishes, i.e. fried chicken, Kung Pao chicken, pan-fried tofu, and stir-fried cabbage were investigated to illustrate the impact of cooking style and material. Fumes of chicken fried with corn, peanut, soybean, and sunflower oils were investigated to demonstrate the influence of cooking oil. A total of 201 chemicals were quantified. Kung Pao chicken emitted more pollutants than other dishes due to its rather intense cooking method. Aromatics and oxygenated compounds were extensively detected among meat-related cooking fumes, while a vegetable-related profile was observed in the emissions of stir-fried cabbage. Ozone formation potential (OFP) was dominated by chemicals in the VOC range. Of the secondary organic aerosol (SOA) estimation, 10.2 %–32.0 % could be explained by S/IVOCs. Pixel-based partial least squares discriminant analysis (PLS-DA) and multiway principal component analysis (MPCA) were utilized for sample classification and component identification. The results indicated that the oil factor explained more variance of chemical compositions than the cooking style factor. MPCA results emphasize the importance of the unsaturated fatty acid-alkadienal-volatile products mechanism (oil autoxidation) accelerated by the cooking and heating procedure.

In this paper, a regional study of mode-2 internal solitary waves (ISWs) at the Pacific coast of Central America is carried out using the seismic reflection method. The observed relationship between the dimensionless propagation speed and the dimensionless amplitude (DA) of the mode-2 ISW is analyzed. When DA < 1.18, the dimensionless propagation speed seems to increase with increasing DA, divided into two parts with different growth rates. When DA > 1.18, the dimensionless propagation speed increases with increasing DA at a relatively small growth rate. We suggest that the influences of seawater depth (submarine topography), pycnocline depth, and pycnocline thickness on the propagation speed of the mode-2 ISW in the study area cause the relationship between dimensionless propagation speed and DA to diversify. The observed relationship between the dimensionless wavelength and the DA of the mode-2 ISW is also analyzed. When DA < 1, the nondimensional wavelengths seem to change from 2.5 to 7 for a fixed nondimensional amplitude. When DA > 1.87, the dimensionless wavelength increases with increasing DA. Additionally, the seawater depth has a great influence on the wavelength of the mode-2 ISW in the study area. Overall, the wavelength increases with increasing seawater depth. As for the vertical structure of the amplitude of the mode-2 ISW in the study area, we find that it is affected by the nonlinearity of the ISW and the pycnocline deviation (especially the downward pycnocline deviation).

Using radiosonde observations at five stations in the tropical western Pacific and reanalysis data for the 15 years from 2005 to 2019, we report an extremely negative anomaly in atmospheric water vapor during the super El Niño winter of 2015/16 and compare the anomaly with that in the other three El Niño winters of the period. A strong specific humidity anomaly is concentrated below 8 km of the troposphere with a peak at 2.5–3.5 km, and a column-integrated water vapor mass anomaly over the five radiosonde sites has a large negative correlation coefficient of −0.63 with the oceanic Niño3.4 index but with a lag of about 2–3 months. In general, the tropical circulation anomaly in the El Niño winter is characterized by divergence (convergence) in the lower troposphere over the tropical western (eastern) Pacific; thus, the water vapor decreases over the tropical western Pacific as upward motion is suppressed. The variability of the Hadley circulation is quite small and has little influence on the observed water vapor anomaly. The anomaly of the Walker circulation makes a considerable contribution to the total anomaly in all four El Niño winters, especially in the 2006/07 and 2015/16 eastern Pacific (EP) El Niño events. The monsoon circulation shows a remarkable change from one event to another, and its anomaly is large in the 2009/10 and 2018/19 central Pacific (CP) El Niño winters and small in the two EP El Niño winters. The observed water vapor anomaly is caused mainly by the Walker circulation anomaly in the super EP event of 2015/16 but is caused by the monsoon circulation anomaly in the strong CP event of 2009/10. The roles of the Hadley, Walker, and monsoon circulations in the EP and CP events are confirmed by the composite EP and CP El Niños based on the reanalysis data for 41 years. Owing to the anomalous decrease in upward transport of water vapor during the El Niño winter, lower cloud amounts and more outgoing longwave radiation over the five stations are clearly presented in satellite observation. In addition, a detailed comparison of water vapor in the reanalysis, radiosonde, and satellite data shows a fine confidence level for the datasets; nevertheless, the reanalysis seems to slightly underestimate the water vapor over the five stations in the 2009/10 winter.

Shoaling internal solitary waves near the Dongsha Atoll in the South China Sea dissipate their energy and enhance diapycnal mixing, which have an important impact on the oceanic environment and primary productivity. The enhanced diapycnal mixing is patchy and instantaneous. Evaluating its spatiotemporal distribution requires comprehensive observation data. Fortunately, seismic oceanography meets the requirements, thanks to its high spatial resolution and large spatial coverage. In this paper, we studied three internal solitary waves in reversing polarity near the Dongsha Atoll and calculated their spatial distribution of diapycnal diffusivity. Our results show that the average diffusivities along three survey lines are 2 orders of magnitude larger than the open-ocean value. The average diffusivity in internal solitary waves with reversing polarity is 3 times that of the non-polarity reversal region. The diapycnal diffusivity is higher at the front of one internal solitary wave and gradually decreases from shallow to deep water in the vertical direction. Our results also indicate that (1) the enhanced diapycnal diffusivity is related to reflection seismic events, (2) convective instability and shear instability may both contribute to the enhanced diapycnal mixing in the polarity-reversing process, and (3) the difference between our results and Richardson-number-dependent turbulence parameterizations is about 2–3 orders of magnitude, but its vertical distribution is almost the same.

The generation and propagation of internal tides in the Andaman Sea are investigated using a three-dimensional high-resolution numerical model. Three categories of experiments, including driving the model with four main semidiurnal tides (M2, S2, N2, and K2), four main diurnal tides (K1, O1, P1, and Q1), and eight main tides (M2, S2, N2, K2, K1, O1, P1, and Q1), are designed to examine the effects of barotropic tides. The results show that the semidiurnal internal tides are dominant in the Andaman Sea, and the inclusion of diurnal barotropic tides negligibly modulates this result. That is partly due to the strength of the diurnal barotropic tides is generally one order smaller than that of the semidiurnal barotropic tides in this region. The sensitivity experiments put this on a firmer footing. In terms of the internal tidal energy, the experiments driven by the diurnal barotropic tides are three orders and one order smaller than those driven by the semidiurnal barotropic tides, respectively, during the spring and neap tides. In addition, the experiments result in total barotropic-to-baroclinic energy conversion rates over the Andaman Sea 29.15 GW (driven by the eight tides), 29.24 GW (driven by the four semidiurnal tides), and 0.05 GW (driven by the fourdiurnal tides) in the spring tidal period and 3.08 GW, 2.56 GW, and 0.31 GW in the neap tidal period, respectively. Four potential generation regions of internal tides are found, three of which are in the Andaman and Nicobar Islands and one in the northeastern Andaman Sea.

The cultural heritage of the Buddhist grottoes on the Gansu section of the Silk Road in China has important historical, cultural and artistic values. In order to better protect these precious cultural sites including through facilitating ongoing assessment of their condition, the Dunhuang Academy has created a pilot inventory based on the Arches open-source software platform. This paper first briefly introduces the most important grottoes and the representative caves in Gansu province. It then, in view of the diversity of the grotto architecture layout, describes the Arches semantic graph data modelling technology and how the flexibility of Arches Designer was used to design the resource models, which not only meet data management needs, but also ensure data interoperability and longevity through incorporation of the CIDOC-CRM ontology. In addition, through design of the grottoes map plugin and resource components within Arches, the simple interaction and intuitive visualization of grottoes information was achieved. The Dunhuang Academy’s deployment of Arches will serve the need for grottoes data management and provide support for the establishment of a more scientific and effective protection and management system for the grottoes on the Silk Road in Gansu province. It is hoped that this project will demonstrate the suitability of the Arches platform to support the management of other grotto sites outside of Gansu province, as well as for heritage management elsewhere in China.

The extended target probability hypothesis density (ET-PHD) filter is a promising approach for multiple extended target tracking. One crucial problem of the ET-PHD filter is partitioning the measurement set. This paper proposes a partitioning algorithm based on clustering by fast search and find density peaks (CFSFDP). Firstly, we adopt CFSFDP algorithm to partition the measurement set and the field theory is introduced to determine the cutoff distance of the CFSFDP algorithm. Then, the cluster center of the CFSFDP algorithm is determined according to solved cutoff distance and measurement rate. Finally, as the CFSFDP algorithm cannot handle the case in which targets are spatially close, an improved sub-partitioning method is implemented. Simulation results show that the proposed algorithm has less computational complexity and stronger robustness than the existing algorithm without losing tracking performance.

An experiment was conducted to investigate the effect of seasonally asymmetric warming on ecosystem respiration (Re), CH4 uptake, and N2O emissions in alpine grassland of the Tianshan of central Asia, from October 2016 to September 2019. The annual means of Re, CH4, and N2O fluxes in growing season were 42.83 mg C m−2 h−1, −41.57 µg C m−2 h−1, and 4.98 µg N m−2 h−1, respectively. Furthermore, warming during the non-growing season increased Re and CH4 uptake by 7.9 % and 10.6 % in the growing season and 10.5 % and 9.2 % in the non-growing season, respectively. However, the increase in N2O emission in the growing season was mainly caused by the warming during the growing season (by 29.7 %). The warming throughout the year and warming during the non-growing season increased N2O emissions by 101.9 % and 192.3 % in the non-growing season, respectively. The Re, CH4 uptake, and N2O emissions were positively correlated with soil temperature. Our results suggested that Re, CH4 uptake, and N2O emissions were regulated by soil temperature, rather than soil moisture, in the case of seasonally asymmetric warming. In addition, the response rate was defined by the changes in greenhouse gas fluxes driven by warming. In our field experiment, we observed the stimulatory effect of warming during the non-growing season on Re and CH4 uptake. In contrast, the response rates of Re and N2O emissions were gradually attenuated by long-term annual warming, and the response rate of Re was also weakened by warming over the growing season. These findings highlight the importance of warming in the non-growing season in regulating greenhouse gas fluxes, a finding which is crucial for improving our understanding of C and N cycles under the scenarios of global warming.

In recent years, the issue of violations or over-exploitation in artificial lake projects have attracted gradual concern. In the initial stage, in order to determine whether there are violations or over-exploitation in the artificial lake project, it is necessary to accurately identify the status of artificial lakes. Because of the relative lack of water resources along the Yellow River, issues such as the contradiction between supply and demand as well as the protection of ecological environment have become increasingly prominent. Besides, the situation of development of artificial lakes along the Yellow River has become serious, increasing the pressure of utilizing water resources along the Yellow River. Therefore, the identification of status of artificial lake is of great significance in ecological supervision of the Yellow River.Under this context, this paper takes the Yellow River as the research area. Based on the land cover data in the multi-temporal geographic national conditions, using spatial analysis, automatic identification and manual differentiation etc., to develop a rapid classification system of artificial lakes to classify the artificial lakes along the Yellow River, tracking and analyzing its annual changes, and finally establishing an artificial lake classification database and an artificial lake change database. Realizing the effective identification and supervision of changes in the artificial lake. Through experiments and analysis of existing data, this paper summarizes a set of indicators for determining the characteristics of artificial lakes, and proposes a set of effective artificial lake identification approaches. It provides an important basis for the follow-up supervision of artificial lake, and also establishes a basic framework for identifying other surface geographic-conditions, which has important theoretical and practical significance in the application of national geographic-conditions data.

To investigate the composition, variation, and sources of nitrated phenols (NPs) in the winter of Beijing, gas-phase NPs were measured by a chemical ionization long time-of-flight mass spectrometer (CI-LToF-MS). A box model was applied to simulate the secondary formation process of NPs. In addition, the primary sources of NPs were resolved by a non-negative matrix factorization (NMF) model. Our results showed that secondary formation contributed 38 %, 9 %, 5 %, 17 %, and almost 100 % of the nitrophenol (NP), methyl-nitrophenol (MNP), dinitrophenol (DNP), methyl-dinitrophenol (MDNP or DNOC), and dimethyl-nitrophenol (DMNP) concentrations. The phenol–OH reaction was the predominant loss pathway (46.7 %) during the heavy pollution episode, which produced the phenoxy radical (C6H5O). The phenoxy radical consequently reacted with NO2 and produced nitrophenol. By estimating the primarily emitted phenol from the ratio of phenol/CO from freshly emitted vehicle exhaust, this study proposed that oxidation of primary phenol contributes much more nitrophenol (37 %) than that from benzene oxidation ( %) in the winter of Beijing. The latter pathway was widely used in models and might lead to great uncertainties. The source apportionment results by NMF indicated the importance of combustion sources (>50 %) to the gas-phase NPs. The industry source contributed 30 % and 9 % to DNP and MDNP, respectively, which was non-negligible. The concentration weighted trajectory (CWT) analysis demonstrated that regional transport from provinces that surround the Yellow and Bohai seas contributed more primary NPs to Beijing. Both primary sources and secondary formation at either local or regional scale should be considered when making control policies of NPs.

Stabilized Criegee intermediates (SCIs) have the potential to oxidize trace species and to produce secondary organic aerosols (SOAs), making them important factors in tropospheric chemistry. This study quantitatively investigates the performance of SCIs in SOA formation at different relative humidity (RH) levels, and the first- and second-generation oxidations of endo- and exocyclic double bonds ozonated in limonene ozonolysis are studied separately. Through regulating SCI scavengers, the yields and rate constants of SCIs in a reaction system were derived, and the quantities of SCIs were calculated. The quantity of SOAs decreased by more than 20 % under low-humidity conditions (10 % RH–50 % RH), compared to that under dry conditions, due to the reactions of SCIs with water, while the inhibitory effect of water on SOA formation was not observed under high-humidity conditions (60 % RH–90 % RH). When using excessive SCI scavengers to exclude SCI reactions, it was found that the effect of water on SOA formation with the presence of SCIs was different from that without the presence of SCIs, suggesting that SCI reactions were relevant to the non-monotonic impact of water. The fractions of the SCI contribution to SOAs were similar between dry and high-humidity conditions, where the SCI reactions accounted for ∼ 63 % and ∼ 73 % in SOA formation in the first- and second-generation oxidation; however, marked differences in SOA formation mechanisms were observed. SOA formation showed a positive correlation with the quantity of SCIs, and the SOA formation potential of SCIs under high-humidity conditions was more significant than that under dry and low-humidity conditions. It was estimated that 20 %–30 % of SCIs could be converted into SOAs under high-humidity conditions, while this value decreased by nearly half under dry and low-humidity conditions. The typical contribution of limonene-derived SCIs to SOA formation is calculated to be (8.21 ± 0.15) × 10−2 µg m−3 h−1 in forest, (4.26 ± 0.46) × 10−2 µg m−3 h−1 in urban areas, and (2.52 ± 0.28) × 10−1 µg m−3 h−1 in indoor areas. Water is an uncertainty in the role SCIs play in SOA formation, and the contribution of SCIs to SOA formation needs consideration even under high RH in the atmosphere.

Soybean is one of the most important economic crops in the world. However, the salinization of soil results in the decrease of soybean yield as it is only a moderately salt-tolerant crop. We treated three soybean cultivars with low and high concentrations of NaCl. The differentially expressed genes between the control group and the salt treatment group were identified by mRNA sequencing and analyzed by gene ontology and Kyoto encyclopedia of genes and genomes annotations. We performed weighted gene co-expression network analysis on all samples and found genes most related to the phenotype. After verifying the results of differentially expressed genes by quantitative PCR, we finally identified Glyma06G01990, Glyma08G22730, Glyma019G05140, and Glyma06G20160 as key genes affecting the soybean growth under salt stress.

A covert communication system is proposed in this study, in which a~ non-zero mean Gaussian sequence is used as a~ random carrier and its mean is modulated by a~ covert binary bit. The aperiodic transmitted signal exhibits the same statistical characteristics as the ambient noise to confuse an~ eavesdropper. The received signal is multiplied with the pseudo-random sequence synchronized with the transmitter to recover these positive and negative mean Gaussian sequence. The sample mean estimator and hard decision are used to determine the covert message, and accordingly, theoretical bit error rate in additive white Gaussian noise channel is also derived. Simulation results are very consistent with the theoretical derivation. The proposed system works in the physical layer with the advantages of simple structure, strong concealment, good BER performance and very suitable for low-cost, resource-limited and low-rate transmission devices.

Estimating area of impervious land cover is the most useful and one of the ecological assessment indexes of urban and regional environment. Global land cover maps are inevitably misclassified, which affects the quality and application of the data. Statistical approach for assessing the accuracy is critical to understand the global change information and area estimation is usually based on sample data with a probability-based estimator. However, research on evaluation of multi-temporal global impervious land cover maps has not been implemented. In this study, spatial characteristics of the data are considered to assess the thematic map accuracy with a two-stage stratified random sampling plan. The first stage of stratification is determined by the global urban ecoregion and the second one is determined by land cover classes. Additionally, sample size of both map stage and pixel stage are calculated using a probability sampling model. A response design is constructed for a per-pixel accuracy assessment and blind interpretation is implemented using sample pixels and its surrounding area. Our method is applied to the multi-temporal global impervious land cover maps between 2000 and 2010 with a time interval of 5 years and the estimated area in different epoch are listed in detail. The main contribution of our research is illustrating the details for calculating the proportion area of impervious land cover and corresponding confidence intervals based on the reference classification. The experimental results show that the increasing area of the impervious surface according to the sample unit shows good agreement with the urbanization and descriptive accuracy assessments by user’s, producer’s and overall accuracy are shown respectively.

Hydrogen peroxide (H2O2) is a vital oxidant in the atmosphere and plays critical roles in the oxidation chemistry of both liquid and aerosol phases. The partitioning of H2O2 between the gas and liquid phases, or the aerosol phase, could affect its abundance in these condensed phases and eventually the formation of secondary components. However, the partitioning processes of H2O2 in gas-liquid and gas-aerosol phases are still unclear, especially in the ambient atmosphere. In this study, field observations of gas-, liquid-, and aerosol-phase H2O2 were carried out in the urban atmosphere of Beijing during the summer and winter of 2018. The effective field-derived mean value of Henry's law constant (HAm, 2.1×105 M atm−1) was 2.5 times of the theoretical value in pure water (HAt, 8.4×104 M atm−1) at 298±2 K. The effective derived gas-aerosol partitioning coefficient (KPm, 3.8×10-3 m3 µg−1) was 4 orders of magnitude higher on average than the theoretical value (KPt, 2.8×10-7 m3 µg−1) at 270±4 K. Beyond following Henry's law or Pankow's absorptive partitioning theory, the partitioning of H2O2 in the gas-liquid and gas-aerosol phases in the ambient atmosphere was also influenced by certain physical and chemical reactions. The average concentration of liquid-phase H2O2 in rainwater during summer was 44.12±26.49 µM. In 69 % of the collected rain samples, the measured level of H2O2 was greater than the predicted value in pure water calculated by Henry's law. In these samples, 41 % of the measured H2O2 was from gas-phase partitioning, while most of the rest may be from residual H2O2 in raindrops. In winter, the level of aerosol-phase H2O2 was 0.093±0.085 ng µg−1, which was much higher than the predicted value based on Pankow's absorptive partitioning theory. The contribution of partitioning of the gas-phase H2O2 to the aerosol-phase H2O2 formation was negligible. The decomposition/hydrolysis rate of aerosol-phase organic peroxides could account for 11 %–74 % of the consumption rate of aerosol-phase H2O2, and the value depended on the composition of organic peroxides in the aerosol particles. Furthermore, the heterogeneous uptake of HO2 and H2O2 on aerosols contributed to 22 % and 2 % of the aerosol-phase H2O2 consumption, respectively.

In the present work, indirect noises generated by compositional disturbances in a non-isotropic convergent nozzle are studied using Large Eddy Simulations (LESs). An in-house compressible LES code, Boundary Fitted Flow Integrator-LESc, is utilized to simulate the noise generation in the system. A non-reflective outlet boundary condition is used to eliminate numerical reflections and to ensure the reproduction of the operating conditions in the experiments. The experiments are designed to feature two configurations with different injection positions, which enable the separation of direct and indirect noises. Different operating conditions are investigated, including different injection gases and air mass flow rates. This present paper compares computational results with the experimental measurements. The results revealed that the processes of direct and indirect noise generation are successfully reproduced in the LES, with the noise magnitudes in good agreement with those in the measurements. Injection of gases with smaller (He) and larger (CO2) molar masses compared to air is found to generate negative and positive indirect noises, respectively, in the LES, which is consistent with the experimental findings. The effect of different air mass flow rates is also investigated and discussed, and the direct noise and indirect noise amplitudes are both found to be closely related to the air mass flow rate.

A pot study was conducted to elucidate morphological and physiological responses of Spatholobus suberectus Dunn accessions to four nitrogen application rates (0, 10, 20, and 30 mmol L-1) and two watering regimes (daily irrigation and irrigation every 10 d for a month). The results showed that water deficit led to marked decreases in leaf water relations and photosynthetic parameters regardless of the N treatments. Under drought conditions, the reductions in net photosynthesis occurred mainly due to stomatal factors. Although the leaf N content and stem mass ratio were positively correlated with N application, other morphological and physiological characteristics were negatively correlated or nonsignificantly correlated with N application under limited water supply. As the root nodules of S. suberectus can provide sufficient N supply, and high N applications made S. suberectus plants more vulnerable to drought stress, we suggest that no N supplement is required for S. suberectus seedlings.

Abstract Background The basidiomycete Phanerochaete carnosa is a white-rot species that has been mainly isolated from coniferous softwood. Given the particular recalcitrance of softwoods to bioconversion, we conducted a comparative transcriptomic analysis of P. carnosa following growth on wood powder from one softwood (spruce; Picea glauca) and one hardwood (aspen; Populus tremuloides). P. carnosa was grown on each substrate for over one month, and mycelia were harvested at five time points for total RNA sequencing. Residual wood powder was also analyzed for total sugar and lignin composition. Results Following a slightly longer lag phase of growth on spruce, radial expansion of the P. carnosa colony was similar on spruce and aspen. Consistent with this observation, the pattern of gene expression by P. carnosa on each substrate converged following the initial adaptation. On both substrates, highest transcript abundances were attributed to genes predicted to encode manganese peroxidases (MnP), along with auxiliary activities from carbohydrate-active enzyme (CAZy) families AA3 and AA5. In addition, a lytic polysaccharide monooxygenase from family AA9 was steadily expressed throughout growth on both substrates. P450 sequences from clans CPY52 and CYP64 accounted for 50% or more of the most highly expressed P450s, which were also the P450 clans that were expanded in the P. carnosa genome relative to other white-rot fungi. Conclusions The inclusion of five growth points and two wood substrates was important to revealing differences in the expression profiles of specific sequences within large glycoside hydrolase families (e.g., GH5 and GH16), and permitted co-expression analyses that identified new targets for study, including non-catalytic proteins and proteins with unknown function.

Volatile organic compounds play an important role in air quality and climate change, largely because they contribute to the formation of oxidizing compounds and secondary organic aerosol (SOA). In this study, a series of products, including peroxides and carbonyl compounds in both gaseous and particulate phases, were simultaneously detected to investigate the oxidation regime and SOA composition in limonene ozonolysis. The roles of different double bonds (DBs), radicals, and water were also examined. In our first investigation, we focused on representative oxidizing compounds produced in limonene ozonolysis, including stabilized Criegee intermediates (SCIs), OH radicals, and peroxides. The dependence of H2O2 and hydroxymethyl hydroperoxide (HMHP) formation on RH demonstrates that the reaction with water is an important reaction pathway for limonene SCIs, and the lower limit SCI yields of endocyclic and exocyclic DBs were estimated to be ∼ 0.24 and ∼ 0.43, respectively. The OH yield was determined by adding sufficient amounts of an OH scavenger, and the OH yields of endocyclic and exocyclic DBs were ∼ 0.65 and ∼ 0.24, respectively. These results indicate that in limonene ozonolysis the endocyclic DB is inclined to generate OH radicals through the hydroperoxide channel, while the exocyclic DB has a higher fraction of forming SCIs. Additionally, other gas-phase and particle-phase peroxides were also studied in this work. The formation of performic acid (PFA) and peracetic acid (PAA) was promoted significantly by increasing RH and the degree of oxidation, and the discrepancy between the experimental and model results suggested some missing formation pathways. Considerable generation of H2O2 from SOA in the aqueous phase was observed, especially at a high [O3] ∕ [limonene] ratio, which was mainly attributed to the hydration and decomposition of unstable peroxides in SOA such as peroxycarboxylic acids and peroxyhemiacetals. Different DBs and OH scavengers had a large impact on the particulate peroxides, and their stability indicated that the types of peroxides in SOA changed under different conditions. As for the contribution of peroxides to SOA, the results demonstrated that the mass fraction of particulate peroxides in limonene SOA was less than 0.2 at a low [O3] ∕ [limonene] ratio, while the mass fraction was 0.4–0.6 at a high [O3] ∕ [limonene] ratio. The partitioning behavior of peroxides showed that multi-generation oxidation helped produce more low-volatility peroxides, which partially explained the higher SOA yield. The partitioning behavior of carbonyls was also examined and the experimental partitioning coefficients (Kp) were found to be typically several orders of magnitude higher than the theoretical values. This study provided new insights into the oxidation regime and SOA composition in limonene ozonolysis, and limonene showed its specificity in many aspects when both endocyclic and exocyclic DBs were ozonated. We suggest that the atmospheric implications of terpenes containing more than one DB and the SOA composition, especially particulate peroxides, need further study.

Abstract Non-polar (11–20) GaN with significantly improved crystal quality has been achieved by means of overgrowth on regularly arrayed micro-rod templates on sapphire in comparison with standard non-polar GaN grown without any patterning processes on sapphire. Our overgrown GaN shows massively reduced linewidth of X-ray rocking curves with typical values of 270 arcsec along the [0001] direction and 380 arcsec along the [1–100] direction, which are among the best reports. Detailed X-ray measurements have been performed in order to investigate strain relaxation and in-plane strain distribution. The study has been compared with the standard non-polar GaN grown without any patterning processes and an extra non-polar GaN sample overgrown on a standard stripe-patterned template. The standard non-polar GaN grown without involving any patterning processes typically exhibits highly anisotropic in-plane strain distribution, while the overgrown GaN on our regularly arrayed micro-rod templates shows a highly isotropic in-plane strain distribution. Between them is the overgrown non-polar GaN on the stripe-patterned template. The results presented demonstrate the major advantages of using our regularly arrayed micro-rod templates for the overgrowth of non-polar GaN, leading to both high crystal quality and isotropic in-plane strain distribution, which is important for the further growth of any device structures.

Disaster preparedness is critical for reducing potential impact. This paper contributes to current knowledge of disaster preparedness using representative national sample data from China, which faces high earthquake risks in many areas of the country. The adoption of earthquake preparedness activities by the general public, including five indicators of material preparedness and five indicators of awareness preparedness, were surveyed and 3245 respondents from all 31 provinces of Mainland China participated in the survey. Linear regression models and logit regression models were used to analyze the effects of potential influencing factors. Overall, the preparedness levels are not satisfied, with a material preparation score of 3.02 (1–5), and awareness preparation score of 2.79 (1–5), nationally. Meanwhile, residents from western China, which has higher earthquake risk, have higher degrees of preparedness. The concern for disaster risk reduction (DRR) and the concern for building safety and participation in public affairs are consistent positive predictors of both material and awareness preparedness. The demographic and socioeconomic variables' effects, such as gender, age, education, income, urban/rural division, and building size, vary according to different preparedness activities. Finally, the paper concludes with a discussion of the theoretical contribution and potential implementation.

Microseismic signals are generally considered to follow the Gauss distribution. A comparison of the dynamic characteristics of sample variance and the symmetry of microseismic signals with the signals which follow α-stable distribution reveals that the microseismic signals have obvious pulse characteristics and that the probability density curve of the microseismic signal is approximately symmetric. Thus, the hypothesis that microseismic signals follow the symmetric α-stable distribution is proposed. On the premise of this hypothesis, the characteristic exponent α of the microseismic signals is obtained by utilizing the fractional low-order statistics, and then a new method of time difference of arrival (TDOA) estimation of microseismic signals based on fractional low-order covariance (FLOC) is proposed. Upon applying this method to the TDOA estimation of Ricker wavelet simulation signals and real microseismic signals, experimental results show that the FLOC method, which is based on the assumption of the symmetric α-stable distribution, leads to enhanced spatial resolution of the TDOA estimation relative to the generalized cross correlation (GCC) method, which is based on the assumption of the Gaussian distribution.

The marine-terminating outlet in Basin 3, Austfonna ice cap, has been accelerating since the mid-1990s. Stepwise multi-annual acceleration associated with seasonal summer speed-up events was observed before the outlet entered the basin-wide surge in autumn 2012. We used multiple numerical models to explore hydrologic activation mechanisms for the surge behaviour. A continuum ice dynamic model was used to invert basal friction coefficient distributions using the control method and observed surface velocity data between April 2012 and July 2014. This has provided input to a discrete element model capable of simulating individual crevasses, with the aim of finding locations where meltwater entered the glacier during the summer and reached the bed. The possible flow paths of surface meltwater reaching the glacier bed as well as those of meltwater produced at the bed were calculated according to the gradient of the hydraulic potential.The inverted friction coefficients show the unplugging of the stagnant ice front and expansion of low-friction regions before the surge reached its peak velocity in January 2013. Crevasse distribution reflects the basal friction pattern to a high degree. The meltwater reaches the bed through the crevasses located above the margins of the subglacial valley and the basal melt that is generated mainly by frictional heating flows either to the fast-flowing units or potentially accumulates in an overdeepened region. Based on these results, the mechanisms facilitated by basal meltwater production, crevasse opening and the routing of meltwater to the bed are discussed for the surge in Basin 3.

For the purpose of extracting productions of some specific branching plants effectively and realizing its 3D reconstruction, Terrestrial LiDAR data was used as extraction source of production, and a 3D reconstruction method based on Terrestrial LiDAR technologies combined with the L-system was proposed in this article. The topology structure of the plant architectures was extracted using the point cloud data of the target plant with space level segmentation mechanism. Subsequently, L-system productions were obtained and the structural parameters and production rules of branches, which fit the given plant, was generated. A three-dimensional simulation model of target plant was established combined with computer visualization algorithm finally. The results suggest that the method can effectively extract a given branching plant topology and describes its production, realizing the extraction of topology structure by the computer algorithm for given branching plant and also simplifying the extraction of branching plant productions which would be complex and time-consuming by L-system. It improves the degree of automation in the L-system extraction of productions of specific branching plants, providing a new way for the extraction of branching plant production rules.

Soil respiration (Rs) is the most important source of carbon dioxide emissions from soil to atmosphere. However, it is unclear what the interactive response of Rs would be to environmental changes such as elevated precipitation, nitrogen (N) deposition and warming, especially in unique temperate desert ecosystems. To investigate this an in situ field experiment was conducted in the Gurbantunggut Desert, northwest China, from September 2014 to October 2016. The results showed that precipitation and N deposition significantly increased Rs, but warming decreased Rs, except in extreme precipitation events, which was mainly through its impact on the variation of soil moisture at 5 cm depth. In addition, the interactive response of Rs to combinations of the factors was much less than that of any single-factor, and the main response was a positive effect, except for the response from the interaction of increased precipitation and high N deposition (60 kg N ha−1 yr−1). Although Rs was found to show a unimodal change pattern with the variation of soil moisture, soil temperature and soil NH4+-N content, and it was significantly positively correlated to soil dissolved organic carbon (DOC) and pH, a structural equation model found that soil temperature was the most important controlling factor. Those results indicated that Rs was mainly interactively controlled by the soil multi-environmental factors and soil nutrients, and was very sensitive to elevated precipitation, N deposition and warming. However, the interactions of multiple factors largely reduced between-year variation of Rs more than any single-factor, suggesting that the carbon cycle in temperate deserts could be profoundly influenced by positive carbon–climate feedback.

By applying 12-year (1998–2009) radiosonde data over a midlatitude station, we studied the vertical wavenumber spectra of three-dimensional wind fluctuations. The horizontal wind spectra in the lower stratosphere coincide well with the well-known universal spectra, with mean spectral slopes of −2.91 ± 0.09 and −2.99 ± 0.09 for the zonal and meridional wind spectra, respectively, while the mean slopes in the troposphere are −2.64 ± 0.07 and −2.70 ± 0.06, respectively, which are systematically less negative than the canonical slope of −3. In both the troposphere and lower stratosphere, the spectral amplitudes (slopes) of the horizontal wind spectra are larger (less negative) in winter, and they are larger (less negative) in the troposphere than in the lower stratosphere. Moreover, we present the first statistical results of vertical wind fluctuation spectra, which revealed a very shallow spectral structure, with mean slopes of −0.58 ± 0.06 and −0.23 ± 0.05 in the troposphere and lower stratosphere, respectively. Such a shallow vertical wind fluctuation spectrum is considerably robust. Different from the horizontal wind spectrum, the slopes of the vertical wind spectra in both the troposphere and lower stratosphere are less negative in summer. The height variation of vertical wind spectrum amplitude is also different from that of the horizontal wind spectrum, with a larger amplitude in the lower stratosphere. These evident differences between the horizontal and vertical wind spectra strongly suggest they should obey different spectral laws. Quantitative comparisons with various theoretical models show that no existing spectral theories can comprehensively explain the observed three-dimensional wind spectra, indicating that the spectral features of atmospheric fluctuations are far from fully understood.

The Miyun Reservoir, the only surface water source for Beijing city, has experienced water supply decline in recent decades. Previous studies suggest that both land use change and climate contribute to the changes of water supply in this critical watershed. However, the specific causes of the decline in the Miyun Reservoir are debatable under a non-stationary climate in the past 4 decades. The central objective of this study was to quantify the separate and collective contributions of land use change and climate variability to the decreasing inflow into the Miyun Reservoir during 1961–2008. Different from previous studies on this watershed, we used a comprehensive approach to quantify the timing of changes in hydrology and associated environmental variables using the long-term historical hydrometeorology and remote-sensing-based land use records. To effectively quantify the different impacts of the climate variation and land use change on streamflow during different sub-periods, an annual water balance model (AWB), the climate elasticity model (CEM), and a rainfall–runoff model (RRM) were employed to conduct attribution analysis synthetically. We found a significant (p p p > 0.1) negative trend in annual precipitation during 1961–2008. We identified two streamflow breakpoints, 1983 and 1999, by the sequential Mann–Kendall test and double-mass curve. Climate variability alone did not explain the decrease in inflow to the Miyun Reservoir. Reduction of water yield was closely related to increase in actual evapotranspiration due to the expansion of forestland and reduction in cropland and grassland, and was likely exacerbated by increased water consumption for domestic and industrial uses in the basin. The contribution to the observed streamflow decline from land use change fell from 64–92 % during 1984–1999 to 36–58 % during 2000–2008, whereas the contribution from climate variation climbed from 8–36 % during the 1984–1999 to 42–64 % during 2000–2008. Model uncertainty analysis further demonstrated that climate warming played a dominant role in streamflow reduction in the most recent decade (i.e., 2000s). We conclude that future climate change and variability will further challenge the water supply capacity of the Miyun Reservoir to meet water demand. A comprehensive watershed management strategy needs to consider the climate variations besides vegetation management in the study basin.

Spatial and temporal variabilities of Kelvin waves during stratospheric sudden warming (SSW) events are investigated by the ERA-Interim reanalysis data, and the results are validated by the COSMIC temperature data. A case study on an exceptionally large SSW event in 2009, and a composite analysis comprising 18 events from 1980 to 2013 are presented. During SSW events, the average temperature increases by 20 K in the polar stratosphere, while the temperature in the tropical stratosphere decreases by about 4 K. Kelvin wave with wave numbers 1 and 2, and periods 10–20 days, clearly appear around the tropical tropopause layer (TTL) during SSWs. The Kelvin wave activity shows obvious coupling with the convection localized in the India Ocean and western Pacific (Indo-Pacific) region. Detailed analysis suggests that the enhanced meridional circulation driven by the extratropical planetary wave forcing during SSW events leads to tropical upwelling, which further produces temperature decrease in the tropical stratosphere. The tropical upwelling and cooling consequently result in enhancement of convection in the equatorial region, which excites the strong Kelvin wave activity. In addition, we investigated the Kelvin wave acceleration to the eastward zonal wind anomalies in the equatorial stratosphere during SSW events. The composite analysis shows that the proportion of Kelvin wave contribution ranges from 5 to 35 % during SSWs, much larger than in the non-SSW mid-winters (less than 5 % in the stratosphere). However, the Kelvin wave alone is insufficient to drive the equatorial eastward zonal wind anomalies during the SSW events, which suggests that the effects of other types of equatorial waves may not be neglected.

By using meteor radar, radiosonde and satellite observations over 20° N and NCEP/NCAR reanalysis data during 81 days from 22 December 2004 to 12 March 2005, a quasi-27-day oscillation propagating from the troposphere to the mesosphere is reported. A pronounced 27-day periodicity is observed in the raw zonal wind from meteor radar. Spectral analysis shows that the oscillation also occurs in the meridional wind and temperature and propagates westward with wavenumber s = 1; thus the oscillation is of Rossby wave type. The oscillation attains a large amplitude of about 12 m s−1 in the eastward wind shear region of the troposphere. When the wind shear reverses, its amplitude rapidly decays, and the background wind gradually evolves to be westward. However, the oscillation can penetrate through the weak westward wind field due to its relatively large phase speed. After this, the oscillation restrengthens with its upward propagation and reaches about 20 m s−1 in the mesosphere. Reanalysis data show that the oscillation can propagate to the mid and high latitudes from the low latitudes and has large amplitudes over there. There is another interesting phenomenon that a quasi-46-day oscillation appears simultaneously in the troposphere, but it cannot penetrate through the westward wind field because of its smaller phase speed. In the observational interval, a quasi-27-day periodicity in outgoing long-wave radiation (OLR) and specific humidity is found in a latitudinal zone of 5–20° N. Thus the quasi-27-day oscillation may be an atmospheric response to forcing due to the convective activity with a period of about 27 days in the tropical region.

The diurnal tide (DT) and its variability in the lower atmosphere over Chongyang (114.14° E, 29.53° N) were studied based on the newly established Wuhan University (WHU) VHF radar observations with the height intervals of 0.145 km (below 9 km) and 0.58 km (above 9 km) in the whole year of 2012. We find that the DT was the dominant tidal component and showed remarkable height and season variations. A prominent seasonally dependent height variability characteristic is that maximum DT amplitude usually occurs around 6 km in the winter and spring months, which might be due to the tidal wave energy concentration arising from the reflections from the strong eastward tropospheric jet around 13 km and the ground surface. Our results suggest that the background wind is a crucial cause for height variability and seasonal variability of DT. In April 2012, a notable strengthening of DT is observed. Meanwhile, the significant higher harmonics of tides, i.e., the semidiurnal, terdiurnal, and even quarterdiurnal tides, can also be observed, which has seldom been reported. Interestingly, these four tidal components displayed consistent short-term variability, implying that they were excited by the same dramatically varying tidal source. In addition, we identified two symptoms of the coupling of DT and planetary waves (PWs), which can also lead to the short-term DT variability. One is the sum and difference interactions between DT and PWs, causing the tidal amplitude short-term variability as a consequence of the energy exchange among the interacting waves. The other one is the modulation of DT by PWs, leading to that the amplitude of DT varies with the periods of the PWs.

Optical properties of aged WS2 monolayers grown by CVD method on Si/SiO2 substrates are studied using temperature dependent photoluminescence and reflectance contrast spectroscopy. Aged WS2 monolayers have a typical surface roughness about 0.5 nm and, in addition, a high density of nanoparticles (nanocaps) with the base diameter about 30 nm and average height of 7 nm. The A-exciton of aged monolayer has a peak position at 1.951 eV while in as-grown monolayer the peak is at about 24 meV higher energy at room temperature. This red-shift is explained using local tensile strain concept, where strain value of 2.1% was calculated for these nanocap regions. Strained nanocaps have lower band gap energy and excitons will funnel into these regions. At T=10K a double exciton and trion peaks were revealed. The separation between double peaks is about 20 meV and the origin of higher energy peaks is related to the optical band gap energy fluctuations caused by random distribution of local tensile strain due to increased surface roughness. In addition, a wide defect related exciton band XD was found at about 1.93 eV in all aged monolayers. It is shown that the theory of localized excitons describes well the temperature dependence of peak position and halfwidth of the A-exciton band. The possible origin of nanocaps is also discussed.

(11-22) semi-polar GaN is expected to exhibit major advantages compared with current c-plane polar GaN in the fabrication of long wavelength such as green and yellow emitters. However, all the III-nitride based semi-/non- polar laser diodes (LDs) reported so far have been achieved exclusively based on homo-epitaxial growth on extremely expensive free-standing GaN substrates with a very limited size. In this paper, we have observed a stimulated emission at room temperature achieved on our semi-polar (11-22) GaN overgrown on a micro-rod arrayed template with an optimized design on m-plane sapphire. This has never been achieved previously on any semi-polar GaN grown on sapphire. Furthermore, an optical gain of 130cm−1 has been measured by means of performing a standard laser stripe-length dependent optical measurement. The values of the threshold and the optical gain obtained are comparable to those of the c-plane GaN reported so far, further validating the satisfactory crystal quality of our overgrown semi-polar (11-22) GaN on sapphire. This represents a major step towards the development of III-nitride semi-polar based LDs on sapphire, especially in the long wavelength regime.

We extended the broad spectral method proposed by Zhang et al. (2013) for the extraction of medium- and high-frequency gravity waves (MHGWs). This method was applied to 11 years (1998–2008) of radiosonde data from 92 stations in the Northern Hemisphere to investigate latitudinal, continuous vertical and seasonal variability of MHGW parameters in the lower atmosphere (2–25 km). The latitudinal and vertical distributions of the wave energy density and horizontal momentum fluxes as well as their seasonal variations exhibit considerable consistency with those of inertial gravity waves. Despite the consistency, the MHGWs have much larger energy density, horizontal momentum fluxes and wave force, indicating the more important role of MHGWs in energy and momentum transportation and acceleration of the background. For the observed MHGWs, the vertical wavelengths are usually larger than 8 km; the horizontal wavelengths peak in the middle troposphere at middle–high latitudes. These characteristics are obviously different from inertial gravity waves. The energy density and horizontal momentum fluxes have similar latitude-dependent seasonality: both of them are dominated by a semiannual variation at low latitudes and an annual variation at middle latitudes; however at high latitudes, they often exhibit more than two peaks per year in the troposphere. Compared with the inertial GWs, the derived intrinsic frequencies are more sensitive to the spatiotemporal variation of the buoyancy frequency, and at all latitudinal regions they are higher in summer. The wavelengths have a weaker seasonal variation; an evident annual cycle can be observed only at middle latitudes.

The outlet glacier of Basin 3 (B3) of Austfonna ice cap, Svalbard, is one of the fastest outlet glaciers in Svalbard, and shows dramatic changes since 1995. In addition to previously observed seasonal summer speed-up associated with the melt season, the winter speed of B3 has accelerated approximately fivefold since 1995. We use the Elmer/Ice full-Stokes model for ice dynamics to infer spatial distributions of basal drag for the winter seasons of 1995, 2008 and 2011. This "inverse" method is based on minimising discrepancy between modelled and observed surface velocities, using satellite remotely sensed velocity fields. We generate steady-state temperature distributions for 1995 and 2011. Frictional heating caused by basal sliding contributes significantly to basal temperatures of the B3 outlet glacier, with heat advection (a longer-timescale process than frictional heating) also being important in the steady state. We present a sensitivity experiment consisting of transient simulations under present-day forcing to demonstrate that using a temporally fixed basal drag field obtained through inversion can lead to thickness change errors of the order of 2 m year−1. Hence it is essential to incorporate the evolution of basal processes in future projections of the evolution of B3. Informed by a combination of our inverse method results and previous studies, we hypothesise a system of processes and feedbacks involving till deformation and basal hydrology to explain both the seasonal accelerations (short residence time pooling of meltwater at the ice–till interface) and the ongoing interannual speed-up (gradual penetration of water into the till, reducing till strength).

The interaction between the climate system and the large polar ice sheet regions is a key process in global environmental change. We carried out dynamic ice simulations of one of the largest drainage systems in East Antarctica: the Lambert Glacier–Amery Ice Shelf system, with an adaptive mesh ice sheet model. The ice sheet model is driven by surface accumulation and basal melt rates computed by the FESOM (Finite-Element Sea-Ice Ocean Model) ocean model and the RACMO2 (Regional Atmospheric Climate Model) and LMDZ4 (Laboratoire de Météorologie Dynamique Zoom) atmosphere models. The change of ice thickness and velocity in the ice shelf is mainly influenced by the basal melt distribution, but, although the ice shelf thins in most of the simulations, there is little grounding line retreat. We find that the Lambert Glacier grounding line can retreat as much as 40 km if there is sufficient thinning of the ice shelf south of Clemence Massif, but the ocean model does not provide sufficiently high melt rates in that region. Overall, the increased accumulation computed by the atmosphere models outweighs ice stream acceleration so that the net contribution to sea level rise is negative.

Starting from a set of fully nonlinear equations, this paper studies that two initial gravity wave packets interact to produce a third substantial packet in a nonisothermal and dissipative atmosphere. The effects of the inhomogeneous temperature and dissipation on interaction are revealed. Numerical experiments indicate that significant energy exchange occurs through the nonlinear interaction in a nonisothermal and dissipative atmosphere. Because of the variability of wavelengths and frequencies of interacting waves, the interaction in an inhomogeneous temperature field is characterised by the nonresonance. The nonresonant three waves mismatch mainly in the vertical wavelengths, but match in the horizontal wavelengths, and their frequencies also tend to match throughout the interaction. Below 80 km, the influence of atmospheric dissipation on the interaction is rather weak due to small diffusivities. With the further propagation of wave above 80 km, the exponentially increasing atmospheric dissipation leads to the remarkable decay and slowly upward propagation of wave energy. Even so, the dissipation below 110 km is not enough to decrease the vertical wavelength of wave. The dissipation seems neither to prevent the interaction occurrence nor to prolong the period of wave energy exchange, which is different from the theoretical prediction based on the linearised equations. The match relationship and wave energy evolution in numerical experiments are helpful in further investigating interaction of gravity waves in the middle atmosphere based on experimental observations.