A Directed-Line-Segment-Based Photogrammetry Method to Survey the Tunneling Working Surface of Mine Roadway

Absolute orientation is a basic technical work in digital image geologic logging of underground coal mine. Traditional control-point-based absolute orientation method requires setting object space control points of the known three-dimensional coordinates, which may lead to low efficiency. Therefore, this paper proposed a point-free close-range photogrammetry absolute orientation algorithm, which utilized direction line segments including plumb line segments and line segments with known directions and lengths to identify the dimensional orientation of a stereoscopic model. Experiment results show that the precision of the orientation results is favorable. σ _X and σ _Y are as high as 0.5 mm, and σ _Z is 0.3 mm. Finally, this paper introduced the application of the proposed algorithm in rapid geological logging of coal mine roadway, which was fast and reliable, convenient and feasible.     

GOCE gravitational gradiometry

GOCE is the first gravitational gradiometry satellite mission. Gravitational gradiometry is the measurement of the second derivatives of the gravitational potential. The nine derivatives form a 3 × 3 matrix, which in geodesy is referred to as Marussi tensor. From the basic properties of the gravitational field, it follows that the matrix is symmetric and trace free. The latter property corresponds to Laplace equation, which gives the theoretical foundation of its representation in terms of spherical harmonic or Fourier series. At the same time, it provides the most powerful quality check of the actual measured gradients. GOCE gradiometry is based on the principle of differential accelerometry. As the satellite carries out a rotational motion in space, the accelerometer differences contain angular effects that must be removed. The GOCE gradiometer provides the components V _xx , V _yy , V _zz and V _xz with high precision, while the components V _xy and V _yz are of low precision, all expressed in the gradiometer reference frame. The best performance is achieved inside the measurement band from 5 × 10^–3 to 0.1 Hz. At lower frequencies, the noise increases with 1/f and is superimposed by cyclic distortions, which are modulated from the orbit and attitude motion into the gradient measurements. Global maps with the individual components show typical patterns related to topographic and tectonic features. The maps are separated into those for ascending and those for descending tracks as the components are expressed in the instrument frame. All results are derived from the measurements of the period from November to December 2009. While the components V _xx and V _yy reach a noise level of about \({10\;\rm{\frac{mE}{\sqrt{Hz}}}}\), that of V _zz and V _xz is about \({20\; \rm{\frac{mE}{\sqrt{Hz}}}}\). The cause of the latter’s higher noise is not yet understood. This is also the reason why the deviation from the Laplace condition is at the \({20 \;\rm{\frac{mE}{\sqrt{Hz}}}}\) level instead of the originally planned \({11\;\rm{\frac{mE}{\sqrt{Hz}}}}\). Each additional measurement cycle will improve the accuracy and to a smaller extent also the resolution of the spherical harmonic coefficients derived from the measured gradients.     

Improved method for estimating the inter-frequency satellite clock bias of triple-frequency GPS

Considering the contribution of the hardware biases to the estimated clock errors, an improved method for estimating the satellite inter-frequency clock bias (IFCB) is presented, i.e., the difference in the satellite clock error as computed from ionospheric-free pseudorange and carrier phase observations using L1/L2 and P1/P2 versus L1/L5 and P1/P5. The IFCB is composed of a constant and a variable part. The constant part is the inter-frequency hardware bias (IFHB). It contains the satellite and receiver hardware delays and can be expressed as a function of the DCBs [DCB (P1 ? P2) and DCB (P1 ? P5)]. When a reference satellite is selected, the satellite IFHB can be computed but is biased by a reference satellite IFHB. This bias will not affect the utilization of IFCB in positioning since it can be absorbed by the receiver clock error. Triple-frequency observations of 30 IGS stations between June 1, 2013, and May 31, 2014, were processed to show the variations of the IFHB. The IFHB values show a long-term variation with time. When a linear and a fourth-order harmonic function are used to model the estimated IFCB, which contains contributions of the hardware delays and clock errors, the results show that 89 % of the IFCB can be corrected given the current five triple-frequency GPS satellites with the averaged fitting RMS of 1.35 cm. Five days of data are processed to test the estimated satellite clock errors using the strategy presented. The residuals of P1/P5 and L1/L5 have a STD of <0.27 m and 0.97 cm, respectively. In addition, most predicted satellite IFCBs reach an accuracy of centimeter level and its mean accuracy of 5 days is better than 7 cm.     

