\section{Introduction}
Millisecond pulsars (MSPs) are a subclass of pulsars that are thought to have experienced a unique evolutionary path, involving the accretion of material from a binary companion. The subsequent transfer of angular momentum from this process results in a rotational frequency up to a factor of a thousand greater than would be expected of an ordinary pulsar,
a reduction in their magnetic field strength down to $\sim 10^8$ G, and makes them far more rotationally stable. Observations of individual members of this subclass have proven incredibly fruitful leading to extensive tests of general relativity \citep{2012MNRAS.423.3328F, 2018Natur.559...73A, 2020A&A...638A..24V}, the discovery of planetary systems hosted by distant stars \citep{Wolszczan_Frail_1992, Wolszczan_planets_1994}, and allowing us to place ever-tightening constraints on nuclear equations of state \citep{Demorest_Nature_2010, Antoniadis1233232, 2021ApJ...915L..12F, 2021arXiv210506979M, 2021arXiv210506980R}. The  precision and stability of these pulsars enables the realisation of sensitive pulsar timing arrays (PTAs), ensembles of pulsars that are consistently observed, and their pulse times of arrival (ToAs) catalogued and measured for up to decades \citep{1978SvA....22...36S, 1979ApJ...234.1100D, 1990ApJ...361..300F}. The primary goal of PTA observations is the detection of nanohertz-frequency gravitational radiation \citep{1990ApJ...361..300F}. The most likely first source of nanohertz-frequency gravitational waves (GWs) is a stochastic gravitational-wave background (SGWB) that is the superposition of many sources \citep{Hellings_Downs_1983}. Several recent searches show potential evidence for a common spectral noise process that could indicate the presence of such a background \cite[][]{NanoGravGWB, 2021arXiv210712112G, 2021MNRAS.508.4970C, 2022MNRAS.510.4873A}, although it remains to be seen if this common red noise process would follow the expected Hellings-Downs (HD) angular correlation. The most likely progenitors of the GW signals in this frequency-band are thought to be inspiralling supermassive black-hole binaries (SMBHBs) \citep{1995ApJ...446..543R, 2003ApJ...583..616J, 2003ApJ...590..691W, 2004ApJ...611..623S, 2011MNRAS.411.1467K, 2012A&A...545A.127R, 2017PhRvL.118r1102T}. However, there may exist additional sources from which a SGWB may arise, such as from cosmic strings \citep{1976JPhA....9.1387K, PhysRevD.81.104028, 2012PhRvD..85l2003S, 2015MNRAS.453.2576L, 2018ApJ...859...47A}, or cosmological phase transitions \citep{1980PhLB...91...99S, 2005PhyU...48.1235G}.

The MeerKAT radio telescope\footnote{Operated by the South African Radio Astronomy Observatory (SARAO).} is a $64$-antenna interferometer located in the Northern Karoo region of South Africa, and is the precursor to the upcoming Square Kilometer Array Mid telescope \citep[SKA;][]{2009IEEEP..97.1482D} demonstrating SKA technology on a SKA site. Millisecond pulsar timing with MeerKAT (The MeerKAT Pulsar Timing Array\footnote{Sometimes referred to as the MeerTime Pulsar Timing Array.}, MPTA), is undertaken as part of the MeerTime Large Survey Project (LSP) \citep{2016mks..confE..11B}. The potential of the MeerKAT PTA project has already been demonstrated through dedicated studies of jitter noise in the MSP population \citep{2021MNRAS.502..407P}, a comprehensive census of the MSPs visible to MeerKAT and their emission properties \citep{spiewak2022}, and the discovery of mode-changing in one of the most precisely timed and often observed millisecond pulsars \citep{2022MNRAS.510.5908M}.

The MPTA is poised to make a significant contribution to the International Pulsar Timing Array \cite[IPTA;][]{2010CQGra..27h4013H}.
The IPTA is currently a collaboration between four individual PTAs (chronologically, the European PTA \citep[EPTA;][]{2008AIPC..983..633J}, the Parkes PTA \citep[PPTA;][]{2008AIPC..983..584M}, the North American Nanohertz Observatory for Gravitational Waves \citep[NANOGrav;][]{2009arXiv0909.1058J}, and the Indian PTA \citep[InPTA;][]{2022arXiv220609289T}) to search for low-frequency gravitational waves. Of the major PTAs that are members of this effort, the PPTA has so far been the only one capable of observing the entire Southern sky. The MPTA can therefore make a key contribution to the observation of Southern MSPs.
MeerKAT has shown itself to be many times more sensitive than the Parkes ({\em Murriyang}) radio telescope, possessing an array gain that is $\sim 4$ times greater, combined with a lower system temperature further bolstering its sensitivity.
In addition to this high sensitivity MeerKAT's ability to quickly slew between sources and the access it has to a large swathe of pulsars (Dec $<+44^{\circ}$), makes the MPTA naturally suited to contribute to the IPTA through high precision timing of a large array of pulsars. Forecasts show that the MeerKAT PTA will be contributing roughly half of the IPTA sensitivity by 2025 with current target lists and cadences \cite[][]{spiewak2022}.

In this paper, we present the initial 2.5-year MPTA data release. These data include both sub-banded and frequency averaged ToAs, the full data archives used to construct this data release, the ephemerides that have been used to perform timing, and band-averaged and sub-banded timing residuals. Also included are both the timing standards (frequency averaged, one dimensional) and frequency resolved standards (portraits), used to calculate each set of ToAs as well as the preliminary noise models constructed to facilitate this work. This is the first such data release from the MPTA suitable for obtaining arrival times. In Section \ref{Section: pulsars, obs, data reduction} we describe the pulsars, observations taken, and the process of data reduction. In Section \ref{Section: data set / data products} we present the final data products. In Section \ref{Section: Discussion} we discuss the performance of the data and explore various features and investigations that are contained in this data release. Finally, in Section \ref{Section: Summary and Conclusions} we conclude and summarise our findings.

\section{Pulsar Selection, Observations, and Data Reduction}
\label{Section: pulsars, obs, data reduction}

This data release comprises observations taken exclusively at the MeerKAT radio telescope, an interferometer consisting of $64$, $13.9$-m dishes, located in the Great Karoo region of South Africa. While MeerKAT currently possesses L-band ($856 - 1712$ MHz) and Ultra-High Frequency UHF ($544 - 1088$ MHz) receivers, all data collected for this project used the L-band system, following its successful commissioning. We have not collected any data with the UHF receiver which was commissioned later. In addition, we expect observations in L-band to be less impacted by noise introduced by interstellar medium (ISM) propagation effects, even though UHF observations potentially provide more accurate measurements of dispersion measure \citep{2010arXiv1010.3785C}. MeerKAT possesses a gain of $2.8$ K\,$\mathrm{Jy}^{-1}$, and a system temperature of $\sim 18$K in this band, which results in a system equivalent flux density of $\sim 7$\,Jy. The data were recorded with the SKA prototype pulsar processors, the pulsar timing user-supplied equipment, collectively known as PTUSE \citep{2020PASA...37...28B}.

\subsection{Selection Criteria and Observing Strategy}
\label{Subsection: selection crit and obs strat}

The data presented here are taken from shortly after the MeerTime project began, spanning from February 2019 through to October 2021 (MJD $58526 - 59510$). This provides approximately $2.5$ years of data for the majority of the pulsars included in the MPTA. The full MPTA includes $88$ pulsars, which have been selected from a census of $189$ MSPs visible to MeerKAT \citep{spiewak2022}. The pulsars were initially chosen to obtain a band averaged timing precision of $<1 \mu\mathrm{s}$ within an integration time of $1000$ seconds. Since the initial selection of the MPTA pulsars, the list has been regularly reviewed, and pulsars have been added and removed from regular timing to maximise the efficiency of the PTA. Similarly, some pulsars have not been included in the MPTA, despite possessing high timing precision, as they are only observable at  declinations that are not ideal for the MeerKAT telescope, so better observed with Northern hemisphere telescopes. The  strategy for the pulsars that are currently included in the MPTA is to regularly observe with a minimum cadence of two weeks per pulsar. There is a minimum integration time of $256$ seconds per observation, if the pulsar is able to achieve sub-microsecond timing precision in that time. Where this cannot be achieved, the pulsar is observed for an integration time required for sub-microsecond precision, if this can be achieved within $2048$ s. While the current MPTA includes $88$ pulsars, this data release includes only the $78$ pulsars that had at least $30$ observations over this observing span. The positions of the pulsars are presented in Figure \ref{fig: skymap}. All of the MSPs experience interstellar scintillation that amplifies and deamplifies their flux densities. The effects of this are more dramatic at low dispersion measures and means that sometimes the pulsars do not reach their target timing residual in their allocated integration times.

\begin{figure*}

    \includegraphics[width=\textwidth]{skymap_wrms_RA_DEC_FIX.pdf}
    \caption[]{MPTA pulsars included in this data release. The color of each pulsar refers to the weighted root-mean-square error of its band-averaged residuals. As most pulsar errors are sub-microsecond, the color is displayed in a log-scale to more easily contrast them. The size of the pulsars refers to their relative median observation times (e.g. larger circles indicate pulsars with typically longer observations).}
    \label{fig: skymap}
\end{figure*}

While most of the observations follow this two week cadence, there are some exceptions. Notably, $38$ of the binary pulsars in the MPTA have been observed by the MeerTime relativistic binary programme \citep{2021MNRAS.504.2094K} which includes longer observations and single orbit observing campaigns to precisely measure Shapiro delays, enabling companion and pulsar masses to be derived along by determining the orbital inclination angle.

A summary of ten of the pulsars that are included in this data release is provided in Table \ref{table: MPTA_summary}, with the full list included in Table \ref{table: MPTA_summary_full_table}. Included in this table are the pulse period ($P$), the dispersion measure (DM), the orbital period where applicable ($\textrm{P}_{\textrm{b}}$), and the median, mean, and standard deviation of the sub-banded ToA uncertainties. Also included is the frequency-averaged effective pulse width ($\mathrm{W_{eff}}$) \citep{2010arXiv1010.3785C, 2021MNRAS.502..407P}, which we discuss in Section \ref{Subsection: Data set performance}. This is presented as it is expected to provide a better prediction of timing precision than more conventional width measurements.

\begin{table*}

