Modeling fish growth and reproduction in the context of the Dynamic Energy Budget theory to predict environmental impact on anchovy spawning duration

Spawning location and timing are critical for understanding fish larval survival. The impact of a changing environment on spawning patterns is, however, poorly understood. A novel approach is to consider the impact of the environment on individual life histories and subsequent spawnings. In the present work, we extend the Dynamic Energy Budget (DEB) theory to investigate how environment variability impacts the spawning timing and duration of a multiple-batch spawning species. The model is successfully applied to reproduce the growth and reproduction of anchovy (Engraulis encrasicolus) in the Bay of Biscay. The model captures realistically the start and ending of the spawning season, including the timing of the spawning events, and the change in egg number per batch. Using a realistic seasonal forcing of temperature and food availability derived from a bio-physical model, our simulation results show that two thirds of the total spawned mass already accumulates before the start of the spawning season and that the condition factor increases with body length. These simulation results are in accordance with previous estimations and observations on growth and reproduction of anchovy. Furthermore, we show how individuals of equal length can differ in reproductive performance according to the environmental conditions they encounter prior to the spawning season. Hatch date turns out to be key for fecundity at age-1 as it partly controls the ability to build up reserves allocated to reproduction. We suggest the model can be used to realistically predict spawning in spatially and temporally varying environments and provide initial conditions for bio-physical models used to predict larval survival.     

Geological structure of the Curonian Spit (of the Baltic Sea) and its evolution history (revised)

The origin and evolution history of the Curonian Spit in the Baltic Sea are discussed based on geological and geomorphologic data. Evidence of the correlation between the spit formation stages and the Holocene sea level oscillations is presented. The ledges of moraine basement and marine accumulative features formed a single barrier. Its subsequent transformation occurred under the influence of the wave-related lithodynamic and eolian processes.     

Meso-beta scale numerical simulations of terrain drag-induced along-stream circulations. Part II: Concentration of potential vorticity within dryline bulges

Summary Hydrostatic and nonhydrostatic simulation models are employed to study the intensification of a terrain drag-induced dryline. The study develops a multi-stage theory for the evolution of the dryline including the concentration of potential vorticity accompanying meso-gamma scale dryline bulges.The numerical simulations indicate three fundamental stages of dryline intensification all of which are either directly or indirectly a result of the terrain-drag on the mid/upper-tropospheric jet stream by the Front Range of the Colorado Rocky Mountains. The first stage involves the downward momentum flux accompanying a large amplitude hydrostatic mountain wave which induces a downslope windstorm along the lee slopes. The surge of momentum (i.e., the dry, warm air associated with the downslope windstorm) propagates down the leeslope and modifies an existing weak dryline boundary. As the downslope windstorm initiates an undular bore along the lee slopes, the high momentum gradient which propagates downstream accompanying the bore, as well as the strong lower tropospheric sinking motions ahead of the bore, contract the scale of the surface moisture boundary between the dry air from above the leeslope and the moist air over the High Plains. This process further strengthens the dryline.The second stage involves the coupling of the terrain drag-induced along-stream ageostrophic front within the midtroposphere to the boundary layer through a thermally-indirect circulation. As the along-stream ageostrophic circulation intensifies within the middle troposphere down-stream from the mountain wave, sinking air parcels originating above 40 kPa descend to below 60 kPa over the High Plains where surface pressures are, only 85 kPa. These descending air parcels within the upstream branch of the along-stream ageostrophic thermally-indirect circulation contain high values of momentum and very low dewpoint values. As the planetary boundary layer (PBL) deepens due to surface warming during the morning hours, momentum and dry air from the midtropospheric along-stream ageostrophic front are entrained into the PBL. This process amplifies the bore-induced hydrostatic dryline bulge via low-level ageostrophic confluence.Finally, regions of low Richardson number (arising from strong vertical shears) within the amplifying midtropospheric along-stream ageostrophic thermally-indirect circulation become preferred regions for the development of non-hydrostatic evanescent internal gravity waves. These waves are embedded within the hydrostatic along-stream front above the low-level dryline and are accomapanied by very significant values of vertical momentum flux which act to focus the meso-gamma scale structure of the dryline into smaller scale bulges where low-level winds and vorticities are very high. This meso-gamma scale process follows the hydrostatic tilting and vortex tube stretching which creates meso-beta scale maxima of mid-lower tropospheric vorticity. The turbulent momentum fluxes accompanying wavebreaking within the nonhydrostatic dryline bulge create very large (i.e., stratospheric values of) potential vorticity near 70 kPa due to the nonconservation of potential vorticity on isentropic surfaces.With 30 Figures     

A comparative study of the performance of various vertical discretization schemes