Proximal Remote Sensing of Herbicide and Drought Stress in Field Grown Colocasia and Sweet Potato Plants by Sunlight-Induced Chlorophyll Fluorescence Imaging

Chlorophyll fluorescence is an indicator of plant photosynthetic activity and has been used to monitor the health status of vegetation. Several studies have exploited the application of red/far-red chlorophyll fluorescence ratio in detecting the impact of various types of stresses in plants. Recently, sunlight-induced chlorophyll fluorescence imaging has been used to detect and discriminate different stages of mosaic virus infection in potted cassava plants with a multi-spectral imaging system (MSIS). In this study, the MSIS is used to investigate the impact of drought and herbicide stress in field grown crop plants. Towards this control and treatment groups of colocasia and sweet potato plants were grown in laterite soil beds and the reflectance images of these crop plants were recorded up to 14-days of treatment at the Fraunhofer lines of O₂ B at 687 nm and O₂ A at 759.5 nm and the off-lines at 684 and 757.5 nm. The recorded images were analyzed using the Fraunhofer Line Discrimination technique to extract the sunlight-induced chlorophyll fluorescence component from the reflectance images of the plant leaves. As compared to the control group, the chlorophyll fluorescence image ratio (F ₆₈₇/F ₇₆₀) in the treatment groups of both the plant varieties shows an increasing trend with increase in the extent of stress. Further, the F ₆₈₇/F ₇₆₀ ratio was found to correlate with the net photosynthetic rate (Pn) and stomatal conductance (g_s) of leaves. The correlation coefficient (R ²) for the relationship of F ₆₈₇/F ₇₆₀ ratio with Pn were found to be 0.78, 0.79 and 0.78, respectively for the control, herbicide treated and drought treated colocasia plants, while these were 0.77, 0.86 and 0.88, respectively for sweet potato plants. The results presented show the potential of proximal remote sensing and the application F ₆₈₇/F ₇₆₀ fluorescence image ratio for effective monitoring of stress-induced changes in field grown plants.     

<Emphasis Type="Italic">C</Emphasis>/<Emphasis Type="Italic">N</Emphasis><Subscript>0</Subscript> estimators for high-sensitivity snapshot GNSS receivers

We address the problem of estimating the carrier-to-noise ratio (C/N₀) in weak signal conditions. There are several environments, such as forested areas, indoor buildings and urban canyons, where high-sensitivity global navigation satellite system (HS-GNSS) receivers are expected to work under these reception conditions. The acquisition of weak signals from the satellites requires the use of post-detection integration (PDI) techniques to accumulate enough energy to detect them. However, due to the attenuation suffered by these signals, estimating their C/N₀ becomes a challenge. Measurements of C/N₀ are important in many applications of HS-GNSS receivers such as the determination of a detection threshold or the mitigation of near-far problems. For this reason, different techniques have been proposed in the literature to estimate the C/N₀, but they only work properly in the high C/N₀ region where the coherent integration is enough to acquire the satellites. We derive four C/N₀ estimators that are specially designed for HS-GNSS snapshot receivers and only use the output of a PDI technique to perform the estimation. We consider four PDI techniques, namely non-coherent PDI, non-quadratic non-coherent PDI, differential PDI and truncated generalized PDI and we obtain the corresponding C/N₀ estimator for each of them. Our performance analysis shows a significant advantage of the proposed estimators with respect to other C/N₀ estimators available in the literature in terms of estimation accuracy and computational resources.     

Application of Remotely-Sensed Data to Estimate a Water Budget for Data-Scarce Endorheic Basins: A Case Study of Lake Urmia basin,Iran

