Kilometer-scale fault-related thermal anomalies in tight gas sandstones

Upper Carboniferous sandstones make one of the most important tight gas reservoirs in Central Europe. This study integrates a variety of geothermometers (chlorite thermometry, fluid inclusion microthermometry and vitrinite reflection measurements) to characterize a thermal anomaly in a reservoir outcrop analog (Piesberg quarry, Lower Saxony Basin), which is assumed responsible for high temperatures of circa 300 °C, deteriorating reservoir quality entirely. The tight gas siliciclastics were overprinted with temperatures approximately 90–120 °C higher compared to outcropping rocks of a similar stratigraphic position some 15 km to the west. The local temperature increase can be explained by circulating hydrothermal fluids along the fault damage zone of a large NNW-SSE striking fault with a displacement of up to 600 m in the east of the quarry, laterally heating up the entire exposed tight gas sandstones. The km-scale lateral extent of this fault-bound thermal anomaly is evidenced by vitrinite reflectance measurements of meta-anthracite coals (VR_rot ∼ 4.66) and the temperature-related diagenetic overprint. Data suggest that this thermal event and the associated highest coalification was reached prior to peak subsidence during Late Jurassic rifting (162 Ma) based on K-Ar dating of the <2 μm fraction of the tight gas sandstones. Associated stable isotope data from fluid inclusions, hosted in a first fracture filling quartz generation (T ∼ 250 °C) close to lithostatic fluid pressure (P ∼ 1000 bars), together with authigenic chlorite growth in mineralized extension fractures, demonstrate that coalification was not subject to significant changes during ongoing burial. This is further evidenced by the biaxial reflectance anisotropy of meta-anthracite coals. A second event of quartz vein formation occurred at lower temperatures (T ∼ 180 °C) and lower (hydrostatic) pressure conditions (P ∼ 400 bars) and can be related to basin inversion. This second quartz generation might be associated with a second event of illite growth and K-Ar ages of 96.5–106.7 Ma derived from the <0.2 μm fraction of the tight gas sandstones.This study demonstrates the exploration risk of fault-bound thermal anomalies by deteriorating entirely the reservoir quality of tight gas sandstones with respect to porosity and permeability due to the cementation with temperature-related authigenic cements. It documents that peak temperatures are not necessarily associated with peak subsidence. Consequently, these phenomena need to be considered in petroleum system models to avoid, for example, overestimates of burial depth and reservoir quality.     

The mechanism and verification analysis of permafrost-associated gas hydrate formation in the Qilian Mountain,Northwest China

Muri Basin in the Qilian Mountain is the only permafrost area in China where gas hydrate samples have been obtained through scientific drilling. Fracture-filling hydrate is the main type of gas hydrate found in the Qilian Mountain permafrost. Most of gas hydrate samples had been found in a thin-layer-like, flake and block group in a fracture of Jurassic mudstone and oil shale, although some pore-filling hydrate was found in porous sandstone. The mechanism for gas hydrate formation in the Qilian Mountain permafrost is as follows: gas generation from source rock was controlled by tectonic subsidence and uplift--gas migration and accumulation was controlled by fault and tight formation--gas hydrate formation and accumulation was controlled by permafrost. Some control factors for gas hydrate formation in the Qilian Mountain permafrost were analyzed and validated through numerical analysis and laboratory experiments. CSMGem was used to estimate the gas hydrate stability zone in the Qilian permafrost at a depth of 100–400 m. This method was used to analyze the gas composition of gas hydrate to determine the gas composition before gas hydrate formation. When the overlying formation of gas accumulation zone had a permeability of 0.05 × 10⁻¹⁵ m² and water saturation of more than 0.8, gas from deep source rocks was sealed up to form the gas accumulation zone. Fracture-filling hydrate was formed in the overlap area of gas hydrate stability zone and gas accumulation zone. The experimental results showed that the lithology of reservoir played a key role in controlling the occurrence and distribution of gas hydrate in the Qilian Mountain permafrost.     

