Unification of New Zealand’s local vertical datums: iterative gravimetric quasigeoid computations

New Zealand uses 13 separate local vertical datums (LVDs) based on geodetic levelling from 12 different tide-gauges. We describe their unification using a regional gravimetric quasigeoid model and GPS-levelling data on each LVD. A novel application of iterative quasigeoid computation is used, where the LVD offsets computed from earlier models are used to apply additional gravity reductions from each LVD to that model. The solution converges after only three iterations yielding LVD offsets ranging from 0.24 to 0.58 m with an average standard deviation of ±0.08 m. The so-computed LVD offsets agree, within expected data errors, with geodetically levelled height differences at common benchmarks between adjacent LVDs. This shows that iterated quasigeoid models have a role in vertical datum unification.     

GPS/IMU辅助机载三线阵影像光束法平差

建立了GPS/IMU观测值与外方位元素之间的数学关系,引入线性多项式、分段多项式和定向片内插3种数学模型用于机载三线阵影像的单航带及区域网光束法平差。ADS40数据的试验结果表明,在少量控制点的支持下,3种数学模型的单航带平差均可取得优于1.1个像素的定位精度,其中定向片内插模型最适合用于三线阵影像的空中三角测量。区域网平差的结果进一步证明,四角布设平高控制点即可得到最优的平差结果。     

基于球面的电子全站仪测高大气竖直折射角严密计算方法研究

大气竖直折射是影响竖直角观测精度的关键因素之一.大气竖直折射的计算一直是困扰测绘界的难题.为了更加合理地计算大气竖直折射,提出了基于球面的电子全站仪测高大气竖直折射计算方法,给出了基于球面的电子全站仪测高严格计算公式.     

Vertical distribution of primary producers and zooplankton grazers during different phases of the Arctic spring bloom

During a circumnavigation of the Svalbard archipelago in May 2006, simultaneous marine environmental (meteorology, heat flux, ocean turbulence, irradiance) and biological (phytoplankton and zooplankton biomass/species) data were sampled at selected stations. The zooplankton data were supplemented by high-resolution, high-speed VPR sampling down to 100 m depth at most stations. We were able to sample different phases of the phytoplankton spring bloom in Arctic as well as in Atlantic waters, and the stations represented different situations with respect to irradiance, turbulence and water-column stability. Phytoplankton growth and depth distribution were physically controlled, while zooplankton distributions were affected by biological parameters and turbulence. Development of the zooplankton followed the phytoplankton bloom phase, which was progressing in a direction from west to east in the waters north of Svalbard, and southwards in the Barents Sea. Our results also showed that the zooplankton did not avoid Phaeocystis pouchetii colonies, which have earlier been described as toxic. Despite an early retreat of the ice this year there was no apparent mismatch between the phytoplankton bloom and the dominant mesozooplankton, Calanus spp.     

Derivation of unified wave overtopping formulas for seawalls with smooth,impermeable surfaces based on selected CLASH datasets

A set of unified formulas for prediction of the mean rate of wave overtopping at coastal structures with smooth, impermeable surfaces have been derived through the analysis of the selected CLASH datasets. The mean wave overtopping rate is expressed as the function of the significant wave height at the structural toe and the relative freeboard. The formulas are applicable for both vertical walls and inclined seawalls with smooth transition between them. The formulas are simple but cover the full range of water depth from the shoreline to deep water. The effects of the toe depth and the seabed slope on wave overtopping rate are duly incorporated in the formulas. Prediction performance of the new formulas is better than the EurOtop formulas for both vertical walls and inclined seawalls.     

Field and laboratory measurements of mean overtopping discharges and spatial distributions at vertical seawalls

Within the CLASH project, wave overtopping at the vertical seawall at Samphire Hoe was measured by HR Wallingford (HRW), and compared laboratory tests in 2 & 3 dimensions carried out at the University of Edinburgh and HRW. At Samphire Hoe, overtopping volumes were captured in three volumetric tanks capable of measuring wave-by-wave and total overtopping volumes. The three tanks were placed progressively farther back from the seawall edge so that the spatial distribution of the overtopping discharges could be determined. The field measurement equipment was successfully deployed on three occasions, and measured overtopping discharges ranged from that barely considered to be hazardous to the public to over q = 3.0 l/s/m. The 2d testing at Edinburgh was modelled at a scale of 1:40, and the 3d model at HRW was modelled at 1:20. For both sets of laboratory tests, a range of conditions, representative of the storm wave conditions and water levels, was reproduced in addition to a set of parametric conditions. The storm conditions allowed a direct comparison between the field and laboratory measurements, and the parametric conditions were used to test the generic overtopping behaviour of the structure. For both sets of laboratory tests, mean overtopping discharges and the spatial distribution were measured separately. Analysis of the distribution data relates the proportion of the discharge that has landed as a function of (x / L_o); where x is the distance behind the crest, and L_o is the offshore wavelength. Analysis of the field, 2d & 3d laboratory data, and empirical prediction methods have not identified any scale effects for overtopping discharges at vertical and near-vertical seawalls.     

Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983

