An Analysis of Turbulent Heat Fluxes and the Energy Balance During the REFLEX Campaign

Three eddy covariance stations were installed at the Barrax experimental farm during the Land-Atmosphere Exchanges (REFLEX) airborne training and measurement campaign to provide ground truth data of energy balance fluxes and vertical temperature and wind profiles. The energy balance closure ratio (EBR) was 105% for a homogeneous camelina site, 86% at a sparse reforestation site, and 73% for a vineyard. We hypothesize that the lower closure in the last site was related to the limited fetch. Incorporating a vertical gradient of soil thermal properties decreased the RMSE of the energy balance at the camelina site by 16 W m^?2. At the camelina site, eddy covariance estimates of sensible and latent heat fluxes could be reproduced well using mean vertical profiles of wind and temperature, provided that the Monin—Obukhov length is known. Measured surface temperature and sensible heat fluxes suggested high excess resistance for heat (kB^?1 = 17).     

Intercomparison of the total storage deficit index (TSDI) over two Canadian Prairie catchments

Retrieval of the terrestrial moisture storage dataset from the Gravity Recovery and Climate Experiment (GRACE) satellite remote sensing system is possible when the catchment of interest is of large spatial scale. These dataset are of paramount importance for the estimation of the total storage deficit index (TSDI), which enables the characterization of a particular drought event from the perspective of the terrestrial moisture storage over that catchment. Incidentally, the GRACE gravity signal over the 13,000 km² Upper Assiniboine River Basin on the drought-prone Canadian Prairie is so poor therefore making the computation of the total storage deficit index for this basin infeasible. Consequently, the estimation of the terrestrial moisture storage from other reliable sources becomes imperative in order to enable the computation of the TSDI over this basin.This study explores the utilization of the Variable Infiltration Capacity (VIC) model, a physically based, spatially distributed hydrologic model to simulate the total moisture storage over the Upper Assiniboine River Basin which was then employed in the estimation of the TSDI over this basin for subsequent characterization of the recent Prairie-wide drought. Interestingly, the temporal patterns in the computed TSDI from the VIC model reveal a strong resemblance with the same drought characterization undertaken over the larger adjacent Saskatchewan River Basin, which was accomplished utilizing terrestrial moisture storage from the GRACE-based approach. Additionally, these independent techniques employed in the characterization of the last Prairie drought over the two adjacently situated basins resulted in similar drought severity classification from the standpoint of the total moisture storage deficits over these basins. This study has therefore shown that in the computation of the total storage deficit index over small-scale catchments during anomalous climatic conditions that propagate extreme dryness through the terrestrial hydrologic systems, simulations of the total water storage from a structurally sound model such as the VIC model could be resourceful for the computation of the monthly total storage deficit index if no constraint is placed on the availability of accurate meteorological forcing.     

全球大震和中国及邻区中强震地震活动(2012年9—10月)

本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出.本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).面波震级MS是对中周期宽频带SK地震仪记录,采用北京     

全球大震和中国及邻区中强震地震活动(2012年11-12月)

本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出. 本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).     

全球大震和中国及邻区中强震地震活动(2012年5-6月)

本目录中的地震参数来自“中国地震台站观测报告”(简称“月报”).其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件.“月报”由中国地震台网中心按月做出. 本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).     

全球大震和中国及邻区中强震地震活动(2012年7-8月)

本目录中的地震参数来自"中国地震台站观测报告"(简称"月报")。其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件。"月报"由中国地震台网中心按月做出。本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC)。震中位置除给出经     

全球大震和中国及邻区中强震地震活动(2012年3-4月)