Modeling free gas content of the Lower Paleozoic shales in the Weiyuan area of the Sichuan Basin,China

Two sets of Lower Paleozoic organic-rich shales develop well in the Weiyuan area of the Sichuan Basin: the Lower Cambrian Jiulaodong shale and the Lower Silurian Longmaxi shale. The Weiyuan area underwent a strong subsidence during the Triassic to Early Cretaceous and followed by an extensive uplifting and erosion after the Late Cretaceous. This has brought about great changes to the temperature and pressure conditions of the shales, which is vitally important for the accumulation and preservation of shale gas. Based on the burial and thermal history, averaged TOC and porosity data, geological and geochemical models for the two sets of shales were established. Within each of the shale units, gas generation was modeled and the evolution of the free gas content was calculated using the PVTSim software. Results show that the free gas content in the Lower Cambrian and Lower Silurian shales in the studied area reached the maxima of 1.98–2.93 m³/t and 3.29–4.91 m³/t, respectively (under a pressure coefficient of 1.0–2.0) at their maximum burial. Subsequently, the free gas content continuously decreased as the shale was uplifted. At present, the free gas content in the two sets of shales is 1.52–2.43 m³/t and 1.94–3.42 m³/t, respectively (under a current pressure coefficient of 1.0–2.0). The results are roughly coincident with the gas content data obtained from in situ measurements in the Weiyuan area. We proposed that the Lower Cambrian and Lower Silurian shales have a shale gas potential, even though they have experienced a strong uplifting.     

基于客观分析的多源卫星叶绿素a浓度产品融合方法研究

本文基于Aqua/MODIS、Terra/MODIS和Envisat/MERIS多源卫星叶绿素a浓度产品,研究了客观分析融合方法,制作了西北太平洋海域(0°~50°N,100°~150°E)叶绿素a浓度融合产品,并从有效数据空间覆盖率和产品精度两个方面对融合方法进行了评价。与单传感器以及欧洲太空局发布的GSM模型业务化融合产品相比,客观分析融合产品空间覆盖率明显提高;与收集的2002-2012年间叶绿素a浓度实测数据比较,GSM模型业务化融合产品的匹配数据点为578个,偏差为-0.20 mg/m3,均方根误差为0.37 mg/m3,客观分析法融合产品的匹配数据点为1432个,偏差为-0.21 mg/m3,均方根误差为0.36 mg/m3。结果表明:本文研究的客观分析融合方法在保证融合产品精度的同时可显著提高产品的空间覆盖率,在海洋水色融合应用前景广阔。     

杭州湾沿海地区地面沉降时序InSAR监测与分析

Time-series InSAR analysis(e.g., permanent scatterers(PSInSAR)) has been proven as an effective technology in monitoring ground deformation over urban areas. However, it is a big challenge to apply this technology in coastal regions due to the lack of man-made targets. An distributed scatterers interferometric synthetic aperture radar(DSInSAR) is developed to solve the problem of insufficient samples and low reliability in monitoring coastal lowland subsidence, by applying a spatially adaptive filter and an eigendecomposition algorithm to estimating the optimal phase of statistically homogeneous distributed scatterers(DSs). Twenty-four scenes of COSMO-Sky Med images acquired between 2013 and 2015 are used to retrieve the land subsidence over the Shangyu District on south coast of the Hangzhou Bay, Zhejiang Province, China. The spatial pattern of the land subsidence obtained by the PS-InSAR and the DSInSAR coincides with each other, but the density of the DSs is three point five times higher than the permanent scatterers(PSs). Validated by precise levelling data over the same period, the DSInSAR method achieves an accuracy of ±5.0 mm/a which is superior to the PS-InSAR with±5.5 mm/a. The land subsidence in the Shangyu District is mainly distributed in the urban areas, industrial towns and land reclamation zones, with a maximum subsidence rate –30.2 mm/a. The analysis of geological data, field investigation and historical reclamation data indicates that human activities and natural compaction of reclamation material are major causes of the detected land subsidence. The results demonstrate that the DSInSAR method has a great potential in monitoring the coastal lowland subsidence and can be used to further investigate subsidence-related environmental issues in coastal regions.     

