Performance evaluation of operational atmospheric correction algorithms over the East China Seas

To acquire high-quality operational data products for Chinese in-orbit and scheduled ocean color sensors, the performances of two operational atmospheric correction(AC) algorithms(ESA MEGS 7.4.1 and NASA Sea DAS 6.1) were evaluated over the East China Seas(ECS) using MERIS data. The spectral remote sensing reflectance R_(rs)(λ), aerosol optical thickness(AOT), and ?ngstr?m exponent(α) retrieved using the two algorithms were validated using in situ measurements obtained between May 2002 and October 2009. Match-ups of R_(rs), AOT, and α between the in situ and MERIS data were obtained through strict exclusion criteria. Statistical analysis of R_(rs)(λ) showed a mean percentage difference(MPD) of 9%–13% in the 490–560 nm spectral range, and significant overestimation was observed at 413 nm(MPD72%). The AOTs were overestimated(MPD32%), and although the ESA algorithm outperformed the NASA algorithm in the blue-green bands, the situation was reversed in the red-near-infrared bands. The value of α was obviously underestimated by the ESA algorithm(MPD=41%) but not by the NASA algorithm(MPD=35%). To clarify why the NASA algorithm performed better in the retrieval of α, scatter plots of the α single scattering albedo(SSA) density were prepared. These α-SSA density scatter plots showed that the applicability of the aerosol models used by the NASA algorithm over the ECS is better than that used by the ESA algorithm, although neither aerosol model is suitable for the ECS region. The results of this study provide a reference to both data users and data agencies regarding the use of operational data products and the investigation into the improvement of current AC schemes over the ECS.     

Evaluation of spectral scale effects in estimation of vegetation leaf area index using spectral indices methods

Spectral index methodology has been widely used in Leaf Area Index (LAI) retrieval at different spatial scales. There are differences in the spectral response of different remote sensors and thus spectral scale effect generated during the use of spectral indices to retrieve LAI. In this study, PROSPECT, leaf optical properties model and Scattering by Arbitrarily Inclined Layers (SAIL) model, were used to simulate canopy spectral reflectance with a bandwidth of 5 nm and a Gaussian spectral response function was employed to simulate the spectral data at six bandwidths ranging from 10 to 35 nm. Additionally, for bandwidths from 5 to 35 nm, the correlation between the spectral index and LAI, and the sensitivities of the spectral index to changes in LAI and bandwidth were analyzed. Finally, the reflectance data at six bandwidths ranging from 40 to 65 nm were used to verify the spectral scale effect generated during the use of the spectral index to retrieve LAI. Results indicate that Vegetation Index of the Universal Pattern Decomposition (VIUPD) had the highest accuracy during LAI retrieval. Followed by Normalized Difference Vegetation Index (NDVI), Modified Simple Ratio Indices (MSRI) and Triangle Vegetation Index (TVI), although the coefficient of determination R ² was higher than 0.96, the retrieved LAI values were less than the actual value and thus lacked validity. Other spectral indices were significantly affected by the spectral scale effect with poor retrieval results. In this study, VIUPD, which exhibited a relatively good correlation and sensitivity to LAI, was less affected by the spectral scale effect and had a relatively good retrieval capability. This conclusion supports a purported feature independent of the sensor of this model and also confirms the great potential of VIUPD for retrieval of physicochemical parameters of vegetation using multi-source remote sensing data.     

Distribution of diatoms and silicoflagellates in surface sediments of the Yellow Sea and offshore from the Changjiang River,China

