The zirconia ceramic restoration (ZCR) is used as substitutes for the metal-ceramic restoration. Clinical studies demonstrating of ZCRs showed a high fracture incidence of veneering layer than metal-ceramic restorations. This attributed to the low bond strength of zirconia to veneering ceramic as a result of lacking of glass content in its matrix. Surface treatment was proposed to improve the bonding strength between zirconia and veneering ceramic. Several studies revealed that some treatment such as airborne particle abrasion (APA) is responsible for generating chipping of veneering ceramic. The study aimed to develop a new zirconia coatings to increase bonding strength between zirconia substrate and veneering porcelain. Three groups of 150 zirconia specimens (discs, rods and crowns) were divided according to the type of surface treatments; Group M1 designated unsintered zirconia specimens coated with a mixture of glaze porcelain powder and partially-sintered zirconia powder prepared at 1100 °C with two particle sizes; A (26 ± 0.3 μm) and size B (47 ± 0.5 μm), group M2 was coated with mixture of ceramic liner paste and same selected particle sizes of partially-sintered zirconia powders. The coated groups (M1 and M2) subdivided according to powder sizes into M1A, M1B, M2A and M2B. APA group (as control) of sintered zirconia specimens abraded with 50 μm aluminium oxide powder. Surface area roughness (Sa) was measured by surface texture analyser. Surface morphology and elemental composition were analysed by scanning electron microscope (SEM). Crystallographic phases were identified by X-ray diffraction (XRD). The coefficient of thermal expansion (CTE) was assessed by a thermomechanical analyser. For shear bond strength (SBS) test was evaluated by a universal testing machine (UTM). The fracture strength (FS) test prepared and measured by Cercon CAD/CAM system and UTM. Fractographic analysis for SBS and FS were examined by a stereomicroscope and SEM. For each test and measurement, 10 specimens were used per group. Finite element analysis (FEA) was used to simulate and predict the iv stress distributions of the static oblique load in coated and non-coated zirconia crown models. The collected data were analysed by one-way ANOVA and Tukey HSD test at (P<0.05). Surface roughness results showed significant differences among all groups (P<0.000). The M1B group exhibited a higher Sa value (10.33 ± 0.59 μm) among the tested groups. XRD analysis showed tetragonal and monoclinic phases in the control group while only tetragonal was detected in the coated groups. The SBS values for coated groups were higher than the control group. The M1B group showed higher and significant differences in SBS value (37.54 ± 4.38 MPa) among other tested groups (P<0.05). The FS test indicated that coated M1B group (647.92 ± 97.33 N) higher than the other groups (P<0.00). FEA showed the coated crown models have lower stress level than the non-coated model. The new coating (M1B) by airbrush spray technique considered as an alternative way to APA treatment to improve the bond strength of zirconia substrate to veneering ceramic and indicating as an applicable surface treatment for improving the clinical performance of the coated ZCR.