Summary Most finite-difference numerical weather prediction models employ vertical discretizations that are staggered, and are low-order (usually second-order) approximations for the important terms such as the derivation of the geopotential from the hydrostatic equation, and the calculation of the vertically integrated divergence. In a sigma-coordinate model the latter is used for computing both the surface pressure change and the vertical velocity. All of the above-mentioned variables can diminish the accuracy of the forecast if they are not calculated accurately, and can have an impact on related quantities such as precipitation.In this study various discretization schemes in the vertical are compared both in theory and in practice. Four different vertical grids are tested: one unstaggered and three staggered (including the widely-used Lorenz grid). The comparison is carried out by assessing the accuracy of the grids using vertical numerics that range from second-order up to sixth-order.The theoretical part of the study examines how faithfully each vertical grid reproduces the vertical modes of the governing equations linearized with a basic state atmosphere. The performance of the grids is evaluated for 2nd, 4th and 6th-order numerical schemes based on Lagrange polynomials, and for a 6th-ordercompact scheme.Our interpretation of the results of the theoretical study is as follows. The most important result is that the order of accuracy employed in the numerics seems to be more significant than the choice of vertical grid. There are differences between the grids at second-order, but these differences effectively vanish as the order of accuracy increases. The sixth-order schemes all produce very accurate results with the grids performing equally well, and with the compact scheme significantly outperforming the Lagrange scheme. A second major result is that for the number of levels typically used in current operational forecast models, second-order schemes (which are used almost universally) all appear to be relatively poor, for other than the lowest modes.The theoretical claims were confirmed in practice using a large number (100) of forecasts with the Australian Bureau of Meteorology Research Centre's operational model. By comparing test model forecasts using the four grids and the different orders of numerics with very high resolution control model forecasts, the results of the theoretical study seem to be corroborated.With 8 Figures     

岩石地球化学图解信息分析辅助系统的研究与开发

在岩石地球化学理论的框架下,用该领域专家在实际工作中的思维方法来描述岩石地球化学数据处理与分析的固有流程,以岩石地球化学参数和图解为应用研究对象,利用岩石地球化学与计算机软件技术交叉的研究方法,概括岩石地球化学数据的特征及其处理方法。通过对双变量图解、三变量图解特征的分析,将数据计算分为主量和微量元素相关参数计算、主量和微量元素标准矿物计算、同位素相关计算,设计合理、有效的算法,实现岩石地球化学图解的数据管理、分类、参数选择与计算、表达式识别、绘制、坐标转换、投点等图解成图全过程的自动化。     

山区谷底沉陷预测模型及其参数反演

我国西部地区地势复杂,沟壑纵横,地下开采极易导致边坡失稳,引发采动滑坡。在地下采动沉降与滑坡体挤压上升的叠加影响下,谷底区域地表沉降值明显小于类似地质采矿条件下的平原地区。为准确预测山区谷底区域地表沉降值,基于简支梁的弹性变形理论,并借助概率密度函数建立了山区谷底区域地表沉陷预计修正模型,明确模型参数物理意义及其取值方法。依据修正模型,以实测值和预测值之差平方和最小为原则构建适应值函数,基于模拟退火粒子群算法提出新的模型参数反演方法,借助MATLAB语言编制了相应的参数反演程序。最后将研究成果应用于山西某矿,得到谷底区域预测结果中误差为73 mm,与实测值基本一致,取得了较好的工程实践效果。     

Flood‐based analysis of high‐magnitude sediment transport using a non‐parametric method

Upland erosion and the resulting reservoir siltation is a serious issue in the Isábena catchment (445 km² Central Spanish Pyrenees). During a three‐month period, water and sediment fluxes have been monitored at the catchment outlet (Capella), two adjacent subcatchments (Villacarli, 41 km²; Cabecera, 145 km²) and the elementary badland catchment Torrelaribera (8 ha). This paper presents the results of the monitoring, a method for the calculation of a sedigraph from intermittent measurements and the derived sediment yields at the monitored locations. The observed suspended sediment concentrations (SSCs) demonstrate the role of badlands as sediment sources: SSCs of up to 280 g l^?1 were encountered for Villacarli, which includes large badland areas. SSCs at the Cabecera catchment, with great areas of woodland, barely exceeded 30 g l^?1. SSCs directly at the sediment source (Torrelaribera) were comparable to those at Villacarli, suggesting a close connection within this subcatchment. At Capella, SSCs of up to 99 g l^?1 were observed. For all sites, SSC displayed only a loose correlation with discharge, inhibiting the application of a simple sediment rating curve. Instead, ancillary variables acting as driving forces or proxies for the processes (rainfall energy, cumulative discharge, rising/falling limb data) were included in a quantile regression forest model to explain the variability in SSC. The variables with most predictive power vary between the sites, suggesting the predominance of different processes. The subsequent flood‐based calculation of sediment yields attests high specific sediment yields for Torrelaribera and Villacarli (6277 and 1971 t km^?2) and medium to high yields for Cabecera and Capella (139 and 410 t km^?2) during the observation period. In all catchments, most of the sediment was exported during intense storms of late summer. Later flood events yield successively less sediment. Relating upland sediment production to yield at the outlet suggests considerable effects of sediment storage within the river channel. Copyright © 2008 John Wiley & Sons, Ltd.     