The spatial distribution of siliceous microfossils (diatoms and silicoflagellates) in the surface sediments was mapped at 113 sites in the Yellow Sea and sea areas adjacent to the Changjiang (Yangtze) River, China. In total, 267 diatom taxa and two silicoflagellate species were identified from the sediments. The spatial variations in abundance and diversity were classified into three distinct geographic patterns using Q mode clustering: a south-north geographic pattern, a coastal-offshore pattern and a unique pattern in the Changjiang River mouth. The south-north geographic pattern was related to the spatial variations in sea temperature. Coscinodiscus oculatus, a warm-water species, indicated these variations by a gradual decrease in abundance from the south to the north. The coastal-offshore pattern was in response to the spatial variations in salinity. Cyclotella stylorum, Actinocyclus ehrenbergii and Dictyocha messanensis, the dominant brackish species in coastal waters, significantly decreased at the isobaths of approximately 30 m, where the salinity was higher than 31. Paralia sulcata and Podosira stelliger indicated the impact of the Yellow Sea Warm Current in the central Yellow Sea. The unique pattern in the Changjiang River mouth showed the highest species diversity but lower abundance, apparently because: freshwater input can significantly increase the proportion of brackish species; nutrients can supply the growth of phytoplankton; and high sedimentation rates can dilute the microfossil abundance in the sediments. Our results show that an integration of environmental factors (e.g., nutrient levels, sedimentation rate, sea temperature, salinity and water depth) determined the spatial characteristics of the siliceous microfossils in the surface sediments.     

Distribution of sediment chloroplastic pigments in the Southern Yellow Sea,China

The distribution of sediment chloroplastic pigments (Chl-a, i.e. chlorophyll a and Pha-a, i.e. phaeophorbide a) in the Southern Yellow Sea of China was studied. Samples were collected from four cruises in January and June 2003, and January and June 2004. The results show that the vertical distribution of Chl-a and Pha-a in the sediment layers 0-2cm, 2-5cm and 5-8cm, follows a stable ratio, 5:3:2. The average ratio of Pha-a to Chl-a in sediment is 2.83. Spearman 2-tailed rank correlation analysis shows that Chl-a and Pha-a contents in each sediment layer have a highly significant correlation. The average contents of Chl-a and Pha-a in the sediment of the 0-8cm layer in the investigated area are 0.31 -0.47μgg-1 and 1.28-1.40 μgg-1 sediment (dry weight), respectively. The average Chl-a and Pha-a contents in sediment are higher in summer than in winter. ANOVA analysis shows that there is a highly significant variation among the Chl-a contents (P = 0.002 <0.01) of the four cruies, but this is not true for the case of Pha-a content (P = 0.766>0.05). The average Chl-a and Pha-a contents in the 2 sediment layers (0-2cm and 2-5cm) have significant or highly significant correlations with organic matter (OM), median diameter (Mdφ), silt plus clay percentage in the January 2003 cruise. In the June 2003 cruise, the average Chl-a content in the 3 sediment layers (0-2cm, 2-5cm, and 5-8cm) has a significant correlation with meiofauna biomass, and Pha-a content has highly significant correlations with water depth, bottom water temperature, OM and Mdφ The contents of Chl-a and Pha-a are lower than those in estuaries and intertidal areas, but close to those in the same area studied previously.     

Effects of Anthropogenic Disturbance on Sediment Organic Carbon Mineralization Under Different Water Conditions in Coastal Wetland of a Subtropical Estuary

The changes in soil organic carbon(C) mineralization as affected by anthropogenic disturbance directly determine the role of soils as C source or sink in the global C budget. The objectives of this study were to investigate the effects of anthropogenic disturbance(aquaculture pond, pollutant discharge and agricultural activity) on soil organic C mineralization under different water conditions in the Minjiang River estuary wetland, Southeast China. The results showed that the organic C mineralization in the wetland soils was significantly affected by human disturbance and water conditions(P 0.001), and the interaction between human disturbance activities and water conditions was also significant(P 0.01). The C mineralization rate and the cumulative mineralized carbon dioxide-carbon(CO_2-C)(at the 49th day) ranked from highest to lowest as follows: Phragmites australis wetland soil aquaculture pond sediment soil near the discharge outlet rice paddy soil. This indicated that human disturbance inhibited the mineralization of C in soils of the Minjiang River estuary wetland, and the inhibition increased with the intensity of human disturbance. The data for cumulative mineralized CO_2-C showed a good fit(R~2 0.91) to the first-order kinetic model C_t = C_0(1 – exp(–kt)). The kinetic parameters C_0, k and C_0 k were significantly affected by human disturbance and water conditions. In addition, the total amount of mineralized C(in 49 d) was positively related to C_0, C_0 k and electrical conductivity of soils. These findings indicated that anthropogenic disturbance suppressed the organic C mineralization potential in subtropical coastal wetland soils, and changes of water pattern as affected by human activities in the future would have a strong influence on C cycling in the subtropical estuarine wetlands.     

