Distinguishing the main components of submarine groundwater and estimating the corresponding fluxes based on radium tracing method—taking the Maowei Sea for example

Submarine groundwater discharge(SGD)is an important part in the land-sea interactions,which mainly contains three components:submarine fresh groundwater discharge(SFGD),tidal flat recirculated saline groundwater discharge(tidal flat RSGD)and subtidal recirculated saline groundwater discharge(subtidal RSGD).In order to make a more accurate assessment of the impact of SGD on coastal ecological environment,it is necessary to distinguish the main components of SGD.In this study,the Maowei Sea,located in the northern part of the Beibu Gulf,was selected as the study area.Based on the radium(Ra)tracing method,we present a new analytical method for distinguishing the three main components of SGD in this area combined with field data.The average daily flow along the coastline of the Maowei Sea of tidal flat RSGD was slightly higher than that of SFGD,and both two were on the magnitude of 1×10^(5)m^(3)/d.The average daily flow for the subtidal RSGD of the entire subtidal zone of the Maowei Sea reached to the magnitude of 1×10^(6)-1×10^(7)m^(3)/d.The long-term variation trend of terrestrial SGD is a valuable information for the study of the influence of terrigenous material on the coastal ecological environment.Based on the results of four sampling periods,it is found that the fluxes of SFGD and tidal flat RSGD in the Maowei Sea had good linear correlation with the net precipitation.As an example,January 2015 to August 2022 were selected as the study periods,and the variation trends of SFGD and tidal flat RSGD were calculated by linear function with net precipitation as the independent variable.The results showed that the flux of tidal flat RSGD was slightly higher than that of SFGD,and the difference between the two is larger in flood season while smaller in dry season.In general,in the coastal range of China,the total SGD flux in the Maowei Sea area is at a high level,and the SFGD flux is at a medium level.

Deviation of the Lagrangian particle tracing method in the evaluation of the Southern Hemisphere annual subduction rate

The classical Lagrangian particle tracing method is widely used in the evaluation of the ocean annual subduction rate.However,our analysis indicates that in addition to neglecting the effect of mixing,there are two possible deviations in the method:one is an overestimation due to not considering that the amount of subducted water at the source location may be inadequate during the late winter of the first year when the mixed layer becomes shallow;the other one is an underestimation due to the neglect of the effective subduction caused by strong vertical pumping.Quantitative analysis shows that these two deviations mainly exist in the low-latitude subduction areas of the South Pacific and South Atlantic.The two deviations have very similar distribution areas and can partially off set each other.However,the overall deviation is still large,and the maximum relative deviation ratio can reach 50%;therefore,it cannot be ignored.

An improved algorithm based on equivalent sound velocity profile method at large incident angle

With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical equivalent sound speed profile(ESSP)method replaces the measured sound velocity profile(SVP)with a simple constant gradient SVP,reducing the computational workload of beam positioning.However,in deep-sea environment,the depth measurement error of this method rapidly increases from the central beam to the edge beam.By analyzing the positioning error of the ESSP method at edge beam,it is discovered that the positioning error increases monotonically with the incident angle,and the relationship between them could be expressed by polynomial function.Therefore,an error correction algorithm based on polynomial fitting is obtained.The simulation experiment conducted on an inclined seafloor shows that the proposed algorithm exhibits comparable efficiency to the original ESSP method,while significantly improving bathymetry accuracy by nearly eight times in the edge beam.

Low-complexity signal detection for massive MIMO systems via trace iterative method

Linear minimum mean square error(MMSE)detection has been shown to achieve near-optimal performance for massive multiple-input multiple-output(MIMO)systems but inevitably involves complicated matrix inversion,which entails high complexity.To avoid the exact matrix inversion,a considerable number of implicit and explicit approximate matrix inversion based detection methods is proposed.By combining the advantages of both the explicit and the implicit matrix inversion,this paper introduces a new low-complexity signal detection algorithm.Firstly,the relationship between implicit and explicit techniques is analyzed.Then,an enhanced Newton iteration method is introduced to realize an approximate MMSE detection for massive MIMO uplink systems.The proposed improved Newton iteration significantly reduces the complexity of conventional Newton iteration.However,its complexity is still high for higher iterations.Thus,it is applied only for first two iterations.For subsequent iterations,we propose a novel trace iterative method(TIM)based low-complexity algorithm,which has significantly lower complexity than higher Newton iterations.Convergence guarantees of the proposed detector are also provided.Numerical simulations verify that the proposed detector exhibits significant performance enhancement over recently reported iterative detectors and achieves close-to-MMSE performance while retaining the low-complexity advantage for systems with hundreds of antennas.

Measurement and verification of concentration-dependent diffusion coefficient:Ray tracing imagery of diffusion process

