集成多视航空影像与LiDAR数据重建3维建筑物模型

集成高分辨率多视航空影像与LiDAR数据,以"轮廓提取—3维轮廓生成—3维模型重建"为主要框架,提出一种多视轮廓与LiDAR数据集成的3维轮廓线生成方法、LiDAR数据支撑下改进的分割-合并-成型算法,形成一套新的、系统的3维建筑物模型重建技术方案。该方案可以有效重建结构复杂的3维建筑物模型,相对单纯基于LiDAR的方案,所重建3维模型的准确性与完整性更好、定位精度更高。

航空遥感中基于高阶INS误差模型的GPS/INS组合定位定向方法

基于随机常数、随机游走和一阶马尔可夫过程组合的高阶IMU误差模型,建立36阶卡尔曼滤波器。通过松散组合模式,最终实现位置标准差±5cm、俯仰/横滚角标准差±0.002°、航向角标准差±0.008°的定位定向精度。

A possible interrelation between Earth rotation and climatic variability at decadal time-scale

Using multichannel singular spectrum analysis （MSSA） we decomposed climatic time se- ries into principal components, and compared them with Earth rotation parameters. The global warming trends were initially subtracted. Similar quasi 60 and 20 year periodic os- cillations have been found in the global mean Earth temperature anomaly （HadCRUT4） and global mean sea level （GMSL）. Similar cycles were also found in Earth rotation variation. Over the last 160 years multi-decadal change of Earth＇s rotation velocity is correlated with the 60-year temperature anomaly, and Chandler wobble envelope reproduces the form of the 60-year oscillation noticed in GMSL. The quasi 20-year oscillation observed in GMSL is correlated with the Chandler wobble excitation. So, we assume that Earth＇s rotation and climate indexes are connected. Despite of all the clues hinting this connection, no sound conclusion can be done as far as ocean circulation modelling is not able to correctly catch angular momentum of the oscillatory modes.

Is the global sea surface temperature rise accelerating?

This is an exploratory investigation to search for the presence of an acceleration in global sea surface temperature rise, which is essential to identify anthropogenic contributions to the climate change during the 20 th century. A weighted statistical model with an acceleration parameter was built progressively to reconstruct the variations in the global sea surface temperature data considering statistically significant confounders and autoregressive disturbances in the process. From the preliminary residual analysis of a weighted regression model, emerged a parsimonious model with first order autoregressive disturbances with a deterministic trend, acceleration and periodicity of 69 yr and its 138 yr subharmonic. The final model solution, selected from 29 alternative combinations of the model parameters using Mallows' s Cp metric, revealed a statistically significant deterministic trend, 0.40 ± 0.03C/c(p < 0.01), and acceleration, 0.67 ± 0.11C/c^2(p < 0.01) explaining 33% of the global sea surface temperature variations. The combined yearly trend and acceleration in global sea surface temperature as predicted by the model,exhibit a strong correlation with the yearly increase in the global CO^2 concentrations observed during the 20th century.

Multichannel singular spectrum analysis of the axial atmospheric angular momentum

Earth＇s variable rotation is mainly produced by the variability of the AAM（atmospheric angular momentum）. In particular, the axial AAM component X_3, which undergoes especially strong variations,induces changes in the Earth＇s rotation rate. In this study we analysed maps of regional input into the effective axial AAM from 1948 through 2011 from NCEP/NCAR reanalysis. Global zonal circulation patterns related to the LOD（length of day） were described. We applied MSSA（Multichannel Singular Spectrum Analysis） jointly to the mass and motion components of AAM, which allowed us to extract annual, semiannual, 4-mo nth, quasi-biennial, 5-year, and low-frequency oscillations. PCs（Principal components） strongly related to ENSO（El Nino southern oscillation） were released. They can be used to study ENSO-induced changes in pressure and wind fields and their coupling to LOD. The PCs describing the trends have captured slow atmospheric circulation changes possibly related to climate variability.

The rigorous adjustment of the global mean sea level budget during 2005-2015

Global mean sea level budget is rigorously adjusted during the period 2005-2015.The emphasis is to provide the best estimates for the linear rates of changes(trends)of the global mean sea level budget components during this period subject to the constraint:Earth’s hydrosphere conserves water.The newly simultaneously adjusted trends of the budget components suggest a larger correction for the global mean sea level trend implicated by the other budget components’trends under the budget constraint.The simultaneous estimation of the linear trends of the budget components subject to the constraint for closure improves their uncertainties and enables a holistic assessment of the global mean sea budget,which has implications for future sea level science studies,including the future Intergovernmental Panel on Climate Change(IPCC)Assessment Reports,and the US Climate Assessment Reports.

Seismic stress perturbation and triggering patterns induced by the 2016 Central Italy earthquake sequences