Predicting Surface Roughness and Moisture of Bare Soils Using Multiband Spectral Reflectance Under Field Conditions

Soil surface roughness, denoted by the root mean square height(RMSH), and soil moisture(SM) are critical factors that affect the accuracy of quantitative remote sensing research due to their combined influence on spectral reflectance(SR). In regards to this issue, three SM levels and four RMSH levels were artificially designed in this study; a total of 12 plots was used, each plot had a size of 3 m × 3 m. Eight spectral observations were conducted from 14 to 30 October 2017 to investigate the correlation between RMSH, SM, and SR. On this basis, 6 commonly used bands of optical satellite sensors were selected in this study, which are red(675 nm), green(555 nm), blue(485 nm), near infrared(845 nm), shortwave infrared 1(1600 nm), and shortwave infrared 2(2200 nm). A negative correlation was found between SR and RMSH, and between SR and SM. The bands with higher coefficient of determination R~2 values were selected for stepwise multiple nonlinear regression analysis. Four characterized bands(i.e., blue, green, near infrared, and shortwave infrared 2) were chosen as the independent variables to estimate SM with R~2 and root mean square error(RMSE) values equal to 0.62 and 2.6%, respectively. Similarly, the four bands(green, red, near infrared, and shortwave infrared 1) were used to estimate RMSH with R~2 and RMSE values equal to 0.48 and 0.69 cm, respectively. These results indicate that the method used is not only suitable for estimating SM but can also be extended to the prediction of RMSH. Finally, the evaluation approach presented in this paper highly restores the real situation of the natural farmland surface on the one hand, and obtains high precision values of SM and RMSH on the other. The method can be further applied to the prediction of farmland SM and RMSH based on satellite and unmanned aerial vehicle(UAV) optical imagery.     

Equality testing for soil grid unit resolutions to polygon unit scales with DNDC modeling of regional SOC pools

Matching soil grid unit resolutions with polygon unit map scales is important to minimize the uncertainty of regional soil organic carbon(SOC) pool simulation due to their strong influences on the modeling.A series of soil grid units at varying cell sizes was derived from soil polygon units at six map scales,namely,1:50 000(C5),1:200 000(D2),1:500 000(P5),1:1 000 000(N1),1:4 000 000(N4) and 1:14 000 000(N14),in the Taihu Region of China.Both soil unit formats were used for regional SOC pool simulation with a De Nitrification-DeC omposition(DNDC) process-based model,which spans the time period from 1982 to 2000 at the six map scales.Four indices,namely,soil type number(STN),area(AREA),average SOC density(ASOCD) and total SOC stocks(SOCS) of surface paddy soils that were simulated by the DNDC,were distinguished from all these soil polygon and grid units.Subjecting to the four index values(IV) from the parent polygon units,the variations in an index value(VIV,%) from the grid units were used to assess its dataset accuracy and redundancy,which reflects the uncertainty in the simulation of SOC pools.Optimal soil grid unit resolutions were generated and suggested for the DNDC simulation of regional SOC pools,matching their respective soil polygon unit map scales.With these optimal raster resolutions,the soil grid units datasets can have the same accuracy as their parent polygon units datasets without any redundancy,when VIV 1% was assumed to be a criterion for all four indices.A quadratic curve regression model,namely,y = – 0.80 × 10~(–6)x~2 + 0.0228 x + 0.0211(R~2 = 0.9994,P 0.05),and a power function model R? = 10.394?~(0.2153)(R~2 = 0.9759,P 0.05) were revealed,which describe the relationship between the optimal soil grid unit resolution(y,km) and soil polygon unit map scale(1:10 000x),the ratio(R?,%) of the optimal soil grid size to average polygon patch size(?,km~2) and the ?,with the highest R~2 among different mathematical regressions,respectively.This knowledge may facilitate the grid partitioning of regions during the investigation and simulation of SOC pool dynamics at a certain map scale,and be referenced to other landscape polygon patches' mesh partition.     

Identifying the driving factors of sediment delivery ratio on individual flood events in a long-term monitoring headwater basin