\caption{MPTA partial parameters and timing statistics for the first 10 pulsars in order of Right Ascension. The full description of these parameters can be found at Table \ref{table: MPTA_summary_full_table}.}
\label{table: MPTA_summary}
\begin{tabular*}{\textwidth}[c]{@{\extracolsep{\fill}} lcccccccc}
\toprule
Pulsar &
  \multicolumn{1}{c}{P} &
  \multicolumn{1}{c}{DM} &
  \multicolumn{1}{c}{$\textrm{P}_{\textrm{b}}$} &
  \multicolumn{1}{c}{$\textrm{W}_{\textrm{eff}}$} &
  \multicolumn{3}{c}{Uncertainty ($\mu$s)} &
  \multicolumn{1}{c}{Span} \\
  \cmidrule(lr){6-8}
  &
  \multicolumn{1}{c}{(ms)} &
  \multicolumn{1}{c}{($\textrm{cm}^{-3}$pc)} &
  \multicolumn{1}{c}{(d)} &
  \multicolumn{1}{c}{($\mu$s)} &
  med &
  $\mu$ &
  $\sigma$ &
  \multicolumn{1}{c}{(yr)} \\ \midrule

J0030+0451 & 4.87 & 4.33 & - & 62.38 & 3.11 & 3.54 & 2.18 & 1.80 \\
J0125-2327 & 3.68 & 9.59 & 7.28 & 66.04 & 0.71 & 0.99 & 0.95 & 2.43 \\
J0437-4715 & 5.76 & 2.65 & 5.74 & 87.51 & 0.13 & 0.13 & 0.01 & 2.50 \\
J0610-2100 & 3.86 & 60.69 & 0.29 & 87.86 & 2.77 & 3.04 & 1.61 & 2.43 \\
J0613-0200 & 3.06 & 38.79 & 1.20 & 38.09 & 1.01 & 1.25 & 0.83 & 2.43 \\
J0614-3329 & 3.15 & 37.05 & 53.58 & 31.29 & 2.43 & 3.15 & 2.92 & 2.43 \\
J0636-3044 & 3.95 & 15.46 & - & 106.92 & 3.91 & 5.32 & 4.26 & 2.37 \\
J0711-6830 & 5.49 & 18.41 & - & 79.99 & 3.39 & 4.67 & 4.15 & 2.44 \\
J0900-3144 & 11.11 & 75.69 & 18.74 & 224.47 & 2.17 & 2.36 & 0.88 & 2.44 \\
J0931-1902 & 4.64 & 41.49 & - & 57.20 & 3.37 & 4.17 & 2.97 & 2.34 \\
...        & ...   & ...   & ...   & ...  & ... & ... & ... & ... \\ \bottomrule
\end{tabular*}%
\end{table*}%

\subsection{Data Reduction and RFI mitigation}
\label{Subsection: data_reduction}