Different styles of metasomatic veining in ultramafic xenoliths from the TUBAF Seamount (Bismarck Microplate,Papua New Guinea)

Petrologic, geochemical and isotopic investigations on two ultramafic xenoliths with metasomatic veins from the TUBAF Seamount in the Bismarck Archipelago NE of Papua New Guinea reveal different styles of metasomatic overprinting. The first xenolith, a clinopyroxene–poor spinel lherzolite, was part of the depleted upper mantle. It contains an orthopyroxene-rich vein that formed by hydrous metasomatism at ~ 980 °C and ~ 1.5 GPa. The second xenolith is a clinopyroxene-dominated spinel olivine websterite that formed as a magmatic cumulate at the transition of the upper mantle to the oceanic crust. The websterite contains a vein with orthopyroxenes and clinopyroxenes, which give evidence for high-temperature crystallization at ~ 1300 °C and < 0.36 GPa. Both xenoliths were transported to the seafloor by a Quaternary trachybasalt in a fore-arc position. The vein minerals show a strong affinity to a supra-subduction zone or island arc setting. The REE pattern of the vein in the clinopyroxene–poor lherzolite strongly resembles the one from the host trachybasalt, with a high enrichment of the LREE and a strong to moderate enrichment of the MREE and HREE. Although broadly similar in shape, the REE pattern of the vein in the websterite shows a much weaker enrichment. The same applies to the trace-element patterns, although there are significant differences in the Eu, Zr, Hf and Nb concentrations. The isotope signatures of both veins suggest a derivation from a subducted slab that had been hydrothermally altered by seawater (high ⁸⁷Sr/⁸⁶Sr values).The contrasting crystallization temperatures of the vein minerals as well as their overall geochemical differences indicate that the metasomatic agents responsible for the vein in the websterite were mobilized from a previously depleted source at a much deeper mantle level than those forming the vein of the clinopyroxene–poor lherzolite. The metasomatic agents may also have been mobilized at different times and from different plates, i.e., the deeply subducted Solomon Sea Microplate (for the veins in the websterite) and the shallow dehydrating Pacific Plate (for the veins in the clinopyroxene–poor lherzolite).Metasomatic agents responsible for similar petrologic phenomena, i.e., modal or cryptic metasomatism, may have distinctly different origins and show contrasting histories. A strongly depleted lherzolite may totally lose its initial geochemical signature by the influence of an enriched metasomatic agent, whereas a primarily enriched ultramafic rock, e.g., a websterite, may strongly obscure the trace-element pattern of a less enriched metasomatic vein. Furthermore, the geochemistry of the ultramafic xenoliths may reflect polyphase cryptic and modal metasomatism related to veining and later transport by the hosting melt to the seafloor.     

The oxygen isotope composition, petrology and geochemistry of mare basalts: Evidence for large-scale compositional variation in the lunar mantle

To investigate the formation and early evolution of the lunar mantle and crust we have analysed the oxygen isotopic composition, titanium content and modal mineralogy of a suite of lunar basalts. Our sample set included eight low-Ti basalts from the Apollo 12 and 15 collections, and 12 high-Ti basalts from Apollo 11 and 17 collections. In addition, we have determined the oxygen isotopic composition of an Apollo 15 KREEP (K - potassium, REE - Rare Earth Element, and P - phosphorus) basalt (sample 15386) and an Apollo 14 feldspathic mare basalt (sample 14053). Our data display a continuum in bulk-rock δ¹⁸O values, from relatively low values in the most Ti-rich samples to higher values in the Ti-poor samples, with the Apollo 11 sample suite partially bridging the gap. Calculation of bulk-rock δ¹⁸O values, using a combination of previously published oxygen isotope data on mineral separates from lunar basalts, and modal mineralogy (determined in this study), match with the measured bulk-rock δ¹⁸O values. This demonstrates that differences in mineral modal assemblage produce differences in mare basalt δ¹⁸O bulk-rock values. Differences between the low- and high-Ti mare basalts appear to be largely a reflection of mantle-source heterogeneities, and in particular, the highly variable distribution of ilmenite within the lunar mantle. Bulk δ¹⁸O variation in mare basalts is also controlled by fractional crystallisation of a few key mineral phases. Thus, ilmenite fractionation is important in the case of high-Ti Apollo 17 samples, whereas olivine plays a more dominant role for the low-Ti Apollo 12 samples.Consistent with the results of previous studies, our data reveal no detectable difference between the Δ¹⁷O of the Earth and Moon. The fact that oxygen three-isotope studies have been unable to detect a measurable difference at such high precisions reinforces doubts about the giant impact hypothesis as presently formulated.