Back to the Roots: Reconstructing Large and Complex Cranial Defects using an Image-based Statistical Shape Model.

Designing implants for large and complex cranial defects is a challenging task, even for professional designers. Current efforts on automating the design process focused mainly on convolutional neural networks (CNN), which have produced state-of-the-art results on reconstructing synthetic defects. However, existing CNN-based methods have been difficult to translate to clinical practice in cranioplasty, as their performance on large and complex cranial defects remains unsatisfactory. In this paper, we present a statistical shape model (SSM) built directly on the segmentation masks of the skulls represented as binary voxel occupancy grids and evaluate it on several cranial implant design datasets. Results show that, while CNN-based approaches outperform the SSM on synthetic defects, they are inferior to SSM when it comes to large, complex and real-world defects. Experienced neurosurgeons evaluate the implants generated by the SSM to be feasible for clinical use after minor manual corrections. Datasets and the SSM model are publicly available at https://github.com/Jianningli/ssm .

A Statistical Approach for the Integration of Multi-Temporal InSAR and GNSS-PPP Ground Deformation Measurements.

Determining and monitoring ground deformations is critical for hazard management studies, especially in megacities, and these studies might help prevent future disaster conditions and save many lives. In recent years, the Golden Horn, located in the southeast of the European part of Istanbul within a UNESCO-protected region, has experienced significant changes and regional deformations linked to rapid population growth, infrastructure work, and tramway construction. In this study, we used Interferometric Synthetic Aperture Radar (InSAR) and Global Navigation Satellite System (GNSS) techniques to investigate the ground deformations along the Golden Horn coastlines. The investigated periods are between 2015 and 2020 and 2017 and 2020 for InSAR and GNSS, respectively. For the InSAR analyses, we used sequences of multi-temporal synthetic aperture radar (SAR) images collected by the Sentinel-1 and ALOS-2 satellites. The ground displacement products (i.e., time series and velocity maps) were then cross-compared with those achievable using the Precise Point Positioning (PPP) technique for the GNSS solutions, which can provide precise positions with a single receiver. In the proposed analysis, we compared the ground displacement velocities obtained by both methods by computing the standard deviations of the difference between the relevant observations considering a weighted least square estimation procedure. Additionally, we identified five circle buffers with different radii ranging between 50 m and 250 m for selecting the most appropriate coherent points to conduct the cross-comparison analysis. Moreover, a vertical displacement rate map was produced. The comparison of the vertical ground velocities derived from PPP and InSAR demonstrates that the PPP technique is valuable. For the coherent stations, the vertical displacement rates vary between -4.86 mm/yr and -23.58 mm/yr and -9.50 and -27.77 mm/yr for InSAR and GNSS, respectively.

Predicting Complicated Mild Traumatic Brain Injury in Adolescent Trauma to Enhance Clinical Decisions in Imaging.

BACKGROUND: The Pediatric Emergency Care Applied Research Network (PECARN) traumatic brain injury algorithm is used to identify children at low risk of clinically significant traumatic brain injuries to reduce computed tomography (CT) exposure. Adapting PECARN rules based on population-specific risk stratification has been suggested to improve diagnostic accuracy. OBJECTIVE: This study sought to identify center-specific patient variables, beyond PECARN rules, that may enhance the identification of patients requiring neuroimaging. METHODS: This single-center, retrospective cohort study was conducted from July 1, 2016, to July 1, 2020, in a Southwestern U.S. Level II pediatric trauma center. The inclusion criteria were adolescents (10-15 years), Glasgow Coma Scale (13-15), with a confirmed mechanical blow to the head. Patients without a head CT were excluded. Logistic regression was performed to identify additional complicated mild traumatic brain injury predictor variables beyond the PECARN. RESULTS: There were 136 patients studied; 21 (15%) presented with a complicated mild traumatic brain injury. Relative to motorcycle collision or all-terrain vehicle trauma (odds ratio [OR] 211.75, 95% confidence interval, CI [4.51, 9931.41], p < .001), an unspecified mechanism (OR 42.0, 95% CI [1.30, 1350.97], p = .03) and consult activation (OR 17.44, 95% CI [1.75, 173.31], p = .01) were significantly associated with complicated mild traumatic brain injury. CONCLUSIONS: We identified additional factors associated with complex mild traumatic brain injury, including motorcycle collision and all-terrain vehicle trauma, unspecified mechanism, and consult activation that are not in the PECARN imaging decision rule. Adding these variables may aid in determining the need for appropriate CT scanning.