Numerous shallow earthquakes, including 24 th August Amatrice, 26 th October Visso, and 30 th October Norcia earthquakes, ruptured the segments of Mount Vettore-Gorzano fault system in the central Apennines(Italy) in 2016. In order to investigate the stress perturbation and triggering patterns among the earthquake sequences, we introduce a more realistic nonplanar coseismic fault geometry model, which improve the rupture model by assimilating relocated aftershocks and the GPS observations. We adopt the seismic slip inversion program of the steepest descent method(SDM) to create the detailed coseismic rupture models and optimize Coulomb Failure Stress model by varying the coefficient of friction and received fault parameters. The results indicate that the nonplanar fault geometry model is more reflective of the deep slip of the coseismic rupture than planar model. As evidenced by the coseismic Coulomb stress changes caused by the three mainshocks at different depth slices, the stress loading mainly distributes on the active fault zones and the stress changes can well explain the spatial distribution of aftershocks. The first large Amatrice mainshock accelerates the occurrence of the Mw 5.9 Visso and Mw 6.6 Norcia earthquakes, with the positive stress changes at the hypocenter exceeding the stress triggering threshold(0.010×10^(6) Pa) and up to 0.015×10^(6) and 0.257×10^(6) Pa, respectively. Furthermore, the Mw 5.9 Visso earthquake as well encourages the occurrence of the Mw 6.6 Norcia event with the increased stress changes of 0.052×10^(6) Pa on the hypocenter. It is concluded that the stress transfer and accumulation play crucial roles on the linkage triggering mechanism among the mainshock-mainshock and mainshockaftershocks. Noteworthily, the cumulative stress changes on the southwest segment of the Norcia Fault(NF), the southeast parts of the Montereale Fault System(MFS) and Mount Gorzano Fault(MGF) of the main regions are up to(1.5~3.5) ×10^(6) Pa. The cumulative stress changes have not been released sufficiently by aftershocks, which may increase the seismic hazard in those regions.

Characterizing receiver clock behaviors onboard Low Earth Orbiters:A case study of GRACE satellites

Accurate estimation of clocks, for example for the Gravity Recovery And Climate Experiment(GRACE)twin-satellites, is a critical part of precise orbit determination(POD) that ensures temporal gravity inversion. Characterizing the periodic variations of the receiver clocks is critical for precise clock modeling and prediction. In this study, the receiver clock is estimated using two different POD procedure: kinematic and reduced-dynamic approaches. Choices and the number of orbital parameters estimated in POD process affect the clock estimates, e.g., there are 8895 and 34,560 total parameters in the reduced-dynamic and kinematic approaches, respectively. In the both cases, the periodic variations of GRACE receiver clock are mainly dominated by the GPS orbit period, as well as once-(1-pr) and twiceper-revolution(2-pr) effects. Here the 1-pr effect is coupled with the relativistic effect, resulting in a difficulty to separate both signals. The clock amplitudes caused by the GPS orbit period, 1-pr and 2-pr are about 0.1, 0.03 and 0.01 ns, respectively. The GPS orbit period is almost one order magnitude larger than the 1-and 2-pr effect. The 0.1-ns amplitude of the 12-h periodic variation is equivalent to a 3-cm error in range. Such a systematic error should be considered in the receiver clock modeling for both the improvement of positioning accuracy and the reduction of number of unknown parameters, if the precise point positioning(PPP) technique is used for the orbit determination of the GRACE.

Photogrammetric Survey of the Historical Church of Nossa Senhora do Rosário and São Benedito in Brazil

In the history of cities and society,architecture plays an important role in the materialization of buildings,in the ways of life and in the constitution of an era.Documenting architectural heritage is an essential task for the conservation,management and collection of knowledge,however it is a time consuming and costly task.Presently,the application of digital aerial and terrestrial photogrammetry has made this process simpler,more agile and of relatively low cost.In this context,the objective of this study was to analyze the visual and geometric quality of the 3D model,of architectural heritage,generated by photogrammetric techniques.To this end,topographic mapping(reference)and aerial and terrestrial photogrammetric surveys of the Church of Nossa Senhora do Rosário and São Benedito,in Estrela do Sul/MG,Brazil,which is a historical and cultural heritage of the country,were carried out.Adopted as a reference,the 3D topographic model was generated by the points collected with a total station.To generate the photogrammetric 3D model,aerial photos were obtained using a UAV(Unmanned Aerial Vehicle)and terrestrial photos,acquired at strategic points.In both cases the PowerShot SX260HS camera was used.The two georeferenced 3D models were analyzed visually,considering texture and completeness,and compared geometrically through coordinates of homologous points.Comparatively,the photogrammetric 3D model presented results with a wealth of details far superior to the topographic ones and a geometry with centimetric discrepancies,showing that the method is suitable for generating 3D models relating to architectural heritage.

Performance of GNSS positioning in PPP mode using MADOCA precise products