Ray tracing method is used to study the propagation of collimated beams in a liquid-core cylindrical lens(LCL),which has dual functions of diffusion cell and image formation.The diffusion images on the focal plane of the used LCL are simulated by establishing and solving both linear and nonlinear ray equations,the calculated results indicate that the complex imaging results of LCL in inhomogeneous media can be treated by the law of ray propagation in homogeneous media under the condition of small refractive index gradient of diffusion solution.Guided by the calculation conditions,the diffusion process of triethylene glycol aqueous solution is experimentally studied at room temperature by using the LCL in this paper.The spatial and temporal concentration profile Ce(z,t)of diffusion solution is obtained by analyzing diffusion image appearing on the focal plane of the LCL;Then,the concentration-dependent diffusion coefficient is assumed to be a polynomial D(C)=D0×(1+α1C+α2C2+α3C3+…).The finite difference method is used to solve the Fick diffusion equation for calculating numerically the concentration profiles Cn(z,t).The D(C)of triethylene glycol aqueous solution is obtained by comparing the Cn(z,t)with Ce(z,t).Finally,the obtained polynomial D(C)is used to calculate the refractive index profiles nn(z,t)s of diffusion solution in the used LCL.Based on the ray propagation law in inhomogeneous media and the calculated n(z,t),the ray tracing method is used again to simulate the dynamic images of the whole experimental diffusion process to varify the correctness of the calculated D(C).The method presented in this work opens up a new way for both measuring and verifying the concentration-dependent liquid diffusion coefficients.

Partial Trace Method for Deriving Density Operators of Light Field

With the help of technique of integration within an ordered product of operators we can fully display the partial trace method to deduce some new density operators of light field. The general form of single-mode density operator is derived in this way from a normalized two-mode generalized squeezed pure state density operator.

A FLUORESCENCE QUENCHING IMNUNOASSAY METHOD FOR DETERMINATION OF TRACE ALBUMIN USING UNLABELED EUROPIUM CHELATE

A fluoroimmunoassay method using unlabeled europium chelate is described.The principle is similar to that of fluoroimmunoassay method using lanthanide chelate as labels.The procedure is simple because labeling process is omitted.The detection limit is about 10^(10) mol/L antigen.The relative standard deviation of immunoassay is less than 10%.The recoveries of human serum albumin and estradiol protein conjugate are 96-105% and 111% respectively.

A QUENCHING FLUORESCENCE IMMUNOASSAY METHOD FOR DETERMINATION OF TRACE ALBUMIN USING UNLABELED TERBIUM CHELATE

A fluoroimmunoassay method using unlabeled Terbium chelate is described.The principle is similar to that of fluoroimmunoassay method using lanthanide chelate as labels.The procedure is simpte because labeling process is unnecessary.The recovery of HSA and albumin in urine is 107% and 95% respectively.The standard deviation is tess than 10%.

A simple and effective model for Jones matrix to evaluate the transmission performance of Savart polariscope

Starting from the wave normal tracing treatment in the Savart polariscope that relates the Jones matrix to its transmission performance, this paper establishes a simple and effective model for the Jones matrix at an arbitrary incidence in the spatial domain. Analytical expressions of all the matrix components are determined with the consideration of all the main impact factors. This model needs only a few parameters hence it is convenient to be employed to evaluate the propagation performance of any birefringent optical system. The simulated results obtained with it demonstrate that this model gives a precise representation of the characteristic of light propagation in the Savart polariscope. This would provide useful suggestions for the design, calibration, and performance improvement of any other birefringent polarisation element and optical system.

Relation between viscoplastic flow, periodic feature and Joule heat validity for electrically assisted friction stir welding

Electrically assisted friction stir welding is a promising way to improve liquidity of viscoplastic metal. However, this local heat source only can take into effect in a certain contact interface status. Semi-analytical model and FEM are adopted to analyze materials flow in electrically assisted FSW. Semi-analysis results show that Joule heating is validity with rotation speed 1 400 r/min and travel speed 400 mm/min. The maximum flow velocity is 0.033 m/s when travel speed is 400 mm/min, and it increases to 0.043 m/s with current assistant. Visco-plasticizing efficiency is a key factor to decide whether Joule heat is significant, and viscoplastic flow can periodically make circuit short. Periodic torque, axial force and onion rings surface show dynamic flow process, which can be used to illustrate the dynamic contact resistance feature and correct the calculation equation of Joule heat. This paper is aim to point out the potential relation between viscoplastic flow, periodic torque and Joule heating validity.

Numerical simulation of friction stir butt-welding of 6061 aluminum alloy

