With the rapid development of the advanced ground-based Global Navigation Satellite Systems (GNSS) tropospheric sounding technique in recent decades, it has offered unprecedented opportunities for atmospheric information acquisition. The variables obtained using this technique, e.g., zenith total delay (ZTD) and tropospheric gradients, are recognized to possess the benefits of high spatiotemporal resolution, long-term stability and all-weather capability. Therefore, in recent years, they have demonstrated their feasibility and great potential in serving the meteorology, hydrology and climate communities. In this data archive, GNSS tropospheric products at 40 ground-based GNSS stations employed in New Zealand over the period January 1-31, 2022 are provided. The whole dataset includes the site coordinates, 5-min ZTD and tropospheric gradients. Spefically, the Bernese GNSS software with the precise point positioning approach was adopted to retrieve ZTD estimates at these stations. In addition, the Vienna Mapping Function 1, the observation window of 27 h, the elevation cut-off angle of 3°, the high quality International GNSS Service Final orbits and clocks were also included in the retrieval process. To further improve the quality of the dataset, during the data processing, the phase wind-up effect, the solid Earth and ocean pole tide effects as well as relativistic delays were all modeled according to the IERS Conventions 2010.

1. All the data files can be viewed through a "Text document reader";2. Naming rules for all the data files are "SITE/ID(4 digits)+Day of Year(3 digits)+0.22zpd", e.g., "auck0010.22zpd". This name indicates it includes the dataset retrieved at the AUCK station on January 1, 2022.3. In the "SITE/ID" part, the latitude, longitude and height for each one of the 40 stations are listed.4. In the "TROP/SOLUTION" part, TROTOT denotes ZTD value; TGNTOT and TGETOT represent the north-south and east-west gradient parameters; the others refer to their respective standard deviations, denoted by "STDDEV". 5. The unit for the "EPOCH" is second. Therefore, the sampling interval of 300 represents the temporal resolution for the whole dataset is 300s, i.e., 5min.