Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Lateral continuity of basement seismic reflections in 15 Ma ultrafast-spreading crust at ODP Site 1256

The Ocean Drilling Program (ODP) initiated drilling at Site 1256D in the Guatemala Basin, about 1,000 km off the East Pacific Rise to penetrate plutonic rocks, anticipated to be relatively shallow in this region, formed at an ultra-fast spreading rate. IODP Expedition E312 successfully drilled into gabbros at ~1,150 m in basement. Multi-channel seismic traces show weak laterally coherent sub-basement reflections at borehole depths. Synthetic reflectivity seismograms were computed using a Ricker wavelet and impedance profiles from borehole sonic logs. These seismograms show significant sub-basement amplitude peaks. A zero-offset vertical seismic profile, shot on E312, was processed to investigate the authenticity of these reflections and their relationship to borehole geology. A dual scheme of the median filtering and F–K dip filtering was used. Tests with synthetic seismograms indicate the approach is effective at reasonable SNR levels. Downgoing energy is clearly identified but negligible upgoing energy is visible over random noise. These results indicate that lava flows and igneous contacts in upper ocean crust have significant topography on lateral scales less than the Fresnel Zone (~300 m) due to igneous and tectonic processes.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Tectonics of southwestern Mexico,isotopic evidence,nuclear Central America,Late Cretaceous break up

In the paleogeographic reconstruction of Mexico and northern Central America, an ever-increasing amount of evidence shows that the entire region is a collage of suspect terranes transported from abroad, whose timing and sense of motion are now beginning to be understood. Among these, the Chortis block (nuclear Central America) and the Baja California Peninsula have been proposed as pieces of continent separated from the Pacific coast of southwestern Mexico, that have moved either southeastward by the Farallon plate or northwestward by the Kula plate. Previous studies mainly confined to the northern margin of the Chortis block, confirmed a left-lateral displacement of 130 km in Neogene time. Further studies made northwestward along the Mexican coast provided a better understanding of magmatic and metamorphic processes in the area, and suggested times of detachment increased to 30 Ma, 40 Ma, and 66 Ma. The pre-detachment westernmost position of the block has changed, depending on the model chosen, from Puerto Vallarta and beyond, to the current position. Here we show that the isotopic mineral ages from coastal granites along the coast from Puerto Vallarta, Jalisco (80 Ma) to Puerto Angel, Oaxaca (11 Ma) record systematic decrease of cooling ages from NW to SE. This pattern is interpreted to result from the progressive uplift of rocks exposed at the present-day coast in that direction, such uplift occurred in response to the development of the Middle America Trench at the newly formed continental margin when the Chortis block was sliding at an average rate of 1.5 cm/year in a sinistral sense to its present position. Our results also constrain the position of the Kula-Farallon spreading axis north of Puerto Vallarta. These observations led us to conclude that several indicators point to this time and region for the onset of strike-slip drifting of the Chortis block toward its current position. Here, we also present several view points in terms of other possilble interpretations to different tectonic, geologic and isotopic data sets published recently by different authors.     

The Upper Middle Rhine Valley as a risk area

The Upper Middle Rhine Valley, granted the status of a World Heritage site, is well known for its unique inner narrow valley of Quaternary age with its historical legacy of numerous medieval castles and old towns. Less known is that this has always been a risk area of floods and gravitative mass movements. Up to the recent past, mainly ice floods caused enormous damage. The inhabitants of the valley were well aware that they lived in a risk area, but they had learned to handle the flood hazard. With the demise of ice floods over the last 40 years, due to climate change and because of the additional heating of the river water by power plants, the awareness of flood hazards has been much diminished, in contrast to that of potential damage by rockfalls and landslides which were also much feared in the past, though at the local level only. Still in the people’s memory is the Kaub catastrophe of March 10, 1876, when 28 persons were killed by a landslide. Nowadays, even minor rockfalls are a major threat, as they will affect the much-used traffic lines on both banks of the river, in particular the railroads. Therefore, since 2002, on behalf of German Rail (Deutsche Bahn, DB), all problematic slopes have been protected by costly steel-ring nets, although they are an aesthetic problem by UNESCO standards. The feeling of absolute safety created among the public is only subjective, though, as planners are well aware of. Moreover, the impact of modern climate change on slope stability is nearly unknown. Therefore, it is still necessary to develop a risk map for the narrow valley, with emphasis on gravitational hazards.     

A Long-Term Decrease of the Mid-Size Segmentation Lengths Observed in the He <Emphasis Type="SmallCaps">ii</Emphasis> (30.4 nm) Solar EUV Emission

