Calcium isotopic signatures of depleted mid-ocean ridge basalts from the northeastern Pacific

A number of high-temperature processes(e.g.,melt-rock reactions,metasomatism,partial melting)can produce signifi cant Ca isotopic fractionation and heterogeneity in the mantle,but the mechanism for such fractionation remains obscure.To investigate the eff ect of mantle partial melting on Ca isotopic fractionation,we reported high-precision Ca isotopic compositions of depleted mid-ocean ridge basalts(MORBs)from the East Pacifi c Rise and Ecuador Rift in the northeastern Pacifi c.Theδ44/40 Ca of these MORB samples exhibit a narrow variation from 0.84‰to 0.88‰with an average of 0.85‰±0.03‰,which are similar to those of reported MORBs(0.83‰±0.11‰)and back-arc basin basalts(BABBs,0.80‰±0.08‰)in literature,but are lower than the estimate value for the bulk silicate Earth(BSE,0.94‰±0.05‰).The lowδ44/40 Ca signatures of MORB samples in this study cannot be caused by fractional crystallization,since intermediate-mafi c diff erentiation has been demonstrated having only limited eff ects on Ca isotopic fractionation.Instead,the off set ofδ44/40 Ca between MORBs and the BSE is most likely produced by mantle partial melting.During this process,the light Ca isotopes are preferentially transferred to the melt,while the heavy ones tend to stay in the residue,which is consistent with the fact thatδ44/40 Ca of melt-depleted peridotites increases with partial melting in literature.The behavior of Ca isotopes during mantle partial melting is closely related to the inter-mineral(Cpx and Opx)Ca isotopic fractionation and melting mineral modes.Mantle partial melting is one of the common processes that can induce lowerδ44/40 Ca values in basalts and Ca isotopic heterogeneity in Earth’s mantle.

Magma Dynamics of Axial Melt Lens at Fast-Spreading Mid-Ocean Ridges

Multichannel seismic studies performed at fastspreading mid-ocean ridges revealed the presence of a thin(tens to hundreds of meters high), narrow(< 1-2 km wide) axial melt lens(AML) in the mid-crust, which is underlain by crystal/melt mush that is in turn laterally surrounded by a transition zone of mostly solidified material. In order to shed light on the complexity of magmatic and metamorphic processes ongoing within and at the roof of axial melt lenses, we have focused on the petrological and geochemical record provided by fossilized AMLs. Of particular significance is Hole 1256D in the equatorial Pacific drilled by the International Ocean Discovery Program(IODP), where for the first time, the transition between sheeted dikes and gabbros in intact fast-spreading crust was penetrated, providing a drill core with a more or less continuous record of the upper part of an AML(Teagle et al., 2006;Koepke et al., 2008). This can be regarded as rosetta stone to answer longstanding questions on the complex magmatic evolution within an AML, as well as on metamorphic and anatectic processes ongoing at the roof of a dynamic AML, rising upward in the midcrust as a consequence of a replenishment event. The plutonic rocks drilled from Hole 1256D consist of quartz-bearing gabbros, diorites and tonalites, which might represent the upper part of a fossilized AML. The gabbros and diorites are consistent with modeled products of MORB fractional crystallization, composed of mixed melt and cumulate in varying ratios. Modeled trace elements support a model in which the tonalites originated from low-degree partial melting of the sheeted dikes overlying the AML, rather than extreme fractional crystallization(Erdmann et al., 2015;Zhang et al., 2017a). Therefore, the upper part of AML, largely composed of low density and high-viscosity felsic magmas, may serve as a barrier to eruptible MORB melts in the lower part of AML. Zoning of apatites from three different lithologies, tonalites, diorites, and gabbros, is common and shows a consistent evolution trend with depletion in Cl and REEs from core to rim. The cores are usually homogenous in composition and interpreted as magmatic origin, whereas zones with lower Cl and REEs are disseminated with heterogeneous concentrations, indicating exchanges with hydrothermal fluids. The high-Cl apatite core indicates assimilation of high-Cl brines at a magmatic stage, which is interpreted as immiscibility product from cycling seawater-derived fluids at a high temperature(Zhang et al., 2017b). The variation of F/Cl and Br/Cl ratios of bull rocks may reflect the mixing between MORB magmas and seawater-derived fluids, crystallization of apatite and amphibole, and/or extraction of magmatic fluids(Zhang et al., 2017c).

