This thesis introduces a comprehensive methodology to improve electromagnetic transient (EMT) modelling of power cables systems. Several improved modelling and validation techniques are proposed at the parameter estimation, time domain simulation and validation stages of the EMT modelling of transmission lines. A novel approach is developed to model sector-shaped cables in electromagnetic transient type programs. First, the applicability of elemental sub-conductor technique is extended to accurately calculate the frequency dependent impedances of sector-shaped cables. The derived admittance and propagation characteristics of the sector-shaped cable are fitted with rational functions using the method of vector fitting in an EMT-type program. The time domain simulations are validated with the numerical inverse Laplace transform method. A novel frequency domain approach is presented to model cascaded transmission systems. The procedure is based on obtaining four composite propagation functions representing the cascaded system. The performance of the technique does not diminish with increased number of cascaded segments and it preserves the intrinsic details of each line segment. This method is capable of modelling cascaded overhead lines or cables with different characteristic admittances and line lengths. This method can be used to validate EMT models of cascaded transmission systems. An improved generalized transmission line model is introduced which is capable of accommodating time steps greater than the travel time of the line. The time step of the conventional EMT models of transmission lines is constrained by the smallest travel time of the line. When the high frequency transients at the line terminations are not of interest, inaccurate nominal π equivalents are used with large time steps to reduce the computational burden. The proposed model not only is more accurate than the π equivalents, but also degenerates to the conventional frequency dependent EMT model when used with time steps smaller than the travel time. Therefore, the proposed model is highly convenient as it can be used for all types of EMT simulations without resorting to nominal π equivalents when the large simulation time steps must be used to reduce computational burden. / February 2017