Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.     

基于EMD与神经网络的机械故障诊断技术

经验模式分解 (EMD)是分析非线性、非平稳信号的有力工具 ,它将信号分解为突出了原信号的不同时间尺度的局部特征信息的内在模函数 (IMF)分量。本文通过将各 IMF分量输入到 BP网络中进行训练学习和故障诊断 ,比直接输入原信号可以提高 BP网络对故障诊断的准确率 ,而且减少了训练时间。     

Transfer-function modelling between environmental variation and mesozooplankton in the Baltic Sea

Time series of freshwater runoff, seawater salinity, temperature and oxygen were used in transfer functions (TF) to model changes of mesozooplankton taxa in the Baltic Sea from the 1960’s to the 1990’s. The models were then compared with long term zooplankton monitoring data from the same period. The TF models for all taxa over the whole Baltic proper and at different depth layers showed statistically significant estimates in t-tests. TF models were further compared using parsimony as a criterion. We present models showing 1) r² > 0.4, 2) the smallest residual standard error with the combination of exploratory variables, 3) the lowest number of parameters and 4) the highest proportional decrease in error term when the TF model residual standard error was compared with those of the univariate ARIMA model of the same response variable. Most often (7 taxa out of a total of 8), zooplankton taxa were dependent on freshwater runoff and/or seawater salinity. Cladocerans and estuarine copepods were more conveniently modelled through the inclusion of seawater temperature and oxygen data as independent variables. Our modelling, however, explains neither the overall increase in zooplankton abundance nor a simultaneous decrease found in the neritic copepod, Temora longicornis. Therefore, biotic controlling agents (e.g. nutrients, primary production and planktivore diets) are suggested as independent variables for further TF modelling. TF modelling enabled us to put the controlling factors in a time frame. It was then possible, despite the inherent multiple correlation among parameters studied to deduce a chain-of-events from the environmental controls and biotic feedback mechanisms to changes in zooplankton species. We suggest that the documented long-term changes in zooplankton could have been driven by climatic regulation only. The control by climate could be mediated to zooplankton through marine chemical and physical factors, as well as biotic factors if all of these were responding to the same external control, such as changes in the freshwater runoff. Increased runoff would explain both the increasing eutrophication, causing the overall increase of zooplankton, and the changes in selective predation, contributing to decline of Temora.     

Shoreline variability via empirical orthogonal function analysis: Part I temporal and spatial characteristics

Empirical orthogonal functions (EOFs) or principal components were used to extract the significant modes of shoreline variability from several data sets collected at three very different locations. Although EOFs have proven to be a valuable tool in the analysis of nearshore data, most applications have focused on the ability of the technique to describe cross-shore or profile variability. Here however, EOFs were used to help identify the dominant modes of longshore shoreline variability at Duck, North Carolina, the Gold Coast, Australia, and at several locations within the Columbia River Littoral Cell in the U.S. Pacific Northwest. In part one of this analysis, characteristic patterns of shoreline variability identified by the EOF analysis are described in detail. At each site, the dominant modes consisting of the first four eigenfunctions were found to describe nearly 95% of the total shoreline variability. At both Duck and the Gold Coast, several interesting longshore periodic features suggestive of sand waves were identified, while boundary effects related to natural headlands and navigational structures/entrances dominated the Pacific Northwest data sets.     

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。     

中国西部地区基于经验模型的土壤水分反演

地表土壤水分含量的时空分布信息是十分重要的,常常作为水文模型、气候模型、生态模型的输入参数,同时,也是干旱预报、农作物估产等工作的重要指标。被动微波遥感是监测土壤含水量最有效的手段之一。相比红外与可见光,它具有波长长,穿透能力强的优势。相比主动微波雷达,被动微波辐射计具有监测面积大、周期短,受粗糙度影响小,对土壤水分更为敏感,算法更为成熟的优势。目前,已研究出许多反演土壤水分的方法.本课题的主要内容是借助AMSR-E土壤水分影像数据、MODIS归一化植被指数(NDVI)影像数据和MODIS分类影像数据,利用ENVI软件进行遥感图像数据处理,运用统计分析方法建立NDVI与土壤水分的经验模型,研究中国西部地区稀疏植被覆盖区土壤水分的反演。     

高位瓦斯抽放钻孔在降低综采工作面瓦斯浓度中的应用

林南仓属于低瓦斯矿井,但存在高瓦斯区域。煤层和采空区是瓦斯的主要来源,尤以采空涌出量大,给煤矿生产和安全带来了极大隐患。通过在1129综采工作面风道施工高位瓦斯孔,把钻孔打到采空区一侧煤层顶板以上冒落裂隙带内,用钻孔进行瓦斯抽放,使采空内的瓦斯通过裂隙带沿钻孔抽出,有效降低综采工作面瓦斯浓度,保证综采工作面正常回采和安全生产。     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

Reconstructing hydrological variability from testate amoebae analysis in Carpathian peatlands

Peatlands offer the potential for high resolution records of water balance over Holocene timescales, yet this potential is under-exploited in many areas of the world. Within Europe, peatlands are mostly confined to areas north of 55° N, but several areas of southern and eastern Europe contain small peatlands which may be suitable for palaeoclimatic reconstruction. In this paper we test the potential of peatlands in the Carpathian region for deriving quantified estimates of water table changes using testate amoebae analysis. A training set for palaeohydrological reconstruction from testate amoebae assemblages was obtained by collecting surface samples from 13 peatlands, including 9 from Hungary and 4 from Transylvania (Romania). Using a simple measure of mean annual water tables estimated from staining of PVC tape, we found that some peatlands were heavily influenced by runoff and groundwater, and were therefore not suitable as modern analogues of ombrotrophic climatically sensitive sites. The relationship between the testate amoebae assemblages in the modern samples and the environmental variables was explored using CCA. The CCA biplot showed that the most important variables are depth to water table and moisture content, confirming that hydrology is a key control on taxon distribution. pH was a secondary gradient. A transfer function for % moisture and depth to water table was established and applied to fossil assemblages from a sequence from Fenyves-tető, Transylvania, Romania. The reconstructed water table shows a number of variations which have parallels with other palaeoclimatic records from Europe and the North Atlantic prominent phases of higher water tables are associated with the periods 8000–8300 cal BP, 3000–2500 cal BP and after 600 cal BP. We suggest that these were periods of particular intensification of westerly airflow which affected eastern Europe as well as western and central Europe.