The Sediment Delivery Ratio (SDR) has multi-fold environmental implications both in evaluating the soil and water losses and the effectiveness of conservation measures in watersheds. Various factors, including hydrological regime and watershed properties, may influence the SDR at interannual timescales. However, the effect of certain important dynamic factors, such as rainfall peak distribution, runoff erosion power and sediment bulk density, on the sediment delivery ratio of single flood events (SDRe) has received little attention. The Qiaogou headwater basin is in the hilly-gully region of the Chinese Loess Plateau, and it encompasses a 0.45 km2 catchment. Three large-scale field runoff plots at different geomorphological positions were chosen to obtain the observation data, and the 20-year period between 1986 and 2005 is presented. The results showed that the SDRe of the Qiaogou headwaters varied from 0.49 to 2.77. Among the numerous influential factors, rainfall and runoff were the driving factors causing slope erosion and sediment transport. The rainfall erosivity had a significant positive relationship with the sediment transport modulus (R²=0.85, P<0.01) but had no significant relationship with SDRe. The rainfall peak coefficient was significantly positively correlated with the SDRe (R²=0.64, P<0.05), indicating the influence of rainfall energy distribution on the SDRe. The runoff erosion power index was not only significantly related to the sediment transport modulus (R²=0.84, P<0.01) but also significantly related to the SDRe (R²=0.57, P<0.01). In addition, the relative bulk density was significantly related to the SDRe, indicating that hyper-concentrated flow characteristics contributed to more transported sediment in the catchment. Thus, the rainfall peak coefficient, runoff erosion power and sediment relative bulk density could be used as dynamic indexes to predict the SDRe in the hilly areas of the Chinese Loess Plateau.     

Changes of Biogenic Elements in <Emphasis Type="Italic">Phragmites australis</Emphasis> and <Emphasis Type="Italic">Suaeda salsa</Emphasis> from Salt Marshes in Yellow River Delta,China

Little information is available on biogenic elements (carbon, nitrogen, phosphorus and sulfur) and the ecological stoichiometric characteristics of plants in coastal wetlands. To investigate the contents of carbon, nitrogen, phosphorus and sulfur of plants, and their ecological stoichiometric characteristics in the Yellow (Huanghe) River Delta, plant samples were collected from two typical salt marshes (Suaeda salsa and Phragmites australis wetlands) during the period of from August to October in 2007, and the ratios of C/N, C/P, N/P, C/N/P and C/N/P/S were calculated. Results showed that during the studying period, plant C, N and P were lower than the global average values, and plant N and P were lower than the China’s average values. Leaf C and S in Suaeda salsa were significantly lower than those in Phragmites australis (P < 0.05), and leaf N and P in Suaeda salsa and Phragmites australis showed no significant differences (P > 0.05). Average C/N ratios were 23.75 in leaf, 73.36 in stem, 65.67 in root of Suaeda salsa, and 33.77 in leaf, 121.68 in stem, 97.13 in root of Phragmites australis. Average C/N ratios of Suaeda salsa and Phragmites australis were all great than 25, indicating the salt marsh in the Yellow River Delta is an N limitation system. Average C/P ratios were 276.78 in leaf, 709.28 in stem and 1031.32 in root of Suaeda salsa, and 536.94 in leaf, 768.13 in stem and 875.22 in root of Phragmites australis. The average N/P ratios of Suaeda salsa were 12.92 in leaf, 10.77 in stem and 10.91 in root, and the average N/P ratios of Phragmites australis were 16.40 in leaf, 7.40 in stem and 6.92 in root, indicating the Suaeda salsa wetlands were N limited and Phragmites australis wetlands were N limited in August and P limited in October in 2007. The average C/N, C/P and C/N/P ratios in Suaeda salsa and Pragmites australis were higher than the global average values, indicating the lower quality of organic matter provided by wetland plants in the Yellow River delta.     

Diel variation in metabolism and ammonia excretion of <Emphasis Type="Italic">Marphysa sanguinea</Emphasis> (Polychaeta: Eunicidae)

