Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Twisted Flux Tube Emergence Evidenced in Longitudinal Magnetograms: Magnetic Tongues

Bipolar active regions (ARs) are thought to be formed by twisted flux tubes, as the presence of such twist is theoretically required for a cohesive rise through the whole convective zone. We use longitudinal magnetograms to demonstrate that a clear signature of a global magnetic twist is present, particularly, during the emergence phase when the AR is forming in a much weaker pre-existing magnetic field environment. The twist is characterised by the presence of elongated polarities, called “magnetic tongues”, which originate from the azimuthal magnetic field component. The tongues first extend in size before retracting when the maximum magnetic flux is reached. This implies an apparent rotation of the magnetic bipole. Using a simple half-torus model of an emerging twisted flux tube having a uniform twist profile, we derive how the direction of the polarity inversion line and the elongation of the tongues depend on the global twist in the flux rope. Using a sample of 40 ARs, we verify that the helicity sign, determined from the magnetic polarity distribution pattern, is consistent with the sign derived from the photospheric helicity flux computed from magnetogram time series, as well as from other proxies such as sheared coronal loops, sigmoids, flare ribbons and/or the associated magnetic cloud observed in situ at 1 AU. The evolution of the tongues observed in emerging ARs is also closely similar to the evolution found in recent MHD numerical simulations. We also found that the elongation of the tongue formed by the leading magnetic polarity is significantly larger than that of the following polarity. This newly discovered asymmetry is consistent with an asymmetric Ω-loop emergence, trailing the solar rotation, which was proposed earlier to explain other asymmetries in bipolar ARs.     

New methods for calculating bare soil land surface temperature over mountainous terrain

Land surface temperature(LST) causes the phase change of water, links to the partitioning of surface water and energy budget, and becomes an important parameter to hydrology, meteorology, ecohydrology, and other researches in the high mountain cold regions. Unlike air temperature, which has common altitudinal lapse rates in the mountainous regions, the influence of terrain leads to complicated estimation for soil LST. This study presents two methods that use air temperature and solar position,to estimate bare LST with high temporal resolution over horizontal sites and mountainous terrain with a random slope azimuth. The data from three horizontal meteorological stations and fourteen LST observation fields with different aspects and slopes were used to test the proposed LST methods. The calculated and measured LST were compared in a range of statistical analysis, and the analysis showed that the average RMSE(root mean square error),MAD(mean absolute deviation), and R~2(correlation coefficient) for three horizontal sites were 5.09℃,3.66℃, 0.92, and 5.03℃, 3.52℃, 0.85 for the fourteen complex terrain sites. The proposed methods showed acceptable accuracy, provide a simple way to estimate LST, and will be helpful for simulating the water and energy cycles in alpine mountainous terrain.     

Groundwater flow system in Bengal Delta,Bangladesh revealed by environmental isotopes

A total number of 328 groundwater samples are analysed to evaluate the groundwater flow systems in Bengal Delta aquifers, Bangladesh using environmental isotope (²H, ¹⁸O, ¹³C, ³H, and ¹⁴C) techniques. A well-defined Local Meteoric Water Line (LMWL) δ²H = 7.7 δ¹⁸O + 10.7 ‰ is constructed applying linear correlation analyses to the monthly weighted rainfall isotopic compositions (δ¹⁸O and δ²H). The δ¹⁸O and δ²H concentrations of all groundwater samples in the study area are plotted more or less over the LMWL, which provides compelling evidence that all groundwaters are derived from rainfall and floodwater with a minor localized evaporation effects for the shallow groundwaters. Tritium concentration is observed in 40 samples out of 41 with values varying between 0.3 and 5.0 TU, which represents an evidence of young water recharge to the shallow and intermediate aquifers. A decreasing trend of ¹⁴C activity is associated with the heavier δ¹³C values, which indicates the presence of geochemical reactions affecting the ¹⁴C concentration along the groundwater flow system. Both vertical and lateral decrease of ¹⁴C activity toward down gradient show the presence of regional groundwater flow commencing from the unconfined aquifers, which discharges along the coastal regions. Finally, shallow, intermediate, and deep groundwater flow dynamics has revealed in the Bengal Delta aquifers, Bangladesh.     

Coastal sea level variability in the US West Coast Ocean Forecast System (WCOFS)

Sea level variability along the US West Coast is analyzed using multi-year time series records from tide gauges and a high-resolution regional ocean model, the base of the West Coast Ocean Forecast System (WCOFS). One of the metrics utilized is the frequency of occurrences when model prediction is within 0.15 m from the observed sea level, F. A target level of F?=?90% is set by an operational agency. A combination of the tidal sea level from a shallow water inverse model, inverted barometer (IB) term computed using surface air pressure from a mesoscale atmospheric model, and low-pass filtered sea level from WCOFS representing the effect of coastal ocean dynamics (DYN) provides the most straightforward approach to reaching levels F>80%. The IB and DYN components each add between 5 and 15% to F. Given the importance of the DYN term bringing F closer to the operational requirement and its role as an indicator of the coastal ocean processes on scales from days to interannual, additional verification of the WCOFS subtidal sea level is provided in terms of the model-data correlation, standard deviation of the band-pass filtered (2–60 days) time series, the annual cycle amplitude, and alongshore sea level coherence in the range of 5–120-day periods. Model-data correlation in sea level increases from south to north along the US coast. The rms amplitude of model sea level variability in the 2–60-day band and its annual amplitude are weaker than observed north of 42 N, in the Pacific Northwest (PNW) coast region. The alongshore coherence amplitude and phase patterns are similar in the model and observations. Availability of the multi-year model solution allows computation and analysis of spatial maps of the coherence amplitude. For a reference location in the Southern California Bight, relatively short-period sea level motions (near 10 days) are incoherent with those north of the Santa Barbara Channel (in part, due to coastal trapped wave scattering and/or dissipation). At a range of periods around 60 days, the coastal sea level in Southern California is coherent with the sea surface height (SSH) variability over the shelf break in Oregon, Washington, and British Columbia, more than with the coastal SSH at the same latitudes.     

