Space weather, geomagnetic, and ionospheric effects are major sources of Global Navigation Satellite System’s (GNSS) Positioning, Navigation, and Timing (PNT) service degradation. Researchers studying the problem often face an obstacle in missing large data sets comprising ionospheric effects variables and space weather ones for their studies. Here a formal assemblage method is proposed for data set assemblage for GNSS ionospheric effects assessment that utilises experimental sets of space weather/geomagnetic observations, and raw GNSS pseudorange observations polluted with uncorrected ionospheric delay effects to estimate Total Electron Content (TEC) as a target variable outlining the ionospheric effects on GNSS. The five Observations Quality Principles (OQP) are proposed for the observations to be used in statistical learning model development to comply to. The proposed assemblage method is demonstrated in a particular case study of an assembled data set for Darwin, NT, Australia, in sub-equatorial zone, using a tailored software developed in the R environment for scientific computing, and a freely available TEC estimation software. The assembled massive data obtained in the demonstration of the proposed method is provided in the open-access manner in support of the international scientific community.