Polychaetes provide an excellent food resource for fish and represent the dominant zoobenthos in marine ecosystems. Diel variation in the rates of metabolism and ammonia-N excretion of Marphysa sanguinea were studied. The worms were grouped according to their wet body weight into small (S; 1.24±0.06 g), medium (M; 4.00±0.30 g), and large (L; 8.54±1.08 g) categories. Their weight-specific metabolic rates, based on aerobic respiration (R), were measured at 16°C (±0.2°C) and classed as either routine (R _R) or standard (R _S) rates. Both respiration types decreased with increasing body weight. Respiration was described by R = a W ^b, where b was -0.400 9 and -0.532 0 for R _R and R _S, respectively. Diurnal changes in R _S for each group was relatively flat, with a slightly increasing trend with time, but was relatively stable as a whole. R _R of the diurnal variation of worms was higher than R S, but both had similar overall trends. The peak values of specific dynamic action (SDA) (R _SDA) in the S, M, and L groups were 2.704, 1.149, and 0.682 mg/(g?h), respectively. The durations of SDA were 13, 6, and 6 h, respectively and the energy expenditures of SDA were 377.98, 117.34, and 74.94 J/g, respectively. These data indicate that the metabolic rates were higher in smaller individuals, which is advantageous for their rapid growth.     

Emissions of biogenic sulfur gases (H<Subscript>2</Subscript>S,COS) from <Emphasis Type="Italic">Phragmites australis</Emphasis> coastal marsh in the Yellow River estuary of China

Emissions of biogenic sulfur gases (hydrogen sulfide (H₂S) and carbonyl sulfide (COS)) from Phragmites australis coastal marsh in the Yellow River estuary of China were determined during April to December in 2014 using static chamber-gas chromatography technique with monthly sampling. The results showed that the fluxes of H₂S and COS both had distinct seasonal and diurnal variations. The H₂S fluxes ranged from 0.09 μg/(m²·h) to 7.65 μg/(m²·h), and the COS fluxes ranged from–1.10 μg/(m²·h) to 3.32 μg/(m²·h). The mean fluxes of H₂S and COS from the P. australis coastal marsh were 2.28 μg/(m²·h), and 1.05 μg/(m²·h), respectively. The P. australis coastal marsh was the emission source of both H₂S and COS over the whole year. Fluxes of H₂S and COS were both higher in plant growing season than in the non-growing season. Temperature had a dramatic effect on the H₂S emission flux, while the correlations between COS flux and the environmental factors were not found during sampling periods. More in-depth and comprehensive research on other related factors, such as vegetation, sediment substrates, and tidal action is needed to discover and further understand the key factors and the release mechanism of sulfur gases.     

Distribution patterns of chlorophyll <Emphasis Type="Italic">a</Emphasis> in spring and autumn in association with hydrological features in the southern Yellow Sea and northern East China Sea

Two field studies were conducted to measure pigments in the Southern Yellow Sea (SYS) and the northern East China Sea (NECS) in April (spring) and September (autumn) to evaluate the distribution pattern of phytoplankton stock (Chl a concentration) and the impact of hydrological features such as water mass, mixing and tidal front on these patterns. The results indicated that the Chl a concentration was 2.43±2.64 (Mean ± SD) mg m^?3 in April (range, 0.35 to 17.02 mg m^?3) and 1.75±3.10 mg m^?3 in September (from 0.07 to 36.54 mg m^?3) in 2003. Additionally, four areas with higher Chl a concentrations were observed in the surface water in April, while two were observed in September, and these areas were located within or near the point at which different water masses converged (temperature front area). The distribution pattern of Chl a was generally consistent between onshore and offshore stations at different depths in April and September. Specifically, higher Chl a concentrations were observed along the coastal line in September, which consisted of a mixing area and a tidal front area, although the distributional pattern of Chl a concentrations varied along transects in April. The maximum Chl a concentration at each station was observed in the surface and subsurface layer (0–10 m) for onshore stations and the thermocline layer (10–30 m) for offshore stations in September, while the greatest concentrations were generally observed in surface and subsurface water (0–10 m) in April. The formation of the Chl a distributional pattern in the SYS and NECS and its relationship with possible influencing factors is also discussed. Although physical forces had a close relationship with Chl a distribution, more data are required to clearly and comprehensively elucidate the spatial pattern dynamics of Chl a in the SYS and NECS.     

Study on the Liquefaction Properties of Silt at Yellow River Delta