Studierfenster: an Open Science Cloud-Based Medical Imaging Analysis Platform.

Imaging modalities such as computed tomography (CT) and magnetic resonance imaging (MRI) are widely used in diagnostics, clinical studies, and treatment planning. Automatic algorithms for image analysis have thus become an invaluable tool in medicine. Examples of this are two- and three-dimensional visualizations, image segmentation, and the registration of all anatomical structure and pathology types. In this context, we introduce Studierfenster ( www.studierfenster.at ): a free, non-commercial open science client-server framework for (bio-)medical image analysis. Studierfenster offers a wide range of capabilities, including the visualization of medical data (CT, MRI, etc.) in two-dimensional (2D) and three-dimensional (3D) space in common web browsers, such as Google Chrome, Mozilla Firefox, Safari, or Microsoft Edge. Other functionalities are the calculation of medical metrics (dice score and Hausdorff distance), manual slice-by-slice outlining of structures in medical images, manual placing of (anatomical) landmarks in medical imaging data, visualization of medical data in virtual reality (VR), and a facial reconstruction and registration of medical data for augmented reality (AR). More sophisticated features include the automatic cranial implant design with a convolutional neural network (CNN), the inpainting of aortic dissections with a generative adversarial network, and a CNN for automatic aortic landmark detection in CT angiography images. A user study with medical and non-medical experts in medical image analysis was performed, to evaluate the usability and the manual functionalities of Studierfenster. When participants were asked about their overall impression of Studierfenster in an ISO standard (ISO-Norm) questionnaire, a mean of 6.3 out of 7.0 possible points were achieved. The evaluation also provided insights into the results achievable with Studierfenster in practice, by comparing these with two ground truth segmentations performed by a physician of the Medical University of Graz in Austria. In this contribution, we presented an online environment for (bio-)medical image analysis. In doing so, we established a client-server-based architecture, which is able to process medical data, especially 3D volumes. Our online environment is not limited to medical applications for humans. Rather, its underlying concept could be interesting for researchers from other fields, in applying the already existing functionalities or future additional implementations of further image processing applications. An example could be the processing of medical acquisitions like CT or MRI from animals [Clinical Pharmacology & Therapeutics, 84(4):448-456, 68], which get more and more common, as veterinary clinics and centers get more and more equipped with such imaging devices. Furthermore, applications in entirely non-medical research in which images/volumes need to be processed are also thinkable, such as those in optical measuring techniques, astronomy, or archaeology.

On the Use of Weighted Least-Squares Approaches for Differential Interferometric SAR Analyses: The Weighted Adaptive Variable-lEngth (WAVE) Technique.

This paper concentrates on the study of the Weighted Least-squares (WLS) approaches for the generation of ground displacement time-series through Differential Interferometric SAR (DInSAR) methods. Usually, within the DInSAR framework, the Weighted Least-squares (WLS) techniques have principally been applied for improving the performance of the phase unwrapping operations as well as for better conveying the inversion of sequences of unwrapped interferograms to generate ground displacement maps. In both cases, the identification of low-coherent areas, where the standard deviation of the phase is high, is requested. In this paper, a WLS method that extends the usability of the Multi-Temporal InSAR (MT-InSAR) Small Baseline Subset (SBAS) algorithm in regions with medium-to-low coherence is presented. In particular, the proposed method relies on the adaptive selection and exploitation, pixel-by-pixel, of the medium-to-high coherent interferograms, only, so as to discard the noisy phase measurements. The selected interferometric phase values are then inverted by solving a WLS optimization problem. Noteworthy, the adopted, pixel-dependent selection of the "good" interferograms to be inverted may lead the available SAR data to be grouped into several disjointed subsets, which are then connected, exploiting the Weighted Singular Value Decomposition (WSVD) method. However, in some critical noisy regions, it may also happen that discarding of the incoherent interferograms may lead to rejecting some SAR acquisitions from the generated ground displacement time-series, at the cost of the reduced temporal sampling of the data measurements. Thus, variable-length ground displacement time-series are generated. The mathematical framework of the developed technique, which is named Weighted Adaptive Variable-lEngth (WAVE), is detailed in the manuscript. The presented experiments have been carried out by applying the WAVE technique to a SAR dataset acquired by the COSMO-SkyMed (CSK) sensors over the Basilicata region, Southern Italy. A cross-comparison analysis between the conventional and the WAVE method has also been provided.