Estimating the water budgets of large basins is a challenge because of the lack of data and information. It becomes more complicated in endorheic basins that consist of separate land and water phases. The application of remotely-sensed data is one solution in this regard. The present study addresses this issue and develops a modeling framework to evaluate a water budget based on remotely-sensed data for endorheic basins. To explore the methodology, Lake Urmia basin was selected as a case study. The lake water level has declined steeply since 1995 and stakeholders have agreed to allocate 3100 MCM of water per year to the lake. This makes it necessary to monitor river inflow into the lake to fulfill the agreement. Gauging stations have been employed around the lake, but they could not account for shortages such as water uptake below the stations. To do this, separate water budgets for the water body and the land were required. More specifically, it was necessary to estimate actual evapotranspiration (ET _a ) from freshwater (E _f ) and saltwater (E _s ) estimated using the SEBAL model. Different methods were applied to estimate soil moisture, groundwater exploitation, and surface-groundwater inflow into the lake. A comparison of the observed and estimated amounts showed good agreement. For instance, the coefficient of determination for the observed/reported and estimated ET _a and E _f were 0.83 and 0.84, respectively. The average annual inflow was estimated to be 2.2 BCM/year for 2002–2008 using the RS model, which is about 84 % of the total inflow from the last recording stations before the lake and shows influence of water exploitation after these stations. Future study should focus on increasing temporal and spatial resolution of the method     

Precise and efficient evaluation of gravimetric quantities at arbitrarily scattered points in space

Gravimetric quantities are commonly represented in terms of high degree surface or solid spherical harmonics. After EGM2008, such expansions routinely extend to spherical harmonic degree 2190, which makes the computation of gravimetric quantities at a large number of arbitrarily scattered points in space using harmonic synthesis, a very computationally demanding process. We present here the development of an algorithm and its associated software for the efficient and precise evaluation of gravimetric quantities, represented in high degree solid spherical harmonics, at arbitrarily scattered points in the space exterior to the surface of the Earth. The new algorithm is based on representation of the quantities of interest in solid ellipsoidal harmonics and application of the tensor product trigonometric needlets. A FORTRAN implementation of this algorithm has been developed and extensively tested. The capabilities of the code are demonstrated using as examples the disturbing potential T, height anomaly \(\zeta \), gravity anomaly \(\Delta g\), gravity disturbance \(\delta g\), north–south deflection of the vertical \(\xi \), east–west deflection of the vertical \(\eta \), and the second radial derivative \(T_{rr}\) of the disturbing potential. After a pre-computational step that takes between 1 and 2 h per quantity, the current version of the software is capable of computing on a standard PC each of these quantities in the range from the surface of the Earth up to 544 km above that surface at speeds between 20,000 and 40,000 point evaluations per second, depending on the gravimetric quantity being evaluated, while the relative error does not exceed \(10^{-6}\) and the memory (RAM) use is 9.3 GB.     

A new GPS/BDS tropospheric delay resolution approach for monitoring deformation in super high-rise buildings

A new approach for deformation monitoring of super high-rise building using GPS/BDS technology is proposed for the case when prior coordinates are known and the baseline is short but has a large height difference. The approach is based on the ambiguity function method (AFM). Considering that the double-differenced (DD) troposphere delay residual error cannot be ignored, the relative zenith tropospheric delay (RZTD) parameter is introduced into the original AFM equation. Thus, the RZTD and 3D coordinate parameters are together obtained through the modified AFM (MAFM). Due to the low computational efficiency of conventional AFM, an improved particle swarm optimization (IPSO) algorithm is used to search the four optimal parameters X/Y/Z/RZTD and replaces the grid search method. In this study, GPS/BDS deformation monitoring data for buildings with approximately 290 m height difference were used to verify the feasibility of the proposed MAFM. Numerical results show a single-epoch average computation time of approximately 0.3 s, which meets the requirements of near-real-time dynamic monitoring. The average accuracy of the GPS single-epoch RZTD solution is better than 1 cm, the combined GPS/BDS MAFM performance outperforms the GPS-only system, and using multi-epoch observations can further improve the accuracy of the RZTD solution. After RZTD correction, GPS/BDS monitoring precision can be improved, particularly the height dimension, whose precision is improved by approximately 6 cm.     

Mapping the Extent of Agroforestry Area in Different Altitudes of Tehri District,North Western Himalaya,India Through GIS and Remote Sensing Data