Deriving fragility functions from bilinearized capacity curves for earthquake scenario modelling using the conditional spectrum

The development of fragility curves to perform seismic scenario-based risk assessment requires a fully probabilistic procedure in order to account for uncertainties at each step of the computation. This is especially true when developing fragility curves conditional on an Intensity Measure that is directly available from a ground-motion prediction equation. In this study, we propose a new derivation method that uses realistic spectra instead of design spectral shapes or uniform hazard spectra and allows one to easily account for the features of the site-specific hazard that influences the fragility, without using non-linear dynamic analysis. The proposed method has been applied to typical school building types in the city of Basel (Switzerland) and the results have been compared to the standard practice in Europe. The results confirm that fragility curves are scenario dependent and are particularly sensitive to the magnitude of the earthquake scenario. The same background theory used for the derivation of the fragility curves has allowed an innovative method to be proposed for the conversion of fragility curves to a common IM (i.e. spectral acceleration or PGA). This conversion is the only way direct comparisons of fragility curves can be made and is useful when inter-period correlation cannot be used in scenario loss assessment. Moreover, such conversion is necessary to compare and verify newly developed curves against those from previous studies. Conversion to macroseismic intensity is also relevant for the comparison between mechanical-based and empirical fragility curves, in order to detect possible biases.     

Remote Sensing of Eupatorium Adenophorum Spreng Based on HJ-A Satellite Data

In the southwest of China, one of the greatest threats to local ecosystem is the area expansion of an invasive species, i.e., Eupatorium adenophorum Spreng (EAS). In this study, the remote-sensing technology was used to detect and map the spatial distribution of EAS in Guizhou Province, China. A series of vegetation indices, including normalized difference vegetation index (NDVI), simple ratio index (SRI) and atmospherically resistant vegetation index (ARVI), were used to identify EAS from HJ-A Chninese satellite data. According to the analysis results of fieldworks from March 21 to 22, 2009, it was found that the vegetation index of {1.9589 ≤ SRI ≤ 4.1095}∩{0.2359 ≤ ARVI ≤ 0.5193} was the optimal remote-sensing parameter for identifying EAS from HJ-A data. According to the spatial distribution of EAS estimated from HJ-A data, it was found that EAS was rather more in southwest of Guizhou Province than in northeast. EAS became sparse from southwest to northeast gradually, and the central Guizhou Province was the ecological corridor linking EAS in southwest to that in northeast. By comparison with validated data collected by the government of Guizhou Province, it was found that the uncertainty of remote-sensing method was 18.52%, 29.31%, 8.77% and 9.46% in grassland, forest, farmland and others respectively, and the mean uncertainty was 13.29%. Owing to the lower height of EAS than many plants in forest, the uncertainty of EAS was the greatest in forest than that in grassland, farmland and so on.     

Analysis of Circular Elastic Plate Resting on Pasternak Foundation by Strain Energy Approach

In this paper, an analytical expression for the deflection of a thin circular elastic plate resting on the Pasternak foundation is derived by adopting the strain energy approach. The parametric study is carried out to observe the location of the plate lift-off and the variation of the deflection profiles for comparing the variation pattern reported in the literature based on an approximate solution technique. It is found that the radial distance of the point of lift-off of the plate decreases with increase in the values of both the shear modulus and the modulus of subgrade reaction of the foundation soil; the reduction being more for their lower values. It is also observed that the variation in deflection of the plate at any radial distance due to change in soil subgrade conditions is significant for lower values of modulus of subgrade reaction.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

Evaluation of the effect of river style framework on water quality: application of geomorphological factors

The Tarwal River basin with an area of 6560.20 km² is located in the eastern part of Iranian Kurdistan Province. This river crosses the Qorveh and Dehgolan plains and joins the Ghezel Ozan River in Zanjan Province. The importance of this river as a source for drinking water and agricultural and industrial uses in the region necessitates the need for research in this field. The main purpose of this study is to identify the natural features of the riverbed from the perspective of river geomorphology and to investigate their impact on water quality and river self-purification capacity. To achieve this, the river style framework was employed. To investigate the effects of each style framework on the river, a total of 20 samples from the entrance and outlet of styles were obtained using Impact Assessment method and sampling standards which were later analyzed for their quality parameters including T, pH, EC, TDS, TSS, Na, Ca, Mg, K, Cl, F, NO₂, NO₃, SO₄, PO₄, DO, COD and BOD. The results indicated that the changes in the styles lead to changes in water quality and the impact of each style is greater on the physical parameters than the chemical parameters. The river self-purification capacity varied depending on the style. The maximum and the minimum self-purifications occurred in fine-grained Anabranching and low-sinuosity fine-grained styles, respectively.