A Phase-Preserving Focusing Technique for TOPS Mode SAR Raw Data Based on Conventional Processing Methods.

We present a new solution for the phase-preserving focusing of synthetic aperture radar (SAR) raw data acquired through the Terrain Observation with Progressive Scan (TOPS) mode. The proposed algorithm consists of a first interpolation stage of the TOPS raw data, which takes into account the Doppler Centroid frequency variations due to the azimuth antenna steering function, and allows us to unfold the azimuth spectra of the TOPS raw data. Subsequently, the interpolated signals are processed by using conventional phase-preserving SAR focusing methods that exploit frequency domain and spectral analyses algorithms, which are extensively used to efficiently process Stripmap and ScanSAR data. Accordingly, the developed focusing approach is easy to implement. In particular, the presented focusing approach exploits one of the available frequency domain Stripmap processing techniques. The only modification is represented by the inclusion, within the 2D frequency domain focusing step, of a spurious azimuth chirp signal with a properly selected azimuthal rate. This allows us to efficiently carry out the TOPS azimuth focusing through the SPECAN method. Furthermore, an important aspect of this algorithm is the possibility to easily achieve a constant and tunable output azimuth pixel size without any additional computing time; this is a remarkable feature with respect to the full-aperture TOPS-mode algorithms available in the existing literature. Moreover, although tailored on Sentinel-1 (S1) raw data, the proposed algorithm can be easily extended to process data collected through the TOPS mode by different radar sensors. The presented experimental results have been obtained by processing real Sentinel-1 raw data and confirm the effectiveness of the proposed algorithm.

A Marker-Less Registration Approach for Mixed Reality-Aided Maxillofacial Surgery: a Pilot Evaluation.

As of common routine in tumor resections, surgeons rely on local examinations of the removed tissues and on the swiftly made microscopy findings of the pathologist, which are based on intraoperatively taken tissue probes. This approach may imply an extended duration of the operation, increased effort for the medical staff, and longer occupancy of the operating room (OR). Mixed reality technologies, and particularly augmented reality, have already been applied in surgical scenarios with positive initial outcomes. Nonetheless, these methods have used manual or marker-based registration. In this work, we design an application for a marker-less registration of PET-CT information for a patient. The algorithm combines facial landmarks extracted from an RGB video stream, and the so-called Spatial-Mapping API provided by the HMD Microsoft HoloLens. The accuracy of the system is compared with a marker-based approach, and the opinions of field specialists have been collected during a demonstration. A survey based on the standard ISO-9241/110 has been designed for this purpose. The measurements show an average positioning error along the three axes of (x, y, z) = (3.3 ± 2.3, - 4.5 ± 2.9, - 9.3 ± 6.1) mm. Compared with the marker-based approach, this shows an increment of the positioning error of approx. 3 mm along two dimensions (x, y), which might be due to the absence of explicit markers. The application has been positively evaluated by the specialists; they have shown interest in continued further work and contributed to the development process with constructive criticism.

Unusual Antioxidant Properties of 26S Proteasome Isolated from Cold-Adapted Organisms.

The oxidative challenge represents an important factor affecting the adaptive strategies in Antarctic fish, but their impact on the protein degradation machinery still remains unclear. The previous analysis of the first 26S proteasome from the Antarctic red-blooded fish Trematomus bernacchii, evidenced improved antioxidant functions necessary to counteract the environmental pro-oxidant conditions. The purpose of this work was to carry out a study on 26S proteasomes from the temperate red-blooded Dicenthrarcus labrax and the icefish Chionodraco hamatus in comparison with the isoform already described from T. bernacchii, to better elucidate the cold-adapted physiological functions of this complex. Therefore, the 26S isoforms were isolated and the complementary DNAs (cDNAs) codifying the catalytic subunits were cloned. The biochemical characterization of Antarctic 26S proteasomes revealed their significantly higher structural stability and resistance to H2O2 with respect to that of the temperate counterpart, as also suggested by a comparative modeling analysis of the catalytic subunits. Moreover, in contrast to that observed in T. bernacchii, the 26S systems from C. hamatus and D. labrax were incapable to hydrolyze oxidized proteins in a ubiquitin-independent manner. Therefore, the 'uncommon' properties displayed by the Antarctic 26S proteasomes can mirror the impact exercised by evolutionary pressure in response to richly oxygenated environments.