Raman spectroscopic measurements of synthetic gas hydrates in the ocean

A Raman spectrometer extensively modified for deep ocean use was used to measure synthetic hydrates formed in an ocean environment. This was the first time hydrates formed in the ocean have been measured in situ using Raman spectroscopy. Gas hydrates were formed in situ in the Monterey Bay by pressurizing a Pyrex cell with various gas mixtures. Raman spectra were obtained for sI methane hydrate and sII methane + ethane hydrate. Gas occlusion resulting from rapid gas growth of methane hydrate was measured immediately after formation. The Raman shift for methane free gas was coincident with that of methane in the small 5¹² hydrate cage. The methane Raman peak widths were used to discriminate between methane in the free gas and hydrate phase. Methane + ethane sII hydrate was formed for 43 days on the seafloor. In this case, gas occlusion was not measured when the gas hydrates were allowed to form over an extended time period. Equivalent Raman spectra were obtained for the in situ and laboratory-formed sII methane + ethane hydrates, under similar p, T, and x conditions. With the Raman spectrometer operating in the ocean, seawater contributes to the Raman spectra obtained. Both the Raman bands for the sulfate ion and water were used to qualitatively determine the distribution of water phases measured (hydrate, seawater) in the Raman spectra.     

Evidence for gas hydrate occurrences in the Canadian Arctic Beaufort Sea within permafrost-associated shelf and deep-water marine environments

The presence of a wedge of offshore permafrost on the shelf of the Canadian Beaufort Sea has been previously recognized and the consequence of a prolonged occurrence of such permafrost is the possibility of an underlying gas hydrate regime. We present the first evidence for wide-spread occurrences of gas hydrates across the shelf in water depths of 60–100 m using 3D and 2D multichannel seismic (MCS) data. A reflection with a polarity opposite to the seafloor was identified ∼1000 m below the seafloor that mimics some of the bottom-simulating reflections (BSRs) in marine gas hydrate regimes. However, the reflection is not truly bottom-simulating, as its depth is controlled by offshore permafrost. The depth of the reflection decreases with increasing water depth, as predicted from thermal modeling of the late Wisconsin transgression. The reflection crosscuts strata and defines a zone of enhanced reflectivity beneath it, which originates from free gas accumulated at the phase boundary over time as permafrost and associated gas hydrate stability zones thin in response to the transgression. The wide-spread gas hydrate occurrence beneath permafrost has implications on the region including drilling hazards associated with the presence of free gas, possible overpressure, lateral migration of fluids and expulsion at the seafloor. In contrast to the permafrost-associated gas hydrates, a deep-water marine BSR was also identified on MCS profiles. The MCS data show a polarity-reversed seismic reflection associated with a low-velocity zone beneath it. The seismic data coverage in the southern Beaufort Sea shows that the deep-water marine BSR is not uniformly present across the entire region. The regional discrepancy of the BSR occurrence between the US Alaska portion and the Mackenzie Delta region may be a result of high sedimentation rates expected for the central Mackenzie delta and high abundance of mass-transport deposits that prohibit gas to accumulate within and beneath the gas hydrate stability zone.     

Measurement of ocean surface winds using synthetic aperture radars

A methodology for retrieving high-resolution ocean surface wind fields from satellite-borne synthetic aperture radar (SAR) data is introduced and validated. The algorithms developed are suited for ocean SAR data, which were acquired at the C band of either vertical (VV) or horizontal (HH) polarization in transmission and reception. Wind directions are extracted from wind-induced streaks that are visible in SAR images of the ocean at horizontal scales greater than 200 m. These wind streaks are very well aligned with the mean surface wind direction. To extract the orientation of these streaks, two algorithms are introduced, which are applied either in the spatial or spectral domain. Ocean surface wind speeds are derived from the normalized radar cross section (NRCS) and image geometry of the calibrated SAR images, together with the local SAR-retrieved wind direction. Therefore, several C-band models (CMOD IFR2, CMOD4, and CMODS) are available, which were developed for VV polarization, and have to be extended for HH polarization. To compare the different algorithms and C-band models as well as demonstrate their applicability, SAR-retrieved wind fields are compared to numerical-model results considering advanced SAR (ASAR) data from Environmental Satellite (ENVISAT), a European satellite.     

