Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Lidar-aided analysis of boreal forest backscatter at Ku band

This paper analyzes the backscatter of the microwave signal in a boreal forest environment based on a Ku -band airborne Frequency-Modulated Continuous Waveform (FMCW) profiling radar—Tomoradar. We selected a half-managed boreal forest in the southern part of Finland for a field test. By decomposing the waveform collected by the Tomoradar, the vertical canopy structure was achieved. Based on the amplitude of the waveform, the Backscattered Energy Ratio of Canopy-to-Total (BERCT) was calculated. Meanwhile, the canopy fraction was derived from the corresponding point cloud recorded by a Velodyne VLP-16 LiDAR mounted on the same platform. Lidar-derived canopy fraction was obtained by counting the number of the first/ the strongest returns versus the total amount of returns. Qualitative and quantitative analysis of radar-derived BERCT on lidar-derived canopy fraction and canopy height are investigated. A fitted model is derived to describe the Ku-band microwave backscatter in the boreal forest to numerically analyze the proportion contributed by four factors: lidar-derived canopy fraction, radar-derived canopy height, the radar-derived distance between trees and radar sensor and other factors, from co-polarization Tomoradar measurements. The Root Mean Squared Error (RMSE) of the proposed model was 0.0958, and the coefficient of determination R2 was 0.912. The fitted model reveals that the correlation coefficient between radar-derived BERCT and lidar-derived canopy fraction is 0.84, which illustrates that lidar surface reflection explains the majority of the profiling /waveform radar response. Thus, vertical canopy structure derived from lidar can be used for the benefit of radar analysis.     

The Upper Devonian Kellwasser event recorded in a regressive sequence from inner shelf to lagoonal pond,Catalan Coastal Ranges,Spain

This study, conducted in the Catalan Coastal Ranges, north‐east Spain, describes the Upper Devonian Kellwasser event in a shallowing‐upward sequence of black shales, siltstones and quartz arenites. This sequence was deposited in a progradational and regressive coastal system where the sedimentary environment evolved from the inner shelf to a lagoonal pond located landward of the shoreline. Three anomalous succeeding steps have been identified by geochemical analysis. The first one, detected on the inner shelf, was characterized by oxygen depletion and high organic productivity. The second, detected in the nearshore, was caused by hydrothermal activity occurring under normal oxic conditions. The third and most intense step was identified in the muds of the lagoonal pond and has been linked to strong anoxic conditions, elevated clastic input derived from changes in the weathering regime at the source area and moderate hydrothermal activity. The Kellwasser event is thus defined in the study area as stepwise and multi‐causal. This is the first time that the Kellwasser event has been identified in a sedimentary environment behind the shoreline. It is also the first time that it has been reported in the Catalan Coastal Ranges.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

WR 7a: a V Sagittae or a qWR star?

The star WR 7a, also known as SPH 2, has a spectrum that resembles that of V Sagittae stars although no O  vi emission has been reported. The Temporal Variance Spectrum – TVS – analysis of our data shows weak but strongly variable emission of O  vi lines which is below the noise level in the intensity spectrum.
Contrary to what is seen in V Sagittae stars, optical photometric monitoring shows very little, if any, flickering. We found evidence of periodic variability. The most likely photometric period is   P _phot= 0.227(±14) d  , while radial velocities suggest a period of   P _spec= 0.204(±13) d  . One-day aliases of these periods can not be ruled out. We call attention to similarities with HD 45166 and DI Cru (= WR 46), where multiple periods are present. They may be associated to the binary motion or to non-radial oscillations.
In contrast to a previous conclusion by Pereira et al., we show that WR 7a contains hydrogen. The spectrum of the primary star seems to be detectable as the N  v 4604 Å  absorption line is visible. If so, it means that the wind is optically thin in the continuum and that it is likely to be a helium main sequence star.
Given the similarity to HD 45166, we suggests that WR 7a may be a qWR – quasi Wolf–Rayet – star. Its classification is WN4h/CE in the Smith, Shara & Moffat three-dimensional classification system.     

Coma Morphology of Three Non-periodic Comets

The coma morphology and short-term evolution was investigated of three non-periodic comets in retrograde orbits, C/2001 Q4 (NEAT), C/2002 T7 (LINEAR), and C/2003 K4 (LINEAR). All three comets display distinct coma features, which were very different from one comet to the next and remained rather constant in shape during the observational period. A single, broad feature perpendicular to the sun-tail direction dominated the coma of C/2003 K4 in all used filters (B,V,R,I), whereas the coma of Comet C/2002 T7 exhibited different features in blue and red filters. C/2001 Q4 showed rather complex coma morphology with clear short-term variability in coma brightness. Therefore, these non-periodic comets neither show a featureless coma nor any similarities of the features detected. The overall distribution of coma material was investigated from the shape of radial coma profiles averaged around the comet nucleus. For C/2001 Q4 and C/2002 T7, the slopes fitted to the linear part of these profiles are flatter in the blue than in the red, which can be explained by the presence of coma gas. For C/2003 K4 no such difference is indicated in the May observations (r = 2.3 AU), while in July (r = 1.7 AU) the profiles in the B-filter are flatter than in V, R, and I, hence gas contamination was relevant at least in the B filter. The R and I filter images were used to determine approximate Afρ values of each comet as a function of time.     

Evidence of climatic warming in the southern Urals region derived from borehole temperatures and meteorological data

Thirty borehole temperature–depth profiles in the central and southern Urals, Russia were scrutinized for evidence of ground surface temperature histories. We explored two inversion schemes: a simple ramp inversion in which solutions are parameterized in terms of an onset time and magnitude of change and a more sophisticated functional space inverse algorithm in which the functional form of the solution is left unspecified. To enhance and potentially identify latitudinal differences in the ground surface temperature signal, we subdivided the data into three groups based on geographic proximity and simultaneously inverted the borehole temperature–depth logs. The simultaneous inversions highlighted 13 temperature–depth logs that could not both fit a common ground surface temperature history and a priori models within reasonable bounds. Our results confirm that this is an effective way to reduce site-specific noise from an ensemble of boreholes. Each inversion scheme gives comparable results indicating locally variable warming on the order of 1°C starting between 1800 and 1900 AD. Similarly surface air temperature records from 12 nearby meteorological stations exhibit locally variable warming also on the order of 1°C of warming during the 20th century. To explore the degree to which borehole temperatures and surface air temperature (SAT) time series are responding to the same signal, we average the SAT data into the same three groups and used these averages as a forcing function at the Earth's surface to generate synthetic transient temperature profiles. Root mean square (RMS) misfits between these synthetic temperature profiles and averaged temperature–depth profiles are low, suggesting that first-order curvature in borehole temperatures and variations in SAT records are correlated.     

Simulation of double cold cores of the 35°N section in the Yellow Sea with a wave-tide-circulation coupled model

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter. This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).