Time required to start multicentre clinical trials within the Italian Medicine Agency programme of support for independent research.

BACKGROUND AND AIM: Time allowed for independent ethics committees (IECs) and administrative offices to assess and activate clinical trials is regulated by law. This study aims to describe time spent activating two multicentre non-profit trials supported by the Italian Medicines Agency (AIFA). Five non-AIFA supported (NAS) trials were used as a benchmark. METHODS: The two AIFA-supported trials were FATA-GIM3 (optimal adjuvant hormonal treatment for breast cancer) and TOSCA (duration of adjuvant FOLFOX in colorectal cancer). The five NAS trials focused on lung or ovarian cancer. The following were measured for all trials: date of submission of trial documentation to peripheral IEC, date of IEC opinion and date trial contracts were signed. Times are reported in months. RESULTS: 106 centres applied to participate in FATA-GIM3 and 137 in TOSCA. An IEC opinion was issued by 100/106 (1 negative opinion) and 137/137 (2 negative opinions) centres, with a median time from submission of 3.6 months (range 0.1-60.2). After a positive IEC opinion, the median time before signing the trial contract was 3.3 months (0.1-59.2). Contracts were signed with 93/99 and 135/135 centres, with a median time from submission of study documentation of 8.4 months (0.5-61.1). Times for NAS trials were not substantially different. CONCLUSIONS: FATA-GIM3 and TOSCA centres were opened after a median of 8 months, consisting of nearly 4 months each for IEC opinion and administrative signature, similar to the NAS trials. The process of trial activation in Italy remains inefficient and takes far longer than legally allowed.

The axillary vein puncture for implantable cardiac defibrillator implantation: 14 years of experience. Analysis of the results.

BACKGROUND: Axillary vein puncture (AVP) is a valid alternative to Subclavan vein puncture for leads insertion in cardiac implantable electronic device implantation, that may reduce acute and delayed complications. Very few data are available about ICD recipients. A simplified AVP technique is described. METHODS: All the patients who consecutively underwent "de novo" ICD implantation, from March 2006 to December 2020 at the University of Verona, were considered. Leads insertion was routinely performed through an AVP, according to a simplified technique. Outcome and complications have been retrospectively analyzed. RESULTS: The study population consisted of 1711 consecutive patients. Out of 1711 patients, 38 (2.2%) were excluded because they were implanted with Medtronic Sprint Fidelis lead. Out of 1673 ICD implantations, 963 (57.6%) were ICD plus cardiac resynchronization therapy, 434 (25.9%) were dual-chamber defibrillators, and 276 (16.5%) were single-chamber defibrillators, for a total of 3879 implanted leads. The AVP success rate was 99.4%. Acute complications occurred in 7/1673 (0.42%) patients. Lead failure (LF) occurred in 20/1673 (1.19%) patients. Comparing the group of patients with lead failure with the group without LF, the presence of three leads inside the vein was significantly associated with LF, and the multivariate analysis confirmed three leads in place as an independent predictor of LF. CONCLUSION: AVP, according to our simplified technique, is safe, effective, has a high success rate, and a very low complication rate. The incidence of LF was exceptionally low. The advantages of AVP are maintained over time in a population of ICD recipients.

Geodetic imaging of magma ascent through a bent and twisted dike during the Tajogaite eruption of 2021 (La Palma, Canary Islands).

On Sept. 19th, 2021, the largest historical eruption on the island of La Palma began, which had a significant scientific, social, and economic impact. The 2021 Tajogaite eruption was characterised by short precursors, lasting only 8 days. The seismicity started on Sept. 11th with a westward and upward migration of hypocenters. Permanent GNSS stations started recording deformation on Sept. 12th on the island's western side, which reached more than 15 cm just before the eruption. After the eruption onset, the ground deformation increased, reaching a maximum on Sept. 22nd and showing a nearly steady deflation trend in the following months. To better understand the dynamics of the eruption, we exploited a joint dataset of GNSS and Sentinel-1 SBAS time series along both ascending and descending orbits. To obtain the geometry of the causative source of the ground deformation, we combined the result of a preliminary non-linear inversion and the precise location of the seismicity. The resulting geometry of the source is that of a twisted dike bending eastward. We performed inverse modelling to obtain the spatiotemporal kinematics of the opening function of the dike. The forward modelling has been realised using a 3D finite-element approach considering the island's topography. Our findings reveal a close correspondence between the magmatic intrusion and pre-eruptive seismicity. The ascent of the magma occurred along two branches, and the rheology of a previously identified ductile layer strongly affected the magma propagation process. Finally, we found evidence of an early shallow deformation, which we interpret as the effect of ascending hydrothermal fluids. Our findings highlight the need for advanced modelling to understand pre-eruptive processes in basaltic volcanoes.