Power spectra of segmentation-cell length (a dominant length scale of EUV emission in the transition region) from full-disk He?ii extreme ultraviolet (EUV) images observed by the Extreme ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SOHO) and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) during periods of quiet-Sun conditions for a time interval from 1996 to 2015 were analyzed. The spatial power as a function of the spatial frequency from about 0.04 to 0.27 (EIT) or up to 0.48 (AIA) Mm^?1 depends on the distribution of the observed segmentation-cell dimensions – a structure of the solar EUV network. The temporal variations of the spatial power reported by Didkovsky and Gurman (Solar Phys. 289, 153, 2014) were suggested as decreases at the mid-spatial frequencies for the compared spectra when the power curves at the highest spatial frequencies of 0.5 pix^?1 were adjusted to match each other. This approach has been extended in this work to compare spectral ratios at high spatial frequencies expressed in the solar spatial frequency units of Mm^?1. A model of EIT and AIA spatial responses allowed us to directly compare spatial spectral ratios at high spatial frequencies for five years of joint operation of EIT and AIA, from 2010 to 2015. Based on this approach, we represent these ratio changes as a long-term network transformation that may be interpreted as a continuous dissipation of mid-size network structures to the smaller-size structures in the transition region. In contrast to expected cycling of the segmentation-cell dimension structures and associated spatial power in the spectra with the solar cycle, the spectra demonstrate a significant and steady change of the EUV network. The temporal trend across these structural spectra is not critically sensitive to any long-term instrumental changes, e.g. degradation of sensitivity, but to the change of the segmentation-cell dimensions of the EUV network structure.     

Climatic and geographic factors affect ecosystem multifunctionality through biodiversity in the Tibetan alpine grasslands

Ecosystem multifunctionality(EMF), the simultaneous provision of multiple ecosystem functions, is often affected by biodiversity and environmental factors. We know little about how the interactions between biodiversity and environmental factors affect EMF. In this case study, a structural equation model was used to clarify climatic and geographic pathways that affect EMF by varying biodiversity in the Tibetan alpine grasslands. In addition to services related to carbon, nitrogen, and water cycling, forage supply, which is related to plantproductivity and palatability, was included in the EMF index. The results showed that 72% of the variation in EMF could be explained by biodiversity and other environmental factors. The ratio of palatable richness to all species richness explained 8.3% of the EMF variation. We found that air temperature, elevation, and latitude all affected EMF, but in different ways. Air temperature and elevation impacted the aboveground parts of the ecosystem, which included plant height, aboveground biomass, richness of palatable species, and ratio of palatable richness to all species richness. Latitude affected EMF by varying both aboveground and belowground parts of the ecosystem, which included palatable speciesrichness and belowground biomass. Our results indicated that there are still uncertainties in the biodiversity–EMF relationships related to the variable components of EMF, and climatic and geographic factors. Clarification of pathways that affect EMF using structural equation modeling techniques could elucidate the mechanisms by which environmental changes affect EMF.     

Artificial light at night: implications for early life stages development in four temperate freshwater fish species

Flora and fauna have both evolved under a natural cycle of light and dark. But especially in urban areas, the night is now increasingly disturbed by artificial light. Many traits and behaviours in fish are triggered by a circadian clock, for example hatching and swim bladder inflation, which predominantly take place at dusk or night. As lighting becomes brighter and extends farther into rural areas, the distinction between day and night becomes increasingly blurred. Therefore, the loss of diurnal trigger by artificial light at night was hypothesized having deleterious effects on these traits and impact fish reproduction. To assess these effects, eggs of four native freshwater fishes, Eurasian perch Perca fluviatilis, roach Rutilus rutilus, bleak Alburnus alburnus and chub Leuciscus cephalus, were incubated under two different light conditions: a photoperiod of 14 h light:10 h darkness (LD) and continuous illumination (LL). The time to hatch and swim bladder inflation was recorded. The species showed inconsistent reactions to the light treatments. In roach and bleak, the time to 50% hatch was longer in LL, whereas continuous lighting had an accelerating effect in chub. Incubation in LL elongated the hatching period in perch and roach and, in perch, the onset of darkness seemed to trigger hatching. The swim bladder inflation was significantly promoted by continuous light in chub and bleak but was not affected in roach. In conclusion, nocturnal artificial illumination could have an effect on hatching and initial swim bladder filling by masking the day–night-change and thereby diminish the trigger effect. However, the reactions were species specific and the increase in variation indicated a lack of diurnal triggering, whilst a general deleterious effect of artificial light at night has not been identified on early life stages.     

海水声速剖面基函数的重构性能分析

海水声速剖面通常使用经验正交函数（Empirical Orthogonal Function,EOF）进行稀疏表示,然而基函数会受到数据完备性和数据测量时间的制约,其代表性误差会导致声速剖面重构精度受限。为了提高声速剖面的重构精度,本文利用模糊C均值聚类对BOA_Argo历史数据集进行聚类分析,探讨不同聚类空间的训练集数据对实测声速剖面重构精度的影响。研究表明,声速剖面具有明显的时空聚集特性,聚类后的历史声速剖面集生成的基函数和平均声速剖面具有最优的重构性能。本文研究结果有助于为历史声速剖面训练集的选取提供实际指导意义,进而提高声速剖面重构精度乃至反演精度。     

岩体中泥化夹层的流变试验及其长期强度的确定

文中首先介绍两种室内试验设备--单剪流变仪和剪应力松弛仪（松弛流变仪）及其试验方法,然后根根岩体中泥化夹层的试验成果,证实陈宗基教授提出的最高屈服值f3--强度的临界值即长期强度值,工程设计要重视它。