本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出. 本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD). 面波震级Ms是对中周期宽频带SK地震仪记录,采用北京台1965年面波震级公式Ms=lg(AH/T)+1.66lg(△)+3.5(1°＜△＜130°)求得.AH是两水平分向最大面波位移的矢量合成位移.Ms7是对763长周期地震仪记录,采用国际上推荐的面波震级公式M7=lg(Av/T)+1.66lg(△)+3.3(20°＜△＜160°)求得.Av是垂直向面波最大地动位移.mb是短周期体波震级,ML是近震震级.为避免混乱,震级之间一律不换算.为方便读者,还给出美国NEIC定出的面波震级Msz和短周期体波震级mb.     

全球大震和中国及邻区中强震地震活动(2012年1-2月)

本目录中的地震参数来自“中国地震台站观测报告”(简称“月报”).其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件.“月报”由中国地震台网中心按月做出. 本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).     

Archaeomagnetic study of ceramics from the Neolithic Los Castillejos multilayer monument (Montefrio,Spain)

The geomagnetic intensity variation from the middle of the 4th millennium BC to the beginning of the 2nd millennium BC is reconstructed from the archaeomagnetic study of ceramic material taken from 24 layers of the Los Castillejos monument (Montefrio, Spain). The general patterns of the geomagnetic intensity variation reconstructed by studying materials from the Cendres Cave and Los Castillejos monuments (Spain) are similar. The intensity level is revised in the time interval including the minimum of the “fundamental” oscillation of the geomagnetic field; the characteristics of geomagnetic intensity variations whose superposition can provide the intensity variation observed in the study time interval are determined. Constraints on the climate humidity variation over the settlement lifetime are obtained from heating-induced variations in the magnetic susceptibility of the ceramics.     

Occurrence of perfluorinated compounds in water, sediment and mussels from the Cantabrian Sea (North Spain)

Perfluorinated compounds (PFCs) have emerged as significant global environmental pollutants with persistent, bioaccumulative and toxic properties. The aim of this study was to determine the occurrence of PFCs in water (wastewater, submarine emissaries and port-waters), sediment and transplanted mussels in estuarine areas of high urban and industrial impact from Northern Spain. Five PFCs of industrial use were studied: perfluorooctanesulfonate, perfluorohexanesulfonate, perfluorobutanesulfonate, perfluorooctanoate acid and perfluorononanoate acid. After selective extraction, samples were analyzed by Ultra Performance Liquid Chromatography coupled to tandem mass spectrometry. ΣPFCs ranged from 0.06 to 10.9 ng/L in water, with higher levels in wastewater treatment plants effluents and port waters than in submarine emissaries. Little accumulation was observed in sediments and mussels with ΣPFCs ranging from 0.01-0.13 ng/g dw and 0.01-0.06 ng/g ww, respectively. Most ubiquitous compounds were PFOS and PFOA. Mass fluxes of PFCs to the Cantabrian Sea are estimated and the impact to the coastal ecosystem is discussed.     

Water vapour mixing ratio profiles over Hornsund, Arctic. Intercomparison of lidar and AIRS results

Since October 2009, a ground-based Raman lidar system has been deployed to perform a regular, night-time, vertical sounding of a water vapour content in the lower and middle troposphere above Polish Polar Station at Hornsund (77.00°N, 15.55°E, 10 m a.s.l.) in the Arctic. The water vapour mixing ratio profiles were obtained for the atmosphere up to 6 km altitude, based on analysis of inelastic Raman backscattering signals from nitrogen molecules (at 387 nm) and water vapour particles (at 407 nm), calibrated with the data from a local Vaisala’s automatic meteorological station. The results obtained for winter seasons in the years 2009–2012 are in a good general agreement with the results obtained from the atmospheric infrared sounder (AIRS) on the Aqua satellite.     

Solar irradiance and total ozone over El Arenosillo (Spain) during the solar eclipse of 3 October 2005

