The debated question on the possible relation between the Earth's magnetic field and climate has been usually focused on direct correlations between different time series representing both systems. ...However, the physical mechanism able to potentially explain this connection is still an open issue. Finding hints about how this connection could work would suppose an important advance in the search of an adequate physical mechanism. Here, we propose an innovative information-theoretic tool, i.e. the transfer entropy, as a good candidate for this scope because is able to determine, not simply the possible existence of a connection, but even the direction in which the link is produced. We have applied this new methodology to two real time series, the South Atlantic Anomaly (SAA) area extent at the Earth's surface (representing the geomagnetic field system) and the Global Sea Level (GSL) rise (for the climate system) for the last 300 years, to measure the possible information flow and sense between them. This connection was previously suggested considering only the long-term trend while now we study this possibility also in shorter scales. The new results seem to support this hypothesis, with more information transferred from the SAA to the GSL time series, with about 90% of confidence level. This result provides new clues on the existence of a link between the geomagnetic field and the Earth's climate in the past and on the physical mechanism involved because, thanks to the application of the transfer entropy, we have determined that the sense of the connection seems to go from the system that produces geomagnetic field to the climate system. Of course, the connection does not mean that the geomagnetic field is fully responsible for the climate changes, rather that it is an important driving component to the variations of the climate.
The study of the preparation phase of large earthquakes is essential to understand the physical processes involved, and potentially useful also to develop a future reliable short-term warning system. ...Here we analyse electron density and magnetic field data measured by Swarm three-satellite constellation for 4.7 years, to look for possible in-situ ionospheric precursors of large earthquakes to study the interactions between the lithosphere and the above atmosphere and ionosphere, in what is called the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC). We define these anomalies statistically in the whole space-time interval of interest and use a Worldwide Statistical Correlation (WSC) analysis through a superposed epoch approach to study the possible relation with the earthquakes. We find some clear concentrations of electron density and magnetic anomalies from more than two months to some days before the earthquake occurrences. Such anomaly clustering is, in general, statistically significant with respect to homogeneous random simulations, supporting a LAIC during the preparation phase of earthquakes. By investigating different earthquake magnitude ranges, not only do we confirm the well-known Rikitake empirical law between ionospheric anomaly precursor time and earthquake magnitude, but we also give more reliability to the seismic source origin for many of the identified anomalies.
Fires and the aerosols that they emit impact air quality, health, and climate, but the abundance and properties of carbonaceous aerosol (both black carbon and organic carbon) from biomass burning ...(BB) remain uncertain and poorly constrained. We aim to explore the uncertainties associated with fire emissions and their air quality and radiative impacts from underlying dry matter consumed and emissions factors. To investigate this, we compare model simulations from a global chemical transport model, GEOS-Chem, driven by a variety of fire emission inventories with surface and airborne observations of black carbon (BC) and organic aerosol (OA) concentrations and satellite-derived aerosol optical depth (AOD). We focus on two fire-detection-based and/or burned-area-based (FD-BA) inventories using burned area and active fire counts, respectively, i.e., the Global Fire Emissions Database version 4 (GFED4s) with small fires and the Fire INventory from NCAR version 1.5 (FINN1.5), and two fire radiative power (FRP)-based approaches, i.e., the Quick Fire Emission Dataset version 2.4 (QFED2.4) and the Global Fire Assimilation System version 1.2 (GFAS1.2). We show that, across the inventories, emissions of BB aerosol (BBA) differ by a factor of 4 to 7 over North America and that dry matter differences, not emissions factors, drive this spread. We find that simulations driven by QFED2.4 generally overestimate BC and, to a lesser extent, OA concentrations observations from two fire-influenced aircraft campaigns in North America (ARCTAS and DC3) and from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network, while simulations driven by FINN1.5 substantially underestimate concentrations. The GFED4s and GFAS1.2-driven simulations provide the best agreement with OA and BC mass concentrations at the surface (IMPROVE), BC observed aloft (DC3 and ARCTAS), and AOD observed by MODIS over North America. We also show that a sensitivity simulation including an enhanced source of secondary organic aerosol (SOA) from fires, based on the NOAA Fire Lab 2016 experiments, produces substantial additional OA; however, the spread in the primary emissions estimates implies that this magnitude of SOA can be neither confirmed nor ruled out when comparing the simulations against the observations explored here. Given the substantial uncertainty in fire emissions, as represented by these four emission inventories, we find a sizeable range in 2012 annual BBA PM2.5 population-weighted exposure over Canada and the contiguous US (0.5 to 1.6 µg/cu. m). We also show that the range in the estimated global direct radiative effect of carbonaceous aerosol from fires (−0.11 to −0.048 W/sq. m) is large and comparable to the direct radiative forcing of OA (−0.09 W/sq. m) estimated in the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC). Our analysis suggests that fire emissions uncertainty challenges our ability to accurately characterize the impact of smoke on air quality and climate.
•New geomagnetic field reconstruction (SHAWQ2k) has been obtained for the last 2 ka.•A novelty weighting scheme based on quality-data is applied in modeling process.•New hints about the past ...dynamical behavior of the South Atlantic Anomaly are shown.•Reversed flux patch linked to the South Atlantic Anomaly is detected since 950 AD.
The South Atlantic Anomaly (SAA) is one of the most outstanding features of the present geomagnetic field. Thus, a good knowledge of the SAA is clue for a better understanding of the dynamical behavior of the geomagnetic field. To achieve this goal, paleomagnetic data are crucial since they provide the unique way to investigate past changes in the Earth's magnetic field. Here, we present a new global geomagnetic field reconstruction, the SHAWQ2k model, which is based on a critical revision of the global archeomagnetic and volcanic dataset. The new model provides an improved description of the geomagnetic field during the last 2 millennia, and yields surprising outcomes about the emergence and development of the SAA. It shows that the reversed flux patch observed at the core-mantle boundary and linked to the SAA, emerged in the Southern Hemisphere from at least 950 AD. This patch moved westward from the equator to southern latitudes, being clearly isolated after 1400 AD. In addition, since 1550 AD a second reversed flux patch moving northeastward is observed in the North Atlantic. The new data now available for the Southern Hemisphere coming from Africa and South America together with the use of an appropriated weighting scheme in the modeling process have improved our understanding of past geomagnetic field behavior and showed new evolutionary aspects of the SAA.
In this work, we present 16 directional and 27 intensity high‐quality values from Iberia. Moreover, we have updated the Iberian archeomagnetic catalogue published more than 10 years ago with a ...considerable increase in the database. This has led to a notable improvement of both temporal and spatial data distribution. A full vector paleosecular variation curve from 1000 BC to 1900 AD has been developed using high‐quality data within a radius of 900 km from Madrid. A hierarchical bootstrap method has been followed for the computation of the curves. The most remarkable feature of the new curves is a notable intensity maximum of about 80 μT around 600 BC, which has not been previously reported for the Iberian Peninsula. We have also analyzed the evolution of the paleofield in Europe for the last three thousand years and conclude that the high maximum intensity values observed around 600 BC in the Iberian Peninsula could respond to the same feature as the Levantine Iron Age Anomaly, after travelling westward through Europe.
Plain language summary
Knowledge of the Earth's magnetic field plays an important role on the understanding of its dynamics. By measuring certain rocks or archeological objects from around the world, we can determine the field's shape and intensity in former times. Knowing its evolution is essential to understand how this field is generated, how it has varied through time and how it may behave in the future. In this work, we present new measurements of the magnetic field from the Iberian Peninsula that provide useful constraints on the magnetic field for archeological times that currently lack information. We have updated the compilation of Iberian data for the last 3,000 years and calculated a new reference curve for the magnetic field for this region. We have found that the magnetic field was particularly intense in the Iberian Peninsula about 2,600 years ago. By comparing this result with data from Europe and the Middle East, we observe that this high intensity has been moving from east to west through southern Europe. This feature is probably related with the rapid intensity change (the geomagnetic spike) recently discovered in the Levantine region.
Key Points
We have obtained 16 archeomagnetic directions and 27 high‐quality archeointensity data
A new full vector paleosecular variation curve for the last three millennia in Iberia is presented
The new curve shows a high‐intensity maximum around 600 BC that might be related to the Levantine Iron Age Anomaly moving through Europe
This work deals with a comprehensive multiparametric and multilayer approach to study earthquake-related processes that occur during the preparation phase of a large earthquake. As a case study, the ...paper investigates the M7.2 Kermadec Islands (New Zealand) large earthquake that occurred on June 15, 2019 as the result of shallow reverse faulting within the Tonga-Kermadec subduction zone. The analyses focused on seismic (earthquake catalogs), atmospheric (climatological archives) and ionospheric data from ground to space (mainly satellite) in order to disclose the possible Lithosphere-Atmosphere-Ionosphere Coupling (LAIC). The ionospheric investigations analysed and compared the Global Navigation Satellite System (GNSS) receiver network with in-situ observations from space thanks to both the European Space Agency (ESA) Swarm constellation and the China National Space Administration (CNSA in partnership with Italian Space Agency, ASI) satellite dedicated to search for possible ionospheric disturbances before medium-large earthquakes, i.e. the China Seismo-Electromagnetic Satellite (CSES-01). An interesting comparison is made with another subsequent earthquake with comparable magnitude (M7.1) that occurred in Ridgecrest, California (USA) on 6 July of the same year but in a different tectonic context. Both earthquakes showed anomalies in several parameters (e.g. aerosol, skin temperature and some ionospheric quantities) that appeared at almost the same times before each earthquake occurrence, evidencing a chain of processes that collectively point to the moment of the corresponding mainshock. In both cases, it is demonstrated that a comprehensive multiparametric and multilayer analysis is fundamental to better understand the LAIC in the occasion of complex phenomena such as earthquakes.
•The 2019 M7.2 Kermadec Islands, New Zealand earthquake was analysed to identify pre-earthquake anomalies from ground to space.•A comprehensive multiparametric and multilayer approach was applied.•A chain of anomalies was found from lithosphere to atmosphere and ionosphere.•Comparable pre-earthquake anomaly times with those of M7.1 Ridgecrest earthquake.•Cumulative number of anomalies points to the earthquake occurrence.
The stimulation of electromagnetic ion cyclotron (EMIC) waves by a magnetospheric compression is perhaps the closest thing to a controlled experiment that is currently possible in magnetospheric ...physics, in that one prominent factor that can increase wave growth acts at a well‐defined time. We present a detailed analysis of EMIC waves observed in the outer dayside magnetosphere by the four Magnetosphere Multiscale (MMS) spacecraft, Van Allen Probe A, and GOES 13 and by four very high latitude ground magnetometer stations in the western hemisphere before, during, and after a modest interplanetary shock on 14 December 2015. Analysis shows several features consistent with current theory, as well as some unexpected features. During the most intense MMS wave burst, which began ~ 1 min after the end of a brief magnetosheath incursion, independent transverse EMIC waves with orthogonal linear polarizations appeared simultaneously at all four spacecraft. He++ band EMIC waves were observed by MMS inside the magnetosphere, whereas almost all previous studies of He++ band EMIC waves observed them only in the magnetosheath and magnetopause boundary layers. Transverse EMIC waves also appeared at Van Allen Probe A and GOES 13 very near the times when the magnetic field compression reached their locations, indicating that the compression lowered the instability threshold to allow for EMIC wave generation throughout the outer dayside magnetosphere. The timing of the EMIC waves at both MMS and Van Allen Probe A was consistent with theoretical expectations for EMIC instabilities based on characteristics of the proton distributions observed by instruments on these spacecraft.
Key Points
MMS observed a burst of independent transverse EMIC waves with orthogonal linear polarizations following the shock
Wave onsets at both MMS and Van Allen Probe A were consistent with theoretical expectations based on particle observations
A rarely observed minimum in wave power at the gyrofrequency of He++ ions was present both before and after the shock
Monitoring levels and trends in premature mortality is crucial to understanding how societies can address prominent sources of early death. The Global Burden of Disease 2016 Study (GBD 2016) provides ...a comprehensive assessment of cause-specific mortality for 264 causes in 195 locations from 1980 to 2016. This assessment includes evaluation of the expected epidemiological transition with changes in development and where local patterns deviate from these trends.
We estimated cause-specific deaths and years of life lost (YLLs) by age, sex, geography, and year. YLLs were calculated from the sum of each death multiplied by the standard life expectancy at each age. We used the GBD cause of death database composed of: vital registration (VR) data corrected for under-registration and garbage coding; national and subnational verbal autopsy (VA) studies corrected for garbage coding; and other sources including surveys and surveillance systems for specific causes such as maternal mortality. To facilitate assessment of quality, we reported on the fraction of deaths assigned to GBD Level 1 or Level 2 causes that cannot be underlying causes of death (major garbage codes) by location and year. Based on completeness, garbage coding, cause list detail, and time periods covered, we provided an overall data quality rating for each location with scores ranging from 0 stars (worst) to 5 stars (best). We used robust statistical methods including the Cause of Death Ensemble model (CODEm) to generate estimates for each location, year, age, and sex. We assessed observed and expected levels and trends of cause-specific deaths in relation to the Socio-demographic Index (SDI), a summary indicator derived from measures of average income per capita, educational attainment, and total fertility, with locations grouped into quintiles by SDI. Relative to GBD 2015, we expanded the GBD cause hierarchy by 18 causes of death for GBD 2016.
The quality of available data varied by location. Data quality in 25 countries rated in the highest category (5 stars), while 48, 30, 21, and 44 countries were rated at each of the succeeding data quality levels. Vital registration or verbal autopsy data were not available in 27 countries, resulting in the assignment of a zero value for data quality. Deaths from non-communicable diseases (NCDs) represented 72·3% (95% uncertainty interval UI 71·2–73·2) of deaths in 2016 with 19·3% (18·5–20·4) of deaths in that year occurring from communicable, maternal, neonatal, and nutritional (CMNN) diseases and a further 8·43% (8·00–8·67) from injuries. Although age-standardised rates of death from NCDs decreased globally between 2006 and 2016, total numbers of these deaths increased; both numbers and age-standardised rates of death from CMNN causes decreased in the decade 2006–16—age-standardised rates of deaths from injuries decreased but total numbers varied little. In 2016, the three leading global causes of death in children under-5 were lower respiratory infections, neonatal preterm birth complications, and neonatal encephalopathy due to birth asphyxia and trauma, combined resulting in 1·80 million deaths (95% UI 1·59 million to 1·89 million). Between 1990 and 2016, a profound shift toward deaths at older ages occurred with a 178% (95% UI 176–181) increase in deaths in ages 90–94 years and a 210% (208–212) increase in deaths older than age 95 years. The ten leading causes by rates of age-standardised YLL significantly decreased from 2006 to 2016 (median annualised rate of change was a decrease of 2·89%); the median annualised rate of change for all other causes was lower (a decrease of 1·59%) during the same interval. Globally, the five leading causes of total YLLs in 2016 were cardiovascular diseases; diarrhoea, lower respiratory infections, and other common infectious diseases; neoplasms; neonatal disorders; and HIV/AIDS and tuberculosis. At a finer level of disaggregation within cause groupings, the ten leading causes of total YLLs in 2016 were ischaemic heart disease, cerebrovascular disease, lower respiratory infections, diarrhoeal diseases, road injuries, malaria, neonatal preterm birth complications, HIV/AIDS, chronic obstructive pulmonary disease, and neonatal encephalopathy due to birth asphyxia and trauma. Ischaemic heart disease was the leading cause of total YLLs in 113 countries for men and 97 countries for women. Comparisons of observed levels of YLLs by countries, relative to the level of YLLs expected on the basis of SDI alone, highlighted distinct regional patterns including the greater than expected level of YLLs from malaria and from HIV/AIDS across sub-Saharan Africa; diabetes mellitus, especially in Oceania; interpersonal violence, notably within Latin America and the Caribbean; and cardiomyopathy and myocarditis, particularly in eastern and central Europe. The level of YLLs from ischaemic heart disease was less than expected in 117 of 195 locations. Other leading causes of YLLs for which YLLs were notably lower than expected included neonatal preterm birth complications in many locations in both south Asia and southeast Asia, and cerebrovascular disease in western Europe.
The past 37 years have featured declining rates of communicable, maternal, neonatal, and nutritional diseases across all quintiles of SDI, with faster than expected gains for many locations relative to their SDI. A global shift towards deaths at older ages suggests success in reducing many causes of early death. YLLs have increased globally for causes such as diabetes mellitus or some neoplasms, and in some locations for causes such as drug use disorders, and conflict and terrorism. Increasing levels of YLLs might reflect outcomes from conditions that required high levels of care but for which effective treatments remain elusive, potentially increasing costs to health systems.
Bill & Melinda Gates Foundation.
Reconstruction of Paleomagnetic Secular Variation (PSV) of the geomagnetic field is fundamental both to assess geodynamo models and to obtain age constraints for rocks, sediments and archaeological ...material. We present refined age-calibrated Holocene PSV and relative paleointensity (RPI) stack curves derived from Arctic marine sediments (Northwestern Barents Sea). The Holocene sections of four sedimentary cores were correlated on the basis of paleomagnetic trends and age models, and stacked. The resultant composite PSV and RPI Holocene records (NBS stack) and the reconstructed Holocene Virtual Geomagnetic Pole (VGP) path were evaluated in comparison with the most recent paleomagnetic stack curves and geomagnetic field models. The data indicate that during the Holocene time, the VGPs moved within the superficial projection of the inner core tangent cylinder, with the exception of short time intervals around 5600 and 3200 cal yr BP when VGPs extended to lower latitudes. These deviations might reflect regional geomagnetic features, such as persistent geomagnetic flux lobes at core-mantle boundary. Our data confirm that the large VGP shift observed around 5600 cal yr BP is the result of an increased radial magnetic field at the core-mantle boundary over North America, whilst the VGP shift around 3200 cal yr BP represents a major swing to middle latitudes toward the Middle East and might be associated to a regional high paleointensity peak, known as Levantine Iron Age Anomaly (LIAA).
•Sedimentary cores from Northwestern Barents Sea were inter-correlated and stacked.•Paleosecular variation and relative paleointensity stacks from Arctic region are presented.•Holocene Virtual Geomagnetic Pole (VGP) path at high latitude was reconstructed.•VGP reached lower latitudes during two intervals of rapid geomagnetic variations.•VGP shifts are related to geomagnetic flux lobes at the core-mantle boundary.
•New 37 archeointensities for the Central Mediterranean (1500 BCE–150 CE).•First evidence of a rapid intensity change between 1070 and 1040 BCE in Greece.•The Levantine Iron Age Anomaly (LIAA) is ...present in the Central Mediterranean.•A global analysis suggests that the LIAA vanishes to the west and to the east.
The magnitude and origin of the Levantine Iron Age geomagnetic Anomaly (LIAA), which spanned the first half of the first millennium before the common era, are not yet well understood. Recent archeomagnetic studies from the Levant and Western Europe suggest a western drift of this feature, stressing the importance of investigating the temporal and spatial behaviour of this event over the Central Mediterranean area. To analyse this issue, we here present 37 new archeointensity data obtained from the archeomagnetic study of 118 ceramics and brick fragments collected in 8 archeological sites in Greece and Italy with ages ranging between 1500 BCE and 150 CE. The samples were analysed using the classical Thellier and Thellier method for paleointensity determination, including the correction for the anisotropy effect of the thermoremanent magnetization (TRM) and for the cooling rate dependence upon TRM acquisition. The results reveal the first evidence of a high-intensity peak in Greece between 1070 and 1040 BCE associated to high virtual axial dipole moment (VADM) values of around 140 ZAm2. A global analysis of available paleointensities suggests that the origin of these high values is the same to the one which produced the maximum VADM of the LIAA in the Levantine region. Our results suggest that the source of the LIAA is located in the Levantine region vanishing to the north, to the west and to the east where lower VADMs are observed. In addition, another high intensity maximum, less pronounced than the one of the LIAA, seems to be present around 500 BCE all over Europe, from the Canary Islands to Turkey showing similar VADM values (around 150 ZAm2) in the different regions. Both events seem to span over a large region at the Earth's surface covering more than 60° of longitude, verifying an Earth's outer core origin for these intensity features.