Type B Aortic Dissection CTA Collection with True and False Lumen Expert Annotations for the Development of AI-based Algorithms.

Aortic dissections (ADs) are serious conditions of the main artery of the human body, where a tear in the inner layer of the aortic wall leads to the formation of a new blood flow channel, named false lumen. ADs affecting the aorta distally to the left subclavian artery are classified as a Stanford type B aortic dissection (type B AD). This is linked to substantial morbidity and mortality, however, the course of the disease for the individual case is often unpredictable. Computed tomography angiography (CTA) is the gold standard for the diagnosis of type B AD. To advance the tools available for the analysis of CTA scans, we provide a CTA collection of 40 type B AD cases from clinical routine with corresponding expert segmentations of the true and false lumina. Segmented CTA scans might aid clinicians in decision making, especially if it is possible to fully automate the process. Therefore, the data collection is meant to be used to develop, train and test algorithms.

Ambulance Response in Eight Minutes or Less: Are Comorbidities a Factor.

A recommended emergency medical services ambulance response time to a medical emergency is within eight minutes for at least 90% of calls. This study aimed to evaluate scene times for rural education and outreach to improve the quality of trauma care. This is a single-center study of Trauma Registry data from July 1, 2016 to February 28, 2022. The inclusion criteria were based upon age (≥18 years). A logistic regression was performed to identify predictor variables on the likelihood that an adult trauma patient will experience scene times greater than eight minutes. 19 321 patients were included in the analysis; 7233 (37%) experienced an elapsed scene time within eight minutes. This research identified an opportunity to improve rural trauma team response time, which is only reaching 37% of the patient population within eight minutes. Prehospital cardiac arrest and unique pre-existing comorbidities may play a role in extended response times by EMS.

Predicting Early Versus Late In-Hospital Mortality in the Trauma Population.

This study aimed to evaluate non-survivors who were admitted to a level I trauma center but later died, in terms of predicting who would expire early vs late. This is a single-center study of Trauma Registry data, from July 3, 2016, to February 24, 2022. The inclusion criteria were based upon age (≥18 years) and in-hospital mortality. 546 patients (mean age 58) were included in the analysis. Trauma patients who may experience an earlier death were those with increasing injury severity scores, activation of massive transfusion protocol, comorbid advanced directive limiting care, COPD, personality disorder, and ED death location. Patients were more likely to experience later in-hospital mortality, including those with increasing ICU stays, and comorbid dementia.

In-Hospital Mortality after Massive Transfusion and Motorcycle Trauma may Predict Highest Injury Severity in a Rural Level I Trauma System.

The leading cause of preventable traumatic death is uncontrolled bleeding. This study aimed to better identify those most likely to experience in-hospital mortality with increasing injury severity scoring (ISS). This is a single-center study of Trauma Registry data, from July 3, 2016, to February 24, 2022. The inclusion criteria were based upon age (≥18 years) and in-hospital mortality. 546 patients (mean age 58) were included in the analysis. There were several significant associations with increasing ISS among those who experienced in-hospital mortality, which included a rising shock index ratio, activation of the massive transfusion protocol, and, most notably, motorcycle trauma. This research reiterates the importance of the "Stop the Bleed" campaign as vital for training laypersons in the life-saving technique for hemorrhage control.

In-Hospital Mortality After Massive Transfusion and Motorcycle Trauma May Predict Highest Injury Severity in a Rural Level I Trauma System.

The leading cause of preventable traumatic death is uncontrolled bleeding. This study aimed to better identify those most likely to experience in-hospital mortality with increasing injury severity scoring (ISS). This is a single-center study of Trauma Registry data, from July 3, 2016, to February 24, 2022. The inclusion criteria were based upon age (≥18 years) and in-hospital mortality. 546 patients (mean age 58) were included in the analysis. There were several significant associations with increasing ISS among those who experienced in-hospital mortality, which included a rising shock index ratio, activation of the massive transfusion protocol, and, most notably, motorcycle trauma. This research reiterates the importance of the "Stop the Bleed" campaign as vital for training laypersons in the life-saving technique for hemorrhage control.