Assessment of area under agroforestry in Tehri district of North Western Himalaya, Uttarakhand, India has been done using GIS and remote sensing technology. The study district characterized by hilly terrain with varying elevations from 288 m to more than 2800 m and generally gentle slopes, valleys, flat land covers and agricultural terraces. High-resolution satellite imageries (spatial resolution 5.8 m) were used in this study for land uses and land covers classification. According to unsupervised classification, highest area was found under forest class (65.22%) followed by cropland (20.41%). Considerable area was also found under snow cover (9.45%) in the district. Area under agroforestry was estimated to be 5572.26 ha (1.53%) by this method, whereas it was estimated to be 7029.06 ha (1.93%) by supervised classification. Estimated cropland area comes out to be about 20.0%. An accuracy of 86.5% was found in this classification for agroforestry class. Highest area under agroforestry of 3707.36 ha was obtained in 1200–2000 m elevations followed by 2231.26 ha in 288–1200 m elevations. Negligible area was found on high elevation zones of more than 2800 m. The major agroforestry systems of dominated by Grewia oppositifolia (Bhimal), Celtis australis (Kharik) and Quercus leucotrichophora (Banj) were identified and mapped and remaining systems were grouped as others class. Estimated area under G. oppositifolia, C. australis and Q. leucotrichophora based systems come out to be 2330.82, 1456.80 and 1129.10 ha, respectively. These systems are multiple usufructs are food, fuelwood, fodder, fiber and small timber. It has been observed from the accuracy assessment that the estimates of area under agroforestry obtained under this study are reliable.     

Monitoring of <Emphasis Type="Italic">Aphis gossypii</Emphasis> Using Greenseeker and SPAD Meter

Cotton aphid (Aphis gossypii) is considered as one of the most important agriculture pest for the cotton production. However, it is generally labor-intensive and time-consuming to obtain some information of Cotton aphid with conventional methods through direct measurement by sampling in the field. This study explores the potential of using a new method to obtain information of the Cotton aphid rapidly. In our study, the cotton canopy spectral indices (NDVI, VI_2, REDrefc, NIRrefc) and chlorophyll concentration, obtained from hand-held high spectrometer GreenSeeker and chlorophyll meter SPAD-502 and Cotton aphid amount derived from the artificial field-based survey were used to uncover the relationship between Cotton aphid amount and canopy spectral index and SPAD value of the cotton in city of Shihezi, China. The results showed that NDVI and NIRrefc were negatively related to Cotton aphid amount. VI_2 content had a significant and positive relationship with its amount. The non-linear three cubic models with alate Aphid amount as independent variables have been established between VI_2 value and alatae Aphid amount, which could explain 92.37 % of VI_2 value variance. SPAD values were also significantly and negatively correlated to the Aphid amount. The non-linear logarithm model with wingless Aphid amount as independent variables was the best for uncovering the relationship between SPAD value and wingless Aphid amount, which could explain 85.48 % of SPAD value variance. The results demonstrate the establishment of the function model provides a theoretical basis and techniques for indirect and rapid monitoring and management of Cotton aphid.     

Global morphology of ionospheric sporadic E layer from the FormoSat-3/COSMIC GPS radio occultation experiment

Ionospheric sporadic-E (Es) activity and global morphology were studied using the 50 Hz signal-to-noise ratio amplitude and excess phase measurements from the FormoSat-3/Constellation Observing System for Meteorology, Ionosphere and Climate (FS3/COSMIC) GPS radio occultation (RO) observations. The results are presented for data collected during the last sunspot cycle from mid-2006 to the end of 2017. The FS3/COSMIC generally performed more than 1000 complete E-region GPS RO observations per day, which were used to retrieve normalized L1-band amplitude standard deviation (SDL1) and relative electron density (N_e) profiles successfully. More or less 31% of those observations were identified as Es events based on SDL1 and peak SDL1 altitude criteria. We found that the peak Es-event i values are approximately proportional to the logarithms of the corresponding peak N_e differences. Five major geographical zones were identified, in which the seasonal and diurnal Es occurrence patterns are markedly different. These five zones include the geomagnetic equatorial zone (??5°?<?magnetic latitude (ML)?<?5°), two extended geomagnetic mid-latitude zones (15°?<?ML?<?55°, and ??55°?<?ML < ??15°), and two auroral zones (70°?<?ML, and ML < ??70°). The Es climatology, namely its variations with each identified zone, altitude, season, and local time has been documented.     