The ground track of the annular eclipse of 3 October 2005 crossed the Iberian Peninsula. The main objective of this work was to analyze the variability of the solar irradiance and the total ozone column during the course of this event at El Arenosillo (Southwestern Spain). For achieving this goal, two Kipp & Zonen broadband radiometers (one for measuring total solar irradiance and other for measuring ultraviolet erythemal solar irradiance), one NILU-UV multi-band instrument and one Brewer spectroradiometer were used in this work. Total irradiance (310–2800 nm), and ultraviolet erythemal radiation (UVER) were recorded at a high frequency of 5 s, showing a strong reduction (higher than 80%) of the irradiance at the maximum solar obscuration which was of 79.6%. The irradiance decrease during the course of the eclipse was positively correlated with the percentage of eclipse obscuration, showing a very high agreement (R²～0.99). The irradiance recorded at selected wavelengths from the NILU-UV instrument shows a more pronounced decrease in the UV irradiance at the lower wavelengths during the solar eclipse. Finally, the evolution of the total ozone column (TOC) derived from Brewer and NILU instruments during the eclipse presented an opposite behavior: while the Brewer derived TOC values increase about 15 DU, the NILU derived TOC values decrease about 11 DU. This opposite behavior is mainly related to an artifact in the spectral irradiances recorded by the two instruments.     

Matrix hydrodynamic properties of carbonate rocks from the Betic Cordillera (Spain)

Laboratory tests performed on 181 rock samples from boreholes drilled in different areas of the Betic Cordillera allow us to calculate their hydraulic conductivity and open porosity values. Higher values are generally associated with Miocene calcareous sandstones, although hydraulic conductivity reaches its highest values in some isolated limestone and dolostone samples. Lower values were found in marly limestones and marbles. Specific yield ranged from 0 to 0·0798, with a mean value of 0·00579. A total of 79 samples did not release water during the specific yield test, while another 11 samples released water for more than 30 min. Such wide ranges of variation show the great diversity of behaviour that the matrix of the carbonate rocks can have, referring to water storage and transfer, and its influence on pollutants spread, for example. A weak relation between interconnected porosity and hydraulic conductivity was found. The relation between interconnected porosity and specific yield is slightly stronger, except in the case of the dolomites, where a high correlation was found. No dependence on depth was found for hydraulic conductivity and interconnected porosity. Copyright © 2004 John Wiley & Sons, Ltd.     

Luminescence and ESR dating of the sedimentary infill from the multi-level cave system of Alkerdi-Zelaieta (Navarre,N Spain)

The finding of Upper Paleolithic engravings in 2016 triggered a multidisciplinary investigation of the Alkerdi cave system (Urdazubi, N Spain). The study of the speleogenetic processes led to the identification of at least 6 paragenetic cave levels with associated sedimentary infill. In order to unravel the timing of changes in the karst dynamics and to get some insights about sediment origin, two sediment samples were collected from cave levels 4 and 1 for numerical dating purpose, using both Optically Stimulated Luminescence (OSL) and Electron Spin Resonance (ESR) methods. One additional modern-age sample was also taken from the entrance of the karst to evaluate the magnitude of the optical bleaching achieved by the quartz grains before entering the cave system. Last, one sample was collected from a flowstone for U-series dating, providing independent age control.OSL measurements were carried out using small quartz multi-grain (MG) and single-grain (SG) Single Aliquot Regenerative-dose (SAR) protocol. Initial Equivalent Dose (D_e) results show evidence of saturated OSL signal for the sample of the higher level (with ∼60% of saturated aliquots). Unlike SA data, SG analyses do not meet the usual quality criteria (e.g., relatively high recuperation ratios, only a few grains with useable signal), suggesting that corresponding dating results should be treated with caution.Hence, both MG and SG TT-OSL and MG ESR methods were subsequently employed, as the corresponding signals are known to have higher saturation levels than OSL ones. ESR analyses were based on the Multiple Centre (MC) approach using the standard multi-grain multi-aliquot additive (MAA) dose method, while TT-OSL measurements were performed following Demuro et al. (2020). The latter, however, did not return any useful results for both samples. MC ESR data show the usual D_e pattern (D_e(Al)>D_e(Ti op. D)>D_e(Ti–H)) indicating that among the three ESR signals analyzed, the Ti–H most likely provides the closest estimate to the true burial age. However, low measurement repeatability and goodness-of-of fit indicate that the reliability of the ESR results may be reasonably questioned.Results derived from this multi-technique dating approach provide the first chronological constraints for the sedimentary infill of the Alkerdi cave system. In particular, they suggest that sediment deposition in Cave Level 4 (ALK-OSL01) occurred around 130 ka, which is in good agreement with the minimum age constraint given by the speleothem (80.5 ± 9.0 ka). In comparison, the young age obtained for sample ALK-OSL02 (5.0 ± 0.9 ka) suggests recent Holocene formation of the lowest cave level 1.Finally, the MG De value obtained for the modern sample is close to zero (<1 Gy), indicating that the OSL signal is almost fully reset at a multi-grain level before entering the cave. Additionally, despite the large De overdispersion measured in this modern sample, the De values are one to two order of magnitude lower than those obtained with the same model in the other two samples. Therefore, the large OSL D_e overdispersion (OD) values of up to 64% obtained for the two samples from Cave Level 4 and 1 are most likely related to re-sedimentation processes inside the cave system.     

