Characterising the joint distribution of extremes of ocean environmental variables such as significant wave height (HS) and spectral peak period (TP) is important for understanding extreme ocean environments and in the design and assessment of marine and coastal structures. Many applications of multivariate extreme value analysis adopt models that assume a particular form of extremal dependence between variables without justification. Models are also typically restricted to joint regions in which all variables are extreme, but regions where only a subset of variables is extreme can be equally important for design. The conditional extremes model of Heffernan and Tawn (2004) provides one approach to overcoming these difficulties.Here, we extend the conditional extremes model to incorporate covariate effects in all of threshold selection, marginal and dependence modelling. Quantile regression is used to select appropriate covariate-dependent extreme value thresholds. Marginal and dependence modelling of extremes is performed within a penalised likelihood framework, using a Fourier parameterisation of marginal and dependence model parameters, with cross-validation to estimate suitable model parameter roughness, and bootstrapping to estimate parameter uncertainty with respect to covariate.We illustrate the approach in application to joint modelling of storm peak HS and TP at a Northern North Sea location with storm direction as covariate. We evaluate the impact of incorporating directional effects on estimates for return values, including those of a structure variable, similar to the structural response of a floating structure. We believe the approach offers the ocean engineer a straightforward procedure, based on sound statistics, to incorporate covariate effects in estimation of joint extreme environmental conditions. © 2013 Elsevier B.V..