Estimation of geometric route distance from its topological distance: application to narrow road networks in Tokyo

The structure of road networks has been investigated in accordance with the development of GIScience. By classifying road networks into wide and narrow ones, we can define the route as the path from the route’s origin (also called the root) on a wide road network to a narrow road segment which consists of the sequence of narrow road segments arranged by ascending order of the number of steps of adjacency to its root. The length of the route can be defined with the following geometric and topological terms: the route distance, measuring the length along the route and the depth, counting the number of road segments on the route. The depth plays the important role of being a substitute for the route distance in modelling road networks as a planar graph. Since road networks clearly exhibit irregular patterns and road segment lengths are non-uniform, it is considered appropriate to adopt a stochastic approach rather than a deterministic one to analyse the route distance. However, the relationship between the route distance and its depth has not been sufficiently investigated stochastically. Therefore, the research question is how can we estimate the route distance from its depth? Based on an empirical study in the Tokyo metropolitan region, it was found that (1) the statistical distribution of the route distance can be formulated as an Erlang distribution whose parameters are its depth and the inverse of the mean length of narrow road segments, and (2) this length is constant and close to 40 m. Therefore, we can estimate the route distance from only one parameter, the depth. Also, as a practical application, accessibility to the kth depth link in terms of firefighting was evaluated because the maximum length of the extension of fire hoses is approximately 200 m. It was found that (1) even if k?≤?5, the probability that the route distance to the kth depth link is equal to or longer than 200 m ranges from 0 to 0.45; and (2) if k?≥?8, the probability is approximately 1. These indicate the limitation of the deterministic approach because, on the basis of complete grid patterns (with intervals of 40 m between intersections), k?=?5 corresponds to a distance of 200 m from wide road networks and the route to the 5th depth link can be covered with fire hoses. Moreover, it was found that the connectivity of wide road networks is higher than that of narrow ones in terms of the smaller ratio of cul-de-sacs and the larger ratio of four-way intersections. These answers contribute substantially not only to constructing a science of cities that provides a simple model and specifies the most important parameter, but also to our understanding of the structure of narrow road networks within several hundred metres of wide road networks.     

Correlation between ROTI and Ionospheric Scintillation Indices using Hong Kong low-latitude GPS data

The correlation between the rate of TEC index (ROTI) and scintillation indices S ₄ and σ _Φ for low-latitude region is analyzed in this study, using data collected from a Global Positioning System (GPS) scintillation monitoring receiver installed at the south of Hong Kong for the periods June–August of 2012 and May 2013 and July–December of 2013. The analysis indicates that the correlation coefficient between ROTI and S ₄/σ _Φ is about 0.6 if data from all GPS satellites are used together. If each individual satellite is considered, the correlation coefficients are above 0.6 on average and sometimes above 0.8. The analysis also shows that the ratio of ROTI and S ₄ varies between 1 and 4. The ratio ROTI/σ _Φ, varies between 2 and 9. In addition, it is also found that there is a good consistency between the temporal variations of ROTI with scintillation activity under different ionospheric conditions. ROTI has a high correlation relationship with scintillation indices on geomagnetically disturbed days or in solar active months. Moreover, the data observed at low elevation angles have weak correlation between ROTI and scintillation indices. These results demonstrate the feasibility of using ROTI derived from GPS observations recorded by common non-scintillation GPS receivers to characterize ionospheric scintillations.     

Multi-Index Soil Moisture Estimation from Satellite Earth Observations: Comparative Evaluation of the Topographic Wetness Index (TWI), the Temperature Vegetation Dryness Index (TVDI) and the Improved TVDI (iTVDI)