Characteristics of observed strong motion accelerograms from the 2011 Lorca (Spain) Earthquake

The Lorca earthquake (southeast Spain) on May 11, 2011; Mw 5.1, and its aftershocks, have provided an important set of accelerograms recorded by the strong motion network of the Instituto Geográfico Nacional of Spain. It is particularly noticeable that the record obtained in Lorca town, very near of the fault rupture, presents a PGA value close to 0.37 g in the N30W component. This paper provides an overview of the strong motion data recorded during the Lorca seismic series, with particular attention to the accelerograms from the mainshock and foreshock and its characteristics. Due to the special circumstances of these two records, a more detailed processing has been required, in which various alternatives to adjust the baseline have been considered and analyzed. Based on this special processing, small residual displacements were obtained and reported in some of these cases. However, given the sensitivity of the process and the small obtained values, these findings should be taken with caution. Besides, response spectra have been analyzed and compared with design spectra proposed by the Spanish Seismic Code (NCSE-02) for the towns of Lorca, Alhama de Murcia and Mula. Large amplitude differences were observed in these spectra if compared to those recorded in Lorca. Also noteworthy is that the design spectra proposed for this town were exceeded by some horizontal directions of the response spectra.     

Local Site Effects in the Town of Benevento (Italy) from Noise Measurements

— The study of ground motion amplification produced by surface geology is extremely interesting in the Benevento area, Southern Italy, as it is characterized by high seismic hazard. The present moderate-to-low seismicity makes the noise method appropriate to estimate the seismic site response in the area. The three components of seismic noise have been recorded in five sites in the Benevento metropolitan area characterized by different surface geology, in order to estimate the seismic site response. In evaluating site amplification effects we used the direct interpretation of amplitude spectra and standard spectral ratio techniques, evaluating sediment-to-bedrock, sediment-to-average and H/V spectral ratios. The temporal evolution of the noise spectra is analysed within one day, in order to assess the stationarity of the noise signal. The noise wavefield properties have been studied through polarization analyses in selected bands of frequency, where spectral peaks are observed to dominate, to better understand the real nature of those peaks. Results give evidence of low amplification levels, missing any correlation between spectral amplitudes and sediment thickness over the basement. We interpret this result as due to the poor impedance contrast between sediments and basement, which is characterized by low values of shear waves velocity. Moreover, sharp amplitude peaks are observed in the raw spectra of the sediment-sites, in the 2–4 Hz frequency band; a numerical simulation interprets this effect as possibly associated with a wide-scale structure, invoking the presence of a sharper impedance contrast at greater depth. At high frequencies the action of ambient noise sources, mainly active on horizontal components of motion, is retained dominant to generate the prominent peaks observed in the H/V spectral ratios; in some cases the presence of a near-surface low-velocity layer can contribute to amplify the seismic motion generated at these frequencies.     