Alteration and multi-stage accumulation of oil and gas in the Ordovician of the Tabei Uplift,Tarim Basin,NW China: Implications for genetic origin of the diverse hydrocarbons

The Ordovician is the most important exploration target in the Tabei Uplift of the Tarim Basin, which contains a range of petroleum types including solid bitumen, heavy oil, light oil, condensate, wet gas and dry gas. The density of the black oils ranges from 0.81 g/cm³ to 1.01 g/cm³ (20 °C) and gas oil ratio (GOR) ranges from 4 m³/m³ to 9300 m³/m³. Oil-source correlations established that most of the oils were derived from the Mid-Upper Ordovician marine shale and carbonate and that the difference in oil properties is mainly attributed to hydrocarbon alteration and multi-stage accumulation. In the Tabei Uplift, there were three main periods of hydrocarbon accumulation in the late Caledonian stage (ca. 450–430 Ma), late Hercynian stage (ca. 293–255 Ma) and the late Himalayan stage (ca. 12–2 Ma). The oil charging events mainly occurred in the late Caledonian and late Hercynian stage, while gas charging occurred in the late Hercynian stage. During the late Caledonian stage, petroleum charged the reservoirs lying east of the uplift. However, due to a crustal uplifting episode in the early Hercynian (ca. 386–372 Ma), most of the hydrocarbons were transformed by processes such as biodegradation, resulting in residual solid bitumen in the fractures of the reservoirs. During the late Hercynian Stage, a major episode of oil charging into Ordovician reservoirs took place. Subsequent crustal uplift and severe alteration by biodegradation in the west-central Basin resulted in heavy oil formation. Since the late Himalayan stage when rapid subsidence of the crust occurred, the oil residing in reservoirs was exposed to high temperature cracking conditions resulting in the production of gas and charged from the southeast further altering the pre-existing oils in the eastern reservoirs. A suite of representative samples of various crude oils including condensates, lights oils and heavy oils have been collected for detailed analysis to investigate the mechanism of formation. Based on the research it was concluded that the diversity of hydrocarbon physical and chemical properties in the Tabei Uplift was mainly attributable to the processes of biodegradation and gas washing. The understanding of the processes is very helpful to predict the spatial distribution of hydrocarbon in the Tabei Uplift and provides a reference case study for other areas.     

The hydroacoustic method for the quantification of the gas flux from a submersed bubble plume

This article presents an inverse hydroacoustic method for the remote quantification of the total gas flux transported from an underwater bubble plume. The method includes the surveying of the bubble plume by a vertically looking echo sounder and the calculation of the flux using the spatial distribution of the ultrasound backscattering at a fixed depth. A simplified parameterization containing only a few parameters is introduced to describe the empirical bubble size distribution. The linear correlation between the backscattering cross section of the bubble stream and the vertical gas flux is found. The calculation procedure takes into account the occurrence of a gas hydrate film at the bubble’s surface. The influence of different parameters on the accuracy of the method is investigated. The resolution volume of the echo sounder corresponding to the fixed distance is considered as a two-dimensional spatial window. The method was applied to quantify the total convective methane flux at the Haakon-Mosby mud volcano (HMMV) depth 1280 m. The calculated values of the total flux near the bottom (100–400 t/year) are in good agreement with the independently estimated flux for the single bubble jet observed from the ROV (70 t/year). These calculations also show significant temporal variability of the flux at the HMMV. The total flux was found to vary by about a factor of 2–3 within time scales of days.     

Seismic velocities on the Nova Scotian margin to estimate gas hydrate and free gas concentrations