Soil moisture estimation from satellite earth observation has emerged effectively advantageous due to the high temporal resolution, spatial resolution, coverage, and processing convenience it affords. In this paper, we present a study carried out to estimate soil moisture level at every location within Enugu State Nigeria from satellite earth observation. Comparative analysis of multiple indices for soil moisture estimation was carried out with a view to evaluating the robustness, correlation, appropriateness and accuracy of the indices in estimating the spatial distribution of soil moisture level in Enugu State. Results were correlated and validated with In-Situ soil moisture observations from multi-sample points. To achieve this, the Topographic Wetness Index (TWI), based on digital elevation data, the Temperature Vegetation Dryness Index (TVDI) and an improved TVDI (iTVDI) incorporating air temperature and a Digital Elevation Model (DEM) were calculated from ASTER global DEM and Landsat images. Possible dependencies of the indices on land cover type, topography, and precipitation were explored. In-Situ soil moisture data were used to validate the derived indices. The results showed that there was a positive significant relationship between iTVDI versus TVDI (R = 0.53, P value < 0.05), while in iTVDI versus TWI (R = 0.00, P value > 0.05) and TVDI versus TWI (R = ?0.01, P value > 0.05) no significant relationship existed. There was a strong relationship between iTVDI and topography, land cover type, and precipitation than other indices (TVDI, TWI). In situ measured soil moisture values showed negative significant relationship with TVDI (R = ?0.52, P value < 0.05) and iTVDI (R = ?0.63, P value < 0.05) but not with TWI (R = ?0.10, P value > 0.05). The iTVDI outperformed the other two index; having a stronger relationship with topography, precipitation, land cover classes and soil moisture. It concludes that although iTVDI outperformed other indices (TVDI, TWI) in soil moisture estimation, the decision of which index to apply is dependent on available data, the intent of usage and spatial scale.     

An Effective Model to Retrieve Soil Moisture from L- and C-Band SAR Data

This study investigated an appropriate method for soil moisture retrieval from radar images and coincident ground measurements acquired over bare soil and sparsely vegetated regions. The adopted approach based on a single scattering integral equation method (IEM) was developed to establish the relationship between backscatter coefficient and surface soil parameters including volumetric soil moisture content and surface roughness. The performance of IEM in 0–7.6 cm is better than that in 0–20 cm. Moreover, IEM can simulate correctly the backscatter coefficients only for the root mean square (RMS) height s < 1.5 cm at C-band and s < 2.5 cm at L-band by using an exponential correlation function and for s > 1.5 cm at C-band and s > 2.5 cm at L-band by using Gaussian function. However, due to the difficulties involved in the parameterization of soil surface roughness, the estimated accuracy is not satisfactory for the inversion of IEM. This paper used a combined roughness parameter and Fresnel reflection coefficient to develop an empirical model. Simulations were performed to support experimental results and to highlight soil moisture content and surface roughness effects in different polarizations. Results showed that a good agreement was found between the IEM simulations and the SAR measurements over a wide range of soil moisture and surface roughness characteristics. The model had a significant operational advantage in soil moisture retrieval. The correlation coefficients were 77.03 % at L-band and 81.45 % at C-band with the RMSEs of 0.515 and 0.4996 dB, respectively. Additionally, this work offered insight into the required application accuracy of soil moisture retrieval at a large area of arid regions.     

Enriching the semantics of the directed polyline–polygon topological relationships: the $$\mathcal {}$$-intersection matrix

In this paper, we define an intersection matrix for enriching the semantics of the topological relationships between a directed polyline and a polygon. In particular, we propose the \(\mathcal {DLP}\)-intersection matrix which enables us to model the origin and destination points, as well as the right- and left-hand sides of the directed polyline. This matrix overcomes the limitation of the well-known DE-9IM, because it allows the representation of the different dimensions of the intersection results at the same time. Accordingly, the geo-operators have been revised and extended in order to address the notions of right- and left-hand sides of a directed polyline, as well as additional notions related to the orientation of the polyline. The \(\mathcal {DLP}\)-intersection matrix has been implemented by extending the Java Topology Suite methods in order to address the new geo-operators based on the notion of orientation.     

Pasture enclosures increase soil carbon dioxide flux rate in Semiarid Rangeland,Kenya

Background

Pasture enclosures play an important role in rehabilitating the degraded soils and vegetation, and may also influence the emission of key greenhouse gasses (GHGs) from the soil. However, no study in East Africa and in Kenya has conducted direct measurements of GHG fluxes following the restoration of degraded communal grazing lands through the establishment of pasture enclosures. A field experiment was conducted in northwestern Kenya to measure the emission of CO₂, CH₄ and N₂O from soil under two pasture restoration systems; grazing dominated enclosure (GDE) and contractual grazing enclosure (CGE), and in the adjacent open grazing rangeland (OGR) as control. Herbaceous vegetation cover, biomass production, and surface (0–10 cm) soil organic carbon (SOC) were also assessed to determine their relationship with the GHG flux rate.