Can Pre-Hospital Medical Management Predict In-Hospital Mortality in Trauma?

The current literature demonstrates an association between both size and presence of TBI and its effects on mortality; however, it does not readily address the morbidity and associated functional outcomes of those who survive. We hypothesize that the likelihood of discharge to home decreases with advancement of age in the presence of TBI. This is a single-center study of trauma registry data, inclusive years July 1, 2016, to October 31, 2021. The inclusion criteria was based upon age (≥40 years), and ICD10 diagnosis of a TBI. Disposition to home without services was the dependent variable. 2031 patients were included in the analysis. We hypothesized correctly that the likelihood of discharge to home decreases (by 6%) with advancement of age (per year) in the presence of intracranial hemorrhage.

Adult Traumatic Brain Injury and Likelihood of Routine Discharge: Do Comorbidities Matter?

The current literature demonstrates an association between both size and presence of traumatic brain injury (TBI) and its effects on mortality, however it does not readily address the morbidity and associated functional outcomes of those who survive. We hypothesize that the likelihood of discharge to home decreases with advancement of age in the presence of TBI. This is a single-center study of Trauma Registry data, inclusive years July 1, 2016 to October 31, 2021. The inclusion criteria was based upon age (≥40 years), and ICD10 diagnosis of a TBI. Disposition to home without services was the dependent variable. 2031 patients were included in the analysis. We hypothesized correctly that the likelihood of discharge to home decreases (by 6%) with advancement of age (per year) in the presence of intracranial hemorrhage.

Traumatic Brain Injury Under Triage Risk: A Rural Trauma System Experience.

BACKGROUND: The process of interfacility transfer may cause a delay in the necessary medical treatment, which could lead to poor outcomes and increased mortality rates. The ACS-COT considers an acceptable under triage rate of <5%. The aim of this research was to identify the likelihood of under triage among transferred-in traumatic brain injury (TBI) patients. METHODS: This is a single-center study of Trauma Registry data, from July 1, 2016, to October 31, 2021. The inclusion criteria were based upon age (≥40 years), ICD10 diagnosis of TBI, and interfacility transfer. Under triage using the Cribari matrix method was the dependent variable. A logistic regression was performed to identify additional predictor variables on the likelihood that an adult TBI trauma patient experienced under triage. RESULTS: 878 patients were included in the analysis; 168 (19%) experienced an under triage. The logistic regression model was statistically significant (N = 837, P < .01). In addition, several significant increases in odds for under triage were identified, which included increasing injury severity score (ISS; OR 1.40, P < .01), increasing AIS head region (OR 6.19, P < .01), and personality disorders (OR 3.61, P = .02). In addition, a reduction in odds in TBI adult trauma under triage is the comorbidity of anticoagulant therapy (OR .25, P < .01). CONCLUSIONS: The likelihood of under triage in the adult TBI trauma population is associated with increasing AIS head injuries and increasing ISS and among those with mental health comorbidities. This evidence and additional protective factors, such as patients on anticoagulant therapy, may aid in education and outreach efforts to reduce under triage among the regional referring centers.

MADIA - Meteorological variables for agriculture: A dataset for the Italian area.

The MADIA gridded dataset provides the dekadal series of the main agro-meteorological variables derived from ERA5 hourly surface data, across Italy for the period 1981-2021, and their respective 1981-2010 and 1991-2020 climate normals, also including absolute minimum and maximum and the main quantiles. Temporal and spatial resolutions are 10-day and 0.25 degrees respectively and the dataset is annotated with standard metadata. The dataset was obtained by: (1) estimating the daily time series of minimum, average and maximum air temperature, minimum and maximum air relative humidity, wind speed, surface solar radiation downwards, precipitation and reference evapotranspiration according to the FAO Penman-Monteith method; (2) summarising them to 10-day series as accumulated values for precipitation and evapotranspiration and mean values for the other variables. The MADIA dataset is provided in both NetCDF and csv format. A complementary vector file is provided which reports for every cell the fractions covered of the total area of each administrative unit considered to derive statistics for Italy on the European Nomenclature of Territorial Units for Statistics levels (NUTS 2 and 3). Other potential dataset reuses are the estimation of bioclimatic indices and statistical downscaling of climate scenarios.