Particulate organic carbon (POC) is vertically transported to the oceanic interior by aggregates and their ballasts, mainly CaCO₃ and biogenic opal, with a smaller role for lithogenic aerosols through the mesopelagic zone. Diel migrating zooplankton communities effect vertical transport and remineralization of POC in the upper layers of the ocean. Below 1.5 km, the presence of zooplankton is reduced and thus the aggregates travel mainly by gravitational transport. We normalized the fluxes of POC, CaCO₃, and biogenic opal from data published on samples collected at 134 globally distributed, bottom-tethered, time-series sediment trap (TS-trap) stations to annual mole fluxes at the mesopelagic/bathypelagic boundary (m/b) at 2 km and defined them as F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio. Using this global data set, we investigated (1) the geographic contrasts of POC export at m/b and (2) the supply rate of ∑CO₂ to the world mesopelagic water column. F_m/bC_org varies from 25 (Pacific Warm Pool) to 605 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; F_m/bC_inorg varies from >8 (high latitude Polar Oceans) or 15 (Pacific Warm Pool) to 459 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; and F_m/bSi_bio, the most spatially/temporally variable flux, ranges from 6 (North Atlantic Drift) to 1118 (Pacific Subarctic Gyre) mmolSi m⁻² yr⁻¹. The oceanic region exhibiting the highest POC flux over a significantly large region is the area of the North Pacific Boreal Gyres where the average F_m/bC_org = 213, F_m/bC_inorg = 126, and F_m/bSi_bio = 578 mmol m⁻² yr⁻¹. F_m/bC_org and F_m/bC_inorg are particularly high in large upwelling margins, including the divergent Arabian Sea and off Cape Verde. One of the data sets showing the lowest flux over a significant region/basin is F_m/bC_org = 39, F_m/bC_inorg = 69, and F_m/bSi_bio = 22 mmol m⁻² yr⁻¹ in the North Pacific subtropical/tropical gyres; Pan-Atlantic average fluxes are similar except F_m/bSi_bio fluxes are even lower. Where C_org/C_inorg and Si_bio/C_inorg are <1 defines the “Carbonate Ocean”, and where these ratios are ?1 defines the “Silica Ocean”. The Carbonate Ocean occupies about 80% of the present world pelagic ocean between the two major oceanographic fronts, the North Pacific Polar Front and the Antarctic Polar Front, and the Silica Ocean is found on the polar sides of these fronts. The total global annual fluxes of F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio at m/b calculated by parameterizations of the export flux data from 134 stations are surprisingly similar; 36.2, 33.8, and 34.6 teramol yr⁻¹ (120, 112, and 114 mmol m⁻² yr⁻¹), respectively, resulting in a near uniform binary ratio between the above three elements of about one. The global ternary % ratios estimated from 152 TS-trap samples of the three elements are 35:32:33. From our global F_m/bC_org and a published model estimate of the global export production, we estimate the regeneration rate of CO₂ through the mesopelagic zone by the biological pump is 441 teramolC yr⁻¹. Based on our global F_m/bC_inorg and recently estimated global primary production of PIC, 36-86 teramolC yr⁻¹ of PIC is assumed to be dissolved within the upper 2 km of the water column.     

Sea ice–ocean coupling using a rescaled vertical coordinate z

Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate “z^*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of “z^*” avoids the difficult issue of vanishing levels under thick ice.Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z^* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z^* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model.     

2017年8月8日九寨沟7.0级地震前地磁垂直强度极化短期异常特征分析

对2017年8月8日九寨沟7.0级地震前后四川及其邻区地磁台站数据进行垂直强度极化方法计算,使用5～100 s周期发现震前地磁极化计算结果存在如下的显著异常特征:①2017年3月震中附近地磁台站同步存在极化超阈值异常,时间上显示为短期异常;②同步异常台站出现最大幅度的时间较一致;③九寨沟地震发生在异常阈值线附近,对该方...     

大气结构特征与雷暴成因

盛夏季节,宜宾雷暴的形成主要与对流层内高温高湿、风及其在对流层的中上部垂直切变小的特点相关联。西宁雷暴的形成主要与上层为深厚的干冷、下层为浅薄的暖湿、对流层的中上层风及其垂直切变大的特点相关联。分析又得到以下结论:西宁雷暴的形成又与高空冷空气加强南移、副热带急流南侧反气旋式切变得到加强和高空出现辐散相关联。