This article provides new constraints on gas hydrate and free gas concentrations in the sediments at the margin off Nova Scotia. Two-dimensional (2-D) velocity models were constructed through simultaneous travel-time inversion of ocean-bottom seismometer (OBS) data and 2-D single-channel seismic (SCS) data acquired in two surveys, in 2004 and 2006. The surveys, separated by ∼5 km, were carried out in regions where the bottom-simulating reflection (BSR) was identified in seismic reflection datasets from earlier studies and address the question of whether the BSR is a good indicator of significant gas hydrate on the Scotian margin. For both datasets, velocity increases by 200–300 m/s at a depth of approximately 220 m below seafloor (mbsf), but the results of the 2006 survey show a smaller velocity decrease (50–80 m/s) at the base of this high-velocity layer (310–330 mbsf) than the results of the 2004 survey (130 m/s). When converted to gas hydrate concentrations using effective medium theory, the 2-D velocity models for both datasets show a gas hydrate layer of ∼100 m thickness above the identified BSR. Gas hydrate concentrations are estimated at approximately 2–10% for the 2006 data and 8–18% for the 2004 survey. The reduction in gas hydrate concentration relative to the distance from the Mohican Channel structure is most likely related to the low porosity within the mud-dominant sediment at the depth of the BSR. Free gas concentrations were calculated to be 1–2% of the sediment pore space for both datasets.     

Structure and oil and gas perspectives of the Californian borderland

The Californian borderland is the largest oil and gas region of the Pacific coast of the United States. Here, a series of large and minor sedimentary basins with established or inferred oil and gas bearing properties are recognized. The Californian borderland is located on a Nevadan-type active continental margin. The bodies of such margins are composed of giant accretion formations and occupy not only their underwater parts but also vast land areas. The accretion formations are dominated by rocks of a deep-water origin. The Californian borderland represents a system of basins and ranges that compose the underwater margin and the Coastal Ranges of California and the basement within the Great Valley. In the depressions of the accretion orogen, which are small in size but feature high subsidence rates, significant thicknesses of deposits are accumulated (up to 6?8 km). In the depressions of the Californian borderland, they are represented by young Neogene (rarer Eocene-Oligocene) formations of mainly terrigenous or siliceous-terrigenous compositions. The active tectonic regime resulted in a sharp reduction of the age range of the oil and gas bearing deposits, which are represented only by Late Cretaceous and Cenozoic rocks. Full-scale oil and gas production here is performed in the Santa Barbara-Ventura and Santa Maria basins. The principal oil and gas bearing unit is the Monterrey Formation of the Middle Miocene. By January 2006, the production on the Californian borderland comprised about 17.8 million tons of oil and over 1 billion m³ of gas. Up to 1.4 billion tons of oil and 200 billion m³ of gas reserves are regarded to be not recovered on the Pacific coast of the United States in the region of the Californian borderland.     

Effects of two offshore gas platforms on soft-bottom benthic communities (northwestern Adriatic Sea,Italy)

The macrozoobenthos living around two offshore gas platforms, Barbara NW (pB) and Calipso (pC) located in the northwestern Adriatic Sea were investigated for three years after their construction to detect eventual effects due to the platforms. The sampling stations were spaced at increasing distance from the platforms up to 1000 m. Both multivariate and univariate analysis showed an initial defaunation (short-term effect) at pB and within a 120 m radius at pC. A general recovery in terms of abundance, species richness and diversity was observed for the benthic communities after one year. During the third monitoring year a mussel mound developed at both the platforms (longer-term effect) extending up to 30 m from pC, whereas similar soft-bottom communities were found at all the other distances. The geographical position plays an important role in the “timing” of the above effects. Moreover, despite the different environmental features, a 3-year monitoring plan can be effective to evaluate the potential impact on benthic communities of offshore gas platforms in the north-central Adriatic Sea.     

Spatial and temporal variations in subsidence due to the natural consolidation and compaction of sediment in the yellow river delta,china

ABSTRACT