Two-dimensional magnetotelluric characterization of the El Hito Anticline (Loranca Basin,Spain)

Magnetotelluric surveys were conducted to investigate the structure in the El Hito Anticline in the southern Loranca Basin, Spain. The structure of this anticline is partly known from nearby data (geological cross sections, geological maps, well-logs) from vintage geophysical and geological basin-scale surveys. Unfortunately, these surveys do not have the appropriate resolution to determine certain characteristics of the anticline accurately, such as the thickness or geometry of geological units. To address this deficiency, magnetotelluric data were acquired at 51 sites along three profiles to image the electrical resistivity of the anticline. To identify a geologically reasonable resistivity model, a conceptual model based on previous geological and geophysical information (cross sections and well-log data) was generated. Several inversions were performed using the conceptual model, which played a key role in the interpretation of the magnetotelluric data, and the construction of a suitable initial model was essential in producing geologically meaningful models. Using these results, we obtained 2-D resistivity models that provide information on the main geological units and structures. We also performed sensitivity tests to understand the morphology and structure of the resistive basement better. The results indicate the presence of basement highs that show different structural styles for the basement and its cover, and an increase in the basement depth towards the south. The final 2-D resistivity models provide new information regarding the structure of the anticline and demonstrate a clear correlation between the main geological units and resistivity, as well as correlations with mapped surface faults.     

Fluxes of fluid and heat from the oceanic crustal reservoir

Recent discoveries define a global scale fluid reservoir residing within the uppermost igneous oceanic crust, a region of seafloor that is both warm and may harbor a substantial biosphere. This hydrothermal fluid reservoir formed initially within volcanic rocks newly erupted at mid-ocean ridges, but extends to the vastly larger and older ridge flanks. Upper oceanic crust is porous and permeable due to the presence of lava drainbacks, fissuring, and inter-unit voids, and this porosity and permeability allows active fluid circulation to advect measurable quantities of lithospheric heat from the crust to an average age of 65 Myr. A compilation of crustal porosities shows that this fluid reservoir contains nearly 2% of the total volume of global seawater. Heat flow and sediment thickness data allow calculation of reservoir temperatures, predicting 40°C mean temperatures in Cretaceous crust. Utilizing these temperature estimates, heat flow measurements and models for the thermal structure and evolution of the oceanic lithosphere, we have computed mean hydrothermal fluxes into the deep ocean as a function of plate age. The total hydrothermal volume flux into the oceans approaches 20% of the total riverine input and may contribute to the global seawater mass balance.     

Seismicity pattern of the Betic Cordillera(southern Spain) derived from the fractal properties of earthquakes and faults

Several studies on earthquake occurrence and associated faulting have demonstrated that both phenomena have a scale-invariant behavior which can be analyzed by means of a set of non-integer dimensions(Dq) describing their fractal properties and the calculation of multi-fractal spectra.It is the case that the behavior of these spectra is asymptotic at the ends of the variation interval of q,which is a real number that enters into the definition of the partition function of the dataset.The difference between the extreme values,called multi-fractal spectrum slope,is used to investigate the heterogeneity of the spatial distribution of earthquakes and fault systems.In this paper we focus on the Betic Cordillera,southeastern Spain,which is commonly considered the contact between the Eurasian and African plates and has an important seismic activity in the context of the Iberian Peninsula.Some of the most conspicuous Iberian earthquakes,such as the 1829 mb6.3 Torrevieja and the 1884 mb6.1 Alhama de Granada earthquakes occurred in this mountain range and both reached intensity X.The present work implies a new analysis based on the slope of multi-fractal spectra and referred to the historical seismicity of the region,specifically b-value(frequency distribution of earthquakes respect to magnitude),epicentral location,seismic energy and faulting.On this basis we propose a seismotectonic zonation that is contrasted with the stress state and the geodynamical evolution of the Betic Cordillera.