The seabed of the Yellow River Delta is formed by the rapid deposition of sediments from the Yellow River. Recent researches have shown that the geological hazards in the Yellow River Delta are mainly related to the liquefaction of silty seabed under cyclic loading. In this paper, based on the theory of Stokes Viscous Principle, a self-design dragging ball apparatus was used to study the fluid characteristics of liquefied and post-liquefied silt more thoroughly. Wave flume and shaking table were used to apply wave loads and vibration loads. The pore pressure, earth pressure, and moving parameters of the dragging ball were recorded during tests. The effect of cyclic loads and excess pore pressure ratio on viscosity of silty soil is discussed. The constitutive model of flow characteristics of liquefied silt is also established. The results indicate that the apparent viscosity of silty soil is much higher than sand. Cyclic loads accelerated the process of silt liquefaction, resulting in the decrease in apparent viscosity. The excess pore pressure ratio (r_u) has effect on silt liquefaction. When the soil had not reached the completely liquefied state (r_u < 1), the apparent viscosity decreased with the increase of excess pore pressure ratio; when the soil had reached the completely liquefied state (r_u = 1), the shear stress decreased with the increase of strain rate. It was determined that silty soil characterized by shear thinning can be viewed as a type of non-Newtonian fluid. The equation of the shear stress and shear strain rates was established by fitting the test data.     

Molecular weight controllable degradation of <Emphasis Type="Italic">Laminaria japonica</Emphasis> polysaccharides and its antioxidant properties

In this study, molecular weight controllable degradation of algal Laminaria japonica polysaccharides（LPS） was investigated by ultrasound combined with hydrogen peroxide. Three main factors, i.e., ultrasonic power（A）, ultrasonic time（B）, and H₂O₂ concentration（C） were chosen for optimizing parameters by employing three-factors, three-levels BBD. The influence of degradation on structure change and antioxidant activities was also investigated. A second-order polynomial equation including molecular weight（Y） of Laminaria japonica polysaccharides and each variable parameter, i.e., ultrasonic power（A）, ultrasonic time（B）, and H₂O₂ concentration（C）, was established: Y=20718.67-4273.13A-4000.38B-1438.75C+2333.25AB+1511.00AC+873.00BC+2838.29A² + 2490.79B²+873.04C². The equation regression coefficient value（R² = 0.969） indicated that this equation was valid. The value of the adjusted determination coefficient（adjusted R² = 0.914） also confirmed that the model was highly significant. The results of selected experimental degradation conditions matched with the predicted value. FT-IR spectra revealed that the structures of LPS before and after degradation were not significantly changed. Antioxidant activities of LPS revealed that low Mws possessed stronger inhibitory than the original polysaccharides. The scavenging effects on superoxide radicals was the highest when IC50 of crude LPS was 4.92 mg mL^-1 and IC50 of Mw 18.576 KDa was 1.02 mg mL^-1, which was fourfold higher than initial polysaccharide.     

Soil Organic Carbon Contents and Stocks in Coastal Salt Marshes with <Emphasis Type="Italic">Spartina alterniflora</Emphasis> Following an Invasion Chronosequence in the Yellow River Delta,China

Plant invasion alters the fundamental structure and function of native ecosystems by affecting the biogeochemical pools and fluxes of materials and energy. Native (Suaeda salsa) and invasive (Spartina alterniflora) salt marshes were selected to study the effects of Spartina alterniflora invasion on soil organic carbon (SOC) contents and stocks in the Yellow River Delta. Results showed that the SOC contents (g/kg) and stocks (kg/m²) were significantly increased (P < 0.05) after Spartina alterniflora invasion of seven years, especially for the surface soil layer (0–20 cm). The SOC contents exhibited an even distribution along the soil profiles in native salt marshes, while the SOC contents were gradually decreased with depth after Spartina alterniflora invasion of seven years. The natural ln response ratios (LnRR) were applied to identify the effects of short-term Spartina alterniflora invasion on the SOC stocks. We also found that Spartina alterniflora invasion might cause soil organic carbon losses in a short-term phase (2–4 years in this study) due to the negative LnRR values, especially for 20–60 cm depth. And the SOCD in surface layer (0–20 cm) do not increase linearly with the invasive age. Spearman correlation analysis revealed that silt + clay content was exponentially related with SOC in surface layer (Adjusted R2 = 0.43, P < 0.001), suggesting that soil texture could play a key role in SOC sequestration of coastal salt marshes.     