All MPTA observations were recorded using the PTUSE backend \citep{2020PASA...37...28B} and then transferred to both the SARAO data archive and the MeerTime data archive hosted on the OzStar supercomputer. The files were produced in \textsc{psrfits} \citep{2004PASA...21..302H} archive format, containing $8$-second sub-integrations, a time resolution of $1024$ bins per folding period, $4$ polarisation products (Stokes I, Q, U $\&$ V), and $1024$ frequency channels dividing the L-band digitized bandwidth ($856 - 1712$ MHz). All of the data are coherently dedispersed at the nominal dispersion measure of the
pulsar and folded at the topocentric period of the pulsar.
Computational limitations with the initial PTUSE machines meant we were not able to record the entire band for earlier observations, instead, we recorded the inner $928$ channels with a total bandwidth of $775.75$ MHz. This has since been remediated, however later data was reduced to keep a consistent band. The system temperature increases in the outer channels so is a relatively small loss in total sensitivity.
Changes in delays and time tagging within the hardware have also introduced known offsets in the pulsar timing data sets. These time offsets have been directly corrected for in this release. The epoch and size of the offsets (often referred to as `jumps') are listed in Table \ref{table: jump_table}.

The custom processing pipeline, \textsc{meerpipe}, was used to further process the archives to a desirable format, and to perform required radio frequency interference (RFI) excision. \textsc{meerpipe} calibrates the full resolution data products in both flux and polarisation. Prior to the $9^\mathrm{th}$ of April 2020, polarisation calibration was performed through the use of the \textsc{psrchive} \footnote{http://psrchive.sourceforge.net/} \citep{2004PASA...21..302H} \texttt{pac} utility. Any data produced after this date are calibrated on-site by the SARAO Science Data Processing pipeline using calibration solutions derived from phase up solutions \citep{2021MNRAS.505.4483S}, prior to its folding in the PTUSE machines. A right-handed circular polarisation convention was adopted, and the pipeline converted the polarisation parameters from the native polarisation (a description of polarisation through the intensity of the two linear components and their cross terms), to the more familiar Stokes parameters. Following this, RFI excision was performed via \textsc{meerguard}, a modified version of the RFI-excision program \textsc{coastguard} \citep{2016MNRAS.458..868L}, which has been optimized for the RFI expected to be present at the MeerKAT telescope location. This software operates by comparing observations of a pulsar to a high-quality, frequency-averaged template profile. The residuals created from subtracting the observation from the profile are then used to identify channels where RFI is present, and the channels are subsequently zero-weighted in the archive. Lastly, the data were decimated to a more suitable format for use in this data release.
The format of these files were frequency averaged from the original $928$ frequency channels, to $16$ sub-bands across the usable frequency range, and were fully time integrated over the entire observation length. These operations are both performed through the \texttt{pam} program contained within \textsc{psrchive}. The nature of the MPTA observing strategy results in most data output by \textsc{meerpipe} to have an integrated length of $\sim 256$ seconds.

\section{The MPTA Data Set}
\label{Section: data set / data products}

The MPTA data release is comprised of the full data archives, fully integrated pulse profiles, the portraits used for the creation of ToAs, the timing ephemerides used in this work, as well as sub-banded and band-averaged ToAs, for all pulsars listed in Table \ref{table: MPTA_summary_full_table}.

\begin{figure*}

    \includegraphics[width=\textwidth]{J0900-3144_grand.dly_4_notebook_comparison.pdf}
    \vspace{-0.5cm}

    \includegraphics[width=\textwidth]{J1909-3744_grand.dly_4_notebook_comparison.pdf}
    \caption[]{Comparison of MPTA portraits built from the full MPTA 2.5 year dataset, via the process described in Section \ref{Subsection: Portrait_creation}, with time-averaged observational data for J0900$-$3144 (Top) and J1909$-$3744 (Bottom). (Left) Direct comparison of \textsc{PulsePortraiture} models (dashed blue) against the observational data (solid red). (Right) The profile residuals from subtracting the model from the observational data. Above each residual, the reduced chi-squared statistic is supplied for both the total pulse profile, and for the off-pulse region. The portraits match the data well enough that it is difficult to discern the differences between the model and the data without the assistance of the residual plot.}
    \label{fig: portrait_example}
\end{figure*}

\subsection{Arrival times}
\label{Subsection: arrival_times}

To produce the arrival times we used the following procedure in the spirit of creating a reproduceable data set.

\begin{enumerate}
    \item From the $16$ sub-banded data output by \textsc{meerpipe}, identify and remove any observations with issues known to be unrecoverable. For clarity, an observation was defined as unrecoverable where there existed phase smearing of $>1$ bin in either the time or frequency space, or where the observation was taken with an incorrect position beyond the tied beam.
    \item Compensate for the system induced, \textit{a priori} known time delays, as per Table \ref{table: jump_table}.
    \item Form initial ToAs from a standard \textsc{psrchive} \texttt{pat} utility, using a Fourier domain monte-carlo algorithm (FDM), and output a tempo2 IPTA format.
    \item Filter for sub-banded observations with signal-to-noise ratio (S/N) $> 10$.
    \item Use \textsc{tempo2} \citep{2006MNRAS.369..655H} to fit for pulsar spin frequency and the frequency time derivative, as well as dispersion measure and its first time derivative. For some pulsars, these initial ephemerides were available from existing PTA data releases \cite[e.g.,][]{2021MNRAS.507.2137R}. For others, we started with ephemerides from the ATNF pulsar catalogue \citep{2005AJ....129.1993M}. Some parameters were of insufficient precision, resulting in significant phase drifts in some of the pulsar timing residuals. To account for this, additional timing parameters were selectively fit for in five pulsars. These parameters included right ascension, declination, orbital parameters, and proper motion.
    \item Inspect outliers and record any that were removed, including a justification for exclusion.
    \item Create an intermediate ephemeris and timing files from \textsc{tempo2}.
    \item Update the data with the ephemeris, and create a corresponding new standard template.
    \item Create a new set of ToAs excluding the outliers. The ToAs were produced as in step (iii).
\end{enumerate}
Following these steps, the pulse profiles were further used to create portraits as described below, and sub-banded ToAs were subsequently calculated using the FDM algorithm and the filtered and updated data. The ephemerides included in this data release are preliminary as we have not fully explored the timing models or undertaken a comprehensive noise analysis. Nevertheless, they are suitable for gravitational wave and other parameter searches.

\subsection{Portrait Creation}
\label{Subsection: Portrait_creation}

Profile evolution, where pulsar-emission characteristics change with frequency, results in the need for frequency dependent average pulse profiles.  This is especially true with the high-sensitivity wide-band observations taken with MeerKAT. We detected profile evolution (of varying degree) in all  pulsars presented in the data set. The portraits used for timing were made with the \textsc{PulsePortraiture} software \citep{2019ApJ...871...34P}, which creates noise-free, high-quality, frequency-dependent templates, that can be used for precise sub-band timing. Models for two pulsars that have been created this way are shown in Figure \ref{fig: portrait_example}.

Producing the pulse portraits was an iterative process. To form the pulse portraits, we first identify high S/N observations. The S/N threshold was varied between pulsars, as each pulsar possessed different median flux values, and many showed frequency dependent scintillation. The full frequency resolution ($928$ channel) archives from these observations were then combined using the \texttt{psradd} utility of \textsc{psrchive} using the pulsar ephemeris to align the observations. \textsc{PulsePortraiture} was used to decompose the profile using principal-component analysis (PCA), with the components allowed to vary in frequency independently.
We found that the process was sensitive to low levels of latent RFI in the data. This is because, in our use of \textsc{PulsePortraiture}, we normalize the data per frequency channel such that the profile evolution becomes the dominant source of variation as opposed to changes in intensity\footnote{We note there are multiple ways to scale the data in \textsc{PulsePortraiture}. However, this method is standard practice in pulsar timing \citep{2019ApJ...871...34P}.}. This normalisation amplified signals in nearly zero-weighted channels affected by RFI. To mitigate this we used \texttt{paz} to median filter, both before and after normalisation. Any remaining channels were manually excised using the \textsc{psrchive} program, \texttt{pazi}. Approximately $150$ channels were removed during this process, which is consistent with the channels that are affected by the RFI as reported in \citet{2020PASA...37...28B}. This was necessary as \textsc{PulsePortraiture} was sensitive to the presence of RFI. Similarly, while the construction of the frequency-dependent model returns a high level of agreement, demonstrated through the small uncertainty distributions in Figure \ref{fig: toa_unc_w_subbands}, elements of structure were still seen to exist in the profile residuals (the difference between the portraits and the data used to produce them), particularly in the lowest and highest sub-bands.
One possible reason for this may be that there exists latent RFI in the archives, that is not significant enough to have been automatically or manually removed, but may influence the components decomposed from the PCA. Alternatively, it could be that there exist subtle variations in the average pulse that were not well modelled through \textsc{PulsePortraiture}. Due to this, in the process of creating the timing residuals, frequency-dependant (FD) parameters \citep{2015ApJ...813...65N} were added to the majority of the pulsar parameter files (discussed in Section \ref{Subsection: Noise modelling}), with the aim to compensate for unmodelled frequency dependent variations. In addition, we note that in the creation of a frequency-resolved portrait, a fiducial DM must be chosen. This could result in differing frequency-resolved portraits, for different choices of DM. When combining MeerKAT data with other IPTA data sets, this may result in both PTA dependent FD parameters and DMs. Alternatively, one could attempt to re-time all the data sets with a common portrait.

\begin{figure}

    \includegraphics[width=\columnwidth]{Uncertainty_inc_subband_density.pdf}
    \caption[]{Arrival time distributions for the MPTA pulsars. The ToA uncertainties for all $78$ pulsars are included, the top panel showing the distribution of the uncertainties of ToAs from all sub-bands. The subsequent 8 panels show the density distribution of the ToAs that make up the total ToA uncertainties, separated by frequency, with the corresponding central frequency of each sub-band labelled. The largest $1\%$ of ToA uncertainties have been excluded from this plot for clarity, and can be treated as non-representative outliers.}
    \label{fig: toa_unc_w_subbands}
\end{figure}

\subsection{Noise Modelling and interpretation}
\label{Subsection: Noise modelling}

Our noise modelling followed approaches taken by other PTAs \citep{2017ApJ...834...35L, 2021MNRAS.502..478G, 2022MNRAS.509.5538C}. We modelled the influence of systematic noise that is uncorrelated through time with the noise parameters EFAC ($\mathrm{E}_\mathrm{{F}}$) and EQUAD ($\mathrm{E}_\mathrm{{Q}}$). $\mathrm{E}_\mathrm{{F}}$ is a parameter used directly as a scale factor of the timing uncertainties associated with the pulsar, and $\mathrm{E}_\mathrm{{Q}}$ influences the timing uncertainties as an addition in quadrature, such that the final timing uncertainties follow the convention \citep{2014MNRAS.437.3004L},

\begin{equation}
\label{eq: white_noise_BM}
    \sigma = \sqrt{\mathrm{E}^{2}_\mathrm{Q} + \mathrm{E}^{2}_\mathrm{F} \times \sigma_\mathrm{ToA}},
\end{equation}

\noindent where $\sigma_\mathrm{ToA}$ are the uncertainties reported by the time-tagging algorithm. We note that the convention of Equation \ref{eq: white_noise_BM} is  different to that used in \textsc{tempo2}, but consistent with that used in \textsc{temponest} \cite[][]{2014MNRAS.437...21M}. For the purposes of our model selection, we refer to this as our \textit{basic model}, or $\mathrm{M}_\mathrm{1}$ noise. In addition to this, another white noise parameter, ECORR ($\mathrm{E}_\mathrm{{C}}$), is included \citep{2015ApJ...813...65N}. This parameter accounts for noise that is correlated between sub-banded ToAs, but uncorrelated between observations. While statistically motivated, this component could account for stochastic pulse-pulse changes in pulse morphology, colloquially referred to as pulse jitter \citep{2014MNRAS.443.1463S, 2019ApJ...872..193L, 2021MNRAS.502..407P}.

In addition to the three white-noise terms, we also searched for evidence of time-correlated noise processes that could be present in the MPTA pulsars. In particular, we modelled frequency-dependent DM variations, and an achromatic red-noise process. Both were modelled as stationary, stochastic processes with temporal variations of the form

\begin{equation}
\label{eq: power_spec_red}
    \mathrm{S}(f) \propto \frac{\mathrm{A}^{2}}{12 \pi^{2}} \left( \frac{f}{f_{c}} \right)^{-\gamma},
\end{equation}

\noindent where A is the signal amplitude, $\gamma$ is the spectral index of the power spectrum, and $f_{\mathrm{c}} =1\,\mathrm{yr}^{-1}$ is the reference frequency.

Prior to determining noise models for each pulsar, a rudimentary significance test was performed on whether frequency-dependent parameters (FD1/FD2), and dispersion measure derivatives (DM1/DM2) were required in the parameter files. These parameters were determined to be significant, and therefore required to be fitted for, if it was found that the value of the parameter was $>3 \sigma$.

To search for the  possible noise sources in each pulsar, we considered eight models created from different possible combinations of the noise processes described here, shown in Table \ref{table: noise_models}. We used the \textsc{temponest} \citep{2014MNRAS.437.3004L} software to perform a Bayesian analysis to appropriately model the noise, whilst also analytically marginalising over the timing model \cite[see Appendix A of ][]{2009MNRAS.395.1005V}). The preferred model was selected by comparing the relative evidences for each model using the following procedure. We first considered the model with the largest evidence. If it was found to have a log-Bayes factor, $\rm{\ln}\mathcal{B}$, of $\rm{\ln}\mathcal{B} \geq 4$ over all other competing models, we selected it as the preferred model. If this was not the case, a model that possessed fewer noise terms that fulfilled the requirement of $\rm{\ln}\mathcal{B} \geq 4$ over other competing models was chosen. Where there were no models that satisfied this requirement, the basic model ($\mathrm{M}_{1}$) was chosen.

Following this model selection process, $31$ of the pulsars present in this data set show evidence for time-correlated DM noise. Significant DM variations in other MSP data sets have been shown \citep{2007MNRAS.378..493Y} to occur on observation spans of equivalent lengths to the data we present here. As such, the presence of time-correlated DM variations is likely to be found at the $\sim 2.5$ yr scale of this release, especially given the precision of the data provided. Generally, the turbulent ISM that is thought to give rise to DM variations in timing data is modelled through a Kolmogorov spectrum \citep{1990ARA&A..28..561R, 2021MNRAS.502..478G} and, assuming this, we would expect a $\gamma_\mathrm{DM} \sim 8/3$. This is relatively inconsistent with our findings, with the majority ($20$) of the pulsars that show evidence for DM noise having spectral indices that are less steep (Table \ref{table: MPTA_results_full}). However, in direct comparison to previous noise studies, e.g. \citet{2021MNRAS.502..478G}, we found that DM noise in many of the MPTA pulsars also found in the PPTA are comparable to within reported uncertainties.

We found evidence for achromatic red noise in $4$ pulsars (Table \ref{table: MPTA_results_full}). The presence of this noise source is of note as it can significantly impact the sensitivity of the pulsar to gravitational waves, particularly if the noise has similar spectral indices to what is expected of a gravitational wave background. Aside from the SGWB, there are a number of possibe sources of red noise, including pulsar intrinsic noise such as interactions between the superfluidic centre of a neutron star and its crust \citep{2014MNRAS.437...21M}, magnetospheric state switching \citep{2010Sci...329..408L}, or small bodies orbiting the pulsar \citep{2013ApJ...766....5S}.
Additionally we note that the values of $\mathrm{E}_\mathrm{{C}}$ obtained both in this data release and elsewhere are historically larger than what is seen in independent studies of noise due to pulse jitter \citep{2021MNRAS.502..407P}, and so it is possible that $\mathrm{E}_\mathrm{{C}}$ includes contributions from another stochastic noise source that has not yet been fully characterised.
We also note that for some pulsars the reduced-$\chi^2$ is less than unity. In many cases, we suspect that this is due to the influence of the $\mathrm{E}_\mathrm{{C}}$ noise parameter presenting at a greater amplitude than expected. An alternative explanation in pulsars possessing fewer observing epochs and binary orbital parameters, is that the generalised least squares fit performed by \textsc{tempo2} is absorbing some of the noise in the deterministic timing parameters.

\begin{table}

\caption{Models compared to determine the most probable representation of the noise processes in each pulsar.}
\label{table: noise_models}
\resizebox{\columnwidth}{!}{%
\begin{tabular}{@{}ccccc@{}}
\toprule
Model & $\mathrm{E}_\mathrm{Q}$ + $\mathrm{E}_\mathrm{F}$ & $\mathrm{E}_\mathrm{C}$ & DM Noise & \begin{tabular}[c]{@{}c@{}}Achromatic \\ Red Noise\end{tabular} \\ \midrule
$\mathrm{M}_{1}$ & \checkmark & -          & -          & -          \\
$\mathrm{M}_{2}$ & \checkmark & \checkmark & -          & -          \\
$\mathrm{M}_{3}$ & \checkmark & -          & -          & \checkmark \\
$\mathrm{M}_{4}$ & \checkmark & -          & \checkmark & -          \\
$\mathrm{M}_{5}$ & \checkmark & \checkmark & -          & \checkmark \\
$\mathrm{M}_{6}$ & \checkmark & \checkmark & \checkmark & -          \\
$\mathrm{M}_{7}$ & \checkmark & \checkmark & \checkmark & \checkmark \\
$\mathrm{M}_{8}$ & \checkmark & -          & \checkmark & \checkmark \\ \bottomrule
\end{tabular}%
}
\end{table}

\subsection{Residuals and noise subtraction}
\label{Subsection: residuals}

We used the timing residuals to assess the quality of the data. Two sets of timing residuals were formed from the sub-banded data: the first with frequency resolution where each sub-banded ToA correlates directly to a timing residual, and the second in which the timing residuals have been averaged together such that a single timing residual is produced for all ToAs that are recorded simultaneously. It was clear through visual inspection that some pulsars possessed evidence for time correlated noise, even on the 2.5-year data release timeline.
We can measure the impact of the noise processes on the data by comparing the residuals before and after subtracting the maximum likelihood realisation of the process. An example of this can be seen in Figure \ref{fig: 1909_noise_comparison}.

In Table \ref{table: MPTA_results_short} we present the MPTA timing results, and the comparisons between the models we have described here for ten pulsars. The full results are provided in Table \ref{table: MPTA_results_full}.

\begin{table*}

\caption{MPTA timing results and indicative noise parameters for the first ten MPTA pulsars in order of Right Ascension. The full results of all timing analyses, and the full indicative noise parameters, can be found in Table \ref{table: MPTA_results_full}.}
\label{table: MPTA_results_short}
\makebox[\textwidth]{%
\begin{tabular}[c]{lrccccccccccc}
\toprule
Pulsar &
  \multicolumn{1}{l|}{$\textrm{N}_{\textrm{ToA}}$} &
  \multicolumn{1}{c|}{$\mathrm{M}_\mathrm{N}$} &
  \multicolumn{4}{c|}{Time-Correlated Noise} &
  \multicolumn{2}{c|}{\begin{tabular}[c]{@{}c@{}}Sub-banded \\ wRMS ($\mu$s)\end{tabular}} &
  \multicolumn{4}{c}{\begin{tabular}[c]{@{}c@{}}Band-averaged \\ wRMS ($\mu$s)\end{tabular}} \\
  \cmidrule(lr){4-7}\cmidrule(lr){8-9}\cmidrule(lr){10-13}
   &
  \multicolumn{1}{l|}{} &
  \multicolumn{1}{c|}{} &
  $\mathrm{log}_{10}\textrm{A}_{\textrm{DM}}$ &
  \multicolumn{1}{c|}{$\gamma_{\textrm{DM}}$} &
  $\mathrm{log}_{10}\textrm{A}_{\text{red}}$ &
  $\gamma_{\textrm{red}}$ &
  Full &
  Whitened &
  Full &
  $\chi^{2}_{\mathrm{r},\mathrm{full}}$ &
  Whitened &
  $\chi^{2}_{\mathrm{r},\mathrm{white}}$ \\ \midrule

J0030+0451 & 747 & 1 & - & - & - & - & 1.895 & - & 0.457 & 1.146 & - & - \\
J0125-2327 & 893 & 6 & -12.3(1) & 4.0(2) & - & - & 0.718 & 0.634 & 0.549 & 2.114 & 0.402 & 1.130 \\
J0437-4715 & 1052 & 6 & -11.8(1) & 1.3(4) & - & - & 0.298 & 0.184 & 0.268 & 2.513 & 0.157 & 0.859 \\
J0610-2100 & 801 & 4 & -11.4(2) & 1.4(8) & - & - & 1.947 & 1.838 & 0.658 & 2.762 & 0.295 & 0.555 \\
J0613-0200 & 687 & 1 & - & - & - & - & 0.832 & - & 0.220 & 1.804 & - & - \\
J0614-3329 & 700 & 1 & - & - & - & - & 1.482 & - & 0.367 & 1.106 & - & - \\
J0636-3044 & 540 & 1 & - & - & - & - & 2.095 & - & 0.411 & 0.697 & - & - \\
J0711-6830 & 775 & 1 & - & - & - & - & 1.475 & - & 0.431 & 1.436 & - & - \\
J0900-3144 & 896 & 6 & -10.9(1) & 1.8(5) & - & - & 3.587 & 2.328 & 2.624 & 5.308 & 0.985 & 0.748 \\
J0931-1902 & 419 & 1 & - & - & - & - & 1.856 & - & 0.602 & 1.679 & - & - \\
... & ... & ... & ... & ... & ... & ... & ... & ... & ... & ... & ... & ... \\ \bottomrule
\end{tabular}%
}
\end{table*}%

\begin{figure}

    \includegraphics[width=\columnwidth]{noise_reduce_steps.pdf}
    \caption[]{Timing residuals for PSR J1909$-$3744. (Top) Full sub-banded timing residuals with no noise reduction. (Middle) Sub-banded timing residuals fully whitened via the removal of all identified red-noise processes. (Bottom) Fully whitened band-averaged timing residuals. Note the change in y-axis scale by a factor of 10 for ease of contrast.}
    \label{fig: 1909_noise_comparison}
\end{figure}

\section{Discussion}
\label{Section: Discussion}
\subsection{Findings from the MPTA Data Set}
\label{Subsection: Data set performance}

The initial timing analysis demonstrates the potential of the MPTA. Of the pulsars included in this data release $26$ are unique to the MPTA, and $67$ are able to achieve sub-microsecond precision using the current observing strategy. In comparison to the major PTAs, the MPTA leads in this performance statistic. Band-averaged, whitened sub-microsecond precision is achieved for $13$ of the PPTA pulsars \citep{2020PASA...37...20K}, for $5$ of the EPTA pulsars \citep{2016MNRAS.458.3341D}, and for $44$ of the pulsars regularly monitored by NANOGrav \citep{2021ApJS..252....4A}.

In the coming years, the primary goal of the MPTA will be the detection of nanohertz-frequency gravitational radiation, collaborating with the IPTA towards the goal of the detection of an angularly correlated signal. It is important to note that while the sensitivity of a PTA to the SGWB increases strongly with the length of the data span, it is also directly proportional to the number of pulsars in the array \citep{2013CQGra..30v4015S}. As noted in \citet{spiewak2022}, the inclusion of the MPTA and its otherwise unobserved pulsars, results in the MPTA potentially becoming the most significant contributor to IPTA searches for the SGWB by 2025, unless cadences or sensitivities at the other PTAs change. In addition, there is work currently underway investigating improvements to the MPTA, potentially resulting in a more aggressive timeline for improved sensitivity (Middleton et al., in prep.).

The sensitivity of the MPTA data set enabled investigations of the relationship between the change in the profile shape through frequency, and the expected timing precision. As noted in Section \ref{Subsection: Portrait_creation}, at least some degree of profile evolution through frequency was noted in all the pulsars in this data release. The effective pulse width,
\begin{equation}
\label{eq: eff_width}
    \mathrm{W_{eff}} = \frac{\Delta \phi}{\Sigma_{\mathrm{i}}[\mathrm{P}(\phi_{\mathrm{i}+1}) - \mathrm{P}(\phi_{\mathrm{i}})]^{2}},
\end{equation}
where $\Delta \phi$ is the resolution of the profile in pulse phase, $\phi$ is defined as the positional pulse phase, and P is the pulse profile flux normalized to unity, is measured for each pulsar's frequency sub-bands. By weighting each sub-band observation by its corresponding S/N, a linear relationship between the $\textrm{W}_\textrm{eff}$ and the timing uncertainty is empirically derived across the MPTA dataset,
\begin{equation}
\label{eq: uncertainty_width}
    \sigma_{\textrm{s}}(\textrm{W}_\textrm{eff}) = \alpha \frac{\textrm{W}_\textrm{eff}}{(\textrm{S}/\textrm{N})_\textrm{obs,s}},
\end{equation}
where $\sigma_\textrm{s}$ is the ToA uncertainty at each sub-band, and $(\textrm{S}/\textrm{N})_\textrm{obs,s}$ is the S/N for each sub-band of the observation. We use the \texttt{polyfit} module of the \textsc{numpy} python package to measure the coefficient of this relationship to be $\alpha = 1.249 \pm 0.002$ and the relative scatter on this approximation is defined as,
\begin{equation}
\label{eq: var_unc_width}
    \delta = \frac{1}{N}\sum\frac{|\sigma_\mathrm{obs} - \sigma_\mathrm{s}(\textrm{W}_\textrm{eff})|}{\sigma_\mathrm{s}(\textrm{W}_\textrm{eff})},
\end{equation}
where N is the number of data points and $\sigma_\mathrm{obs}$ is the uncertainty at each observation. Evaluating this expression results in $\delta = 0.201$, implying a fractional uncertainty of $\sim 20\%$ of the data relating to this model.

We also investigated how the effective width varies with frequency. Naively, one would expect the effective width to vary monotonically with frequency. For example, if the pulse profile is affected by scattering from a thin screen in the ISM we would expect the relationship between the thin-screen scattering timescale ($\tau_\textrm{s}$) and frequency ($f$) to evolve as \citep{1977ARA&A..15..479R},

\begin{equation}
\label{eq: scattering_timescale}
    \tau_\textrm{s} \propto f^{-4}.
\end{equation}

It is also thought that higher-frequency emission mechanisms occur closer to the surface of the pulsar, resulting in narrower pulse widths from the relative radius of the emission cone \citep{1978ApJ...222.1006C, 1999ApJ...526..957K, 2015MNRAS.449.3223D}.

While we found $44$ of the MPTA pulsars follow power-law relationships of effective width and frequency, others demonstrated far more complicated relationships better described by polynomials of order $>1$ (Figure \ref{fig: effective_width_unc}). These relationships are counter-intuitive to how pulsar emission is thought to occur, as well as how it is thought to be affected by the ISM, with recent studies of non-recycled pulsars identifying an expected, largely monotonic, width evolution through frequency with both positive and negative slopes \citep{2021MNRAS.508.4249P}. By combining the effective width measurements with the flux density it will be possible to further optimise observing strategies, and assess the utility of observations in other bands at MeerKAT, or with other telescopes such as the SKA.

\begin{figure}

    \includegraphics[width=\columnwidth]{eff_width_unc.pdf}
    \caption[]{Two pulsars (Left: J1909-3744; Right: J1327-0755) that possess decidedly different relationships between effective width ($W_\text{eff}$) and frequency. The relationship exhibited by J1909$-$3744 is seen more often, whereas J1327$-$0755 exhibits nonmonotonic, and  almost parabolic, evolution through frequency.}
    \label{fig: effective_width_unc}
\end{figure}

\subsection{Clock Correction Recovery}
\label{Subsection: Clock recovery}

The main purpose of a PTA is to detect signals common to the ensemble of pulsars. To demonstrate the potential of the MPTA, we use the data release to recover the MeerKAT observatory time standard.

The MeerKAT station clock uses two Hydrogen masers, two Rubidium atomic clocks, and a quartz crystal. This station clock produces a time standard (Karoo Telescope Time, KTT). The time standard drifts by a few nanoseconds per day relative to the Coordinate Universal Time (UTC) as reported by the Global Positioning Service (GPS). Every few months the masers are adjusted to keep the time within 500~ns of UTC (GPS) \citep{2020PASA...37...28B}.

For accurate pulsar timing, the drifting offset of KTT relative to UTC  is measured through a system described in detail in \citet{2022JATIS...8a1013B}, and results in the production of a clock correction\footnote{https://archive-gw-1.kat.ac.za/public/tfr/mk2utc.clk}. This correction is applied using pulsar timing software (in our case \textsc{tempo2}), which corrects for the resultant offset in the pulse arrival times to a precision of less than $4.925$ ns \citep{2022JATIS...8a1013B}.

Here we derive the KTT from the pulsar timing measurements. To do this, we use the approach described in \citet{2020MNRAS.491.5951H}, in order to constrain the properties of the clock signal through a Bayesian methodology. This was chosen as opposed to the frequentist methods, which is also presented in \citet{2020MNRAS.491.5951H}, as the Bayesian approach is now adopted by most PTA collaborations.

We modelled the clock signal as a power-law, red noise process of the form
\begin{equation}
\label{eq: clock_power_spectrum}
    S_{\mathrm{clk}}(f) = \frac{{A^{2}_{\mathrm{clk}}}}{12\pi^{2}}\left(\frac{f}{f_c}\right)^{-\gamma_\mathrm{clk}},
\end{equation}
where the coefficients are of the same form as Equation \ref{eq: power_spec_red}.

To perform parameter estimation we employed a Parallel Tempered Markov Chain Monte Carlo (PTMCMC) \citep{justin_ellis_2017_1037579} sampler via the \textsc{enterprise} software \citep{2019ascl.soft12015E}, and allowed \textsc{enterprise} to simultaneously model individual pulsar red noise processes (isolated to a general achromatic red noise process and a dispersion measure noise process) as well as a common red noise process of the form given in Equation \ref{eq: clock_power_spectrum} across all pulsars in the sample. This parameter estimation resulted in a common signal possessing the noise characteristics $\log_{10}\mathrm{A}_\mathrm{clk} = -13.35$ and $\gamma_\mathrm{clk} = 2.1$.

Following this, we constructed noise covariance matrices for each pulsar of the form
\begin{equation}
\label{eq: cov_total}
    \textbf{C}_\mathrm{PSR} = \textbf{C}_\mathrm{white} + \textbf{C}_\mathrm{red} + \textbf{C}_\mathrm{DM},
\end{equation}
where $\textbf{C}_\mathrm{white}$ models the white noise components, including $\mathrm{E}_\mathrm{C}$ (noise that is correlated between sub-banded ToAs) where required. Without the addition of $\mathrm{E}_\mathrm{C}$, this is a diagonal matrix populated with the arrival time uncertainties, influenced by the white noise terms as per Equation \ref{eq: white_noise_BM}. With the addition of $\mathrm{E}_\mathrm{C}$, the matrix is no longer diagonal, but accounts for a covariance value between arrival times from the same observation.

$\textbf{C}_\mathrm{DM}$ and $\textbf{C}_\mathrm{red}$ are the covariance matrices relating to the DM noise and the achromatic red noise processes, accounting for their presence where required. We followed the definitions of the covariance matrices presented in \citet{2014MNRAS.441.2831L}. The covariance between the ToAs $i$ and $j$ are

\begin{equation}
\label{eq: covariance_DM}
    \textbf{C}_\mathrm{DM,ij} = \frac{\kappa^{2}\int_{1/T}^{\infty}\mathrm{S}_\mathrm{DM}(f)\cos{(2\pi f t_{ij})}\,df}{\nu^{2}_\mathrm{i} \nu^{2}_\mathrm{j}},
\end{equation}
and
\begin{equation}
\label{eq: covariance_red}
    \textbf{C}_\mathrm{red,ij} = \int_{1/T}^{\infty}\mathrm{S}_\mathrm{red}(f)\cos{(2\pi f t_{ij})}\,df,
\end{equation}
where T is the total time span of the data, $f$ is the fluctuation frequency (as opposed to observing frequency $\nu$), $\mathrm{t}_\mathrm{ij}$ is the time lag between arrival times i and j, and $\kappa$ is the dispersion constant $\kappa = 2.410 \times 10^{-16}$ $\mathrm{Hz}^{-2}\mathrm{cm}^{-3}\mathrm{pc}\mathrm{s}^{-1}$.
$\mathrm{S}_\mathrm{red}(f)$ and $\mathrm{S}_\mathrm{DM}(f)$ are the power spectral densities of the achromatic red noise and DM noise respectively and are defined as,
\begin{equation}
\label{eq: red_power_spectrum}
    S_{\mathrm{red}}(f) = \frac{{A^{2}_{\mathrm{red}}}}{12\pi^{2}}\left(\frac{f}{f_c}\right)^{-\gamma_\mathrm{red}},
\end{equation}
and
\begin{equation}
\label{eq: DM_power_spectrum}
    S_{\mathrm{DM}}(f) = {A^{2}_\mathrm{DM}}\left(\frac{f}{f_c}\right)^{-\gamma_\mathrm{DM}},
\end{equation}
noting the lack of the $1/12\pi^{2}$ normalisation factor in Equation \ref{eq: DM_power_spectrum}, following the convention of \citet{2014MNRAS.437.3004L}.

Using this methodology we could reconstruct the clock signal over an arbitrary date range, that can extend beyond our data span.
We did this by the use of ten of the most precise pulsars presented in this data release, combining the covariance matrices we describe above via
\begin{equation}
\label{eq: mega_matrix}
    \textbf{C}_\mathrm{PTA,10} = \begin{bmatrix}
\textbf{C}_{1} & 0 & 0 & ... & 0 \\
0 & \textbf{C}_{2} & 0 & ... & 0 \\
0 & 0 & \textbf{C}_{3} & ... & 0 \\
... & ... & ... & ... & ... \\
0 & 0 & 0 & ... & \textbf{C}_{10}
\end{bmatrix},
\end{equation}

where $\textbf{C}_\mathrm{N}$ represents each of the 10 pulsars' respective $\textbf{C}_\mathrm{PSR}$ per Equation \ref{eq: cov_total}.

In order to evaluate the clock signal, we defined $N_\mathrm{rec}=500$ epochs at which we wanted to reconstruct it. We chose the initial epoch to be when KTT was first published at MJD 58484, and we chose to extend the reconstruction beyond our observations to MJD 59611, to investigate the ability to extrapolate the clock signal.

To do this we then extended $\textbf{C}_\mathrm{PTA,10}$ as per \citet{2012MNRAS.424..244D}, by $N_\mathrm{rec}$ rows and columns to create a matrix of the form

\begin{equation}
\label{eq: final_matrix}
    \textbf{C}_\mathrm{Final} = \begin{bmatrix}
\textbf{C}_\mathrm{PTA,10} & 0 \\
0 &  0
\end{bmatrix},
\end{equation}
populated by zeroes.

The clock covariance is then evaluated at not only our observations, but at the epochs of our reconstructed signal

\begin{equation}
    \textbf{C}_\mathrm{Total} = \textbf{C}_\mathrm{Final} + \textbf{C}_\mathrm{clk},
\end{equation}

where the covariance of the clock signal is defined as
\begin{equation}
\label{eq: covariance_clock}
    \textbf{C}_\mathrm{clk,ij} = \int_{1/T_\mathrm{clk}}^{\infty}\mathrm{S}_\mathrm{clk}(f)\cos{(2\pi f t_{ij})}\,df,
\end{equation}
with equation parameters as described in Equations \ref{eq: covariance_DM} and \ref{eq: covariance_red}. We note that the clock covariance is calculated simultaneously for both the epochs of the observations, and the epochs for the reconstructed signal. Further, $T_\mathrm{clk}$ is the span over which the clock signal is evaluated.

We defined the total residuals $\textbf{t}_\mathrm{res}$ at which this signal is evaluated, as the vector containing the observed timing residuals for the 10 pulsars $\textbf{t}_\mathrm{pf,PTA} = (\textbf{t}_\mathrm{pf,1}, ..., \textbf{t}_\mathrm{pf,10})$, and the perfect residuals we assigned to the epochs of our reconstructed signal $\textbf{t}_\mathrm{ideal}=\textbf{0}$.

This final covariance matrix was used to measure the recovered clock waveform \citep{2014MNRAS.441.2831L, 2020MNRAS.491.5951H},

\begin{equation}
\label{eq: clock_waveform}
    \textbf{t}_\mathrm{clk} = \textbf{C}_\mathrm{clk}\textbf{C}^{-1}_\mathrm{Total}\textbf{t}_\mathrm{res}.
\end{equation}

The covariance of the clock signal is
\begin{equation}
\label{eq: clock_uncertainties}
    \mathbf{\sigma}_\mathrm{clk} = \textbf{C}_\mathrm{clk} - \textbf{C}_\mathrm{clk}\textbf{C}^{-1}_\mathrm{Total}\textbf{C}_\mathrm{clk},
\end{equation}
where the uncertainty of the clock signal is the square root of the diagonal of $\mathbf{\sigma}_\mathrm{clk}$.

The recovered clock signal is presented in Figure \ref{fig: clock_recovery}. The precision to which the MPTA is able to extract an injected signal is evident, recovering the signal to within $50$ ns for the majority of the data span. There is evidence of a $\approx 150$\,n s deviation in the signal between MJDs $\sim58700$ and $\sim58900$, the cause of which is currently uncertain.
While the beginning of this apparent  offset is close to one of the known timing offsets in the data, the end is not associated with one.

First, we considered if the deviation was caused by an individual pulsar with poorly specified noise. To test this, we removed one of the pulsars from the method and reformed the clock signal. We found the feature persisted in all cases, suggesting an effect common to all of the pulsars. Possible other causes of the deviation could be unidentified systematic errors in the clock correction, an additional noise process that is common to the pulsars, or even the influence of a SGWB.

Owing to the unknown pulsar spin frequency and frequency derivatives, pulsars are not sensitive to linear or quadratic variations in time standards. We attempted to account for a possible quadratic drift in KTT, and found that the residual of the expected and the recovered clock signal is not able to be well modelled by a quadratic. We conclude that the deviation we observe is not expected to be caused by a quadratic. We also note that the signal we reconstructed is not well-described by a Gaussian process. Arising from this misspecification, it is possible that the reconstructed clock signal will not well recover the sharp (high fluctuation frequency) features in KTT. This also may contribute to the larger discrepancy between MJDs $\sim58700$ and $\sim58900$.

Uncertainties that exist in a telescope time standard must be kept far lower than the amplitude of the expected SGWB signal. Typically, observatory time standards are assumed to have precision much lower than the GW signal, however, this is the first time that this has been demonstrated through pulsar timing \citep{2016MNRAS.455.4339T}. The difference between the recovered clock corrections and the corresponding time standard can then be used as a measure of the accuracy of a timing array. We propose that the root-mean-square (rms) residual and its spectral characteristics form a quantitative measure of the quality of the timing array. These can then be compared with the amplitudes of any purported gravitational wave background to assess array performance. For the MPTA the rms residual is 62 ns over 2.5 years. Assuming a SGWB amplitude of $2 \times 10^{-15}$ (similar to the amplitude observed in the purported common red signal), the expected induced rms variations in five years of observational data is expected to be to be $\sim 127$ ns \cite[][]{2013MNRAS.428.1147V}.

\begin{figure*}

    \includegraphics[width=\textwidth]{clock_recovery_residual.pdf}
    \caption[]{(Top) The recovered clock correction signal and corresponding uncertainty as defined in Section \ref{Subsection: Clock recovery} (green, shaded), juxtaposed against the clock correction signal that was to be extracted from the data and corresponding uncertainty (blue, shaded). Note the uncertainty rapidly increasing following the end of the data span, as well as the marked deviation between MJDs $\sim58700$ and $\sim58900$. (Middle) The residual difference (teal) between the recovered clock signal and the KTT-UTC(GPS) clock correction signal. The teal shading corresponds to the uncertainties of both the recovered signal and KTT-UTC(GPS) summed together in quadrature. The dashed grey line and corresponding shading represent a $50$ ns window around $0$ ns, to demonstrate the precision of the recovered signal. (Bottom) The MJDs corresponding to observations used for this reconstruction.}
    \label{fig: clock_recovery}
\end{figure*}

\section{Summary and Conclusions}
\label{Section: Summary and Conclusions}

We have presented here the first data release from the MPTA, including the observations of $78$ pulsars spanning February 2019 to October 2021. As well as the recorded sub-banded and band-averaged arrival times for each pulsar, we also provide initial timing ephemerides, timing portraits, sub-banded and band-averaged timing residuals, and preliminary noise models. Of these pulsars, $67$ report sub-microsecond timing precision, with $2$ possessing a precision of less than $100$ ns from whitened, band-averaged data. We detail the observing strategy of the MPTA, and selection criteria for a pulsar to be deemed suitable for inclusion into it. Additionally, the data reduction pipeline used for this data release, and the nature of the data products output by it is outlined. We also describe how the accompanying portraits are created through the use of \textsc{PulsePortraiture}, and the process of noise model selection and the consequent subtraction of noise models from the data in order to provide more accurate timing residuals. All data products described here are open to the scientific community, and can be accessed via the Australian gravitational wave data portal\footnote{http://dx.doi.org/10.26185/6392814b27073}. The collaboration expects to release subsequent data releases with both the four and five year MPTA data sets when they are mature.

We have also investigated the relationship of the effective-width metric ($\mathrm{W}_\mathrm{eff}$) and observing frequency, and find somewhat counter-intuitive results that may influence the observing priority in frequency bands outside of the L-band range we report on. We provide initial noise models, showing that DM variations are detected in a substantial sub-sample of our pulsars. These noise models can be used as the starting point to develop more sophisticated models.

Lastly, we are able to demonstrate the precision of the MPTA data set by using a subset of the full MPTA to recover the observatory time standard, Karoo Telescope Time (KTT). Through most of the data we see good agreement, with the recovered signal within $50$ ns of KTT. We find a modest ($150$ ns) disagreement for $\sim 150$ days that appears to be common to all the pulsars. Work is ongoing to see if we can explain this discrepancy. We test if this deviation could be attributed to any individual pulsar we include in the signal recovery.  As we finding that it cannot, we briefly explore possible other causes.

The MPTA data set includes the largest number of regularly timed pulsars of any current PTA. In this data release, we do not include an analysis of the merits of wide-band timing methods as opposed to the sub-banded arrival times we calculate, and would encourage any interested parties to fully utilise the components of this data release to explore this technique. Individually, the data released here has the potential for a wide range of studies, many of which are already well underway by members of the MeerTime collaboration. This includes the study of individual pulsars including updated measurements of orbital parameters. We can also study the population as a whole. This includes an analysis of the astrometry in MPTA pulsars, aiming to improve MSP population-synthesis models; an investigation of the flux density variations in the MPTA population, both to study the ISM through diffractive and refractive scintillation, and to optimise pulsar timing observations via dynamic scheduling; and, naturally, searches for common noise and correlated signals to identify the influence of gravitational waves on the MPTA pulsars, as well as other signals that are common in the array.

Finally, the inclusion of the MPTA data set to the IPTA data streams is expected to make a significant impact in the coming years \citep{spiewak2022}, both through the inclusion of high-quality timing data of pulsars that are already observed, and through the inclusion of data from an additional 26 pulsars that are not yet observed in any other PTA.

\section*{Acknowledgements}

We thank T. Pennucci for discussions.
The MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation.
MTM, RMS, MB, and DR acknowledge support through the ARC centre of Excellence grant CE17010004 (OzGrav).
RMS acknowledges support through Australian Council Future Fellowship FT190100155.
VVK acknowledges continuing support from the Max Planck society.
GT acknowledges financial support from Agence Nationale de la Recherche (ANR-18-CE31-0015), France.
AK acknowledges funding from the STFC consolidated grant to Oxford Astrophysics, code ST/000488.
PTUSE was developed with support from the Australian SKA Office and Swinburne University of Technology, with financial contributions from the MeerTime collaboration members.
This work used the OzSTAR national facility at Swinburne University of Technology.
OzSTAR is funded by Swinburne University of Technology and the National Collaborative Research Infrastructure Strategy (NCRIS).

\section*{Data Availability}

All data used in this release available at the Australian Gravitational Wave Data Centre, at \url{http://dx.doi.org/10.26185/6392814b27073}.
The data provided includes both sub-banded and frequency averaged ToAs, the full data archives used to construct this data release, the ephemerides that have been used to perform timing, and band-averaged and sub-banded timing residuals. Also included are both the timing standards (frequency averaged, one dimensional) and frequency resolved portraits, used to calculate each set of ToAs as well as the preliminary noise models constructed to facilitate this work.

The archives, standards, and portraits are in \textsc{psrfits} file format. The ephemerides and residuals are in a standard \texttt{ascii} text file format, and the arrival times are supplied as IFF data. Information regarding the preliminary noise models can be found contained within the ephemerides.
Data from four pulsars (PSRs J0614$-$3329, J1036$-$8317, J1756$-$2251, and J1902$-$5105) will be available after the 18-month SARAO embargo is lifted. An updated DOI will be provided at this time.

\bibliographystyle{mnras}

\section{Introduction}

The journal \textit{Monthly Notices of the Royal Astronomical Society} (MNRAS) encourages authors to prepare their papers using \LaTeX.
The style file \verb'mnras.cls' can be used to approximate the final appearance of the journal, and provides numerous features to simplify the preparation of papers.
This document, \verb'mnras_guide.tex', provides guidance on using that style file and the features it enables.

This is not a general guide on how to use \LaTeX, of which many excellent examples already exist.
We particularly recommend \textit{Wikibooks \LaTeX}\footnote{\url{https://en.wikibooks.org/wiki/LaTeX}}, a collaborative online textbook which is of use to both beginners and experts.
Alternatively there are several other online resources, and most academic libraries also hold suitable beginner's guides.

For guidance on the contents of papers, journal style, and how to submit a paper, see the MNRAS Instructions to Authors\footnote{\label{foot:itas}\url{http://www.oxfordjournals.org/our_journals/mnras/for_authors/}}.
Only technical issues with the \LaTeX\ class are considered here.

\section{Obtaining and installing the MNRAS package}
Some \LaTeX\ distributions come with the MNRAS package by default.
If yours does not, you can either install it using your distribution's package manager, or download it from the Comprehensive \TeX\ Archive Network\footnote{\url{http://www.ctan.org/tex-archive/macros/latex/contrib/mnras}} (CTAN).

The files can either be installed permanently by placing them in the appropriate directory (consult the documentation for your \LaTeX\ distribution), or used temporarily by placing them in the working directory for your paper.

To use the MNRAS package, simply specify \verb'mnras' as the document class at the start of a \verb'.tex' file:

\begin{verbatim}
\documentclass{mnras}
\end{verbatim}
Then compile \LaTeX\ (and if necessary \bibtex) in the usual way.

\section{Preparing and submitting a paper}
We recommend that you start with a copy of the \texttt{mnras\_template.tex} file.
Rename the file, update the information on the title page, and then work on the text of your paper.
Guidelines for content, style etc. are given in the instructions to authors on the journal's website$^{\ref{foot:itas}}$.
Note that this document does not follow all the aspects of MNRAS journal style (e.g. it has a table of contents).

If a paper is accepted, it is professionally typeset and copyedited by the publishers.
It is therefore likely that minor changes to presentation will occur.
For this reason, we ask authors to ignore minor details such as slightly long lines, extra blank spaces, or misplaced figures, because these details will be dealt with during the production process.

Papers must be submitted electronically via the online submission system; paper submissions are not permitted.
For full guidance on how to submit a paper, see the instructions to authors.

\section{Class options}
\label{sec:options}
There are several options which can be added to the document class line like this:

\begin{verbatim}
\documentclass[option1,option2]{mnras}
\end{verbatim}
The available options are:
\begin{itemize}
\item \verb'letters' -- used for papers in the journal's Letters section.
\item \verb'onecolumn' -- single column, instead of the default two columns. This should be used {\it only} if necessary for the display of numerous very long equations.
\item \verb'doublespacing' -- text has double line spacing. Please don't submit papers in this format.
\item \verb'referee' -- \textit{(deprecated)} single column, double spaced, larger text, bigger margins. Please don't submit papers in this format.
\item \verb'galley' -- \textit{(deprecated)} no running headers, no attempt to align the bottom of columns.
\item \verb'landscape' -- \textit{(deprecated)} sets the whole document on landscape paper.
\item \verb"usenatbib" -- \textit{(all papers should use this)} this uses Patrick Daly's \verb"natbib.sty" package for citations.
\item \verb"usegraphicx" -- \textit{(most papers will need this)} includes the \verb'graphicx' package, for inclusion of figures and images.
\item \verb'useAMS' -- adds support for upright Greek characters \verb'\upi', \verb'\umu' and \verb'\upartial' ($\upi$, $\umu$ and $\upartial$). Only these three are included, if you require other symbols you will need to include the \verb'amsmath' or \verb'amsymb' packages (see section~\ref{sec:packages}).
\item \verb"usedcolumn" -- includes the package \verb"dcolumn", which includes two new types of column alignment for use in tables.
\end{itemize}

Some of these options are deprecated and retained for backwards compatibility only.
Others are used in almost all papers, but again are retained as options to ensure that papers written decades ago will continue to compile without problems.
If you want to include any other packages, see section~\ref{sec:packages}.

\section{Title page}

If you are using \texttt{mnras\_template.tex} the necessary code for generating the title page, headers and footers is already present.
Simply edit the title, author list, institutions, abstract and keywords as described below.

\subsection{Title}
There are two forms of the title: the full version used on the first page, and a short version which is used in the header of other odd-numbered pages (the `running head').
Enter them with \verb'\title[]{}' like this:
\begin{verbatim}
\title[Running head]{Full title of the paper}
\end{verbatim}
The full title can be multiple lines (use \verb'\\' to start a new line) and may be as long as necessary, although we encourage authors to use concise titles. The running head must be $\le~45$ characters on a single line.

See appendix~\ref{sec:advanced} for more complicated examples.

\subsection{Authors and institutions}

Like the title, there are two forms of author list: the full version which appears on the title page, and a short form which appears in the header of the even-numbered pages. Enter them using the \verb'\author[]{}' command.

If the author list is more than one line long, start a new line using \verb'\newauthor'. Use \verb'\\' to start the institution list. Affiliations for each author should be indicated with a superscript number, and correspond to the list of institutions below the author list.

For example, if I were to write a paper with two coauthors at another institution, one of whom also works at a third location:
\begin{verbatim}
\author[K. T. Smith et al.]{
Keith T. Smith,$^{1}$
A. N. Other,$^{2}$
and Third Author$^{2,3}$
\\
$^{1}$Affiliation 1\\
$^{2}$Affiliation 2\\
$^{3}$Affiliation 3}
\end{verbatim}
Affiliations should be in the format `Department, Institution, Street Address, City and Postal Code, Country'.

Email addresses can be inserted with the \verb'\thanks{}' command which adds a title page footnote.
If you want to list more than one email, put them all in the same \verb'\thanks' and use \verb'\footnotemark[]' to refer to the same footnote multiple times.
Present addresses (if different to those where the work was performed) can also be added with a \verb'\thanks' command.

\subsection{Abstract and keywords}

The abstract is entered in an \verb'abstract' environment:
\begin{verbatim}
\begin{abstract}
The abstract of the paper.
\end{abstract}
\end{verbatim}
\noindent Note that there is a word limit on the length of abstracts.
For the current word limit, see the journal instructions to authors$^{\ref{foot:itas}}$.

Immediately following the abstract, a set of keywords is entered in a \verb'keywords' environment:
\begin{verbatim}
\begin{keywords}
keyword 1 -- keyword 2 -- keyword 3
\end{keywords}
\end{verbatim}
\noindent There is a list of permitted keywords, which is agreed between all the major astronomy journals and revised every few years.
Do \emph{not} make up new keywords!
For the current list of allowed keywords, see the journal's instructions to authors$^{\ref{foot:itas}}$.

\section{Sections and lists}

Sections and lists are generally the same as in the standard \LaTeX\ classes.

\subsection{Sections}
\label{sec:sections}
Sections are entered in the usual way, using \verb'\section{}' and its variants. It is possible to nest up to four section levels:
\begin{verbatim}
\section{Main section}
 \subsection{Subsection}
  \subsubsection{Subsubsection}
   \paragraph{Lowest level section}
\end{verbatim}
\noindent The other \LaTeX\ sectioning commands \verb'\part', \verb'\chapter' and \verb'\subparagraph{}' are deprecated and should not be used.

Some sections are not numbered as part of journal style (e.g. the Acknowledgements).
To insert an unnumbered section use the `starred' version of the command: \verb'\section*{}'.

See appendix~\ref{sec:advanced} for more complicated examples.

\subsection{Lists}

Two forms of lists can be used in MNRAS -- numbered and unnumbered.

For a numbered list, use the \verb'enumerate' environment:
\begin{verbatim}
\begin{enumerate}
 \item First item
 \item Second item
 \item etc.
\end{enumerate}
\end{verbatim}
\noindent which produces
\begin{enumerate}
 \item First item
 \item Second item
 \item etc.
\end{enumerate}
Note that the list uses lowercase Roman numerals, rather than the \LaTeX\ default Arabic numerals.

For an unnumbered list, use the \verb'description' environment without the optional argument:
\begin{verbatim}
\begin{description}
 \item First item
 \item Second item
 \item etc.
\end{description}
\end{verbatim}
\noindent which produces
\begin{description}
 \item First item
 \item Second item
 \item etc.
\end{description}

Bulleted lists using the \verb'itemize' environment should not be used in MNRAS; it is retained for backwards compatibility only.

\section{Mathematics and symbols}

The MNRAS class mostly adopts standard \LaTeX\ handling of mathematics, which is briefly summarised here.
See also section~\ref{sec:packages} for packages that support more advanced mathematics.

Mathematics can be inserted into the running text using the syntax \verb'$1+1=2$', which produces $1+1=2$.
Use this only for short expressions or when referring to mathematical quantities; equations should be entered as described below.

\subsection{Equations}
Equations should be entered using the \verb'equation' environment, which automatically numbers them:

\begin{verbatim}
\begin{equation}
 a^2=b^2+c^2
\end{equation}
\end{verbatim}
\noindent which produces
\begin{equation}
 a^2=b^2+c^2
\end{equation}

By default, the equations are numbered sequentially throughout the whole paper. If a paper has a large number of equations, it may be better to number them by section (2.1, 2.2 etc.). To do this, add the command \verb'\numberwithin{equation}{section}' to the preamble.

It is also possible to produce un-numbered equations by using the \LaTeX\ built-in \verb'\['\textellipsis\verb'\]' and \verb'$$'\textellipsis\verb'$$' commands; however MNRAS requires that all equations are numbered, so these commands should be avoided.

\subsection{Special symbols}

\begin{table}
 \caption{Additional commands for special symbols commonly used in astronomy. These can be used anywhere.}
 \label{tab:anysymbols}
 \begin{tabular*}{\columnwidth}{@{}l@{\hspace*{50pt}}l@{\hspace*{50pt}}l@{}}
  \hline
  Command & Output & Meaning\\
  \hline
  \verb'\sun' & \sun & Sun, solar\\[2pt]
  \verb'\earth' & \earth & Earth, terrestrial\\[2pt]
  \verb'\micron' & \micron & microns\\[2pt]
  \verb'\degr' & \degr & degrees\\[2pt]
  \verb'\arcmin' & \arcmin & arcminutes\\[2pt]
  \verb'\arcsec' & \arcsec & arcseconds\\[2pt]
  \verb'\fdg' & \fdg & fraction of a degree\\[2pt]
  \verb'\farcm' & \farcm & fraction of an arcminute\\[2pt]
  \verb'\farcs' & \farcs & fraction of an arcsecond\\[2pt]
  \verb'\fd' & \fd & fraction of a day\\[2pt]
  \verb'\fh' & \fh & fraction of an hour\\[2pt]
  \verb'\fm' & \fm & fraction of a minute\\[2pt]
  \verb'\fs' & \fs & fraction of a second\\[2pt]
  \verb'\fp' & \fp & fraction of a period\\[2pt]
  \verb'\diameter' & \diameter & diameter\\[2pt]
  \verb'\sq' & \sq & square, Q.E.D.\\[2pt]
  \hline
 \end{tabular*}
\end{table}

\begin{table}
 \caption{Additional commands for mathematical symbols. These can only be used in maths mode.}
 \label{tab:mathssymbols}
 \begin{tabular*}{\columnwidth}{l@{\hspace*{40pt}}l@{\hspace*{40pt}}l}
  \hline
  Command & Output & Meaning\\
  \hline
  \verb'\upi' & $\upi$ & upright pi\\[2pt]
  \verb'\umu' & $\umu$ & upright mu\\[2pt]
  \verb'\upartial' & $\upartial$ & upright partial derivative\\[2pt]
  \verb'\lid' & $\lid$ & less than or equal to\\[2pt]
  \verb'\gid' & $\gid$ & greater than or equal to\\[2pt]
  \verb'\la' & $\la$ & less than of order\\[2pt]
  \verb'\ga' & $\ga$ & greater than of order\\[2pt]
  \verb'\loa' & $\loa$ & less than approximately\\[2pt]
  \verb'\goa' & $\goa$ & greater than approximately\\[2pt]
  \verb'\cor' & $\cor$ & corresponds to\\[2pt]
  \verb'\sol' & $\sol$ & similar to or less than\\[2pt]
  \verb'\sog' & $\sog$ & similar to or greater than\\[2pt]
  \verb'\lse' & $\lse$ & less than or homotopic to \\[2pt]
  \verb'\gse' & $\gse$ & greater than or homotopic to\\[2pt]
  \verb'\getsto' & $\getsto$ & from over to\\[2pt]
  \verb'\grole' & $\grole$ & greater over less\\[2pt]
  \verb'\leogr' & $\leogr$ & less over greater\\
  \hline
 \end{tabular*}
\end{table}

Some additional symbols of common use in astronomy have been added in the MNRAS class. These are shown in tables~\ref{tab:anysymbols}--\ref{tab:mathssymbols}. The command names are -- as far as possible -- the same as those used in other major astronomy journals.

Many other mathematical symbols are also available, either built into \LaTeX\ or via additional packages. If you want to insert a specific symbol but don't know the \LaTeX\ command, we recommend using the Detexify website\footnote{\url{http://detexify.kirelabs.org}}.

Sometimes font or coding limitations mean a symbol may not get smaller when used in sub- or superscripts, and will therefore be displayed at the wrong size. There is no need to worry about this as it will be corrected by the typesetter during production.

To produce bold symbols in mathematics, use \verb'\bmath' for simple variables, and the \verb'bm' package for more complex symbols (see section~\ref{sec:packages}). Vectors are set in bold italic, using \verb'\mathbfit{}'.

For matrices, use \verb'\mathbfss{}' to produce a bold sans-serif font e.g. \mathbfss{H}; this works even outside maths mode, but not all symbols are available (e.g. Greek). For $\nabla$ (del, used in gradients, divergence etc.) use \verb'$\nabla$'.

\subsection{Ions}

A new \verb'\ion{}{}' command has been added to the class file, for the correct typesetting of ionisation states.
For example, to typeset singly ionised calcium use \verb'\ion{Ca}{ii}', which produces \ion{Ca}{ii}.

\section{Figures and tables}
\label{sec:fig_table}
Figures and tables (collectively called `floats') are mostly the same as built into \LaTeX.

\subsection{Basic examples}
\begin{figure}
 \includegraphics[width=\columnwidth]{example}
 \caption{An example figure.}
 \label{fig:example}
\end{figure}
Figures are inserted in the usual way using a \verb'figure' environment and \verb'\includegraphics'. The example Figure~\ref{fig:example} was generated using the code:
\begin{verbatim}
\begin{figure}
 \includegraphics[width=\columnwidth]{example}
 \caption{An example figure.}
 \label{fig:example}
\end{figure}
\end{verbatim}

\begin{table}
 \caption{An example table.}
 \label{tab:example}
 \begin{tabular}{lcc}
  \hline
  Star & Mass & Luminosity\\
   & $M_{\sun}$ & $L_{\sun}$\\
  \hline
  Sun & 1.00 & 1.00\\
  $\alpha$~Cen~A & 1.10 & 1.52\\
  $\epsilon$~Eri & 0.82 & 0.34\\
  \hline
 \end{tabular}
\end{table}
The example Table~\ref{tab:example} was generated using the code:
\begin{verbatim}
\begin{table}
 \caption{An example table.}
 \label{tab:example}
 \begin{tabular}{lcc}
  \hline
  Star & Mass & Luminosity\\
   & $M_{\sun}$ & $L_{\sun}$\\
  \hline
  Sun & 1.00 & 1.00\\
  $\alpha$~Cen~A & 1.10 & 1.52\\
  $\epsilon$~Eri & 0.82 & 0.34\\
  \hline
 \end{tabular}
\end{table}
\end{verbatim}

\subsection{Captions and placement}
Captions go \emph{above} tables but \emph{below} figures, as in the examples above.

The \LaTeX\ float placement commands \verb'[htbp]' are intentionally disabled.
Layout of figures and tables will be adjusted by the publisher during the production process, so authors should not concern themselves with placement to avoid disappointment and wasted effort.
Simply place the \LaTeX\ code close to where the figure or table is first mentioned in the text and leave exact placement to the publishers.

By default a figure or table will occupy one column of the page.
To produce a wider version which covers both columns, use the \verb'figure*' or \verb'table*'  environment.

If a figure or table is too long to fit on a single page it can be split it into several parts.
Create an additional figure or table which uses \verb'\contcaption{}' instead of \verb'\caption{}'.
This will automatically correct the numbering and add `\emph{continued}' at the start of the caption.
\begin{table}
 \contcaption{A table continued from the previous one.}
 \label{tab:continued}
 \begin{tabular}{lcc}
  \hline
  Star & Mass & Luminosity\\
   & $M_{\sun}$ & $L_{\sun}$\\
  \hline
  $\tau$~Cet & 0.78 & 0.52\\
  $\delta$~Pav & 0.99 & 1.22\\
  $\sigma$~Dra & 0.87 & 0.43\\
  \hline
 \end{tabular}
\end{table}
Table~\ref{tab:continued} was generated using the code:

\begin{verbatim}
\begin{table}
 \contcaption{A table continued from the previous one.}
 \label{tab:continued}
 \begin{tabular}{lcc}
  \hline
  Star & Mass & Luminosity\\
   & $M_{\sun}$ & $L_{\sun}$\\
  \hline
  $\tau$~Cet & 0.78 & 0.52\\
  $\delta$~Pav & 0.99 & 1.22\\
  $\sigma$~Dra & 0.87 & 0.43\\
  \hline
 \end{tabular}
\end{table}
\end{verbatim}

To produce a landscape figure or table, use the \verb'pdflscape' package and the \verb'landscape' environment.
The landscape Table~\ref{tab:landscape} was produced using the code:
\begin{verbatim}
\begin{landscape}
 \begin{table}
  \caption{An example landscape table.}
  \label{tab:landscape}
  \begin{tabular}{cccccccccc}
    \hline
    Header & Header & ...\\
    Unit & Unit & ...\\
    \hline
    Data & Data & ...\\
    Data & Data & ...\\
    ...\\
    \hline
  \end{tabular}
 \end{table}
\end{landscape}
\end{verbatim}
Unfortunately this method will force a page break before the table appears.
More complicated solutions are possible, but authors shouldn't worry about this.

\begin{landscape}
 \begin{table}
  \caption{An example landscape table.}
  \label{tab:landscape}
  \begin{tabular}{cccccccccc}
    \hline
    Header & Header & Header & Header & Header & Header & Header & Header & Header & Header\\
    Unit & Unit & Unit & Unit & Unit & Unit & Unit & Unit & Unit & Unit \\
    \hline
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    Data & Data & Data & Data & Data & Data & Data & Data & Data & Data\\
    \hline
  \end{tabular}
 \end{table}
\end{landscape}

\section{References and citations}

\subsection{Cross-referencing}

The usual \LaTeX\ commands \verb'\label{}' and \verb'\ref{}' can be used for cross-referencing within the same paper.
We recommend that you use these whenever relevant, rather than writing out the section or figure numbers explicitly.
This ensures that cross-references are updated whenever the numbering changes (e.g. during revision) and provides clickable links (if available in your compiler).

It is best to give each section, figure and table a logical label.
For example, Table~\ref{tab:mathssymbols} has the label \verb'tab:mathssymbols', whilst section~\ref{sec:packages} has the label \verb'sec:packages'.
Add the label \emph{after} the section or caption command, as in the examples in sections~\ref{sec:sections} and \ref{sec:fig_table}.
Enter the cross-reference with a non-breaking space between the type of object and the number, like this: \verb'see Figure~\ref{fig:example}'.

The \verb'\autoref{}' command can be used to automatically fill out the type of object, saving on typing.
It also causes the link to cover the whole phrase rather than just the number, but for that reason is only suitable for single cross-references rather than ranges.
For example, \verb'\autoref{tab:journal_abbr}' produces \autoref{tab:journal_abbr}.

\subsection{Citations}
\label{sec:cite}

MNRAS uses the Harvard -- author (year) -- citation style, e.g. \citet{author2013}.
This is implemented in \LaTeX\ via the \verb'natbib' package, which in turn is included via the \verb'usenatbib' package option (see section~\ref{sec:options}), which should be used in all papers.

Each entry in the reference list has a `key' (see section~\ref{sec:ref_list}) which is used to generate citations.
There are two basic \verb'natbib' commands:
\begin{description}
 \item \verb'\citet{key}' produces an in-text citation: \citet{author2013}
 \item \verb'\citep{key}' produces a bracketed (parenthetical) citation: \citep{author2013}
\end{description}
Citations will include clickable links to the relevant entry in the reference list, if supported by your \LaTeX\ compiler.

\defcitealias{smith2014}{Paper~I}
\begin{table*}
 \caption{Common citation commands, provided by the \texttt{natbib} package.}
 \label{tab:natbib}
 \begin{tabular}{lll}
  \hline
  Command & Ouput & Note\\
  \hline
  \verb'\citet{key}' & \citet{smith2014} & \\
  \verb'\citep{key}' & \citep{smith2014} & \\
  \verb'\citep{key,key2}' & \citep{smith2014,jones2015} & Multiple papers\\
  \verb'\citet[table 4]{key}' & \citet[table 4]{smith2014} & \\
  \verb'\citep[see][figure 7]{key}' & \citep[see][figure 7]{smith2014} & \\
  \verb'\citealt{key}' & \citealt{smith2014} & For use with manual brackets\\
  \verb'\citeauthor{key}' & \citeauthor{smith2014} & If already cited in close proximity\\
  \verb'\defcitealias{key}{Paper~I}' &  & Define an alias (doesn't work in floats)\\
  \verb'\citetalias{key}' & \citetalias{smith2014} & \\
  \verb'\citepalias{key}' & \citepalias{smith2014} & \\
  \hline
 \end{tabular}
\end{table*}

There are a number of other \verb'natbib' commands which can be used for more complicated citations.
The most commonly used ones are listed in Table~\ref{tab:natbib}.
For full guidance on their use, consult the \verb'natbib' documentation\footnote{\url{http://www.ctan.org/pkg/natbib}}.

If a reference has several authors, \verb'natbib' will automatically use `et al.' if there are more than two authors. However, if a paper has exactly three authors, MNRAS style is to list all three on the first citation and use `et al.' thereafter. If you are using \bibtex\ (see section~\ref{sec:ref_list}) then this is handled automatically. If not, the \verb'\citet*{}' and \verb'\citep*{}' commands can be used at the first citation to include all of the authors.

\subsection{The list of references}
\label{sec:ref_list}

It is possible to enter references manually using the usual \LaTeX\ commands, but we strongly encourage authors to use \bibtex\ instead.
\bibtex\ ensures that the reference list is updated automatically as references are added or removed from the paper, puts them in the correct format, saves on typing, and the same reference file can be used for many different papers -- saving time hunting down reference details.
An MNRAS \bibtex\ style file, \verb'mnras.bst', is distributed as part of this package.
The rest of this section will assume you are using \bibtex.

References are entered into a separate \verb'.bib' file in standard \bibtex\ formatting.
This can be done manually, or there are several software packages which make editing the \verb'.bib' file much easier.
We particularly recommend \textsc{JabRef}\footnote{\url{http://jabref.sourceforge.net/}}, which works on all major operating systems.
\bibtex\ entries can be obtained from the NASA Astrophysics Data System\footnote{\label{foot:ads}\url{http://adsabs.harvard.edu}} (ADS) by clicking on `Bibtex entry for this abstract' on any entry.
Simply copy this into your \verb'.bib' file or into the `BibTeX source' tab in \textsc{JabRef}.

Each entry in the \verb'.bib' file must specify a unique `key' to identify the paper, the format of which is up to the author.
Simply cite it in the usual way, as described in section~\ref{sec:cite}, using the specified key.
Compile the paper as usual, but add an extra step to run the \texttt{bibtex} command.
Consult the documentation for your compiler or latex distribution.

Correct formatting of the reference list will be handled by \bibtex\ in almost all cases, provided that the correct information was entered into the \verb'.bib' file.
Note that ADS entries are not always correct, particularly for older papers and conference proceedings, so may need to be edited.
If in doubt, or if you are producing the reference list manually, see the MNRAS instructions to authors$^{\ref{foot:itas}}$ for the current guidelines on how to format the list of references.

\section{Appendices and online material}

To start an appendix, simply place the \verb'