Results

Vegetation cover was higher enclosure systems and ranged from 20.7% in OGR to 40.2% in GDE while aboveground biomass increased from 72.0 kg DM ha^?1 in OGR to 483.1 and 560.4 kg DM ha^?1 in CGE and GDE respectively. The SOC concentration in GDE and CGE increased by an average of 27% relative to OGR and ranged between 4.4 g kg^?1 and 6.6 g kg^?1. The mean emission rates across the grazing systems were 18.6 μg N m^?2 h^?1, 50.1 μg C m^?2 h^?1 and 199.7 mg C m^?2 h^?1 for N₂O, CH₄, and CO₂, respectively. Soil CO₂ emission was considerably higher in GDE and CGE systems than in OGR (P?<?0.001). However, non-significantly higher CH₄ and N₂O emissions were observed in GDE and CGE compared to OGR (P?=?0.33 and 0.53 for CH₄ and N₂O, respectively). Soil moisture exhibited a significant positive relationship with CO₂, CH₄, and N₂O, implying that it is the key factor influencing the flux rate of GHGs in the area.

Conclusions

The results demonstrated that the establishment of enclosures in tropical rangelands is a valuable intervention for improving pasture production and restoration of surface soil properties. However, a long-term study is required to evaluate the patterns in annual CO₂, N₂O, CH₄ fluxes from soils and determine the ecosystem carbon balance across the pastoral landscape.

     

A first modeling of dynamic and static crustal strain field from near-field dilatation measurements: example of the 2013 $$M_w$$ 6.2 Ruisui earthquake,Taiwan

We analyze the high-resolution dilatation data for the October 2013 \(M_w\) 6.2 Ruisui, Taiwan, earthquake, which occurred at a distance of 15–20 km away from a Sacks–Evertson dilatometer network. Based on well-constrained source parameters (\(\hbox {strike}=217^\circ \), \(\hbox {dip}=48^\circ \), \(\hbox {rake}=49^\circ \)), we propose a simple rupture model that explains the permanent static deformation and the dynamic vibrations at short period (\(\sim \)3.5–4.5 s) for most of the four sites with less than 20 % of discrepancies. This study represents a first attempt of modeling simultaneously the dynamic and static crustal strain using dilatation data. The results illustrate the potential for strain recordings of high-frequency seismic waves in the near-field of an earthquake to add constraints on the properties of seismic sources.     

An a priori solar radiation pressure model for the QZSS Michibiki satellite

It has been noted that the satellite laser ranging (SLR) residuals of the Quasi-Zenith Satellite System (QZSS) Michibiki satellite orbits show very marked dependence on the elevation angle of the Sun above the orbital plane (i.e., the \(\beta \) angle). It is well recognized that the systematic error is caused by mismodeling of the solar radiation pressure (SRP). Although the error can be reduced by the updated ECOM SRP model, the orbit error is still very large when the satellite switches to orbit-normal (ON) orientation. In this study, an a priori SRP model was established for the QZSS Michibiki satellite to enhance the ECOM model. This model is expressed in ECOM’s D, Y, and B axes (DYB) using seven parameters for the yaw-steering (YS) mode, and additional three parameters are used to compensate the remaining modeling deficiencies, particularly the perturbations in the Y axis, based on a redefined DYB for the ON mode. With the proposed a priori model, QZSS Michibiki’s precise orbits over 21 months were determined. SLR validation indicated that the systematic \(\beta \)-angle-dependent error was reduced when the satellite was in the YS mode, and better than an 8-cm root mean square (RMS) was achieved. More importantly, the orbit quality was also improved significantly when the satellite was in the ON mode. Relative to ECOM and adjustable box-wing model, the proposed SRP model showed the best performance in the ON mode, and the RMS of the SLR residuals was better than 15 cm, which was a two times improvement over the ECOM without a priori model used, but was still two times worse than the YS mode.