Effects of vegetation type on surface elevation change in Liaohe River Delta wetlands facing accelerated sea level rise

Rising sea levels threaten the sustainability of coastal wetlands around the globe. The ability of coastal marshes to maintain their position in the intertidal zone depends on the accumulation of both organic and inorganic materials, and vegetation is important in these processes. To study the effects of vegetation type on surface elevation change, we measured surface accretion and elevation change from 2011 to 2016 using rod surface elevation table and feldspar marker horizon method(RSET-MH) in two Phragmites and two Suaeda marshes in the Liaohe River Delta. The Phragmites marshes exhibited higher rates of surface accretion and elevation change than the Suaeda marshes. The two Phragmites marsh sites had average surface elevation change rates at 8.78 mm/yr and 9.26 mm/yr and surface accretion rates at 17.56 mm/yr and 17.88 mm/yr, respectively. At the same time, the two Suaeda marsh sites had average surface elevation change rates at 5.77 mm/yr and 5.91 mm/yr and surface accretion rates at 13.42 mm/yr and 14.38 mm/yr, respectively. The elevation change rates in both the Phragmites marshes and the Suaeda marshes in the Liaohe River Delta could keep pace and even continue to gain elevation relative to averaged sea level rise in the Bohai Sea reported by the 2016 State Oceanic Administration, People's Republic of China projection(2.4–5.5 mm/yr) in current situations. Our data suggest that vegetation is important in the accretionary processes and vegetation type could regulate the wetland surface elevation. However, the vulnerability of coastal wetlands in the Liaohe River Delta need further assessment considering the accelerated sea level rise, the high rate of subsidence, and the declining sediment delivery, especially for the Suaeda marshes.     

Global water vapor content decreases from 2003 to 2012: An analysis based on MODIS data

Water vapor in the earth′s upper atmosphere plays a crucial role in the radiative balance, hydrological process, and climate change. Based on the latest moderate-resolution imaging spectroradiometer(MODIS) data, this study probes the spatio-temporal variations of global water vapor content in the past decade. It is found that overall the global water vapor content declined from 2003 to 2012(slope b = –0.0149, R = 0.893, P = 0.0005). The decreasing trend over the ocean surface(b = –0.0170, R = 0.908, P = 0.0003) is more explicit than that over terrestrial surface(b = –0.0100, R = 0.782, P = 0.0070), more significant over the Northern Hemisphere(b = –0.0175, R = 0.923, P = 0.0001) than that over the Southern Hemisphere(b = –0.0123, R = 0.826, P = 0.0030). In addition, the analytical results indicate that water vapor content are decreasing obviously between latitude of 36°N and 36°S(b = 0.0224, R = 0.892, P = 0.0005), especially between latitude of 0°N and 36°N(b = 0.0263, R = 0.931, P = 0.0001), while the water vapor concentrations are increasing slightly in the Arctic regions(b = 0.0028, R = 0.612, P = 0.0590). The decreasing and spatial variation of water vapor content regulates the effects of carbon dioxide which is the main reason of the trend in global surface temperatures becoming nearly flat since the late 1990 s. The spatio-temporal variations of water vapor content also affect the growth and spatial distribution of global vegetation which also regulates the global surface temperature change, and the climate change is mainly caused by the earth's orbit position in the solar and galaxy system. A big data model based on gravitational-magmatic change with the solar or the galactic system is proposed to be built for analyzing how the earth's orbit position in the solar and galaxy system affects spatio-temporal variations of global water vapor content, vegetation and temperature at large spatio-temporal scale. This comprehensive examination of water vapor changes promises a holistic understanding of the global climate change and potential underlying mechanisms.     

Spatial distribution of snow depth based on geographically weighted regression kriging in the Bayanbulak Basin of the Tianshan Mountains,China