A two-dimensional computational fluid dynamics model was established to simulate the friction stir butt-welding of 6061 aluminum alloy. The dynamic mesh method was applied in this model to make the tool move forward and rotate in a manner similar to a real tool, and the calculated volumetric source of energy was loaded to establish a similar thermal environment to that used in the experiment. Besides, a small piece of zinc stock was embedded into the workpiece as a trace element. Temperature fields and vector plots were determined using a finite volume method, which was indirectly verified by traditional metallography. The simulation result indicated that the temperature distribution was asymmetric but had a similar tendency on the two sides of the welding line. The maximum temperature on the advancing side was approximately 10 K higher than that on the retreating side. Furthermore, the precise process of material flow behavior in combination with streamtraces was demonstrated by contour maps of the phases. Under the shearing force and forward extrusion pressure, material located in front of the tool tended to move along the tangent direction of the rotating tool. Notably, three whirlpools formed under a special pressure environment around the tool, resulting in a uniform composition distribution.

Multi-point Contact of the High-speed Vehicle-turnout System Dynamics

The wheel-rail contact problems, such as the number, location and the track of contact patches, are very important for optimizing the spatial structure of the rails and lowering the vehicle-turnout system dynamics. However, the above problems are not well solved currently because of having the difficulties in how to determine the multi-contact, to preciously present the changeable profiles of the rails and to establish an accurate spatial turnout system dynamics model. Based on a high-speed vehicle-turnout coupled model in which the track is modeled as flexible with rails and sleepers represented by beams, the line tracing extreme point method is introduced to investigate the wheel-rail multiple contact conditions and the key sections of the blade rail, longer nose rail, shorter rail in the switch and nose rail area are discretized to represent the varying profiles of rails in the turnout. The dynamic interaction between the vehicle and turnout is simulated for cases of the vehicle divergently passing the turnout and the multi-point contact is obtained. The tracks of the contact patches on the top of the rails are presented and the wheel-rail impact forces are offered in comparison with the contact patches transference on the rails. The numerical simulation results indicate that the length of two-point contact occurrence of a worn wheel profile and rails is longer than that of the new wheel profile and rails; The two-point contact definitely occurs in the switch and crossing area. Generally, three-point contact doesn’t occur for the new rail profile, which is testified by the wheel-rails interpolation distance and the first order derivative function of the tracing line extreme points. The presented research is not only helpful to optimize the structure of the turnout, but also useful to lower the dynamics of the high speed vehicle-turnout system.

INVESTIGATIONS ON MOISTURE TRANSPORTS, BUDGETS AND SOURCES RESPONSIBLE FOR THE DECADAL VARIABILITY OF PRECIPITATION IN SOUTHERN CHINA

In the context of global warming,apparent decdal-interdecdal variabilities can be detected in summer precipitation in southern China.Especially around the 1990 s,precipitation in South China experienced a phase transition from a deficiency regime to an abundance regime in the early 1990 s,while the Yangtze River Valley witnessed a phase shift of summer precipitation from abundance to deficiency in the late 1990 s.Pertinent analyses reveal a close relationship between such decadal precipitation shifts and moisture budgets,which is mainly modulated by the meridional component.This relationship can be attributed to large-scale moisture transport anomalies.Further,the HYSPLIT model is utilized to quantitatively evaluate relative moisture contributions from diverse sources during different regimes.It can be found that during the period with abundant precipitation in South China,the moisture contribution from the source of Indochina Peninsula-South China Sea increased significantly,while during the deficient precipitation regime in the Yangtze River Valley,moisture from local source,western Pacific and Indochina Peninsula-South China Sea contributed less to precipitation.It means some new features of relative moisture contributions from diverse sources to precipitation anomaly in southern China took shape after 1990 s.

An efficient adaptive space partitioning algorithm for electromagnetic scattering calculation of complex 3D models

The space partitioning algorithm based on the rounding and addressing operations has been proved to be an efficient space partitioning algorithm with the potential for real-time calculation.An improvement on this kind of space partitioning algorithms for solving complex 3D models is presented.Numerical examples show that the efficiency of the improved algorithm is better than that of the original method.When the size of most target elements is smaller than the size of spatial grids,the efficiency of the improved method can be more than four times of that of the original method.An adaptive method of space partitioning based on the improved algorithm is developed by taking the surface element density or the curvature as the threshold for deep partitioning and conducting the deep partitioning using the octree method.A computer program implementation for applying the method in some typical applications is discussed,and the performance in terms of the efficiency,reliability,and resource use is evaluated.Application testing shows that the results of the adaptive spacing partitioning are more convenient for the follow-up use than that of the basic uniform space partitioning.Furthermore,when it is used to calculate the electromagnetic scattering of complex targets by the ray tracing(RT)method,the adaptive space partitioning algorithm can reduce the calculation time of the RT process by more than 40%compared with the uniform space segmentation algorithm.

Thermal vacuum state for the two-coupled-oscillator model at finite temperature: Derivation and application