The Global Navigation Satellite System (GNSS) is widely utilized for accurate positioning.One commonly applied method to obtain precise coordinate estimates is by implementing the relative positioning in network mode.However,this approach can be complex and challenging.Fortunately,The Japan Aerospace Exploration Agency (JAXA) offers freely available satellite orbit and clock correction products called Multi-GNSS Advanced Demonstration Tool for Orbit and Clock Analysis (MADOCA),which can enhance positioning accuracy through the precise point positioning (PPP) method.This study focuses on evaluating PPP static mode positioning using MADOCA products and comparing the results with the highly precise relative positioning method.By analyzing a network of 20 GNSS stations in Indonesia,we found that the PPP method using MADOCA products provided favorable positioning estimates.The median discrepancies and the corresponding median absolute deviation (MAD) for easting,northing,and up components were estimated as 9±18 mm,10±9 mm,and 3±40 mm,respectively.These results indicate that PPP with MADOCA products can be a reliable alternative for establishing Indonesia's horizontal control networks,particularly for orders 0,1,2,and 3,and for a broad spectrum of geoscience monitoring activities.However,considerations such as epoch transformations and seismic activities should be taken into account for accurate positioning applications that comply with the definition of the national reference framework.

基于等几何分析的太阳帆板热结构耦合随机动力响应

航天器周期性地出入日光区和背光区,会受到太阳辐射而引起周期性热应力的影响.基于等几何分析法,研究了具有随机场参数的太阳能帆板在热结构耦合下的随机动力响应.首先基于非均匀有理B样条构造了帆板的几何模型,得到太阳能帆板的热传导微分方程和振动微分方程;利用Newmark-β法求解了热-结构耦合下的微分方程,得到结构温度场和随机动力响应;最后,研究了不同随机场参数对结构动力响应的影响.

Effects of Gaussian filter in processing GRACE data: Gravity rate of change at Lhasa,southern Tibet

In this paper, the spatial gravity distribution over Tibetan Plateau and the gravity rate of change at Lhasa for different Gaussian filter radii are computed using GRACE data. Results show that the estimate of the gravity rate of change is spatialradius-dependent of the Ganssian filter. The GRACE-estimated gravity rate of change agrees well with the surface measured one. In other words, the GRACE-estimated gravity rate of change has a limited value as that obtained by surface measurement when the spatial filter radius reaches zero. Then numerical simulations are made for different spatial radii of the Gaussian filter to investigate its behaviors when applied to surface signals. Results show that the estimate of a physical signal is filter-radius dependent. If the computing area is equal to or less than the mass area, especially for a uniformly distributed mass, the estimate gives an almost correct result, no matter what filter radius is used. The estimate has large error because of the signal leakage caused by harmonic truncation if the computing area is much bigger than the mass distribution （or inverse for a small mass anomaly）. If a mass anomaly is too small, it is difficult to recover it from space observation unless the filter radius is extremely small. If the computing point （or area） is outside the mass distribution, the estimated result is almost zero, particularly for small filter radii. These properties of the Gaussian filter are helpful in applying GRACE data in different geophysical problems with different spatial position and geometrical size. We further discuss physical sources causing the scalar gravity change at Lhasa. Discussions indicate that the gravity rate of change at Lhasa is not caused by the present-day ice melting （PDIM） （or Little Ice Age, LIA） effect because no ice melting occurs in Lhasa city and nearby. The gravity rate of change is attributable mainly to tectonic deformation associated with the Indian Plate collision. Simultaneous surface displacement, surface denudation, and GIA effects are not negligible.

Antarctic-wide annual ice flow maps from Landsat 8 imagery between 2013 and 2019

Ice velocity constitutes a key parameter for quantifying ice-sheet discharge rates and is thus crucial for improving the coupled models of the Antarctic ice sheet towards accurately predict its contribution to future global sea-level rise.However,in Antarctica,high-resolution and continuous ice velocity estimates remain elusive,which is key to unravel Antarctica’s present-day ice mass balance processes.Here,we present a suite of newly estimated Antarctic-wide,annually-sampled ice velocity products at 105-m grid-spacing observed by Landsat 8 optical images data.We first describe a procedure that can automatically calibrate and integrate ice displacement maps to generate Antarcticwide seamless ice velocity products.The annual ice velocity mosaics are assembled using a total of 250,000 displacement maps inferred from more than 80,000 Landsat 8 images acquired between December 2013 and April 2019.The new annual Antarctic ice velocity data product exhibits an improved quantification of near-decadal Antarctic-wide ice flow,and an opportunity to investigate ice dynamics at a higher spatial resolution and annual sampling,as compared to existing data products.Validation studies confirmed improved accuracy and consistency of this new data product,when compared with independently estimated optical and radar ice velocity data products,as well as in situ data.

Comparative analysis of sol-air temperature in typical open and semi-closed courtyard spaces

The building cooling load is calculated in most standards using the sol-air temperature(T_(sol-air)),which depends on the surface properties and the building configuration.This study presents an experimental and numerical analysis to find a relationship between the aspect ratio of the courtyard and the T_(sol-air) at different height levels.Besides,the influence of the surface reflection factor on the T_(sol-air) is analyzed.The results show that the T_(sol-air) inside the courtyard differs compared to the usual outside facade.This difference is increased by increasing the aspect ratio of the courtyard and can reach up to 60°C at low levels of the deep courtyards.It is also found that the courtyard geometry is much more important than the thermal behavior of the materials and the T_(sol-air) within the deep courtyard is less affected by the reflection factor of materials.