Snow depth is a general input variable in many models of agriculture,hydrology,climate and ecology.This study makes use of observational data of snow depth and explanatory variables to compare the accuracy and effect of geographically weighted regression kriging(GWRK)and regression kriging(RK)in a spatial interpolation of regional snow depth.The auxiliary variables are analyzed using correlation coefficients and the variance inflation factor(VIF).Three variables,Height,topographic ruggedness index(TRI),and land surface temperature(LST),are used as explanatory variables to establish a regression model for snow depth.The estimated spatial distribution of snow depth in the Bayanbulak Basin of the Tianshan Mountains in China with a spatial resolution of 1 km is obtained.The results indicate that 1)the result of GWRK's accuracy is slightly higher than that of RK(R~2=0.55 vs.R~2=0.50,RMSE(root mean square error)=0.102 m vs.RMSE=0.077 m);2)for the subareas,GWRK and RK exhibit similar estimation results of snow depth.Areas in the Bayanbulak Basin with a snow depth greater than 0.15m are mainly distributed in an elevation range of 2632.00–3269.00 m and the snow in this area comprises 45.00–46.00% of the total amount of snow in this basin.However,the GWRK resulted in more detailed information on snow depth distribution than the RK.The final conclusion is that GWRK is better suited for estimating regional snow depth distribution.     

Trophic Interaction in a <Emphasis Type="Italic">Portunus rituberculatus</Emphasis> Polyculture Ecosystem Based on Carbon and Nitrogen Stable Isotope Analysis

Trophic interaction among various biomass groups in a swimming crab Portunus trituberculatus polyculture pond was investigated using carbon and nitrogen stable isotope analysis. The polycultured animal species also included white shrimp Litopenaeus vannamei, short-necked clam Ruditapes philippinarum, and redlip mullet Liza haematochila. The mean δ¹³C value for all the biomass groups in polyculture ecosystem ranged from ?25.61‰ to ?16.60‰, and the mean δ¹⁵N value ranged from 6.80‰ to 13.09‰. Significant difference in the δ¹³C value was found between particulate organic matter (POM) and sediment organic matter (SOM) (P < 0.05), indicating that these two organic matter pools have different material sources. Assuming that a ¹³C-enrichment factor of 1.00‰ and a ¹⁵N-enrichment factor of 2.70‰ existed between consumer and prey, diets of the four cultured animals were estimated using a stable isotope mixing model. The estimated model results indicated that P. trituberculatus mainly feed on Aloidis laevis; L. vannamei mainly feed on shrimp feed; while A. laevis, R. philippinarum and L. haematochelia mainly feed on POM. Shrimp feed was also an important food source of R. philippinarum and L. haematochelia. The diets of P. trituberculatus, L. vannamei, R. philippinarum, and L. haematochila showed complementary effects in this polyculture ecosystem. Our finding indicated that the polyculture of these four organisms with suitable farming density could make an effective use of most of the food sources, which can make a highly efficient polyculture ecosystem.     

Two new species of the genus <Emphasis Type="Italic">Cobitis</Emphasis> Linnaeus (Teleostei: Cobitidae) from southern China

Two new species of the genus Cobitis from southern China, C. hereromacula from the Luohe River in Guangdong Province and C. baishagensis from the Nandujiang River in Hainan Province, are described and illustrated here. C. hereromacula can be distinguished from its congeners by possessing the following combination of characteristics: absence of the second and third pigmentary zones of Gambetta; 13–16 oval blotches on the dorsum and 10–13 vertical, elongated triangular blotches below the midlateral line with more than 20 vertical dark brown bars between them; 6–7 narrow rows of dark spots on the caudal fin; a vertical oval spot smaller than the eye diameter on the upper part of the caudal peduncle; pointed mental lobes of the lower lip pointed with a slightly filiform tip; one slender and long needle-shaped lamina circularis at the base of the first branched ray of the male pectoral fins. C. baishagensis can be distinguished from its congeners by the fourth Gambetta zone being covered by 10–12 transverse elongated blotches; 4–5 narrow rows of dark spots on the caudal fin; a vertical blotch smaller than the eye diameter on the upper part of the caudal peduncle; males with a slender and long needle-shaped lamina circularis at the second branched pectoral fin ray in males; large scales with a slightly large focal zone; undeveloped mental lobes with a lower lip that does not end posteriorly in a filiform tip.