Following the spirit of thermo field dynamics initiated by Takahashi and Umezawa, we employ the technique of integration within an ordered product of operators to derive the thermal vacuum state （TVS） for the Hamiltonian H of the two-coupled-oscillator model. The ensemble averages of the system are derived conveniently by using the TVS. In addition, the entropy for this system is discussed based on the relation between the generalized Hellmann-Feynman theorem and the entroy variation in the context of the TVS.

A Tracing Algorithm for Surface—Surface Intersections on Surface Boundaries

In this paper we present an algorithm with a new trace-terminating conditionfor tracing along surface-surface intersection curves on surface boundaries, while several tracingmethods and embedding methods that include tracing scheme may cause false termination with atraditional trace-terminating condition: tracing stops when the surface-domain's boundary isreached. And we also suggest another iterative method to calculate intersection points on surfaceboundaries with parallel surface normal. Some numerical examples with these two ideas andcomparisons to 'DESIGNBASE', 'ACIS', and 'Parasolid' are included to demonstrate the effectivenessof our algorithm.

Uniform Deterministic Discrete Method for Three Dimensional Systems

For radiative direct exchange areas in three dimensional system, the Uniform Deterministic Discrete Method (UDDM) was adopted. The spherical surface dividing method for sending area element and the regular icosahedron for sending volume element can meet with the direct exchange area computation of any kind of zone pairs. The numerical examples of direct exchange area in three dimensional system with nonhomogeneous attenuation coefficients indicated that the UDDM can give very high numerical accuracy.

Simple and cost-effective methods for precise analysis of trace element abundances in geological materials with ICP-MS

Inductively coupled plasma mass spectrometry （ICP-MS） is the most commonly used technique to deter- mine the abundances of trace elements in a wide range of geological materials. However, incomplete sample digestion, isobaric interferences and instrumental drift remain obvious problems that must be overcome in order to obtain precise and accurate results, For this reason, we have done many experi- ments and developed a set of simple, cost-effective and practical methods widely applicable for precise and rapid determination of trace element abundances in geological materials using ICP-MS. Commonly used high-pressure digestion technique is indeed effective in decomposing refractory phases, but this inevitably produces fluoride complexes that create new problems. We demonstrate that the fluoride complexes formed during high-pressure digestion can be readily re-dissolved using high-pressure vessel at 190 ℃ for only 2 h for 50 mg sample. In the case of isobaric interferences, although oxide （e.g., MO^＋/M^＋） and hydroxide （e.g., MO^＋/M^＋） productivity is variable between runs, the （MO^＋/M^＋）/（CeO^＋/Ce^＋） and （MOH^＋/M^＋）/（CeO^＋/Ce^＋） ratios remain constant, making isobaric interference correction for all other elements of interest straightforward, for which we provide an easy-to-use off-line procedure. We also show that mass-time-intensity drift curve is smooth as recognized previously, for which the correction can be readily done by analyzing a quality-control （QC） solution and using off-line Excel VBA procedure without internal standards. With these methods, we can produce data in reasonable agreement with rec- ommended values of international rock reference standards with a relative error of 〈8% and precision generally better than 5%. Importantly, compared to the widely used analytical practice, we can effectively save 〉60% of time （e.g., 〈24 h vs. 〉60 h）.

TRACE AVERAGING DOMAIN DECOMPOSITION METHOD WITH NONCONFORMING FINITE ELEMENTS

We consider, in this paper, the trace averaging domain decomposition method for the second order self-adjoint elliptic problems discretized by a class of nonconforming finite elements, which is only continuous at the nodes of the quasi-uniform mesh. We show its geometric convergence and present the dependence of the convergence factor on the relaxation factor, the subdomain diameter H and the mesh parameter h. In essence;, this method is equivalent to the simple iterative method for the preconditioned capacitance equation. The preconditioner implied in this iteration is easily invertible and can be applied to preconditioning the capacitance matrix with the condition number no more than O((1 + In H/h)max(1 + H-2, 1 + In H/h)).

Thermal Emission of a Disc Body of SemitransparentMaterial

By introducing the concept of radiosity intensity to diffuse surfaces, the ray tracing method is improved to analyze the thermal emission of a disc body of gray semitransparent material. The two plane sur-faces of the disc body are both specularly reflecting, and the fiank surface is either diffusely reflecting or specularly reflecting. The apparent thermal emission from one plane sllrface is investigated with considering the infiuences of the characteristic optical thickness, the dimensionless radius, the refrac-tive index of the material and the reflecting characteristics of the flank surface. The directional and hemispherical emissions show considerable differences under different refiecting characteristics of the flank surface. Moreover, in some cases, the emission not only varies with the viewing direction but also with the apparent emitting position on the plane surface. Some interesting results are presented and discussed.