Model Forming of Mesoarchean Gridino Mafic Dyke Swarm during Subduction "Mid-Ocean Ridge – Continent"

The high-pressure metamorphosed Gridino dyke swarm comprises a major group of Mesoarchean 2.87-2.82 Ga mafic dykes intruded within the Mesoarchean continental crust of the Kola craton(the Belomorian tectonic

Implications of the melting depth and temperature of the Atlantic mid-ocean ridge basalts

Mid-ocean ridge basalts(MORBs) are characterized by large variations in trace element compositions and isotopic ratios, which are difficult to be interpreted solely by using magmatic process such as partial melting of a peridotitic mantle and subsequently fractional crystallization. Geochemical diversity of MORBs have been attributed to large-scale heterogeneity within the underlying mantle, and the heterogeneity might have been caused by addition of recycled crustal component, subcontinental lithosphere, metasomatized lithosphere and outer core contribution. In this study, we investigated the MORBs along the Mid-Atlantic Ridge(MAR) by estimating the temperature and pressure of partial melting, and comprehensively comparing trace element and isotope ratios. The data for MORBs from areas close to mantle plumes show large variations. Mantle plumes can affect mid-oceanic ridges 1 400 km away, but plume effects did not cover all of the ridge segments, and those segments without plume effects did not have any abnormalities in temperature, trace element or isotope ratios.We ascribed the above phenomena to result from the shapes of the plume flow, which we categorized as "pipelike channels" and "pancake-like channels". The "pancake-like channels" plumes affected the ambient mantle nondirectionally, but the range of the mantle affected by the "pipe-like channels" plumes were selective. Element ratios of MORBs reveal that the mantle source of the MORBs along the MAR is highly heterogeneous. We suggest that most of source heterogeneities of the MORBs may be due to the presence of subducted slab and delaminated lower crust in the source. In addition, the plume that carried materials from the core-mantle boundary may affect some of the segments.

A Quantitative Method for Active Fault Migration Distance Assessment on both Sides of Mid-Ocean Ridges——Based on Multi-Beam Data

Fracture-fissure systems found at mid-ocean ridges are dominating conduits for the circulation of metallogenic fluid.Ascertaining the distribution area of active faults on both sides of mid-ocean ridges will provide a useful tool in the search for potential hydrothermal vents,thus guiding the exploration of modern seafloor sulfides.Considering the MidAtlantic Ridge 20°N–24°N(NMAR)and North Chile Rise(NCR)as examples,fault elements such as Fault Spacing(?S)and Fault Heave(?X)can be identified and quantitatively measured.The methods used include Fourier filtering of the multi-beam bathymetry data,in combination with measurements of the topographic slope,curvature,and slope aspect patterns.According to the Sequential Faulting Model of mid-ocean ridges,the maximal migration distance of an active fault on either side of mid-ocean ridges—that is,the distribution range of active faults—can be measured.Results show that the maximal migration distance of active faults at the NMAR is 0.76–1.01 km(the distance is larger at the center than at the ends of this segment),and at the NCR,the distribution range of active faults is 0.38–1.6 km.The migration distance of active faults on the two study areas is positively related to the axial variation of magma supply.In the NCR study area,where there is an abundant magma input,the number of faults within a certain distance is mainly affected by the variation of lithospheric thickness.Here a large range of faulting clearly corresponds to a high proportion of magmatism to seafloor spreading near mid-ocean ridges(M)value,and in the study area of the NMAR,there is insufficient magmatism,and the number of faults may be controlled by both lithospheric thickness and magma supply,leading to a less obvious positive correlation between the distribution range of active faults and M.

Molybdenum isotope composition of the upper mantle and its origin:insight from mid-ocean ridge basalt

The molybdenum(Mo)isotope system is pivotal in reconstructing marine redox changes throughout Earth’s history and has emerged as a promising tracer for igneous and metamorphic processes.Understanding its composition and variation across major geochemical reservoirs is essential for its application in investigating high-temperature processes.However,there is debate regarding theδ^(98/95)Mo value of the Earth’s mantle,with estimates ranging from sub-chondritic to super-chondritic values.Recent analyses of global mid-ocean ridge basalt(MORB)glasses revealed significantδ^(98/95)Mo variations attributed to mantle heterogeneity,proposing a two-component mixing model to explain the observed variation.Complementary studies confirmed the sub-chondriticδ^(98/95)Mo of the depleted upper mantle,suggesting remixing of subduction-modified oceanic crust as a plausible mechanism.These findings underscore the role of Mo isotopes as effective tracers for understanding dynamic processes associated with mantle-crustal recycling.

Singularity of lithosphere mass density over the mid-ocean ridges and implication on floor depth and heat flow

The relation of heat flow and floor depth across the mid-ocean ridges versus lithosphere age can be described by linear functions of square root of age according to plate thermal conductive Half Space Models(HSM).However,one of the long-standing problems of these classical models is the discrepancies between predicted and observed heat flow and floor depth for very young and very old lithosphere.There have been several recent attempts to overcome this problem:one model incorporates temperature-and pressure-dependent parameters and the second model includes an additional low-conductivity crustal layer or magma rich mantle layer(MRM).Alternatively,in the current paper,the ordinary density of lithosphere in the plate conductive models is substituted with a reduction of lithosphere density towards axis that features the irregularity and nonlinearity of plates across the mid-ocean ridges.A new model is formulated incorporating the new form of density for predicting both peak heat flow and floor depth.Simple solutions of power-law forms derived from the model can significantly improve the predicting results of heat flow and floor depth over the mid-ocean ridges.Several datasets in the literature were reutilized for model validation and comparison.These datasets include both earlier datasets used for original model calibration and the more recently compiled high-quality datasets with both sedimentary and crustal loading corrections.The results indicate that both the heat flow and the slope(first orderderivative)of sea floor approach infinity(undifferentiability or singularities)around the mid-ocean ridges.These singularities are partially due to the boundary condition as it has been already known in the literature and partially to the reduction of density of lithosphere as discovered for the first time in the current research.

Hydrographic Characteristics and Oceanic Heat Flux in the Upper Arctic Ocean over the Alpha Ridge Observed by the DTOP Platform in 2018 and 2021

In 2018 and 2021,the Drift-Towing Ocean Profilers(DTOP)provided extensive temperature and salinity data on the upper 120m ocean through their drifts over the Alpha Ridge north of the Canada Basin.The thickness and temperature maximum of Alaska Coastal Water(ACW)ranged from 20m to 40m and-1.5℃to-0.8℃,respectively,and the salinity generally maintained from 30.2 to 32.5.Comparison with World Ocean Atlas 2018’s climatology manifested a 40m-thick and warm ACW roughly ex-ceeding the temperature maximum by 0.4–0.5℃in June–August 2021.This anomalously warm ACW was highly related to the ex-pansion of the Beaufort Gyre in the negative Arctic Oscillation phase.During summer,the under-ice oceanic heat flux F_(w)^(OHF)was elevated,with a maximum value of above 25Wm^(-2).F_(w)^(OHF)was typically low in the freezing season,with an average value of 1.2Wm^(-2).The estimates of upward heat flux contributed by ACW to the sea ice bottom F_(w)^(OHF)were in the range of 3–4Wm^(-2)in June–August 2021,when ACW contained a heat content of more than 80MJm^(-2).The heat loss over this period was driven by a weak stratification upon the ACW layer associated with a surface mixed layer(SML)approaching the ACW core.After autumn,F_(w)^(OHF)was reduced(<2 Wm^(-2))except during rare events when it elevated F_(w)^(OHF)slightly.In addition,the intensive and widespread Ekman suction,which created a violent upwelling north of the Canada Basin,was largely responsible for the substantial cooling and thinning of the ACW layer in the summer of 2021.

Nuclear charge radius predictions by kernel ridge regression with odd-even effects

The extended kernel ridge regression(EKRR)method with odd-even effects was adopted to improve the description of the nuclear charge radius using five commonly used nuclear models.These are:(i)the isospin-dependent A^(1∕3) formula,(ii)relativistic continuum Hartree-Bogoliubov(RCHB)theory,(iii)Hartree-Fock-Bogoliubov(HFB)model HFB25,(iv)the Weizsacker-Skyrme(WS)model WS*,and(v)HFB25*model.In the last two models,the charge radii were calculated using a five-parameter formula with the nuclear shell corrections and deformations obtained from the WS and HFB25 models,respectively.For each model,the resultant root-mean-square deviation for the 1014 nuclei with proton number Z≥8 can be significantly reduced to 0.009-0.013 fm after considering the modification with the EKRR method.The best among them was the RCHB model,with a root-mean-square deviation of 0.0092 fm.The extrapolation abilities of the KRR and EKRR methods for the neutron-rich region were examined,and it was found that after considering the odd-even effects,the extrapolation power was improved compared with that of the original KRR method.The strong odd-even staggering of nuclear charge radii of Ca and Cu isotopes and the abrupt kinks across the neutron N=126 and 82 shell closures were also calculated and could be reproduced quite well by calculations using the EKRR method.

Gravity Fault Subsidence and Beach Ridges Progradation in Quinta-Cassino (RS) Coastal Plain, Brazil

Ground penetrating radar (GPR) surveys have being applied to investigate very near-surface stratification of sedimentary units in coastal plains and to define their depositional conditions. This paper presents, however, low-frequency GPR survey to investigate fault-related depositional systems at greater depths. The Quinta-Cassino area in the Rio Grande do Sul Coastal Plain (RGSCP, Brazil) shows a wide strandplain that is made off by very long, continuous, and linear geomorphic features (beach ridges). This strandplain extends for ~70 km southward. The beach ridges show low-angle truncations against the Quinta escarpment, and also truncations in the strandplain. The traditional approach points that RGSCP was developed by juxtaposition of four lagoons/barrier systems as consequence of sea level changes;previous model assumes that no deformational episode occurred in RGSCP. The geophysical and geological surveys carried out in this area showed the existence of listric fault controlling the beach ridges in the escarpments and hanging-wall blocks. The radargrams could distinguish Pleistocene basement unit anticlockwise rotation, thickening of beach ridges radarfacies close to listric normal faults, and horst structures. These deformational features indicate that the extensional zone of a large-scale gravity-driven structure controlled the mechanical subsidence, the Holocene sedimentation and its stratigraphic and geomorphic features in the Quinta-Cassino area to build up an asymmetric delta. The results point to a new approach in dealing with RGSCP Holocene evolution.

基于Bayesian-Ridge模型的煤炭企业净资产收益率影响因素

为评估各因素对净资产收益的影响程度,以2002—2022年我国上市煤炭企业的相关数据为研究样本,基于机器学习中的Bayesian-Ridge模型进行训练,分析各自变量对因变量的影响程度,并与传统OLS模型进行对比。实验结果表明:在影响因素权重大小方面,Bayesian-Ridge模型与OLS模型均显示经营负债比率对净资产收益率影响最大,其次是国有持股比例和营业利润率;在方法层面,随样本点数据增加,Bayesian-Ridge模型可降低先验分布的影响;Bayesian-Ridge模型决定系数R~2为0.002 3,高于OLS模型的0.001 8,在该问题研究中Bayesian-Ridge模型优于OLS模型。

Review on the fully mulched ridge–furrow system for sustainable maize production on the semi-arid Loess Plateau

The fully mulched ridge–furrow(FMRF) system has been widely used on the semi-arid Loess Plateau of China due to its high maize(Zea mays L.) productivity and rainfall use efficiency. However, high outputs under this system led to a depletion of soil moisture and soil nutrients, which reduces its sustainability in the long run. Therefore, it is necessary to optimize the system for the sustainable development of agriculture. The development, yield-increasing mechanisms,negative impacts, optimization, and their relations in the FMRF system are reviewed in this paper. We suggest using grain and forage maize varieties instead of regular maize;mulching plastic film in autumn or leaving the mulch after maize harvesting until the next spring, and then removing the old film and mulching new film;combining reduced/notillage with straw return;utilizing crop rotation or intercropping with winter canola(Brassica campestris L.), millet(Setaria italica), or oilseed flax(Linum usitatissimum L.);reducing nitrogen fertilizer and partially replacing chemical fertilizer with organic fertilizer;using biodegradable or weather-resistant film;and implementing mechanized production. These integrations help to establish an environmentally friendly, high quality, and sustainable agricultural system, promote highquality development of dryland farming, and create new opportunities for agricultural development in the semi-arid Loess Plateau.

Genesis of ^(230)Th excess in basalts from mid-ocean ridges and ocean islands:Constraints from the global U-series isotope database and major and rare earth element geochemistry

Based on 230Th-238U disequilibrium and major element data from mid-ocean ridge basalts(MORBs) and ocean island basalts(OIBs),this study calculates mantle melting parameters,and thereby investigates the origin of 230Th excess.(230Th/238U) in global MORBs shows a positive correlation with Fe8,Po,Na8,and Fmelt(Fe8 and Na8 are FeO and Na2O contents respectively after correction for crustal fractionation relative to MgO = 8 wt%,Po=pressure of initial melting and Fmelt=degree of melt),while 230Th excess in OIBs has no obvious correlation with either initial mantle melting depth or the average degree of mantle melting.Furthermore,compared with the MORBs,higher(230Th/238U) in OIBs actually corresponds to a lower melting degree.This suggests that the 230Th excess in MORBs is controlled by mantle melting conditions,while the 230Th excess in OIBs is more likely related to the deep garnet control.The vast majority of calculated initial melting pressures of MORBs with excess 230Th are between 1.0 and 2.5 GPa,which is consistent with the conclusion from experiments in recent years that DU>DTh for Al-clinopyroxene at pressures of >1.0 GPa.The initial melting pressure of OIBs is 2.2-3.5 GPa(around the spinel-garnet transition zone),with their low excess 226Ra compared to MORBs also suggesting a deeper mantle source.Accordingly,excess 230Th in MORBs and OIBs may be formed respectively in the spinel and garnet stability field.In addition,there is no obvious correlation of K2O/TiO2 with(230Th/238U) and initial melting pressure(Po) of MORBs,so it is proposed that the melting depth producing excess 230Th does not tap the spinel-garnet transition zone.OIBs and MORBs in both(230Th/238U) vs.K2O/TiO2 and(230Th/238U) vs.Po plots fall in two distinct areas,indicating that the mineral phases which dominate their excess 230Th are different.Ce/Yb-Ce curves of fast and slow ridge MORBs are similar,while,in comparison,the Ce/Yb-Ce curve for OIBs shows more influence from garnet.The mechanisms generating excess 230Th in MORBs and OIBs are significantly different,with formation of excess 230Th in the garnet zone only being suitable for OIBs.

Origin of submarine canyon-channel systems along the middle segment of West Mariana Ridge,Philippine Sea

Submarine canyon-channel systems have been documented in the Parece Vela Basin,West Mariana Ridge;however,the mechanism of the formation and controlling factors remain poorly understood.Based on high-resolution multibeam bathymetric data and two-dimensional(2D)seismic profiles,we identified and mapped the submarine canyon-channel system along the middle segment of West Mariana Ridge in the Philippine Sea.These submarine canyon-channels show a main W-E orientation at depth of 2000–4500 m.They are approximately 72–128 km in length and 1.3–15 km in width,and their canyon heads are adjacent to the seamounts with several branches.The upper reaches of submarine canyon-channels are characterized by deeply incised,narrow,V-shaped thalwegs,suggesting the powerful erosion of gravity flows.The distinguished sediment waves are suggested to be resulted from the interaction of turbidity currents and seafloor.Our observations demonstrate that gravity flows originated from the collapses of seamount flanks plays a vital role in developing the submarine canyonchannel system along the West Mariana Ridge.This work provides better understanding of erosion,transport,and deposition of sediments from subducting ridges to deep-water basins,and also new insights into the origin and evolution of submarine canyon-channel systems along subducting ridges.

Export of Greenland Sea water across the Mohn Ridge induced by summer storms

The Mohn Ridge separates the Greenland Sea and the Lofoten Basin.Previous studies identified the export across the Mohn Ridge(EMR)from the Greenland Basin into the Lofoten Basin using water mass analysis and the tracer diffusion method,but there is still lack of direct current measurements.A surface mooring with four current meters was deployed on the Mohn Ridge from June 5 to June 18 in 2015,when three cyclones passed in the adjacent area.In the absence of cyclones,the flow on the Mohn Ridge was northeastward,parallel to the ridge.When cyclones appeared,the EMR occurred to transport Greenland Sea water into the Lofoten Basin.The probable mechanism is the sea level height variation caused by the perturbation of the low pressure of cyclones,which yields an outward pressure gradient force to drive the outflow.Our results suggest that the outflow is intermittent and only happens during cyclone activities.The annual mean volume flux of EMR was roughly estimated by the limited data,which is about 3.8×10^(6)m^(3)/s,a little bit smaller than the estimation based on volume conservation.The results indicate that the contribution of the cyclonic storms to EMR is a greatly important mechanism that potentially influences the global thermohaline circulation through the Greenland-Scotland Ridge overflow.

Critical thresholds for stage division of water erosion process in different ridge systems in mollisol region of Northeast China

Extreme rainfall events on a slope under ridge tillage systems cause concentrated stream soil loss.To analyse the critical thresholds for different stages of water erosion process of ridge systems,simulated rainfall-erosion experiments for the contour wide ridge(CWR),contour narrow ridge(CNR),longitudinal wide ridge(LWR),and longitudinal narrow ridge(LNR)were conducted under four rainfall intensities,with slope gradients of 3°and 5°.For the runoff event,the runoff depth order was LNR>LWR>CWR>CNR;the soil loss order was CNR>LNR>CWR>LWR.The product of slope factor(S)and rainfall erosivity(R)or runoff depth(D),can be adopted as critical thresholds for different stages of runoff and soil erosion process.For the longitudinal ridge systems,R values were provided for LWR and LNR and were the beginning of sheet flow,whereas the product of rainfall erosivity and slope factor(RS)values were provided for LWR and LNR as the beginning of the accelerated concentrated flow.For the contour ridge systems,R values were provided for CWR and CNR as critical thresholds for the beginning of overflow.The product of runoff depth and slope factor(DS)values were 9.98 and 7.73 mm for CWR and CNR,respectively,and were critical thresholds for the beginning of ridge failure;the DS values were 18.45 and 12.75 mm for CWR and CNR,respectively,and were critical thresholds for the beginning of the formation of ephemeral gully erosion.The critical thresholds can distinguish different stages of soil erosion process modelling.

Geological context and vents morphology in the ultramafic-hosted Tianxiu field, Carlsberg Ridge

The Tianxiu hydrothermal field(TXHF) located on Carlsberg Ridge is one of the few active ultramafic-hosted venting systems known in the Indian Ocean. Despite numerous investigations, there is limited understanding of its sulfide structure morphology, and the factors controlling the formation of TXHF are poorly understood. In this study, we conducted detailed seafloor mapping using visual data obtained by dives using the human-occupied vehicle(HOV) Jiaolong. The TXHF is found to be an active, off-axis, ultramafic-hosted, high-temperature hydrothermal area in which serpentine peridotite is exposed. Two main hydrothermal sites were identified, i.e., P and Y, both of which feature a complex of chimneys and beehive diffusers constituting a “chimney jungle” and isolated large steep-sided structures developed on flat-lying sulfide mounds. In addition, some sporadic inactive chimneys and outcrops of hydrothermal deposits were noted. The chimneys are rich in Fe and Zn sulfide, and lack the central fluid channel formed by focused high-temperature fluid flow. Hydrothermal venting at TXHF is likely related to low-angle detachment faults that focus and transport hydrothermal fluids away from a heat source along the valley wall. Our results complement and expand upon previous works concerning sulfide chimney morphology and their corresponding mineral paragenesis in ultramafic-hosted hydrothermal systems in the Indian Ocean and further our understanding of modern seafloor hydrothermal systems.

Ridge regression energy levels calculation of neutral ytterbium(Z=70)

In view of the difficulty in calculating the atomic structure parameters of high-Z elements,the Hartree–Fock with relativistic corrections(HFR)theory in combination with the ridge regression(RR)algorithm rather than the Cowan code’s least squares fitting(LSF)method is proposed and applied.By analyzing the energy level structure parameters of the HFR theory and using the fitting experimental energy level extrapolation method,some excited state energy levels of the Yb I(Z=70)atom including the 4f open shell are calculated.The advantages of the ridge regression algorithm are demonstrated by comparing it with Cowan code’s LSF results.In addition,the results obtained by the new method are compared with the experimental results and other theoretical results to demonstrate the reliability and accuracy of our approach.

Analysis of Ship Motion in An Ice Ridge Field Based on Empirical Approach

The estimation of ship speed in ice ridge fields is important for both route planning and prediction of emergency response time.An analytical method for estimating ship motion in first-year ice ridge is developed based on ice resistance models and ship motion equations,in which the effect of ship speed on ridge resistance is taken into account.Two model tests in level ice and one model test in ice ridge for an icebreaking multipurpose vessel are used to validate and benchmark the presented method.The predicted results including level ice resistances,net thrust and ship motion in the ice ridge field are compared with the model test data.The comparisons show that the presented method can generate reasonable results.The effects of input parameters on ship speed,penetration depth and number of necessary rams to transit ridge have been studied.Based on the calibrated model,insights into the ice resistance and the ship motion are obtained.It is found that the energy consumption of the keel obtained by integral calculation of the keel resistance at the penetration distance is with the same magnitude as the result of the maximum keel resistance multiplied by the ridge length.In addition,the effect of ridge width and keel depth on keel resistance and average transit speed is investigated.

Formation and modification of wrinkle ridges in the central Tharsis region of Mars as constrained by detailed geomorphological mapping and landsystem analysis

Wrinkle ridges are common landforms documented on all rocky planets and the Moon in the inner solar system.Despite the long research history,their formation mechanisms remain debated.A key unresolved issue is whether the wrinkle-ridge formation is related to igneous processes.This is because wrinkle ridges are mostly associated in space and possibly in time with the occurrence of flood-basalt volcanism in all cases in the inner solar system.To address this issue,we conducted geomorphological mapping,a topographic-data analysis,and a detailed landform and landsystem analysis of satellite images at a resolution of 25 cm/pixel to 6 m/pixel in the central Tharsis region of Mars.The main results of this work are in the form of(1)a regional geomorphological map at a resolution of 6 m/pixel and(2)a local geomorphological map at a resolution of 50 cm/pixel.Our work suggests the following older-to-younger sequence of geological events in the study area:(1)formation of a northeast-trending mountain range(i.e.,the Thaumasia plateau)along the eastern margin of the Tharsis rise that was created by the Himalayan-style crustal-scale thrusting;(2)coeval volcanic-plateau construction west of the thrusting-induced rising mountain range;(3)eastward-flowing lavas that were sourced from a volcanic plateau to the west terminated at the rising Thaumasia plateau to the east;(4)wrinkle-ridge development by decollement folding of recently emplaced warm,ductile volcanic-lava piles;(5)emplacement of a regionally extensive ice sheet over the central Tharsis region that produced extensive boulder-bearing materials,striated surfaces,and boulder-bearing dendritic-ridge networks possibly representing subglacial eskers;and(6)local deposition of highly concentrated glacial flours resulted in the formation of mantled terrain on plains between wrinkle ridges.Our work supports the early suggestion that the Tharsis wrinkle ridges were created by horizontal shortening induced by crustal-scale tectonic processes.In detail,however,the occurrence of flow-front-like fold margins associated with many mapped wrinkle ridges suggests the involvement of ductile-flow deformation during ridge formation.We attribute the flow-like fold fronts to ductile deformation of thermally weakened lava piles that were emplaced during or immediately before the folding event.Our compression-induced wrinkle-ridge model also differs from the early hypotheses in that the thin-skinned folding is associated with basement subduction,which explains the lack of coeval and parallel folding and extensional faulting associated with wrinkle ridge formation in the study area.Post-folding glacial modification means that the present wrinkle-ridge morphologies may differ significantly from the original fold shapes,which prevents the utility of using topographic profiles across wrinkle ridges for inverting the underlying thrust geometries.