We present a Virtual Research Environment (VRE), developed in the form of a Jupyter Notebook, devoted to a multidisciplinary and multiparametric analysis of an earthquake. In particular, the VRE retrieves and analyses lithospheric, atmospheric and ionospheric parameters from various sources. Among them, we selected the Open Data of the European Plate Observing System (EPOS) for the earthquake catalogue, integrated with the atmospheric climatological archive MERRA-2 provided by NASA, and the satellite magnetic data of the Swarm mission produced by ESA. The open-source code and friendly environment of Jupyter Notebook allow future users to personalise the research parameters (e.g., investigated time, area of interest) and apply the same VRE to other earthquakes. The VRE is designed for the 2016 seismic sequence that affected Central Italy, characterised by the Mw = 6.0 earthquake near Accumoli and Amatrice on August 24, 2016, and the Mw = 6.5 earthquake near Norcia on October 30, 2016. The Jupyter Notebook performs: 1) seismological investigations such as the calculus of cumulated Benioff strain and released energy weighted for the distance of the target event; 2) extraction of typical atmospheric values and identification of possible anomalies and 3) searching for small oscillations of magnetic field not likely related to geomagnetic activity. Finally, a graphical comparison of all the investigations is provided to search for possible interactions among the different geo-layers. In the context of studying the preparation phase of the 2016–2017 Italian Seismic Sequence, we provide a tool for exploiting and integrating multiparametric EPOS data with other datasets.