The Yellow River Delta, which is the second-largest delta in China, has experienced varying degrees of land subsidence since the late 1970s. Although recent studies have identified the natural consolidation and compaction of sediment among the most important contributors to geologic processes, their processes have rarely been quantified. We estimated the sediment compaction over different time ranges to determine the temporal evolution of subsidence parameters (i.e., cumulative compaction). Estimates of primary consolidation, secondary consolidation, and the degree of consolidation in 152 boreholes revealed the spatial–temporal characteristics of sediment compaction and consolidation using geotechnical parameters collected from 152 boreholes, soil mechanics equations and the Kriging interplolation method. In addition, these estimates were partially constrained and cross-validated using the interferometric synthetic aperture radar (InSAR) results from early 2007 to late 2010 which were provided in a previous study. By performing a comparison analysis between theoretical evaluations of compaction for borehole data and InSAR observations, we were able to quantify subsidence due to sediment compaction. The comparison results suggest that the theoretical solutions agreed well with the measurements recorded by the well-validated, advanced InSAR method and that the deviations between the InSAR technique and geotechnical evaluations ranged from ?22 to 3?mm. The results reveal that the land subsidence of the chosen borehole sites during the investigative period was dominated by the primary consolidation and compaction of sediment. The underprediction of subsidence may be explained by fluid withdrawal, oil exploitation and engineering construction. To speculate, more geological disasters may occur if the current subsidence condition extends into the future.     

Evaluation of ENVISAT ASAR data for sea surface wind retrieval in Hong Kong coastal waters of China

The C-band wind speed retrieval models, CMOD4, CMOD - IFR2, and CMOD5 were applied to retrieval of sea surface wind speeds from ENVISAT （European environmental satellite） ASAR （advanced synthetic aperture radar） data in the coastal waters near Hong Kong during a period from October 2005 to July 2007. The retrieved wind speeds are evaluated by comparing with buoy measurements and the QuikSCAT （quick scatterometer） wind products. The results show that the CMOD4 model gives the best performance at wind speeds lower than 15 m/s. The correlation coefficients with buoy and QuikSCAT winds are 0.781 and 0.896, respectively. The root mean square errors are the same 1.74 m/s. Namely, the CMOD4 model is the best one for sea surface wind speed retrieval from ASAR data in the coastal waters near Hong Kong.     

Differential selection for vegetative structure between juvenile brown shrimp (Penaeus aztecus) and white shrimp (P. setiferus), and implications in predator-prey relationships

The distributions of juvenile brown shrimp, Penaeus aztecus, and white shrimp, P. setiferus, (57–72 mm TL) were monitored in partially vegetated laboratory enclosures to examine selection for vegetative structure. Simulated Spartina culms, composed of green straws, were used as structure in the experiments. Brown shrimp selected for this structure during the day but not at night. White shrimp showed no strong selective preference for structure during the day or night.Selection in the field was measured by examining densities of the two species of shrimp (40–80 mm TL) in a Galveston Bay, Texas, salt marsh. Brown shrimp were more abundant in vegetated habitats (Spartina alterniflora) than in adjacent non-vegetated areas. Differences in white shrimp densities between the two habitats were variable, and no consistent pattern was observed.laboratory predation experiments were conducted in non-vegetated and vegetated cages using Atlantic croaker (Micropogonias undulatus) as predators. In non-vegetated cages there was no significant selection for either species of shrimp by the fish. In partially and completely vegetated cages, fish selected white shrimp over brown shrimp. These data suggest that the presence of vegetation alters prey selection by Atlantic croaker, and this interaction may result in greater mortality of white shrimp compared with brown shrimp in salt marshes.     

Spatial Variation of Sea Level Trend Along the Bangladesh Coast

The Bangladesh coast is threatened by rising sea level due to various factors. The results based on the analysis of past 22 years of tidal data of the Bangladesh coast reveal that the annual mean tidal level in the eastern Bangladesh coast is rising at an alarmingly high rate of 7.8 mm/year, which is almost twice the observed rate in the western region. This type of sea level trend seems to be the result of changing local conditions like increased precipitation and land subsidence during the recent decades. It seems that the higher rate of land subsidence in the eastern Bangladesh coast is the main causative factor for the steeper sea level trends there. The differential sea level trends show that the subsidence component in the sea level rise may be as high as 4 mm/year in the eastern Bangladesh coast. However, this needs to be verified with actual geological observations.     

Spatial Variation of Sea Level Trend Along the Bangladesh Coast

The Bangladesh coast is threatened by rising sea level due to various factors. The results based on the analysis of past 22 years of tidal data of the Bangladesh coast reveal that the annual mean tidal level in the eastern Bangladesh coast is rising at an alarmingly high rate of 7.8 mm/year, which is almost twice the observed rate in the western region. This type of sea level trend seems to be the result of changing local conditions like increased precipitation and land subsidence during the recent decades. It seems that the higher rate of land subsidence in the eastern Bangladesh coast is the main causative factor for the steeper sea level trends there. The differential sea level trends show that the subsidence component in the sea level rise may be as high as 4 mm/year in the eastern Bangladesh coast. However, this needs to be verified with actual geological observations.     

Recent sea level trends and accelerations: Comparison of tide gauge and satellite results

A comprehensive set of 456 monthly tide gauge records is analyzed for trend and acceleration over the same period that satellite altimetry was analyzed (1993 to 2011). Additionally, a 90 tide gauge record subset is analyzed for which GPS data are available. The selection criterion for the tide gauge data is 85% data completion. All measurements are adjusted for vertical land motion. Results from 456 pairs of tide gauges, adjusted for Global Isostatic Adjustment, and satellite recordings located within 1° root-mean-square latitude and longitude separation differences are compared. The tide gauge trends and accelerations are adjusted for spatial bias using the more globally dense satellite data.The average trends of the 456 and 90 gauge sets (3.26 and 2.68 mm/year, respectively) agree reasonably well with the global trend average of the satellite data (3.09 mm/year). Average trends for the 456 tide gauges are also in good agreement (within 95% confidence limits) with trends based on satellite data within the 1° satellite proximity criterion (3.26 and 3.31 mm/year, respectively). The trends for the 90 gauges with GPS nearby and qualifying satellite locations are 2.68 and 2.74 mm/year, respectively. For all datasets analyzed, the accelerations are quite strongly negative but the uncertainty is relatively large. Adjustment of the tide gauge trends for spatial bias modified both trends and accelerations significantly and decreased trend differences between the 456 and 90 gauge datasets. The spatially adjusted tide gauge trends (2.95 and 2.72 mm/year, respectively for the 456 and 90 tide gauges sets) are somewhat less than the 1° spatially adjusted satellite data (3.09 mm/year). Whether the increased sea level trend of approximately 3 mm/year measured by the satellites since the 1990's is a long-term increase from the 20th Century value of approximately 1.7 mm/year or part of a cycle will require longer records; however, the negative accelerations support some cyclic character.     

海?气界面气体交换速率与二氧化碳气体通量的估算

海?气界面CO₂通量的估算采用块体公式,其等于气体交换速率、CO₂溶解度以及海水与大气的CO₂分压差的乘积,其中的气体交换速率通常与风速相联系,不同作者提出了气体交换速率为风速不同幂次多项式的参数化方案。本文对比了气体交换速率为风速函数的主要研究结果,发现与风速多项式的依赖关系相比,观测数据所基于的观测方法对于气体交换速率的影响更大。在此基础上,本文用多种不同的气体交换速率参数化公式计算了1982?2018年全球的CO₂通量,海洋整体上是大气CO₂的汇,赤道海区是源,南北半球40°附近的海域构成沿纬向的强吸收带。37 a间,海洋CO₂通量的年平均值(以碳计)为(?1.53±0.15) Pg/a, 1999年前,海洋吸收量逐年减小,1999年达到最小值,之后海洋吸收量开始增大,海洋吸收量的增大主要发生在南大洋。